Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools

You are here: metafor.ltas.ulg.ac.be » doc » Documentation » elements » General Points » Volume element
Trace:
doc:user:elements:volumes:volumeelement

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
Last revisionBoth sides next revision
doc:user:elements:volumes:volumeelement [2015/10/13 21:03] – [Volume[2|3]DElement] wauteletdoc:user:elements:volumes:volumeelement [2020/12/28 09:46] – [QuadTetraVolume3DElement] tanaka
Line 3: Line 3:
 ===== Introduction ===== ===== Introduction =====
  
-In this section, Metafor volume element are described. To them is associated a ''[[doc:user:elements:volumes:volumeinteraction|FieldApplicator]]'' interaction.+In this section, Metafor volume element are described. A ''[[doc:user:elements:volumes:volumeinteraction|FieldApplicator]]'' interaction is associated to them.
  
 ===== Volume[2|3]DElement ===== ===== Volume[2|3]DElement =====
Line 9: Line 9:
 === Description === === Description ===
  
-''Volume2DElement'' and ''Volume3DElement'' are basic elements in Metafor. ''Volume2DElement'' is a quadrangle with 4 nodes in 2D, when ''Volume3DElement''is a hexaedron with 8 nodes in 3D. Each element has respectively 8 and 24 degrees of freedom.+''Volume2DElement'' and ''Volume3DElement'' are basic elements in Metafor. ''Volume2DElement'' is a quadrangle with 4 nodes in 2D, while ''Volume3DElement'' is a hexaedron with 8 nodes in 3D. Each element has respectively 8 and 24 degrees of freedom.
  
 By default, if $n$ is the degree of mechanical interpolation and $d$ the element dimension, the number of integration points in the deviatoric part of the stress field is $(n+1)^{d}$, unless stated otherwise with the definition of ''NIPDKSI'', ''NIPDETA'' or ''NIPDZETA'' (see [[#Parameters ]]). \\ For example, the interpolation of ''Volume2DElement'' is of the first degree ($n=1$), and it is a 2D element ($d=2$), leading to 4 integration points. By default, if $n$ is the degree of mechanical interpolation and $d$ the element dimension, the number of integration points in the deviatoric part of the stress field is $(n+1)^{d}$, unless stated otherwise with the definition of ''NIPDKSI'', ''NIPDETA'' or ''NIPDZETA'' (see [[#Parameters ]]). \\ For example, the interpolation of ''Volume2DElement'' is of the first degree ($n=1$), and it is a 2D element ($d=2$), leading to 4 integration points.
Line 16: Line 16:
  
 === Remarks === === Remarks ===
 +
 == Centrifugal forces == == Centrifugal forces ==
 +
   * The angular velocity ''OMEGA'' is used when loading comes from centrifugal forces (whatever the integration scheme).   * The angular velocity ''OMEGA'' is used when loading comes from centrifugal forces (whatever the integration scheme).
   * Both start and end points of the rotation axis can be moved (see fixation and loading).   * Both start and end points of the rotation axis can be moved (see fixation and loading).
 +
 == Initial Rotative Balancing ==  == Initial Rotative Balancing == 
  
-= Old Metafor Version <= 2422 = +== Old Metafor Version <= 2422 == 
- * When ''MDE_IQSI=1'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, they are replaced by the real structural rotation.+ 
 +  * When ''MDE_IQSI=1'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, they are replaced by the real structural rotation.
   * If ''MDE_IQSI=2'' and ''MDE_NDYN=2'': The quasi-static equilibrium phase is applied, but the frame of the blade is still used.    * If ''MDE_IQSI=2'' and ''MDE_NDYN=2'': The quasi-static equilibrium phase is applied, but the frame of the blade is still used. 
-= New Metafor Version > 2422 =+ 
 +== New Metafor Version > 2422 == 
 + 
 + * ''MDE_IQSI=0'' :  
 +<code> 
 +ti.setUseInitialRotationBalancing(False)  
 +</code> 
 + 
 + * ''MDE_IQSI=1'' :  
 +<code> 
 +ti.setUseInitialRotationBalancing(True)  
 +ti.setShiftToFixedFrame(True)  
 +</code> 
 + * ''MDE_IQSI=2'' : 
 +<code> 
 +ti.setUseInitialRotationBalancing(True)  
 +ti.setShiftToFixedFrame(False)  
 +</code> 
 + * ''MDR_DTMR'' : 
 +<code> 
 +ti.setInitialRotationFactorIncrement(_initialRotationFactorIncrement) 
 +</code> 
 + 
 +ti is a reference towards the current dynamic time integration scheme. 
  
 == Integration Method  == == Integration Method  ==
 +
   * :!: In axisymmetric modeling, [[doc:user:elements:volumes:elements_formulation#Selective Reduced Integration with Pressure Report]] ''VES_CMVIM_SRIPR'' must be used.   * :!: In axisymmetric modeling, [[doc:user:elements:volumes:elements_formulation#Selective Reduced Integration with Pressure Report]] ''VES_CMVIM_SRIPR'' must be used.
   * :!: [[doc:user:elements:volumes:elements_formulation#EAS Formulation]] is not valid in axisymmetric modeling, some EAS modes are missing along $\theta$.   * :!: [[doc:user:elements:volumes:elements_formulation#EAS Formulation]] is not valid in axisymmetric modeling, some EAS modes are missing along $\theta$.
Line 73: Line 101:
 Quadratic, 6-nodes triangular element in 2D, 12 dofs. The standard formulation must be used (''CAUCHYMECHVOLINTMETH = VES_CMVIM_STD''). Quadratic, 6-nodes triangular element in 2D, 12 dofs. The standard formulation must be used (''CAUCHYMECHVOLINTMETH = VES_CMVIM_STD'').
  
-By default, stresses are integrate over 3 integration points in the deviatoric part.+By default, stresses are integrated over 3 integration points in the deviatoric part.
  
 __Remark__: To use these higher order elements, the mesh must be also defined as second or third degree.  __Remark__: To use these higher order elements, the mesh must be also defined as second or third degree. 
Line 81: Line 109:
 Quadratic, 10-nodes tetrahedron in 3D, 30 dofs. As with quadratic elements in 2D, the standard formulation must be used (''CAUCHYMECHVOLINTMETH = VES_CMVIM_STD''). Quadratic, 10-nodes tetrahedron in 3D, 30 dofs. As with quadratic elements in 2D, the standard formulation must be used (''CAUCHYMECHVOLINTMETH = VES_CMVIM_STD'').
  
-By default, stresses are integrate over 4 integration points in the deviatoric part.+By default, stresses are integrated over 4 integration points in the deviatoric part.
  
 __Remark__: To use these higher order elements, the mesh must be also defined as second or third degree.  __Remark__: To use these higher order elements, the mesh must be also defined as second or third degree. 
Line 101: Line 129:
 ^      Code Metafor        ^          Description                                      ^  Dependency  ^ ^      Code Metafor        ^          Description                                      ^  Dependency  ^
 | ''MATERIAL''             | Number of the volume material to consider   |      -    | | ''MATERIAL''             | Number of the volume material to consider   |      -    |
-| ''STIFFMETHOD''          | Method used to compute the stiffness matrix \\ = ''STIFF_ANALYTIC'' : analytic matrix (default) \\ = ''STIFF_NUMERIC''  : numerical matrix                    |     -       |                      +| ''STIFFMETHOD''          | Method used to compute the stiffness matrix \\ = ''STIFF_ANALYTIC'' : analytic matrix (default) \\ = ''STIFF_NUMERIC''  : numerical matrix                    |     -       |     
 +| ''DAMPSTIFF''            | Part of the Stiffness matrix used to build damping matrix used by [[doc:user:integration:scheme:dynimpl#damped_alpha-generalized_family|Damped Time Integration scheme]] \\ usual value : 1.0e-7 - 1.0e-5  |     Time (need DMUPERSTEP or DMUPERSTAGE)      | 
 +| ''DAMPMASS''             | Part of the Mass matrix used tobuild damping matrix used by [[doc:user:integration:scheme:dynimpl#damped_alpha-generalized_family|Damped Time Integration scheme]] \\ usual value : 1.0e3 - 1.0e5    |     Time  (need DMUPERSTEP or DMUPERSTAGE)     |
 | ''OMEGA''                | Angular speed (°/s) to simulate a rotating frame (calculation of centrifugal and Coriolis forces). \\ If ''MDE_IQSI=1'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, they are replaced by the real structural rotation. \\ If ''MDE_IQSI=2'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, the simulation is carried in the moving frame, so centrifugal and Coriolis forces are taken into account in the model.  |  time | | ''OMEGA''                | Angular speed (°/s) to simulate a rotating frame (calculation of centrifugal and Coriolis forces). \\ If ''MDE_IQSI=1'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, they are replaced by the real structural rotation. \\ If ''MDE_IQSI=2'' and ''MDE_NDYN=2'', centrifugal forces are computed during a quasi-static equilibrium phase. Then, the simulation is carried in the moving frame, so centrifugal and Coriolis forces are taken into account in the model.  |  time |
 | ''OMEGA_PT1''            | Number of the first point which defines the rotation axis      |          -           | | ''OMEGA_PT1''            | Number of the first point which defines the rotation axis      |          -           |
 | ''OMEGA_PT2''            | Number of the second point which defines the rotation axis    |          -           | | ''OMEGA_PT2''            | Number of the second point which defines the rotation axis    |          -           |
 | ''CORIOLIS''             | Boolean to take into account Coriolis forces in a moving frame, during a dynamic simulation.\\ True (default) / False   |   -   | | ''CORIOLIS''             | Boolean to take into account Coriolis forces in a moving frame, during a dynamic simulation.\\ True (default) / False   |   -   |
-| ''GRAVITY_X'', ''GRAVITY_Y'', ''GRAVITY_Z''  | Gravity       |        time |+| ''GRAVITY_X'', ''GRAVITY_Y'', ''GRAVITY_Z''  | Gravity       |        time   |
 | ''CAUCHYMECHVOLINTMETH'' | [[doc:user:elements:volumes:elements_formulation|Method]] used to integrate the mechanical part of volume elements (stresses) \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_standard|VES_CMVIM_SRI]]'' \\ = ''[[doc:user:elements:volumes:elements_formulation#sous_integration_selective_avec_report_de_pression|VES_CMVIM_SRIPR]]'' (default) \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_standard|VES_CMVIM_STD]]'' \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_eas|VES_CMVIM_EAS]]'' |  -  | | ''CAUCHYMECHVOLINTMETH'' | [[doc:user:elements:volumes:elements_formulation|Method]] used to integrate the mechanical part of volume elements (stresses) \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_standard|VES_CMVIM_SRI]]'' \\ = ''[[doc:user:elements:volumes:elements_formulation#sous_integration_selective_avec_report_de_pression|VES_CMVIM_SRIPR]]'' (default) \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_standard|VES_CMVIM_STD]]'' \\ = ''[[doc:user:elements:volumes:elements_formulation#formulation_eas|VES_CMVIM_EAS]]'' |  -  |
 | ''CONSERVINGMETHOD''     | Use of ''[[doc:user:integration:scheme:dynimpl#conservative algorithm]]'' \\ = ''VES_WITHOUTCORRECTION'' : plastic correction tensors neglected (default) \\ = ''VES_WITHCORRECTION'' : plastic correction tensors considered |  -  | | ''CONSERVINGMETHOD''     | Use of ''[[doc:user:integration:scheme:dynimpl#conservative algorithm]]'' \\ = ''VES_WITHOUTCORRECTION'' : plastic correction tensors neglected (default) \\ = ''VES_WITHCORRECTION'' : plastic correction tensors considered |  -  |
Line 117: Line 147:
 | ''TEAS'' | Transformation of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes from the isoparametric space \\  =0 : Simo-Armero transformation (default) \\ =1 : Glaser-Armero transformation [do not use, still under development] |  -  | | ''TEAS'' | Transformation of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes from the isoparametric space \\  =0 : Simo-Armero transformation (default) \\ =1 : Glaser-Armero transformation [do not use, still under development] |  -  |
 | ''EEAS'' | Extrapolation of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes \\ =0 : modes set to 0 for each time step (safe but slow). \\ =1 : classical extrapolation (default) \\ =2 : initialization to the value corresponding to the previous step. \\ =3 : set to 0 for each iteration |  -  | | ''EEAS'' | Extrapolation of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes \\ =0 : modes set to 0 for each time step (safe but slow). \\ =1 : classical extrapolation (default) \\ =2 : initialization to the value corresponding to the previous step. \\ =3 : set to 0 for each iteration |  -  |
-| ''PEAS'' | Accuracy of resolution of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes (default: 1.0e-8)|  -  |+| ''PEAS'' | Accuracy of resolution of [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] modes (default: 1.0e-8) unfortunately NOT adimensional!!! \\ => = 1.0e-8 for "small tests" in mm \\ => = 1.0e-6 for "real tests" in mm \\ => = 1.0e-9 for "real tests" in m |  -  |
 | ''VERBOSE'' | (bool) Debug information concerning resolution of  [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] (default: false)|  -  | | ''VERBOSE'' | (bool) Debug information concerning resolution of  [[doc:user:elements:volumes:elements_formulation#formulation_eas|EAS]] (default: false)|  -  |
 +
doc/user/elements/volumes/volumeelement.txt · Last modified: 2020/12/29 18:27 by tanaka
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki