$$$$ MODE NOTICE FANDEUR 14/01/22 21:15:10 7891 DATE 14/01/22 Operateur MODE (MODELISER) Voir aussi : OPTI -------------------------- MATE CARA MODL1 = MODE GEO1 FOR1 MAT1 ( ...MATn ) ... | ... ('EPSI' MOT1 ) ( ELEM1 ... ELEMn ) ... | | ... ('INCO' MOT2 (MOT3) ) ... | | ( | 'NUME_LOI' ILOI1 | ) ( 'PARA_LOI' LMOTS1 ) | 'NOM_LOI' CLOI16 | ( 'C_MATERIAU' LMOTS2 ) ( 'C_VARINTER' LMOTS3 ) ('INTEGRE' N1 ) ( 'CONS' MOT4 ) ( 'DPGE' P1) ('PHAS' MOT5 ) ; | MODL2 ( ... MODLn) ; Description : ___________ The object of the MODE (MODELISER) operator is to associate a formulation, a model of material behavior, a type of finite element to be used, and possibly a name of constituent, with a mesh.The number of integration points for DKT shell element may be given ('INTEGRE' N1) and the point support of additional unknowns for generalized plane strains must be given ( 'DPGE' P1). In the case of an external model of non linear mechanical behavior, conceived by the user, additionnal informations have to be given : number or name attributed to the behavior law, list of the external parameters, list of the material components, list of the internal variables of the law. The MELANGE PARALLELE modelling enables to associate elementary modellings labelled with phase names and add the effects affected linearly by the phase parameters proportions. The different datas and keywords must be ordered as shown in the general syntax above. Contents : _________ GEO1 : geometry (MAILLAGE type) FOR1 : formulation defined by one or several words chosen from: - simple formulations : | 'MECANIQUE' | 'LIQUIDE' | 'POREUX' | 'CONTACT' | 'THERMIQUE' | 'DARCY' | 'MAGNETODYNAMIQUE' | 'NAVIER_STOKES' | 'EULER' | 'MELANGE' | 'FISSURE' | 'THERMOHYDRIQUE' | 'LIAISON' | 'DIFFUSION' - coupled formulation : 'LIQUIDE' 'MECANIQUE' MAT1 (...MATn ) : type of material with as many words as necessary (MOT type). The possible types are listed further down. (EPSI MOT1 ) : This data allow to change the way of computing strains. By default the way is using 'QUADRATIQUE' terms. This way may be also defineded by OPTION EPSI ... Possibiliies are for MOT1: LINEAIRE (Linear) QUADRATIQUE (Quadratic) TRUESDELL JAUMANN UTILISATEUR (user) In any case the code gives to the COMP operator (material behaviour) the gradient of displacements from initial state at the begining of the time step and at the end of the time step. In case of UTILISATEUR the strains are linear (LINEAIRE) but computed with the semi-step configuration ( ELEM1...ELEMn ): specific finite element to be used (MOT type). By default, we will use the finite element whose name is the same as that of the geometrical support. This list is required for shells and joints. The possible types are listed further down. for the 'NAVIER_STOKES' model we can use the generic names LINE (linear) MACRO (iso-P2 iso-P1) or QUAF (quadratique for the fluids). ( INCO MOT2 (MOT3) ) : only for the DIFFUSION formulation: names of the nodal (primal) unknown MOT2 and its dual MOT3. ATTENTION : the name of the nodal unknown is limited to 2 letters (CO by default) to be able to name the components of the gradient (CO,X...). MOT3 is optionnal. By default, the name of the dual is 'QXX', if 'XX' is the name of the primal. Group of 4 datas defining an external behavior's model, conceived by the user, for the 'MECANIQUE' formulation : => 'NON_LINEAIRE' 'UTILISATEUR' law, which can be associated to - 'ELASTIQUE' 'ISOTROPE', - 'ELASTIQUE' 'ORTHOTROPE', - 'ELASTIQUE' 'ANISOTROPE', - 'ELASTIQUE' 'UNIDIRECTIONNEL'. => behavior's laws of the group 'VISCO_EXTERNE', which can be associated only to 'ELASTIQUE' 'ISOTROPE'. ( ILOI1 ) : Number attributed to the behavior's law. ENTIER type :1 .LE. ILOI1 .LE. 999999 ( CLOI16 ) : Name attributed to the behavior's law. MOT type with at the most 16 characters ILOI1 OR CLOI16 = Compulsory data for all external model. ( LMOTS1 ) : Object of type LISTMOTS, giving the list of the names of the external parameters of the behavior's law. Optional data. If the temperature 'T ' is one of the external parameters of the model, it has to be declared first of the list. Redundant parameters are forbidden. ( LMOTS2 ) : Object of type LISTMOTS, giving the list of the names of the material components of the behavior's law. Compulsory data for a 'NON_LINEAIRE' 'UTILISATEUR' law : all the components of the behavior (linear and non linear domains) must be declared. Redundant components are forbidden. Ignored data in other cases : all the laws of the 'VISCO_EXTERNE' group have the same material components as the 'ELASTIQUE' 'ISOTROPE' behavior. ( LMOTS3 ) : Object of type LISTMOTS, giving the list of the names of the internal variables of the behavior's law. Optional data. Redundant internal variables are forbidden. Laws of the 'VISCO_EXTERNE' group have 4 pre-defined internal variables : 'EC0 ', 'ESW0', 'P ' and 'QTLD' In this case, the LMOTS3 object gives the names of the additionnal internal variables, if there are. Notes : => A 'NON_LINEAIRE' 'UTILISATEUR' law can be applied to all type of finite element, excepted shells whithout any integration point in the thickness. For other elements than tridimensionnal, the usefull geometrical properties can be declared among the external parameters of the law. => Laws of the 'VISCO_EXTERNE' group can be applied only to tridimensionnal elements. => Redundant names between material components and internal variables are forbidden. => If the 'NON_LINEAIRE' 'UTILISATEUR' external model computes inelastic strains, these can be passed out if they have been declared among the internal variables. (N1 ) : number of integration points (between 1 and 15) in the thickness when using the shell elements with numerical integration in the thickness . It is recommended that an odd number be used. ( MOT4 ) : name of constituent (MOT type, 16 characters at the most). This possibility will be used when the association of several models with a mesh object in the same calculation is requested. The MOT4 name will then enable the user to identify the constituents. By default, an unique constituent is generated for each model. ( MOT5 ) : name of phase (MOT type, 8 characters at the most) This possibility will be used when the association of several models with a mesh object in the same calculation is requested in order to add the effects like in phases mix. MOT5 is thus a phase name. The default value is ' '. MODL2 ( .. MODn) : model object labelled with a phase name (MMODEL type) MODL1 : generated model object (MMODEL type) ----------------------------- | Material name in MECANIQUE | ----------------------------- First of all, it is required that the linear behavior be defined. Then, if the case arises, the non-linear behavior will be defined according to the following instructions. The following data, for instance : ELASTIQUE ISOTROPE PLASTIQUE ISOTROPE correspond to a material whose linear behavior is both elastic and isotropic, and whose non-linear behavior is plastic, according to Von Mises' isotropic stress hardening model. - Linear behaviors : ------------------ * 'ELASTIQUE' | ('ISOTROPE') (option by default) | 'ORTHOTROPE' | 'ANISOTROPE' | 'UNIDIRECTIONNEL' | 'HOMOGENEISE' | 'SECTION' | 'ARMATURE' | 'MODAL' | 'STATIQUE' The case SECTION can be used for Timoschenko beam elements to describe the linear and non linear behaviours of the cross section.If the beam element is bidimensional (plain stresses or plain strain), the Oy local axe of the beam element corresponds to the Ox axe of the section. The case ARMATURE can be used for the (active or passive) tendons of the reinforced concrete ( finite element : BARRE, geometrical support : SEG2). The 'MODAL' material label enables to describe the displacements on the eigenmodes basis (see VIBR operator). The associated geometry is a single point mesh. Likely the 'STATIQUE' material label enables to add to the set of eigenmodes the elastic solution obtained from a local loading problem. - Non-linear behaviors : ---------------------- * 'NON_LINEAIRE' : non linear elasticity or external non linear behavior. No default option. Must be followed by a word chosen from : |---------------------------------------------------------------- | 'EQUIPLAS' : for loading case the behaviour is strictly the | same as the PLASTIQUE ISOTROPE with isotropic | stress hardening. | |---------------------------------------------------------------- | 'UTILISATEUR' : external non linear behavior's model, conceived | by the user. | | A number ILOI1 or a CLOI16 name has to be attributed to the | external behavior's law. | The list of the external parameters of the law is given by the | LMOTS1 object, if such parameters exist. | The list of the material components is given by LMOTS2. This | list covers the linear and non linear domains of the behavior. | The list of the internal variables of the law is given by the | LMOTS3 object, if such variables exist. | | The 'NON_LINEAIRE' 'UTILISATEUR' model is integrated by a | formulation coded by the user in the external subroutine UMAT. | This subroutine is called by COMP operator to compute the | solution at each integration point of the model. | With the following given informations : | - an initial state, at time t0 : stresses, total strains, | internal variables ; | - a time increment dt ; | - an increment of total strains, imposed from t0 to t0+dt ; | - the values of the external parameters at t0 and t0+dt ; | the UMAT subroutine computes the final state at t0+dt : new | stresses and internal variables (inelastic strains might be | among them). | | Remarks : | | Some laws defined by its number are pre-existing in the code. | We can list : | 5 for PLASTIQUE ISOTROPE | 7 for VISCOPLASTIQUE' 'GATT_MONERIE' | 10 for ELASTICITE ISOTROPE ISOTHERME | 11 for ELASTICITE ISOTROPE NON ISOTHERME | 12 for ELASTICITE ORTHOTROPE ISOTHERME | 21,22 for 'FLUAGE' 'POLYNOMIAL' (massif) | 31 for model Mooney-Rivlin(2D Plan CONT, 3D massive) | 32 for model neo-hookien (hyperelastic) | 33 for model GD (hyperelastic) | 34 for model Hart-Smith (hyperelastic) | 35 for model Biderman (hyperelastic) | 36 for model 8chaines (hyperelastic) | | exemples of use of these laws can be found in the .dgibi files | | Operating rules : | --------------- | 1 - Code the UMAT subroutine and its dependancies : | - Take the example UMAT subroutine integrated in CAST3M. | It is necessary to keep the first card of the example | subroutine, in order to compile successfully. | - Keep the calling interface of the routine. | 2 - Compile the external subroutine UMAT and its dependancies, | then edit links with the rest of the code. | 3 - Use the MODE and COMP operators, the PASAPAS procedure | with respect to their syntax in Gibiane language. | | Interface of the external subroutine UMAT : | ----------------------------------------- | SUBROUTINE UMAT ( STRESS, STATEV, DDSDDE, SSE, SPD, SCD, | & RPL, DDSDDT, DRPLDE, DRPLDT, | & STRAN, DSTRAN, TIME, DTIME, | & TEMP, DTEMP, PREDEF, DPRED, | & CMNAME, NDI, NSHR, NTENS, NSTATV, | & PROPS, NPROPS, COORDS, | & DROT, PNEWDT, CELENT, DFGRD0, DFGRD1, | & NOEL, NPT, LAYER, KSPT, KSTEP, KINC ) | IMPLICIT REAL*8(A-H,O-Z) | | CHARACTER*16 CMNAME | | INTEGER NDI, NSHR, NTENS, NSTATV, NPROPS, | & NOEL, NPT, LAYER, KSPT, KSTEP, KINC | | REAL*8 STRESS(NTENS), STATEV(*), | & DDSDDE(NTENS,NTENS), | & SSE, SPD, SCD, | & RPL, DDSDDT(NTENS), DRPLDE(NTENS), DRPLDT, | & STRAN(NTENS), DSTRAN(NTENS), | & TIME(2), DTIME, | & TEMP, DTEMP, PREDEF(*), DPRED(*), | & PROPS(NPROPS), | & COORDS(3), | & DROT(3,3), | & PNEWDT, | & CELENT, | & DFGRD0(3,3), DFGRD1(3,3) | | | IN/OUT : STRESS : REAL*8(NTENS), stress tensor | On input : stresses at t0 | On output : stresses at t0+dt | IN/OUT : STATEV : REAL*8(*), internal variables | On input : internal variables at t0 | On output : internal variables at t0+dt | | | OUT : DDSDDE : REAL*8(NTENS,NTENS), jacobian matrix of the | model at t0+dt | Optional output, not used further by CAST3M | | IN/OUT : SSE : REAL*8, specific elastic strain energy | SPD : REAL*8, plastic dissipation | SCD : REAL*8, "creep" dissipation | Values at t0 on input, at t0+dt on output | Optional input/output, not used further by | CAST3M | | | OUT : RPL : REAL*8, volumetric heat generation per unit | time at t0+dt, caused by mechanical working | of the material | DDSDDT : REAL*8(NTENS), variation of the stress | increment with respect to the temperature | DRPLDE : REAL*8(NTENS), variation of RPL with respect | to the strain increments | DRPLDT : REAL*8, variation of RPL with respect to the | temperature | Optional outputs, not used further by CAST3M | These outputs are usefull for a fully | coupled temperature-displacement analysis | | | IN : STRAN : REAL*8(NTENS), total strains tensor at t0 | DSTRAN : REAL*8(NTENS), total strain increments | tensor imposed from t0 to t0+dt | | | IN : TIME : REAL*8(2), TIME(1) = 0 | TIME(2) = t0 | DTIME : REAL*8, DTIME = dt | t0 : time at last reached balance | dt : new time increment proposed by PASAPAS | to reach the mechanical balance with | the imposed increments of total strains | (DSTRAN) | | | IN : TEMP : REAL*8, temperature at t0 | DTEMP : REAL*8, increment of temperature at t0+dt | PREDEF : REAL*8(*), values of the external parameters | of the behavior's law at t0 | DPRED : REAL*8(*), increments of external parameters | at t0+dt | | | IN : CMNAME : CHARACTER*16, identifier of the behavior law | 1) The law is identified by the number | which has been attributed to it : ILOI1, | behind the keyword 'NUME_LOI'. | The number is encoded in the last four | characters of the string CMNAME, and has to | be recovered in UMAT by an instruction such | as K4ILOI = CMNAME(13:16) | where K4ILOI is a local variable of type | CHARACTER*4 | 2) The law is identified by the name which has | been attributed to it : CLOI16, behind the | keyword 'NOM_LOI'. | The name can be directly recovered in UMAT, | because CMNAME = CLOI16 (16 characters) | | | IN : NDI : INTEGER defining the CAST3M type of analysis | | | IN : NSHR : INTEGER | Inactive input | | | IN : NTENS : INTEGER, number of components of the stress | tensor | NSTATV : INTEGER, number of internal variables | | | IN : PROPS : REAL*8(NPROPS), values of the material | constants | NPROPS : INTEGER, number of material constants | | | IN : COORDS : REAL*8(3), coordinates of the current | integration point | | | IN : DROT : REAL*8(3,3), transfer matrix from the local | solid finite element base to the general | mesh base | Inactive input for other types of finite | elements | | | OUT : PNEWDT : REAL*8, ratio between the new suggested time | increment and the time increment on input | (NEWDT = PNEWDT * DTIME) | | | IN : CELENT : REAL*8, characteristic element length | Computed as the maximum distance between two | nodes of the element | | | IN : DFGRD0 : REAL*8(3,3), gradient of displacement | tensor at t0 | DFGRD1 : REAL*8(3,3), gradient of displacement | tensor at t0+dt | | | IN : NOEL : INTEGER, number of the current element | Note : it is a local number (number of the | element in the current region of the model) | NPT : INTEGER, number of the current integration | point | | | IN : LAYER : INTEGER | KSPT : INTEGER | Inactive input | | | IN : KSTEP : INTEGER | KINC : INTEGER | These inputs have no sense when | linking to CAST3M | | Programming rule : the code of the user may use KINC as a | return flag of the UMAT subroutine, with respect to the two | rules below : | 1 - No re-initialization of KINC when entering into UMAT. | KINC has been initialized to 1 before calling UMAT. | 2 - On error, KINC is affected by a value not equal to 1 | |---------------------------------------------------------------- * 'PLASTIQUE' followed by a word chosen from : | ('ISOTROPE') : Von Mises' model with isotropic stress | hardening (option by default) | | 'PARFAIT' : Von Mises' model without stress | hardening | | 'CINEMATIQUE' : Von Mises' model with linear kinematic | stress hardening | | 'DRUCKER_PARFAIT' : Drucker-Prager's model with | associated yield and without | stress hardening | | 'DRUCKER_PRAGER' : Drucker-Prager's model with non- | associated yield and with stress | hardening | | 'BETON' : Concrete model (in plane stresses only) | | 'CHABOCHE1' : Chaboche's model with a single centre | (and isotropic stress hardening) | | 'CHABOCHE2' : Chaboche's model with two centres | (and isotropic stress hardening) | | 'TUYAU_FISSURE' : Plasticity model (ideal)/(stress | hardening) for the TUYAU FISSURE | element | | 'ENDOMMAGEABLE' : Model of elastoplastic material that | could be damaged (Lemaitre Chaboche) | | 'GAUVAIN' : Gauvain's model | | 'BILIN_MOMY' : bilinear bending or shear beam model for | 'BILIN_EFFZ' POUT and TIMO elements | | 'TAKEDA_MOMY' : takeda bending or shear beam model for | 'TAKEMO_EFFZ' POUT and TIMO elements | | 'BA1D' : Global model for Reinforced concrete columns | for POUT finite elements | | 'LINESPRING' : Plasticity model (ideal)/(stress | hardening) for the LINESPRING element | | 'UBIQUITOUS' : Ubiquitous model | | 'GLOBAL' : Plasticity model for reinforced concrete | elements which enable the global behavior | laws which differ according to the type | of loading to be taken into account | | 'CAM_CLAY' : Model with non linear elastic behaviour | in volume and plastic in volume and | distortion | | 'HUJEUX' : Model for soils and some granular media, | with plastic behavior in volume and | distortion | | 'GURSON' : Plastic flow model based on the Gurson | yield surface which is pressure dependent. | Porosity is an internal variable. This model | is used for ductile fracture by void growth. | It works only for the 3D and axisym. options. | | 'JOINT_DILATANT' : Joint model based on a Mohr-Coulomb | criterion and non associated flow rule | (only in 2D plane stress, plane strain) | | 'JOINT_SOFT' : Joint model based on a Mohr-Coulomb | criterion with softening in shear and tension | and without dilatancy. | (only in 2D plane stress, plane strain) | | 'JOINT_COAT' : Joint model in shear based on a plastic | criterion with softening and damage. | (only in 2D plane stress, plane strain) | | 'ANCRAGE_ELIGEHAUSEN' : Joint model with bond slip uniaxial law for shear | The monotonic behaviour is based on the Eligehausen model | Joint behaviour remains elastic in tension. | (only in 2D plane stress, plane strain) | | 'COULOMB' : Dilatant joint model based on a Mohr-Coulomb | criterion and associated flow rule | (works in 3D isotropic case, 2D axisym., | plane stress, plane strain) | | 'AMADEI' : Non linear incremental model with strain | softening behaviour in shear | ( en 3D isotrope, 2D axi, defo. planes) | | 'ACIER_UNI': Menegotto-Pinto steel model for elements | describing the cross section of a beam | | 'ACIER_ANCRAGE': Bond slip uniaxial law for anchorages or lap splices. | This law is based on the steel law ACIER_UNI and | the bond slip law ANCRAGE_ELIGEHAUSEN | It is valid for rod elements and beam fiber model | | 'BETON_UNI': Hognestad concrete model for elements | describing the cross section of a beam | | 'BETON_BAEL': BAEL concrete model for elements | describing the cross section of a beam | | 'MAZARS' : Scalar damage model for concrete | for elements describing the cross section | of a beam (well suited to monotonic loadings) | (model in its 3D formulation and not uniaxial) | | 'UNILATERAL' : Damage model with 2 scalar variables | for concrete for elements describing | the cross section of a beam | (well suited to cyclic loadings) | | 'FRAGILE_UNI': Damage model with softening under tension | and compression for elements | describing the cross section of a beam | | 'STRUT_UNI': Strut and Tie non-linear shear model for elements | describing the cross section of a beam | | 'CISAIL_NL' : Simplified non linear shear beam model | with softening for Timoshenko beam element | and for elements describing the cross section | of a beam | (in the Oz direction of the local beam axis) | | 'INFILL_UNI': Simplified non linear strut model without strength | in traction. | Behaviour in compression is based on damage | and plasticity. | It is unilateral in traction/compression. | This law can be used with 2 strut elements | as global model of masonry infill panel. | | 'PARFAIT_UNI': Elastic perfectly plastic uniaxial law | for elements describing the cross section of a beam | | 'PARFAIT_ANCRAGE': Bond slip uniaxial law for anchorages or lap splices. | This law is based on the steel law PARFAIT_UNI and | the bond slip law ANCRAGE_ELIGEHAUSEN | It is valid for rod elements and beam fiber model | | 'OTTOSEN' : According to the fictitious crack model approach, | OTTOSEN model is well adapted for brittle material | such as concrete or ceramic. When the damage is | localized, a size parameter is taking into account | in order to provide objective result to mesh | refinement. | | 'BETOCYCL' : Plastic model with compressive strength | degradation under cyclic loading. | The model is defined by two plasticity surfaces | (interior and exterior surfaces) with a tensile | and a compressive Rankyne-type mechanisms. | It has been developped in order to study | the behaviour of masonry structures under | in-plane cyclic loading | (only in 2D plane stress) | | 'STEINBERG': Steinberg-Cochran-Guinan steel model for | shocks ( 3D , 2D axi. , plane strain ). | | 'ZERILLI' : Zerilli-Armstrong model proposing a | mathematical formulation of the | Von Mises's yielding stress based on | the theory of the dislocations. | | 'PRESTON' : Preston-Tonks-Wallace model proposing a | mathematical formulation of the | Von Mises's yielding stress based on | the theory of the dislocations and the | dimensional analysis. | | 'HINTE' : Hinte model is a join model with degradation till | rupture. This dammaging model is not time | dependant. The finite element support is the | RACCOR element (JOI2) | | | 'J2' : Von Mises' model with exponential/linear | isotropic hardening | Consistent tangent matrix available | (Option for calculation: plan defo, axis) | | 'RH_COULOMB': Rounded Hyperbolic Mohr-Coulomb model | Perfect and associated plasticity | Consistent tangent matrix available | Substepping available (see SSTE) | (Option for calculation: plan defo, axis) | | 'MRS_LADE' : MRS-Lade model for granular noncohesive materials | Cap-model, non-associated hardening/softening cone | Consitent tangent matrix available | Substepping available (see SSTE) | (Option for calculation: plan defo, axis) | | 'VMT_FEFP' : Von Mises - Tresca model with exponential/linear | isotropic hardening, hyperelasticity and finite | deformations. | Consistent tangent matrix available (see ECFE) | (Option for calculation: plan defo, axis) | | 'RHMC_FEFP' : Rounded Hyperbolic Mohr-Coulomb plastic model | with hyperelasticity and finite deformations. | Perfect and associated plasticity. | Consistent tangent matrix available (see ECFE) | (Option for calculation: plan defo, axis) | | 'POWDER_FEFP': Elliptic density dependent plastic model with | hyperelasticity and finite deformations. | Used in cold powder compaction simulations. | Consistent tangent matrix available (see ECFE) | (Option for calculation: plan defo, axis) | | 'POWDERCAP_FEFP': Elliptic density dependent plastic model with | hyperelasticity and finite deformations. | Used in cold powder compaction simulations. | Consistent tangent matrix available (see ECFE) | (Option for calculation: plan defo, axis) | | 'DRUCKER_PRAGER_2' : Plastic model with softening for concrete. | It is formulated with two threshold | surfaces: a Rankine criterion in tension | and a Drucker-Prager criterion on compression. | | 'RICBET_UNI' : Damage model for concrete for elements describing | the cross section of a beam (well suited for cyclic | and seismic loadings). The original feature of this | model is that it accounts for local hysteretic | effects. * 'ENDOMMAGEMENT' followed by a word chosen from : | 'MAZARS' : scalar damage model for concrete | (well suited to monotonic loadings) | | | 'UNILATERAL' : damage model with 2 scalar variables | for concrete (well suited to cyclic | loadings) | | 'ROTATING_CRACK' : damage model for concrete | | 'SIC_SIC' : damage model for SiC/SiC ceramic matrix | composites | 'MVM' : Modified Von Mises damage model. | nonlocal scalar damage model for | quasibrittle materials. | | 'SICSCAL' : 'Scalar' damage model for woven SiCf/SiC composite | developped at ONERA with 3 damage variables. | Damage variables d1 and d2 correspond to cracks planes | normal to the 2 fibres directions. Damage variable d3 | corresponds to cracks in the ply plane. | | 'SICTENS' : 'Pseudo-tensorial' damage model for woven SiCf/SiC | composite developped at ONERA with 5 damage variables. | Damage variables d1 and d2 correspond to cracks planes | normal to the 2 fibres directions. Damage variable d3 | corresponds to cracks in the ply plane. variables d4 | and d5 correspond to cracks planes normal to the axes | orientated at + and - 45° from the fibres directions | (rotation in the ply plane). | | 'DAMAGE_TC' : Damage model with two scalar variables: one in tension | and one in compression. This model can be used for | cyclic applications and is regularized by the fracture | energy in ension. | | 'DESMORAT' : Anisotropic damage model for concrete (well suited to | monotonic loadings). | Can be used with the non-local theory | | 'FATSIN' : Isotropic scalar damage model for material tested under | sinusoidal fatigue (local or nonlocal implementation) | | 'RICRAG' : Isotropic continuum damage mechanics based model for con | It is recommended not only for monotonic loadings but al | cyclic ones with medium level. It has been implemented f | 2D/3D finite elements. | | 'GLRC_DM' : Isotropic continuum damage model aiming at describing th | of a representative volume element (RVE) of reinforced c | subjected to cyclic loadings. This model is expressed in | generalized stresses. It is implemented in shell element | notice that there is not softening branch. Therefore, no | localization can appear what ensures the mesh independen | | 'RICBET' : Isotropic continuum damage model aiming at describing th | of plain concrete. It has been developped for cyclic and | applications. Among their features, one can point out th | ones : disymmetry between tension and compression, perma | in tension and compression, hysteretic effects in tensio | a full linear unilateral effect. It can be used with the | approach. | | 'RICCOQ' : Simplified isotropic damage for cyclic and seismic appli | It is regularized by fracture energy technique. Simple a | * 'FLUAGE' followed by a word chosen from : | 'NORTON' : NORTON's creep model | | 'BLACKBURN' : Blackburn's creep model | | 'RCCMR_316' : creep model RCC-MR for concrete | 316-SS | | 'RCCMR_304' : creep model RCC-MR for concrete | 304-SS | | 'LEMAITRE' : Lemaitre's creep model | | 'POLYNOMIAL' : Polynomial creep model | | 'CERAMIQUE' : Ceramic's creep model | Above a transition temperature the material | creeps folowing Norton's law. | Below this temperature it behaves folowing the | Ottosen's model well adapted for brittle materials | such as concrete or ceramic. | This model can also take into account | a complete damage of the material above a given | creep deformation * 'PLASTIQUE_ENDOM' followed by a word chosen from : | 'PSURY' : the P/Y fracture model | | 'ROUSSELIER': Rousselier's ductile fracture model | | 'GURSON2' : Gurson's ductile fracture model modified | by Needleman and Tvergaard | | 'DRAGON' : Dragon's damage model | | 'BETON_URGC' : concrete model | | 'BETON_INSA' : concrete model from INSA (Lyon) | | 'BETON_DYNAR_LMT' : concrete model for fast dynamic analysis * 'VISCOPLASTIQUE' followed by a word chosen from : | 'GUIONNET' : Guionnet's viscoplastic model | | 'CHABOCHE' : Chaboche's unified viscoplastic model | | 'OHNO' : Ohno's viscoplastic model | | 'VISCODOMMAGE' : Lemaitre and Chaboche's damage | viscoplastic model | | 'PARFAIT' : viscoplasticity without hardening | | 'POUDRE_A' : Viscoplastic model of ABOUAF for powder | compaction | | 'DDI' : Two inelastic strain model. | | 'KOCKS' : KOCKS visco plastic model | | 'NOUAILLA_A' : Nouaillas' viscoplastique model | | 'NOUAILLA_B' : Chaboche's viscoplastique model | | 'VISK2' : A viscoplastic model with 2 internal variables | built on plasticity with kinematic hardening | | 'VISCOHINTE' : This model is a time dependant dammaging model. | Dammage will occur til rupture. The finite | element support for this model is the RACCOR | type element JOI2. | | 'MISTRAL' : General Strain Model (MaterIal STRAin Laws). | This model requires an orthotropic elastic behaviour. | | 'GATT_MONERIE' : GATT_MONERIE's viscoplastic model describing | the behaviour of UO2 fuel | | 'UO2': Model coupling the viscoplasticity behaviour of UO2 fuel | (based on the GATT_MONERIE model) with its brittle | behaviour under tensile loading (based on the | fictitious crack model approach proposed by OTTOSEN) | | 'VISCODD': Viscoplastic model with isotropic hardening | and two damage mechanisms : one for smooth isotropic damage | and the second for anisotropic creep damage | | 'SYCO1' : standard 'SYMONDS & COWPER' law (cf. MATE) | 'SYCO2' : modified 'SYMONDS & COWPER' law (cf. MATE) | available for 3D and 2D plane stress and plane strain | formulations and with massive and XFEM elements. * 'VISCO_EXTERNE' : external viscoplastic behavior's model. No default option. Must be followed by a word chosen from : | 'GENERAL' : this model evaluates an inelastic strain as the | combination of a swelling term with a creep term, | using a general criterion on stresses | (Von Mises criterion) | | 'DRUCKER_PRAGER' : Option not available for calculation | | 'COHESION' : Option not available for calculation | | 'CONSOLIDATION' : Option not available for calculation | | 'UTILISATEUR' : Option not available for calculation | | | Contents : | -------- | A number ILOI1 or a CLOI16 name has to be attributed to any | behavior's law of thhe 'VISCO_EXTERNE' group. | The list of the external parameters of the law is given by the | LMOTS1 object, if such parameters exist. | The material components are the same as these of the | 'ELASTIQUE' 'ISOTROPE' behavior. | Four internal variables are pre-defined for any law of the | 'VISCO_EXTERNE' group : | - 'EC0 ' : uniaxial equivalent creep strain, | - 'ESW0' : uniaxial equivalent swelling strain, | - 'P ' : equivalent pressure stress = 1/3 trace (SIGMA), | - 'QTLD' : equivalent deviatoric stress (Von Mises' criterion) | = SQRT(3/2 S:S), | where S is the deviatoric stress tensor. | The list of the additional internal variables of the law is | given by the LMOTS3 object, if such variables exist. | | The 'VISCO_EXTERNE' laws are integrated by a general second | order Runge-Kutta method. The integration scheme calls the | user-defined subroutine CREEP for the computation of creep | and swelling equivalent strain rates. | With the following given informations : | - the uniaxial equivalent creep and swelling strains EC0 and | ESW0, at the start of the time increment ; | - the equivalent pressure and deviatoric stresses P and QTLD, | at the start OR at the end of the time increment ; | - if applicable, the additional internal variables at the | start of the time increment ; | - a time increment dt ; | - a start/end of time increment flag ; | the CREEP subroutine computes : | - two increments of equivalent creep and swelling strains, | defined as the products of the equivalent strain rates | (calculated at the start OR at the end of the time increment) | multiplied by the time increment ; | - the values of the additional internal variables at the end | of time increment, if the routine has been called with the | flag "end of time increment". | | Operating rules : | --------------- | 1. Code the CREEP subroutine and its dependancies : | - Take the example CREEP subroutine integrated in CAST3M. | It is necessary to keep the first card of the example | subroutine, in order to compile successfully. | - Keep the calling interface of the routine. | 2. Compile the external subroutine CREEP and its dependancies, | then edit links with the rest of the code. | 3. Use the MODE and COMP operators, the PASAPAS procedure | with respect to their syntax in Gibiane language. | |---------------------------------------------------------------- | Interface of the external subroutine CREEP : | ------------------------------------------ | | SUBROUTINE CREEP (DECRA,DESWA,STATEV,serd,EC0,ESW0,P,QTILD, | & TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME, | & CMNAME,leximp,LEND,COORDS,NSTATV,NOEL,NPT, | & layer,kspt,KSTEPC,KINC) | IMPLICIT REAL*8(A-H,O-Z) | | | CHARACTER*16 CMNAME | | INTEGER leximp, LEND, NSTATV, NOEL, NPT, layer, kspt, | & KSTEPC, KINC | | REAL*8 DECRA(5), DESWA(5), STATEV(*), | & serd, EC0, ESW0, P, QTILD, | & TEMP,DTEMP, PREDEF(*),DPRED(*), TIME(2),DTIME, | & COORDS(*) | |---------------------------------------------------------------- | OUT : DECRA : REAL*8(5) | DECRA(1) : scalar increment of equivalent | creep strain, calculated at the start of | time increment if LEND=0, at the end of | time increment if LEND=1. | Components 2 to 5 are not used when linking | to CAST3M. | | OUT : DESWA : REAL*8(5) | DESWA(1) : scalar increment of equivalent | swelling strain, calculated at the start of | time increment if LEND=0, at the end of time | increment if LEND=1. | Components 2 to 5 are not used when linking | to CAST3M. | | Note : the CREEP subroutine evaluates the strain RATES (creep | and swelling strain rates) at the start or at the end of the | time increment, depending on the value of the flag LEND. | The strain increments DECRA(1) and DESWA(1) are defined as the | products of the strain rates by the time increment DTIME. | | IN/OUT : STATEV : REAL*8(*), additional internal variables. | In addition to the four pre-defined internal | variables 'EC0 ', 'ESW0', 'P ' and 'QTLD'. | Values at the start of time increment. | Might be updated by the values at the end of | time increment if LEND=1. | | IN/OUT : serd : REAL*8, magnitude of the strain energy rate | density (contracted product of stress tensor | by the inelastic strain rates tensor). | Optional input/output, not used further by | CAST3M | | IN : EC0 : REAL*8, uniaxial equivalent creep strain. | 1st of the 4 pre-defined internal variables. | Value at the start of time increment. | | IN : ESW0 : REAL*8, uniaxial equivalent swelling strain. | 2nd of the 4 pre-defined internal variables. | Value at the start of time increment. | | IN : P : REAL*8, equivalent pressure stress | = 1/3 trace(SIGMA) | 3rd of the 4 pre-defined internal variables. | Value at the start of time increment if | LEND=0, at the end if LEND=1. | | IN : QTILD : REAL*8, equivalent deviatoric stress | (Von Mises' criterion) = SQRT(3/2 S:S), | where S is deviatoric stress tensor. | 4th of the 4 pre-defined internal variables. | Value at the start of time increment if | LEND=0, at the end if LEND=1. | | IN : TEMP : REAL*8, temperature at the end | of time increment. | DTEMP : REAL*8, increment of temperature during | the time increment. | PREDEF : REAL*8(*), external parameters of the | behavior's law, values at the end of | time increment. | DPRED : REAL*8(*), increments of external parameters | during the time increment. | | IN : TIME : REAL*8(2) | TIME(1) : value of step time at the end of | the increment. The step time is defined as | the cumulation of the increments during the | integration. | TIME(2) : value of total time at the end of | the increment. | DTIME : REAL*8, time increment. | | Note : the total time at the last point of convergence is | TIME(2)-DTIME. | | IN : CMNAME : CHARACTER*16, identifier of the behavior law | 1) The law is identified by the number | which has been attributed to it : ILOI1, | behind the keyword 'NUME_LOI'. | The number is encoded in the last four | characters of the string CMNAME, and has to | be recovered in CREEP by an instruction such | as K4ILOI = CMNAME(13:16) | where K4ILOI is a local variable of type | CHARACTER*4 | 2) The law is identified by the name which has | been attributed to it : CLOI16, behind the | keyword 'NOM_LOI'. | The name can be directly recovered in CREEP | because CMNAME = CLOI16 (16 characters) | | IN : leximp : INTEGER | Inactive input | | IN : LEND : INTEGER, start/end of time increment flag. | LEND=0 : the inputs P and QTILD are defined | at the start of time increment ; the outputs | DECRA(1) and DESWA(1) are calculated at the | start of time increment. | LEND=1 : the inputs P and QTILD are defined | at the end of time increment ; the outputs | DECRA(1) and DESWA(1) are calculated at the | end of time increment. | | IN : COORDS : REAL*8(*), coordinates of the current | integration point. | | IN NSTATV : INTEGER, number of additional internal | variables (in addition to the 4 pre-defined) | | IN : NOEL : INTEGER, number of the current element | Note : it is a local number (number of the | element in the current region of the model). | NPT : INTEGER, number of the current integration | point. | | IN : layer : INTEGER | kspt : INTEGER | Inactive input | | IN/OUT : KSTEPC : INTEGER | This input has no sense when linking | to CAST3M | | Programming rule : the code of the user may use KSTEPC as a | return flag of the CREEP subroutine, with respect to the two | rules below : | 1. No re-initialization of KSTEPC when entering into CREEP. | KSTEPC has been initialized to 1 before calling CREEP. | 2. On error, KSTEPC is affected by a value not equal to 1 | | IN : KINC : INTEGER, counter of time increments. | This counter is updated by the calling | subroutine. | |---------------------------------------------------------------- Notes : * To date, it is not possible to use the non-linear ________ behaviors simultaneously * To date, only the isotropic and elastic linear behaviors may be used simultaneously with a non-linear behavior, except for the MISTRAL model which requires an orthotropic elastic behaviour. -------------------------------------------- | Material name for the LIQUIDE formulation| -------------------------------------------- Only the formulation data is required. ------------------------------------------------------- | Material name for the LIQUIDE MECANIQUE formulation | ------------------------------------------------------- Only the formulation data is required. -------------------------------------------- | Material name for the POREUX formulation | -------------------------------------------- First of all, it is required that the linear behavior be defined. Then, if the case arises, the non-linear behavior will be defined according to the following instructions : - Linear behaviors : ------------------ * 'ELASTIQUE' ('ISOTROPE') - Non-linear behaviors : ---------------------- * 'PLASTIQUE' followed by a word chosen from : | ('ISOTROPE') : Von Mises' model with isotropic | stress hardening | | 'PARFAIT' : Von Mises' model without stress | hardening | | 'CINEMATIQUE' : Von Mises' model with linear | kinematic stress hardening | | 'DRUCKER_PRAGER' : Drucker-Prager's model with non- | associated yield | | 'DRUCKER_PARFAIT' : Drucker-Prager's simplified model | (ideal) --------------------------------------------- | Material name for the CONTACT formulation | --------------------------------------------- The given mesh must be made by the 'IMPO' operator. - No friction : ------------ * ('UNILATERAL') (no other possibility for the time being) - Friction : --------- * 'FROTTANT' followed by one of the following : | 'COULOMB' : friction with coulomb law between two lines (2D) | or two surfaces (3D). | | 'FROCABLE' : friction of tendons in concrete with french | regulatory laws from BPEL99. This word must be | followed by the MODEL object of tendons and by | the mesh of the volume inside which tendons are | embedded. ------------------------------- | Material name in THERMIQUE | ------------------------------- | ('ISOTROPE') | 'ORTHOTROPE' | 'ANISOTROPE' * Behaviours | 'CONDUCTION' It is the default behaviour | | 'CONVECTION' Convection between two parts of the structure or | between the structure and the environnement. | For shell only one of the words : | 'SUPERIEURE' | 'INFERIEURE' if the exchange is through upper or | under layer of the shell | | 'RAYONNEMENT' Exchange is done by radiation | | - 'INFINI' radiation going to infinity | - 'FAC_A_FAC' Radiation between very close surfaces | GEO1 surface 1 | GEO2 surface 2 | GEO3 mesh of seg2 linking GEO1 and GEO | MOD1 "thermique convection" model made | with racc element joining the | two surfaces. | - 'CAVITE' radiation of the cavity with itself. | 'CONVEXE' if the cavity is convexe | 'FERME' if it is a closed cavity. | 'SYMETRIQUE' if only a symetric part of | the cavity is represented. | 'FAC_FORME' if the user wants to use the | method with computation of | the forms factors (more | expensive) | 'ABSO' if the inside of the cavity is | absorbing | | Remark : a thermique rayonnement model can be defined in several | parts if they share the same CONSTITUANT | | | 'PHASE' In case of phase changing. This is a non linear part | of the conduction behaviour | 'ADVECTION' Represents an exchange with a fluid moving inside a | pipe. The temperature of the fluid is constant in | a cross section of the pipe.This is a non linear part | of the conduction behaviour --------------------------- | Material name in DARCY | --------------------------- | ('ISOTROPE') | 'ORTHOTROPE' | 'ANISOTROPE' ---------------------------------------------------- | Material name for the MAGNETODYNAMIQUE formulation | ---------------------------------------------------- The one presently available formulation is the thin shell formulation 'POTENTIEL_VECTEUR' (default option). - syntax : -------- * 'POTENTIEL_VECTEUR' | ('ISOTROPE') (default option) | 'ORTHOTROPE' -------------------------------------------------- | Material name for the NAVIER_STOKES formulation | -------------------------------------------------- Only the formulation data is required -------------------------------------------------- | Material name for the EULER formulation | -------------------------------------------------- Only the formulation data is required ------------------------------------------------------- | Material name for the MELANGE formulation | ------------------------------------------------------- 'CEREM' : phase transformation modelling for 16MND5 steel. 'PARALLELE' : allows a parametric addition of elementary modelling provided with PHASE names. 'ZTMAX' : a modelling of state change with a transition. In the default case, the state variable is temperature T. ------------------------------------------------------- | Material name for the FISSURE formulation | ------------------------------------------------------- | ('PARF') | ('MASS') | ('POISEU_BLASIUS') | ('REEL') | ('FILM') | ('POISEU_COLEBROOK') | ('FROTTEMENT1') | ('FROTTEMENT2') | ('FROTTEMENT3') | ('FROTTEMENT4') * choice of a law of steam | 'PARF' : perfect gas | 'REEL' : real gas * choice of a condensation model | 'MASS' : bulk condensation | 'FILM' : film condensation * choice of a law for the friction coefficient (the various laws are describ in the information on FISS operator) | ('POISEU_BLASIUS') : flow in a smooth rectangular duct | ('POISEU_COLEBROOK') : flow in a rough rectangular duct | ('FROTTEMENTi') : laws with coefficients chosen by the user * default values are : PARF MASS POISEU_BLASIUS ------------------------------------------------------- | Material name for the THERMOHYDRIQUE formulation | ------------------------------------------------------- 'SCHREFLER' : non iso-thermal air-flow and water-flow in a deforming porous medium model. -------------------------------------------------- | Material name for the LIAISON formulation | -------------------------------------------------- The A basis is a set of displacements fields (eigen-modes or static solutions for instance) ; the B basis is the standard finite elements one. - Liaisons available with PASAPAS and DYNE (please see the DYNE operator notice) 'POINT_PLAN' 'FLUIDE' : writes either in A basis or B basis 'POINT_PLAN' 'FROTTEMENT' : B basis 'POINT_PLAN' : A basis or B basis 'POINT_POINT' 'FROTTEMENT' : B basis 'POINT_POINT' 'DEPLACEMENT_PLASTIQUE' : B basis 'POINT_POINT' 'ROTATION_PLASTIQUE' : B basis 'POINT_POINT' : B basis 'POINT_CERCLE' 'MOBILE' : B basis 'POINT_CERCLE' 'FROTTEMENT' : B basis 'POINT_CERCLE' : B basis 'CERCLE_PLAN' 'FROTTEMENT' : B basis 'CERCLE_CERCLE' 'FROTTEMENT' : B basis 'PROFIL_PROFIL' 'INTERIEUR' : B basis 'PROFIL_PROFIL' 'EXTERIEUR' : B basis 'LIGNE_LIGNE' 'FROTTEMENT' : B basis 'LIGNE_CERCLE' 'FROTTEMENT' : B basis 'PALIER_FLUIDE' 'RHODE_LI' : B basis 'COUPLAGE' 'DEPLACEMENT' : A basis 'COUPLAGE' 'VITESSE' : A basis 'POLYNOMIALE' : A basis - Liaisons available only with PASAPAS 'NEWMARK' 'MODAL' : A basis ------------------------------------------------ | Material name for the DIFFUSION formulation | ------------------------------------------------ 'FICK' : Fick's law (J = -D.gradC), model by default. Note : _______ The names of formulation and material behavior will be given in words. The following tables indicate, for each formulation, which finite elements associated with a given geometrical support are available, as well as the degree(s) of their interpolation functions, the options for the calculations in which they may be used (see OPTI), and the corresponding nodal unknowns. -------------------------------------------------------------------- | Finite elements in MECANIQUE formulation | -------------------------------------------------------------------- |Geometrical | Finite | Description|Degree| Option for| Nodal | | support | element | | |calculation| unknowns | -------------------------------------------------------------------- | POI1 | POI1 | point | | PLAN GENE | UX UY | | | | | |see Remark | UZ RX RY | -------------------------------------------------------------------- | POI1 | CERC | circle | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT | -------------------------------------------------------------------- | SEG2 | BARR | bar | 1 | PLAN CONT | UX UY | | | | | | PLAN DEFO | UX UY | | | | tendon | | TRID | UX UY UZ | | |------------------------------------------------------ | | BAEX | bar | 1 | PLAN CONT | UX UY | | | | excentred | | PLAN DEFO | UX UY | | | | | | TRID | UX UY UZ | | |-----------------------------------------------------| | | (M1D2) | solid (1D) | 1 | UNID PLAN | UX | | | | 2 nodes | | UNID AXIS | UR | | | | | | UNID SPHE | UR | | |------------------------------------------------------ | | COQ2 | thin shell | 1 | PLAN CONT | UX UY RZ | | | | | | PLAN DEFO | UX UY RZ | | | | | and | PLAN GENE | UX UY RZ | | | | | |see Remark | UZ RX RY | | | | | 3 | AXIS | UR UZ RT | | | | | | FOUR | UR UZ UT | | | | | | | RT | | |------------------------------------------------------ | | POUT | beam | 1 | TRID | UX UY UZ | | | | | and | | RX RY RZ | | | | | 3 | PLAN CONT | UX UY RZ | | | | | | PLAN DEFO | UX UY RZ | | |------------------------------------------------------ | | TIMO |Timoschenko's 1 | TRID | UX UY UZ | | | | beam | | | RX RY RZ | | | | | | PLAN CONT | UX UY RZ | | | | | | PLAN DEFO | UX UY RZ | | |------------------------------------------------------ | | TUYA | right pipe | 1 | TRID | UX UY UZ | | | | and | and | | RX RY RZ | | | | elbow | 3 | | | | |------------------------------------------------------ | | TUFI | cracked | | TRID | UX UY UZ | | | | pipe | | | RX RY RZ | -------------------------------------------------------------------- | SEG3 | BAR3 | barre | 2 | TRID | UX UY UZ | | |------------------------------------------------------ | | (M1D3) | solid (1D) | 1 | UNID PLAN | UX | | | | 3 nodes | | UNID AXIS | UR | | | | | | UNID SPHE | UR | -------------------------------------------------------------------- | TRI3 | (TRI3) | solid | 1 | PLAN CONT | UX UY | | | | triangle | | PLAN DEFO | UX UY | | | | with 3 | | PLAN GENE | UX UY | | | | nodes | |see Remark | UZ RX RY | | | | | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT |* | |------------------------------------------------------ | | ICT3 |incompressib.| 1 | PLAN CONT | UX UY | | | | triangle | | PLAN DEFO | UX UY | | | | with 3 | | AXIS | UR UZ | | | | nodes | | | | | |------------------------------------------------------ | | COQ3 | thin shell | 1 | TRID | UX UY UZ | | | | unfinished | and | | RX RY RZ | | | | polynomial | 3 | | | | | | in flexion | | | | | |------------------------------------------------------ | | DKT | thin shell | 1 | TRID | UX UY UZ | | | | Kirchhoff's| and | | RX RY RZ | | | | discrete | 3 | | | | | | hypothesis | | | | | |------------------------------------------------------ | | DST | shell with | 1 | TRID | UX UY UZ | | | | transverse | and | | RX RY RZ | | | | shear | 3 | | | | |------------------------------------------------------ | | TRIH |homogeneized| | PLAN CONT | P PI UX UY| | | | element | | PLAN DEFO | P PI UX UY| | | | | | AXIS | P PI UR RT| | | | | | FOUR | P PI UR UT| | | | | | | RR RT | | |------------------------------------------------------ | | TRIS |beam cross- | 1 | PLAN | | | | | -section | | | | | | | element | | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | solid | 1 | PLAN CONT | UX UY | | | | quadrangle | | PLAN DEFO | UX UY | | | | with 4 | | PLAN GENE | UX UY | | | | nodes | |see Remark | UZ RX RY | | | | | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT | | |------------------------------------------------------ | | Q4RI | solid | 1 | PLAN CONT | UX UY | | | | quadrangle | | PLAN DEFO | UX UY | | | | with 4 | | PLAN GENE | UX UY | | | | nodes and | |see Remark | UZ RX RY | | | | 1x1 Gauss | | AXIS | UR UZ | | | | points | | FOUR | UR UZ UT | | |-----------------------------------------------------| | | ICQ4 |incompressib.| 1 | PLAN CONT | UX UY | | | | quadrangle | | PLAN DEFO | UX UY | | | | with 4 | | AXIS | UR UZ | | | | nodes | | | | | |------------------------------------------------------ | | COQ4 | thin shell | 1 | TRID | UX UY UZ | | | | with | | | RX RY RZ | | | | transverse | | | | | | | shear | | | | | |------------------------------------------------------ | | QUAS |beam cross- | 1 | PLAN | | | | | -section | | | | | | | element | | | | | |------------------------------------------------------ | | QUAH | homogenized | | AXIS | P PI UR RT| | | | element | | FOUR | P PI UR UT| | | | | | | RR RT | -------------------------------------------------------------------- | TRI6 | (TRI6) | solid | 2 | PLAN CONT | UX UY | | | | triangle | | PLAN DEFO | UX UY | | | | with 6 | | PLAN GENE | UX UY | | | | nodes | |see Remark | UZ RX RY | | | | | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT | | |------------------------------------------------------ | | ICT6 |incompressib.| 2 | PLAN CONT | UX UY | | | | triangle | | PLAN DEFO | UX UY | | | | with 6 | | AXIS | UR UZ | | | | nodes | | | RX RY RZ | | |------------------------------------------------------ | | COQ6 | thick | 2 | TRID | UX UY UZ | | | | shell | | | RX RY RZ | | |------------------------------------------------------ | | TRH6 | homogenized | 2 | PLAN CONT | P PI UX UY| | | | element | | PLAN DEFO | P PI UX UY| -------------------------------------------------------------------- | QUA8 | (QUA8) | solid | 2 | PLAN CONT | UX UY | | | | quadrangle| | PLAN DEFO | UX UY | | | | with 8 | | PLAN GENE | UX UY | | | | nodes | |see Remark | UZ RX RY | | | | | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT | | |------------------------------------------------------ | | Q8RI | solid | 2 | PLAN CONT | UX UY | | | | quadrangle | | PLAN DEFO | UX UY | | | | with 8 | | PLAN GENE | UX UY | | | | nodes and | |see Remark | UZ RX RY | | | | 2x2 Gauss | | AXIS | UR UZ | | | | points | | FOUR | UR UZ UT | | |-----------------------------------------------------| | | ICQ8 |incompressib.| 2 | PLAN CONT | UX UY | | | | quadrangle | | PLAN DEFO | UX UY | | | | with 8 | | AXIS | UR UZ | | | | nodes | | | RX RY RZ | | |------------------------------------------------------ | | COQ8 | thick | 2 | TRID | UX UY UZ | | | | shell | | | RX RY RZ | -------------------------------------------------------------------- | POLY | POLY |solid Polygon| 1 | PLAN CONT | UX UY | | | | number of | | PLAN DEFO | UX UY | | | | nodes = | | PLAN GENE | UX UY | | | | number of | |see Remark | UZ RX RY | | | | sides | | AXIS | UR UZ | | | | | | FOUR | UR UZ UT | -------------------------------------------------------------------- | CUB8 | (CUB8) | solid | 1 | TRID | UX UY UZ | | | | cube with | | | | | | | 8 nodes | | | | | |------------------------------------------------------ | | CUBH | homogenised | | TRID | P PI UX RY| | | | element | | | UY RX| | | | | | | | -------------------------------------------------------------------- | | SHB8 | Shell | | TRID | UX UY UZ | | | | element | | | | | | | | | | | -------------------------------------------------------------------- | TET4 | (TET4) | solid | 1 | TRID | UX UY UZ | | | | tetrahedron| | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | PRI6 | (PRI6) | solid | 1 | TRID | UX UY UZ | | | | prism | | | | | | |with 6 nodes| | | | -------------------------------------------------------------------- | PYR5 | (PYR5) | solid | 1 | TRID | UX UY UZ | | | | pyramid | | | | | | |with 5 nodes| | | | -------------------------------------------------------------------- | CU20 | (CU20) | solid | 2 | TRID | UX UY UZ | | | | cube | | | | | | |with 20 nodes | | | -------------------------------------------------------------------- | TE10 | (TE10) | solid | 2 | TRID | UX UY UZ | | | | tetrahedron| | | | | | |with 10 nodes | | | -------------------------------------------------------------------- | PR15 | (PR15) | solid | 2 | TRID | UX UY UZ | | | | prism | | | | | | |with 15 nodes | | | -------------------------------------------------------------------- | PY13 | (PY13) | solid | 2 | TRID | UX UY UZ | | | | pyramid | | | | | | |with 13 nodes | | | -------------------------------------------------------------------- | RAC2 | LISP | linespring | | TRID | UX UY UZ | | | | element | | | RX RY RZ | | |------------------------------------------------------ | | LISM | modified | | TRID | UX UY UZ | | | | linespring | | | RX RY RZ | | | | element | | | | | |------------------------------------------------------ | | JOI2 | joint | 1 | PLAN DEFO | UX UY | | | | element | | PLAN CONT | UX UY | | | |with 4 nodes | | AXIS | UR UZ | -------------------------------------------------------------------- | RAC3 | JOI3 | joint | 2 | PLAN DEFO | UX UY | | | | element | | PLAN CONT | UX UY | | | |with 6 nodes | | AXIS | UR UZ | -------------------------------------------------------------------- | LIA3 | JOT3 | joint | 1 | TRID | UX UY UZ | | | | element | | | | | | |with 6 nodes | | | | -------------------------------------------------------------------- | LIA4 | JOI4 | joint | 1 | TRID | UX UY UZ | | | | element | | | | | | |with 8 nodes | | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in FLUIDE formulation | -------------------------------------------------------------------- |Geometrical | Finite | Description|Degree| Option for| Nodal | | support | element | | |calculation| unknowns | -------------------------------------------------------------------- | SEG2 | LSE2 | pipe | 1 | TRID | P PI | | | | element | | | | | | | pure | | | | | | | acoustics | | | | | |------------------------------------------------------ | | LSU2 | free | 1 | PLAN CONT | P PI UZ | | | | surface | | PLAN DEFO | P PI UZ | | | | element | | AXIS | P PI UZ | | | | | | FOUR | P PI UZ | -------------------------------------------------------------------- | TRI3 | LTR3 | solid | 1 | PLAN CONT | P PI | | | | triangle | | PLAN DEFO | P PI | | | | with 3 | | AXIS | P PI | | | | nodes | | FOUR | P PI | | |------------------------------------------------------ | | LSU3 | free | 1 | TRID | P PI UZ | | | | surface | | | | | | | element | | | | -------------------------------------------------------------------- | QUA4 | LQU4 | solid | 1 | PLAN CONT | P PI | | | | quadrangle | | PLAN DEFO | P PI | | | | with 4 | | AXIS | P PI | | | | nodes | | FOUR | P PI | | |------------------------------------------------------ | | LSU4 | free | 1 | TRID | P PI UZ | | | | surface | | | | | | | element | | | | -------------------------------------------------------------------- | CUB8 | LCU8 | solid | 1 | TRID | P PI | | | | cube | | | | | | |with 8 nodes| | | | -------------------------------------------------------------------- | TET4 | LTE4 | solid | 1 | TRID | P PI | | | |tetrahedron | | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | PRI6 | LPR6 | solid | 1 | TRID | P PI | | | | prism | | | | | | |with 6 nodes| | | | -------------------------------------------------------------------- | PYR5 | LPY5 | solid | 1 | TRID | P PI | | | | pyramid | | | | | | |with 5 nodes| | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in FLUIDE MECANIQUE formulation | -------------------------------------------------------------------- |Geometrical | Finite | Description |Degree|Calculation| Nodal | | support | element | | | option | unknowns | -------------------------------------------------------------------- | SEG2 | LITU | liquid | 1 | TRID | P PI UX UY| | | | pipe | | | UZ | | | | transition | | | | -------------------------------------------------------------------- | RAC2 | (RAC2) | liquid | 1 | PLAN CONT | P PI UX UY| | | | solid | | PLAN DEFO | P PI UX UY| | | | transition | | AXIS | P PI UR UZ| | | | | | FOUR | P PI UR UZ| | |------------------------------------------------------ | | RACO | liquid | 3 | PLAN CONT | P PI UX UY| | | | shell | | | RZ | | | | transition | | PLAN DEFO | P PI UX UY| | | | | | | RZ | | | | | | AXIS | P PI UR UZ| | | | | | | RT | | | | | | FOUR | P PI UR UZ| | | | | | | RT | | LIA3 | (LIA3) | liquid | 1 | TRID | P PI | | | | solid | | | UX UY UZ | | | | transition | | | | | |------------------------------------------------------ | | LICO | liquid | 3 | TRID | P PI | | | | shell | | | UX UY UZ | | | | transition | | | RX RY RZ | -------------------------------------------------------------------- | LIA4 | (LIA4) | liquid | 1 | TRID | P PI | | | | solid | | | UX UY UZ | | | | transition | | | | | |------------------------------------------------------ | | LIC4 | liquid | 3 | TRID | P PI | | | | shell | | | UX UY UZ | | | | transition | | | RX RY RZ | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in POREUX formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree| Option for| Nodal | | support | element | | |calculation| unkowns | -------------------------------------------------------------------- | TRI6 | (TRIP) | solid | 2 | PLAN CONT | UX UY P | | | | triangle | and | PLAN DEFO | UX UY P | | | |with 6 nodes| 1 | AXIS | UR UZ P | -------------------------------------------------------------------- | QUA8 | (QUAP) | solid | 2 | PLAN CONT | UX UY P | | | | quadrangle| and | PLAN DEFO | UX UY P | | | |with 8 nodes| 1 | AXIS | UR UZ P | -------------------------------------------------------------------- | CU20 | (CUBP) | solid | 2 | TRID | UX UY UZ | | | | cube | and | | P | | | |with 20 nodes 1 | | | -------------------------------------------------------------------- | TE10 | (TETP) | solid | 2 | TRID | UX UY UZ | | | | tetrahedron| and | | P | | | |with 10 nodes 1 | | | -------------------------------------------------------------------- | PR15 | (PRIP) | solid | 2 | TRID | UX UY UZ | | | | prism | and | | P | | | |with 15 nodes 1 | | | -------------------------------------------------------------------- | RAP3 | (JOP3) | 8 nodes | 2 | PLAN CONT | UX UY P | | | | joint | and | PLAN DEFO | UX UY P | | | | element | 1 | AXIS | UR UZ P | -------------------------------------------------------------------- | LIP6 | (JOP6) | 15 nodes | 2 | TRID | UX UY UZ | | | | joint | and | | P | | | | element | 1 | | | -------------------------------------------------------------------- | LIP8 | (JOP8) | 20 nodes | 2 | TRID | UX UY UZ | | | | joint | and | | P | | | | element | 1 | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in THERMIQUE CONDUCTION PHASE formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree| Option for| Nodal | | support | element | | |calculation| unknowns | -------------------------------------------------------------------- | SEG2 | BARR | bar | 1 | PLAN | T | | | | | | TRID | T | | |------------------------------------------------------ | | (M1D2) | solid (1D) | 1 | UNID PLAN | T | | | | 2 nodes | | UNID AXIS | T | | | | | | UNID SPHE | T | | |------------------------------------------------------ | | COQ2 | thin shell | 1 | PLAN | T TSUP | | | | parabolic | | | TINF | | | | variation | | AXIS | T TSUP | | | | in | | | TINF | | | | thickness | | | | -------------------------------------------------------------------- | SEG3 | (T1D3) | solid (1D) | 1 | UNID PLAN | T | | | | 3 nodes | | UNID AXIS | T | | | | | | UNID SPHE | T | -------------------------------------------------------------------- | TRI3 | (TRI3) | solid | 1 | PLAN | T | | | | triangle | | AXIS | T | | | |with 3 nodes| | | | | |------------------------------------------------------ | | COQ3 | thin shell | 1 | TRID | T TSUP | | | | parabolic | | | TINF | | | | variation | | | | | | | in | | | | | | | thickness | | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | solid | 1 | PLAN | T | | | | quadrangle | | AXIS | T | | | |with 4 nodes| | | | | |------------------------------------------------------ | | COQ4 | thin shell | 1 | TRID | T TSUP | | | | parabolic | | | TINF | | | | variation | | | | | | | in | | | | | | | thickness | | | | -------------------------------------------------------------------- | TRI6 | (TRI6) | solid | 2 | PLAN | T | | | | triangle | | AXIS | T | | | |with 6 nodes| | | | | |------------------------------------------------------ | | COQ6 | thick | 2 | TRID | T TSUP | | | | shell | | | TINF | | | | parabolic | | | | | | | variation | | | | | | | in | | | | | | | thickness | | | | -------------------------------------------------------------------- | QUA8 | (QUA8) | solid | 2 | PLAN | T | | | | quadrangle | | AXIS | T | | | |with 8 nodes| | | | | |------------------------------------------------------ | | COQ8 | thick | 2 | TRID | T TSUP | | | | shell | | | TINF | | | | parabolic | | | | | | | variation | | | | | | | in | | | | | | | thickness | | | | -------------------------------------------------------------------- | CUB8 | (CUB8) | solid | 1 | TRID | T | | | | cube | | | | | | |with 8 nodes| | | | -------------------------------------------------------------------- | TET4 | (TET4) | solid | 1 | TRID | T | | | |tetrahedron | | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | PRI6 | (PRI6) | solid | 1 | TRID | T | | | | prism | | | | | | |with 6 nodes| | | | -------------------------------------------------------------------- | PYR5 | (PYR5) | solid | 1 | TRID | T | | | | pyramid | | | | | | |with 5 nodes| | | | -------------------------------------------------------------------- | CU20 | (CU20) | solid | 2 | TRID | T | | | | cube | | | | | | |with 20 nodes | | | -------------------------------------------------------------------- | TE10 | (TE10) | solid | 2 | TRID | T | | | |tetrahedron | | | | | | |with 10 nodes | | | -------------------------------------------------------------------- | PR15 | (PR15) | solid | 2 | TRID | T | | | | prism | | | | | | |with 15 nodes | | | -------------------------------------------------------------------- | PY13 | (PY13) | solid | 2 | TRID | T | | | | pyramid | | | | | | |with 13 nodes | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in THERMIQUE CONVECTION formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | POI1 | (POI1) | exchange | 1 | UNID PLAN | T | | | | element 1D | | UNID AXIS | T | | | | for solid | | UNID SPHE | T | -------------------------------------------------------------------- | SEG2 | (SEG2) | exchange | 1 | PLAN | T | | | | element | | AXIS | T | | | | for solid | | | | | |------------------------------------------------------ | | COQ2 | exchange | 1 | AXIS | TSUP | | | | element | | |ou TINF | | | | for shell | | | | | |------------------------------------------------------ | | (SEG2) | exchange | 1 | UNID PLAN | T | | | | element 1D | | UNID AXIS | T | | | | opposite | | UNID SPHE | T | | | | 1x1 nodes | | | | -------------------------------------------------------------------- | SEG3 | (SEG3) | exchange | 2 | PLAN | T | | | | element | | AXIS | T | | | | for solid | | | | -------------------------------------------------------------------- | TRI3 | (TRI3) | exchange | 1 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ3 | exchange | 1 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | exchange | 1 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ4 | exchange | 1 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | TRI6 | (TRI6) | exchange | 2 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ6 | exchange | 2 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | QUA8 | (QUA8) | exchange | 2 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ8 | exchange | 2 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | RAC2 | (RAC2) | exchange | 1 | PLAN | T | | | | element | | AXIS | | | | | opposite | | | | | | | 2x2 nodes | | | | -------------------------------------------------------------------- | RAC3 | (RAC3) | exchange | 2 | PLAN | T | | | | element | | AXIS | | | | | opposite | | | | | | | 2x3 nodes | | | | -------------------------------------------------------------------- | LIA3 | (LIA3) | exchange | 1 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x3 nodes | | | | -------------------------------------------------------------------- | LIA4 | (LIA4) | exchange | 1 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x4 nodes | | | | -------------------------------------------------------------------- | LIA6 | (LIA6) | exchange | 2 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x6 nodes | | | | -------------------------------------------------------------------- | LIA8 | (LIA8) | exchange | 2 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x8 nodes | | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in THERMIQUE RAYONNEMENT formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | SEG2 | (SEG2) | radiation | 1 | PLAN | T | | | | element | | AXIS | T | | | | for solid | | | | | |------------------------------------------------------ | | COQ2 | radiation | 1 | AXIS | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | SEG3 | (SEG3) | radiation | 2 | PLAN | T | | | | element | | AXIS | T | | | | for solid | | | | -------------------------------------------------------------------- | TRI3 | (TRI3) | radiation | 1 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ3 | radiation | 1 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | radiation | 1 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ4 | radiation | 1 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | TRI6 | (TRI6) | radiation | 2 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ6 | radiation | 2 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | QUA8 | (QUA8) | radiation | 2 | TRID | T | | | | element | | | | | | | for solid | | | | | |------------------------------------------------------ | | COQ8 | radiation | 2 | TRID | TSUP | | | | element | | |ou TINF | | | | for shell | | | | -------------------------------------------------------------------- | RAC2 | (RAC2) | radiation | 1 | PLAN | T | | | | element | | AXIS | | | | | opposite | | | | | | | 2x2 nodes | | | | -------------------------------------------------------------------- | RAC3 | (RAC3) | radiation | 2 | PLAN | T | | | | element | | AXIS | | | | | opposite | | | | | | | 2x3 nodes | | | | -------------------------------------------------------------------- | LIA3 | (LIA3) | radiation | 1 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x3 nodes | | | | -------------------------------------------------------------------- | LIA4 | (LIA4) | radiation | 1 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x4 nodes | | | | -------------------------------------------------------------------- | LIA6 | (LIA6) | radiation | 2 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x6 nodes | | | | -------------------------------------------------------------------- | LIA8 | (LIA8) | radiation | 2 | TRID | T | | | | element | | | | | | | opposite | | | | | | | 2x8 nodes | | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in THERMIQUE ADVECTION formuulation | -------------------------------------------------------------------- | Support | Element | Description |Degre| Option | Inconnues | | geometrique| fini | | | de calcul | nodales | -------------------------------------------------------------------- | SEG2 | (SEG2) | element de | 1 | PLAN | T | | | | tuyauterie | | TRID | | | | | | | | | |------------|------------------------------------------------------ | | (SEG3) | element de | 2 | PLAN | | | SEG3 | | tuyauterie | | TRID | | | | | | | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in DARCY formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | TRI7 | (HYT3) | hybrid | | PLAN | TH | | | | triangle | | | | | | |with 3 nodes| | | | -------------------------------------------------------------------- | QUA9 | (HYQ4) | hybrid | | PLAN | TH | | | | quadrangle | | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | CU27 | (HYC8) | hybrid | | TRID | TH | | | | cube | | | | | | |with 8 nodes| | | | -------------------------------------------------------------------- | TE15 | (HYT4) | hybrid | | TRID | TH | | | |tetrahedron | | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | PR21 | (HYP6) | hybrid | | TRID | TH | | | | prism | | | | | | |with 6 nodes| | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements in FROTTEMENT formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description|Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | MULT | (FRO3) | element | | PLAN | UX UY LX | | | | with | | AXIS | UR UZ LX | | | | 3 + 2 nodes| | | | | MULT | (FRO4) | 2+ n noeuds| 1 | 3D | UX UY UZ | | | | | | | LX | | | | | | | | | MULT | | 2 noeuds | 1 | 2D-3D | Ui LX | | | | | | | | -------------------------------------------------------------------- -------------------------------------------------------------------- --------------------------------------------------------------------- | Finite elements for the MAGNETODYNAMIQUE formulation | --------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation | unknowns | --------------------------------------------------------------------- | TRI3 | (ROT3) | thin shell | 1 | TRID | FC | | | | element for | | | | | | | magneto- | | | | | | | dynamic | | | | | | | formulation | | | | --------------------------------------------------------------------- --------------------------------------------------------------------- | Finite elements for the NAVIER_STOKES formulation | --------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | QUAF LINE and MACR are generic | | | | | names for elements family listed | | | | | below in the corresponding column. | | | | -------------------------------------------------------------------- | QUAF | LINE | | 1 U | | | | | | | 0 P | | | | TRI7 QUA9 |TRF3 QUF4| TRI3 QUA4 | | PLAN AXI | UX UY | | CU27 PR21 |CUF8 PRF6| CUB8 PRI6 | | TRID | UX UY UZ | | TE15 PR19 |TEF4 PYF5| TET4 PYR5 | | TRID | UX UY UZ | | | | Pression nc | | | | | | | P0 | | | | --------------------------------------------------------------------- | QUAF | MACR | | 1 U | | | | | | | 0 P | | | | TRI7 QUA9 |MTR6 MQU9|4xTRI3 4xQUA4| | PLAN AXI | UX UY | | CU27 PR21 |MC27 MP18|8xCUB8 8xPRI6| | TRID | UX UY UZ | | TE15 PR19 |MT10 MP14|8xTET4 | | TRID | UX UY UZ | | | | Pression nc | | | | | | | iso P1 | | | | --------------------------------------------------------------------- | QUAF | QUAF | | 2 U | | | | | | | 1 P | | | | TRI7 QUA9 |TRF7 QUF9| TRI7 QUA9 | | PLAN AXI | UX UY | | CU27 PR21 |CF27 PF21| CU27 PR21 | | TRID | UX UY UZ | | TE15 PR19 |TF15 PF19| TE15 PR19 | | TRID | UX UY UZ | | | | Pression nc | | | | | | | P1 | | | | --------------------------------------------------------------------- --------------------------------------------------------------------- | Finite volumes for the EULER formulation | --------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation | unknowns | --------------------------------------------------------------------- | TRI3 QUA4 | | | | PLAN | | | CUB8 PRI6 | | | | TRID | | | TET4 PYR5 | | | | TRID | | --------------------------------------------------------------------- --------------------------------------------------------------------- | Finite elements for the FISSURE formulation | --------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation | unknowns | --------------------------------------------------------------------- | SEG2 | | | | MONOD | | | | | | | PLAN | | | | | | | TRID | | --------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements for the THERMOHYDRIQUE formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation | unknowns | -------------------------------------------------------------------- | TRI3 | (TRI3) | solid | 1 | PLAN | PG,PC,T | | | | triangle | | AXIS | PG,PC,T | | | |with 3 nodes| | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | solid | 1 | PLAN | PG,PC,T | | | | quadrangle | | AXIS | PG,PC,T | | | |with 4 nodes| | | | -------------------------------------------------------------------- | TRI6 | (TRI6) | solid | 2 | PLAN | PG,PC,T | | | | triangle | | AXIS | PG,PC,T | | | |with 6 nodes| | | | -------------------------------------------------------------------- | QUA8 | (QUA8) | solid | 2 | PLAN | PG,PC,T | | | | quadrangle | | AXIS | PG,PC,T | | | |with 8 nodes| | | | -------------------------------------------------------------------- | CUB8 | (CUB8) | solid | 1 | TRID | PG,PC,T | | | | cube | | | | | | |with 8 nodes| | | | -------------------------------------------------------------------- | TET4 | (TET4) | solid | 1 | TRID | PG,PC,T | | | |tetrahedron | | | | | | |with 4 nodes| | | | -------------------------------------------------------------------- | PRI6 | (PRI6) | solid | 1 | TRID | PG,PC,T | | | | prism | | | | | | |with 6 nodes| | | | -------------------------------------------------------------------- | PYR5 | (PYR5) | solid | 1 | TRID | PG,PC,T | | | | pyramid | | | | | | |with 5 nodes| | | | -------------------------------------------------------------------- | CU20 | (CU20) | solid | 2 | TRID | PG,PC,T | | | | cube | | | | | | |with 20 nodes | | | -------------------------------------------------------------------- | TE10 | (TE10) | solid | 2 | TRID | PG,PC,T | | | |tetrahedron | | | | | | |with 10 nodes | | | -------------------------------------------------------------------- | PR15 | (PR15) | solid | 2 | TRID | PG,PC,T | | | | prism | | | | | | |with 15 nodes | | | -------------------------------------------------------------------- | PY13 | (PY13) | solid | 2 | TRID | PG,PC,T | | | | pyramid | | | | | | |with 13 nodes | | | -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite element for the LIAISON formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation| unknowns | ------------------------------------------------------------------- | POI1 | (POI1) | point | 1 | PLAN | ALFA BETA| | | | | | TRID |in A basis| -------------------------------------------------------------------- -------------------------------------------------------------------- | Finite elements for the DIFFUSION formulation | -------------------------------------------------------------------- | Geometrical| Finite | Description |Degree|Option for | Nodal | | support | element | | |calculation| unknowns | -------------------------------------------------------------------- | SEG2 | (M1D2) | solid (1D) | 1 | UNID PLAN | (CO) | | | | at 2 nodes | | UNID AXIS | (CO) | | | | | | UNID SPHE | (CO) | -------------------------------------------------------------------- | SEG3 | (M1D3) | solid (1D) | 1 | UNID PLAN | (CO) | | | | at 3 nodes | | UNID AXIS | (CO) | | | | | | UNID SPHE | (CO) | -------------------------------------------------------------------- | TRI3 | (TRI3) | solid | 1 | PLAN | (CO) | | | | triangle | | AXIS | (CO) | | | | at 3 nodes | | | | -------------------------------------------------------------------- | QUA4 | (QUA4) | solid | 1 | PLAN | (CO) | | | | quadrangle | | AXIS | (CO) | | | | at 4 nodes | | | | -------------------------------------------------------------------- | TRI6 | (TRI6) | solid | 2 | PLAN | (CO) | | | | triangle | | AXIS | (CO) | | | | at 6 nodes | | | | -------------------------------------------------------------------- | QUA8 | (QUA8) | solid | 2 | PLAN | (CO) | | | | quadrangle | | AXIS | (CO) | | | | at 8 nodes | | | | -------------------------------------------------------------------- | CUB8 | (CUB8) | solid cube | 1 | TRID | (CO) | | | | at 8 nodes | | | | -------------------------------------------------------------------- | TET4 | (TET4) | solid | 1 | TRID | (CO) | | | | tetrahedron | | | | | | | at 4 nodes | | | | -------------------------------------------------------------------- | PRI6 | (PRI6) | solid prism | 1 | TRID | (CO) | | | | at 6 nodes | | | | -------------------------------------------------------------------- | PYR5 | (PYR5) | solid | 1 | TRID | (CO) | | | | pyramid | | | | | | | at 5 nodes | | | | -------------------------------------------------------------------- | CU20 | (CU20) | solid cube | 2 | TRID | (CO) | | | | at 20 nodes | | | | -------------------------------------------------------------------- | TE10 | (TE10) | solid | 2 | TRID | (CO) | | | | tetrahedron | | | | | | | at 10 nodes | | | | -------------------------------------------------------------------- | PR15 | (PR15) | solid prism | 2 | TRID | (CO) | | | | at 15 nodes | | | | -------------------------------------------------------------------- | PY13 | (PY13) | solid | 2 | TRID | (CO) | | | | pyramid | | | | | | | at 13 nodes | | | | -------------------------------------------------------------------- Remark: Correspondence between the nodal (primal) unknowns (P) and ------- their duals (D): P : UX UY UZ UT RX RY RZ RT RR P PI T TSUP TINF LX TH FC PG PC D : FX FY FZ FT MX MY MZ MT MR FP FPI Q QSUP QINF FLX FLUX ED QG QC P : ALFA BETA CO D : FALF FBET QCO The nodal unknowns of the generalized plane strains (UZ,RX,RY) and their duals (FZ,MX,MY) are supported by the point that has been defined for the MODE PLAN GENE option. For the 1D-calculus modes, the nodal unknows linked to the generalized plane strain(s) (UZ,UY) and their duals (FZ,FY) are supported by the point definid for the MODE UNID PLAN CYGZ/DYGZ/GYCZ/GYDZ/GYGZ option or the MODE UNID AXIS AXGZ option. Note : The behaviour ELASTIQUE UNIDIRECTIONNEL does not work ----- in a tridimensionnal solid case. $$$$
© Cast3M 2003 - All rights reserved.
Disclaimer