User Tools

Site Tools


prepro

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
prepro [2019/06/13 10:00]
helene [16.2 Description of the foundation segments]
prepro [2023/11/23 17:38] (current)
gilles [3. Nodes coordinates]
Line 37: Line 37:
 |:::|= 6 forging simulation, different from Dyduch| |:::|= 6 forging simulation, different from Dyduch|
 |:::|= 7 proportional fatigue computation (see fatigue calculation in Hill_3D_KI law), you must generate the file 61 (.ntfdam) [[prepro#​(0) NTFDAM file structure|(0)]]| |:::|= 7 proportional fatigue computation (see fatigue calculation in Hill_3D_KI law), you must generate the file 61 (.ntfdam) [[prepro#​(0) NTFDAM file structure|(0)]]|
 +|IFEM2<​sup>​3</​sup>​|Type of multi-scale problem:|
 +|:::| =0 - Only macro-scale problem|
 +|:::| =1 - Micro hydro-mechanical law with interfaces ([[laws:​fe2wg|FE2WG]])|
 +|:::| =2 - Micro hydraulic law ([[laws:​hmic|HMIC]])|
 +|:::| =3 - Micro hydraulic + pollutants law ([[laws:​hypofe2|HYPOFE2]])|
 ^2.2 (15I5)^^ ^2.2 (15I5)^^
 |IRENU<​sup>​5</​sup>​| Index for the DOF renumbering \\ = 0 nil \\ = 1 reading of the RENUM card and renumbering of the DOF| |IRENU<​sup>​5</​sup>​| Index for the DOF renumbering \\ = 0 nil \\ = 1 reading of the RENUM card and renumbering of the DOF|
Line 45: Line 50:
 |MSIGP<​sup>​30</​sup>​| Maximum number of stresses per point. [[prepro#​(3)|(3)]] \\ Default value : 18| |MSIGP<​sup>​30</​sup>​| Maximum number of stresses per point. [[prepro#​(3)|(3)]] \\ Default value : 18|
 |MPARI<​sup>​35</​sup>​| Maximum number of integer parameters for any constitutive law. [[prepro#​(3)|(3)]] \\ Default value : 20 (Be careful, often not sufficient for thermo-mechanical-metallurgical coupling)| |MPARI<​sup>​35</​sup>​| Maximum number of integer parameters for any constitutive law. [[prepro#​(3)|(3)]] \\ Default value : 20 (Be careful, often not sufficient for thermo-mechanical-metallurgical coupling)|
-|IMDIS<​sup>​40</​sup>​| = 1 displacements can be printed (the CONEC dimension increases) (compulsory for BEM2D) \\ = 0 displacements are not printed NCOL = 1 for static analysis, NCOL = 3 for dynamic analysis. \\ = 2 NCOL = 2 for static analysis useful for OCASFO routine (see data appendix 11) \\ ≥10 presence of a dead load, the NCOL value is incremented by 1. (see Wang)| +|IMDIS<​sup>​40</​sup>​| = 1 displacements can be printed (the CONEC dimension increases) (compulsory for BEM2D) \\ = 0 displacements are not printed NCOL = 1 for static analysis, NCOL = 3 for dynamic analysis. \\ = 2 NCOL = 2 for static analysis useful for OCASFO routine (see data [[appendices:​a11|appendix 11]]) \\ ≥10 presence of a dead load, the NCOL value is incremented by 1. (see Wang)| 
-|IPIFO<​sup>​45</​sup>​| = 0 nil \\ = 1 force pilotage by the routine OCASFO (see appendix)|+|IPIFO<​sup>​45</​sup>​| = 0 nil \\ = 1 force pilotage by the routine OCASFO (see [[appendices:​a11|appendix ​11]])|
 |ICEVO<​sup>​50</​sup>​| = 0 nil \\ = 1 checking of the total volume of the piece computed by the Jacobian matrix determinant (ready for BLZ3D)| |ICEVO<​sup>​50</​sup>​| = 0 nil \\ = 1 checking of the total volume of the piece computed by the Jacobian matrix determinant (ready for BLZ3D)|
 |ITREE<​sup>​55</​sup>​| = 0 nil \\ = 1 search of the 3D contact by arborescence;​ \\ ITREE = number of foundations using the sub-tree method [[prepro#​(4)|(4)]]| |ITREE<​sup>​55</​sup>​| = 0 nil \\ = 1 search of the 3D contact by arborescence;​ \\ ITREE = number of foundations using the sub-tree method [[prepro#​(4)|(4)]]|
Line 100: Line 105:
 \\ \\
 A negative value of NTANA is used for dynamic analysis. \\ A negative value of NTANA is used for dynamic analysis. \\
-|NTANA |(2 translations + water pressure + air pressure + temperature)| +^NTANA^(2 translations + water pressure + air pressure + temperature)^ 
-|:::| = 1 Plane mechanical analysis \\ NDOFN = NSPAC = 3 : X,Y,T \\ (no equation for DOF 3)| + = 1  ​| ​Plane mechanical analysis \\ NDOFN = NSPAC = 3 : X,Y,T \\ (no equation for DOF 3)| 
-|:::|= 2 3D mechanical analysis \\ NDOFN = NSPAC = 4 : X, Y, Z, T \\ (no equation for DOF 4)| + = 2  |3D mechanical analysis \\ NDOFN = NSPAC = 4 : X, Y, Z, T \\ (no equation for DOF 4)| 
-|:::|= 3 Plane thermal analysis \\ NDOFN = NSPAC = 3 : X, Y, T \\ (no equation for DOF 1 and 2)| + = 3  |Plane thermal analysis \\ NDOFN = NSPAC = 3 : X, Y, T \\ (no equation for DOF 1 and 2)| 
-|:::|= 4 Plane thermomech. analysis \\ NDOFN = NSPAC = 3 : X, Y, T| + = 4  |Plane thermomech. analysis \\ NDOFN = NSPAC = 3 : X, Y, T| 
-|:::|= 5 2D coupled mechanical water – thermal – air flow analysis \\ NDOFN = NSPAC = 5 : X, Y, PW, Pa, T| + = 5  |2D coupled mechanical water – thermal – air flow analysis \\ NDOFN = NSPAC = 5 : X, Y, PW, Pa, T| 
-|:::|= 6 3D thermal. analysis \\ NDOFN = NSPAC = 4 : X,Y,Z,T| + = 6  |3D thermal. analysis \\ NDOFN = NSPAC = 4 : X,Y,Z,T| 
-|:::|= 7 Plane mechanical analysis with shells \\ NDOFN = NSPAC = 3 \\ (2 translations,​ 1 rotation) X,Y,ZZ | + = 7  |Plane mechanical analysis with shells \\ NDOFN = NSPAC = 3 \\ (2 translations,​ 1 rotation) X,Y,ZZ | 
-|:::|= 8 3D mechanical analysis with shells \\ NDOFN = NSPAC = 6 : X, Y, Z, ψ1, ψ2, ψ3 \\ (3 translations + 3 coordinates of rotation vector)| + = 8  |3D mechanical analysis with shells \\ NDOFN = NSPAC = 6 : X, Y, Z, ψ1, ψ2, ψ3 \\ (3 translations + 3 coordinates of rotation vector)| 
-|:::|= 9 3D thermomechanical. analysis \\ NDOFN = NSPAC = 4 : X, Y, Z, T| + = 9  |3D thermomechanical. analysis \\ NDOFN = NSPAC = 4 : X, Y, Z, T| 
-|:::|= 10 Plane mechanical analysis, second gradient method (Ouafa ElHammoumi) \\ NDOFN = NSPAC = 5 : X, Y, λ, $\frac{d\lambda}{dX}$,​$\frac{d\lambda}{dY}$| + = 10  |Plane mechanical analysis, second gradient method (Ouafa ElHammoumi) \\ NDOFN = NSPAC = 5 : X, Y, λ, $\frac{d\lambda}{dX}$,​$\frac{d\lambda}{dY}$| 
-|:::|= 13 3D coupled mechanical water – thermal – air flow analysis \\ NDOFN = NSPAC = 6 : X, Y, Z, PW, Pa, T| + = 13  |3D coupled mechanical water – thermal – air flow analysis \\ NDOFN = NSPAC = 6 : X, Y, Z, PW, Pa, T| 
-|:::|= 14 Plane mechanical analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 6 : X, Y, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| + = 14  |Plane mechanical analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 6 : X, Y, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| 
-|:::|= 15 2D coupled mechanical water – thermal – air flow – chemical analysis \\ NDOFN = NSPAC = 6 : X, Y, PW, Pa, T, c| + = 15  |2D coupled mechanical water – thermal – air flow – chemical analysis \\ NDOFN = NSPAC = 6 : X, Y, PW, Pa, T, c| 
-|:::|= 16 Plane coupled mechanical water flow analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 7 : X, Y, PW, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| + = 16  |Plane coupled mechanical water flow analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 7 : X, Y, PW, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| 
-|:::|= 17 Plane coupled mechanical water – thermal – air flow analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 9 : X, Y, PW, Pa, T, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| + = 17  |Plane coupled mechanical water – thermal – air flow analysis, second gradient method (Grenoble) \\ NDOFN = NSPAC = 9 : X, Y, PW, Pa, T, V11, V12, V21, V22 (ou λ11, λ12, λ21, λ22 for central node)| 
-|:::|= 18 2D coupled mechanical water – thermal – air flow – chemical analysis \\ NDOFN = NSPAC = 8 : X, Y, PW, Pa, T, c1, c2, c3|+ = 18  |2D coupled mechanical water – thermal – air flow – chemical analysis \\ NDOFN = NSPAC = 8 : X, Y, PW, Pa, T, c1, c2, c3
 +|  = 19  |21 Dofs (3 disp + 18 GND)|
  
 ==== (2) Remark on IANA=5 (generalized plane state) ==== ==== (2) Remark on IANA=5 (generalized plane state) ====
Line 142: Line 148:
 |INODE|Node number| |INODE|Node number|
 |XYZ(INODE,​J)J=1:​NSPAC|Node coordinates| |XYZ(INODE,​J)J=1:​NSPAC|Node coordinates|
 +^3.3 Micro-nodes definition (I5, 2G10.0) if IFEM2 $\neq$ 0^^
 +|TITLE|"​Microstructure"​ written in the first 25 columns|
 +|INODE2|Micro-node number|
 +|XY(INODE2,​J)J=1:​2|Micro-node coordinates|
 __Remark__ : convention \\ __Remark__ : convention \\
 The plane or axisymmetric state axis system is direct (right-handed). Conventionally the Y axis is directed toward the top. In the axisymmetric state, the Y axis is the same as the symmetry axis. All the nodes, elements, and local axis lists definitions must be given in the same axis system, implicit in the whole of this handbook. The plane or axisymmetric state axis system is direct (right-handed). Conventionally the Y axis is directed toward the top. In the axisymmetric state, the Y axis is the same as the symmetry axis. All the nodes, elements, and local axis lists definitions must be given in the same axis system, implicit in the whole of this handbook.
Line 192: Line 202:
 |TITLE|"​DUPLI"​ written in the first 5 columns| |TITLE|"​DUPLI"​ written in the first 5 columns|
 ^6.2 (2I5) or (2I10) if I99999>​1^^ ^6.2 (2I5) or (2I10) if I99999>​1^^
-|IDUPL(I,​1),​ IDUPL(I,2)|I=1, NDUPL \\ Nodes list with IDUPL(I,​1)<​IDUPL(I,​2),​ \\ and IDUPL(I,2) ≠ IDUPL(II,1) \\ and IDUPL(I,​2) ​ ≠ IDUPL(II,​2),​ II=1, I-1|+|IDUPL(I,​1),​ IDUPL(I,​2) ​\\ I=1, NDUPLNodes list with IDUPL(I,​1)<​IDUPL(I,​2),​ \\ and IDUPL(I,2) ≠ IDUPL(II,1) \\ and IDUPL(I,​2) ​ ≠ IDUPL(II,​2),​ II=1, I-1|
  
 ---- ----
Line 221: Line 231:
 ===== 9. Setting of the DOF to a non-zero value ===== ===== 9. Setting of the DOF to a non-zero value =====
  
-This section only exists if NDISP ≠ 0. +-This section only exists if NDISP ≠ 0.  
-After a sub-title card, the program reads the cards 9.2 until it founds NDISP imposed displacements to non zero values.+  
 +-After a sub-title card, the program reads the cards 9.2 until it founds NDISP imposed displacements to non zero values.
 ^9.1 Sub-title (A5)^^ ^9.1 Sub-title (A5)^^
 |TITLE|"​DISPL"​ written in the first 5 columns| |TITLE|"​DISPL"​ written in the first 5 columns|
-^9.2 Imposed displacements (I5, 6G10.0) or (I10, 6G10.0) if I99999=1^^ +^9.2 Imposed displacements (I5, 6G10.0) or (I10, 6G10.0) if I99999=1^^   
-|INODE|Node number|+|INODE|Node number|  
 |DISPL(I) \\ I=1,NDOFN| Imposed displacements at these node's DOF. \\ An imposed to zero displacement is ignored.| |DISPL(I) \\ I=1,NDOFN| Imposed displacements at these node's DOF. \\ An imposed to zero displacement is ignored.|
  
Line 311: Line 322:
  
 ^16.1 Sub-title (A5)^^ ^16.1 Sub-title (A5)^^
-|TITLE| ​"FOUND written in the first five columns|+|TITLE| FOUND written in the first five columns|
 ==== 16.2 Description of the foundation segments ==== ==== 16.2 Description of the foundation segments ====
 This section must be repeated NFOUN times. This section must be repeated NFOUN times.
Line 351: Line 362:
 === Remarks === === Remarks ===
  
-  * For ICODE = 6, the first node represents the circle center, the second node has the (R, φ) coordinates where φ represents the circle rotation (initially φ = 0). \\  +  * __For ICODE = 6__, the first node represents the circle center, the second node has the (R, φ) coordinates where φ represents the circle rotation (initially φ = 0). \\  
-  * For ICODE = 7, the first two nodes represent the cylinder axis limits. The last one has the (R, φ, 0) coordinates where φ represent the cylinder rotation (initially φ = 0). \\ +  * __For ICODE = 7__, the first two nodes represent the cylinder axis limits. The last one has the (R, φ, 0) coordinates where φ represent the cylinder rotation (initially φ = 0). \\ 
-  * For ICODE = 8, the first two nodes represent the cone frustrum axis limits. The last one has the (R<​sub>​1</​sub>,​ R<​sub>​2</​sub>,​ φ) coordinates where R<​sub>​1</​sub>​ and R<​sub>​2</​sub>​ are the radius at nodes 1 and 2, and φ represents the cone frustum rotation (initially φ = 0). +  * __For ICODE = 8__, the first two nodes represent the cone frustrum axis limits. The last one has the (R<​sub>​1</​sub>,​ R<​sub>​2</​sub>,​ φ) coordinates where R<​sub>​1</​sub>​ and R<​sub>​2</​sub>​ are the radius at nodes 1 and 2, and φ represents the cone frustum rotation (initially φ = 0). 
-  * For ICODE = ±9, the first two nodes represent the cone frustum axis limits (N1 and N2). The third node represents the middle of the arc in the plane perpendicular to the cone frustum axis, coming from N1. The fourth node has the (R1, R2, θ) coordinates where R1 and R2 are the radius at the nodes 1 and 2, and θ is the angle (in radian) of the cone frustum (see <imgref fig3>) – only one case has been tested presently: the cone frustum axis is parallel to the z axis.+  * __For ICODE = ±9__, the first two nodes represent the cone frustum axis limits (N1 and N2). The third node represents the middle of the arc in the plane perpendicular to the cone frustum axis, coming from N1. The fourth node has the (R1, R2, θ) coordinates where R1 and R2 are the radius at the nodes 1 and 2, and θ is the angle (in radian) of the cone frustum (see <imgref fig3>) – only one case has been tested presently: the cone frustum axis is parallel to the z axis.
 <​imgcaption fig3|>{{ :​prepro_fig3.png?​400 |}}</​imgcaption>​ <​imgcaption fig3|>{{ :​prepro_fig3.png?​400 |}}</​imgcaption>​
-  * For ICODE = 10, the first node represents the sphere center, the second node has the (R, F) coordinates where F allows to see only one part of the sphere. \\ F=0.abcdef \\ a=1 : positive X part visible \\ b=1 : negative X part visible \\ c=1 : positive Y part visible \\ d=1 : negative Y part visible \\ e=1 : positive Z part visible \\ f=0 : negative Z part visible. \\ Example: F=0.101010 draws the X, Y, Z positive parts. +  * __For ICODE = 10__, the first node represents the sphere center, the second node has the (R, F) coordinates where F allows to see only one part of the sphere. \\ F=0.abcdef \\ a=1 : positive X part visible \\ b=1 : negative X part visible \\ c=1 : positive Y part visible \\ d=1 : negative Y part visible \\ e=1 : positive Z part visible \\ f=0 : negative Z part visible. \\ Example: F=0.101010 draws the X, Y, Z positive parts. 
-  * For ICODE = 98 and 99, this segment is composed by only one node which DOF pilot the supposed rigid foundation. \\ This segment type is unique for a given foundation and __MUST__ be placed in the first position. The pilot node DOF can be free or fixed (fixed or imposed displacement),​ but __all the other foundation nodes must be fixed__. The pilot node DOF must be consistent with the analysis type. In case of rotations, it means finite rotations. The foundation follows its pilot, its displacement can be defined by DMULT or *.DEP. In the case ICODE = 98, the foundation can have rotations. +  * __For ICODE = 98 and 99__, this segment is composed by only one node which DOF pilot the supposed rigid foundation. \\ This segment type is unique for a given foundation and __MUST__ be placed in the first position. The pilot node DOF can be free or fixed (fixed or imposed displacement),​ but __all the other foundation nodes must be fixed__. The pilot node DOF must be consistent with the analysis type. In case of rotations, it means finite rotations. The foundation follows its pilot, its displacement can be defined by DMULT or *.DEP. In the case ICODE = 98, the foundation can have rotations. 
-  * Case ICODE = 98, for the 2D case, the DOF 1 represents the rotation around the perpendicular to the plan (2), on the other hand for the 3D case, the DOF 1 represents the rotation around x, the DOF 2 represents the rotation around y and the DOF 3 the rotation around z. All these rotations are given in radian.+  * __For ICODE = 98__, for the 2D case, the DOF 1 represents the rotation around the perpendicular to the plan (2), on the other hand for the 3D case, the DOF 1 represents the rotation around x, the DOF 2 represents the rotation around y and the DOF 3 the rotation around z. All these rotations are given in radian.
  
 ---- ----
prepro.1560412854.txt.gz · Last modified: 2020/08/25 15:36 (external edit)