FEDEASLab

Utilities > PreProcessing > Sections > Create_MPMesh4Circ

Create_MPMesh4Circ

Purpose

generate the mesh for a circular disc or annulus

Synopsis

[R1,R2,thfib,dth] = Create_MPMesh4Circ (Ri,Ro,nrfib,nthfib,MeshOpt)

Description

CREATE_MPMESH4CIRC generate the mesh for a circular disc or annulus
  [R1,R2,THFIB,DTH] = CREATE_MPMESH4CIRC (RI,RO,NRFIB,NTHFIB,MESHOPT)
  the function generates the mesh for a circular disc or annulus with inner radius RI and
  outer radius RO; there are NRFIB radial and NTHFIB circumferential subdivisions;
  the function returns the radial subdivision in vectors R1 and R2 and the circumferential
  subdivisions for each ring in the cellarray THFIB;
  DTH is a vector with the circumfererential angle increment for each ring;
  the discretization is controlled by the optional data structure MESHOPT with the following fields:
  MESHOPT.BASIC = true (for NRFIB x NTHFIB mesh) or false (for NRFIB x variable NTHFIB mesh)
         .CP    = true (with IP at center)       or false (without IP at center)
         .UW    = true (for uniform IP weight)   or false (for variable IP weight)
         .ThInc = increment in number of circumferential fibers/quadrant
         .Phi   = angle of first integration point in circumferential direction relative to y-axis
                  if Phi='var' each ring is offset by half the angle between successive IPs
         .Nthst = initial number of circumferential fibers/quadrant (default = 4*Ri/Ro)
  Default values for MESHOPT fields:
         BASIC = true, CP = false, UW = false, ThInc = 1, Phi = 0
  the function returns

Cross-Reference Information

This function calls:

This function is called by:
  • Plot_Mesh4Circ plots mesh for a circular disc or annulus in current window
  • MP2dInt4Circ integration points and weights for 2d-midpoint rule of circular disc or annulus