Make Geometry with Gmsh and define Physical groups |
cl1 = 5; Point(1) = {10, 0, 0, cl1}; Point(2) = {50, 0, 0, cl1}; Point(3) = {100, 0, 0, cl1}; Line(1) = {1, 2}; Line(2) = {2, 3}; Extrude {0, 10, 0} { Line{1, 2}; } Physical Surface(181) = {6}; Physical Surface(182) = {10}; Physical Line(183) = {4}; Physical Line(184) = {7}; |
|
Mesh with 2D-2order |
Mesh 2D set Order 2 |
|
Mesh save |
Save as gmsh.inp (Mesh -- Abaqus INP) |
|
Mesh Convert to Calculix |
>python gi2ci2.py gmsh.inp calc.inp 2 CPS |
1) gmsh.inp ---- input file (gmsh's abaqus format ) 2) calc.inp ---- output file (calculix format , line elements are converted to node group ) 3) 2 ---- 2D 4) CPS --- PreWord 'CPS' will make Element Type 'CPS6' |
Check and make BC with calculix |
>cgx -c calc.inp |
|
prnt set comp CURVE184 down plot e SURFACE181 g view elem on view edge off frame plus e SURFACE182 b plus n CURVE183 plus f CURVE184 send all abq names send CURVE184 abq pres 1.0 |
||
make solver file (sol.inp) |
*INCLUDE, INPUT=calc.inp *INCLUDE, INPUT=all.nam *MATERIAL, NAME=STEEL *ELASTIC 200000 , 0.3 *MATERIAL, NAME=PLA *ELASTIC 1000 , 0.3 *SOLID SECTION , Elset=SURFACE181 , Material=STEEL 20 *SOLID SECTION , Elset=SURFACE182 , Material=PLA 20 *STEP *STATIC *BOUNDARY CURVE183,1,3,0 *DLOAD *include,input=CURVE184.dlo *NODE PRINT,NSET=Nall U *EL PRINT,ELSET=Eall S *NODE FILE U *EL FILE S *END STEP |
|
solver start |
>ccx sol |
|
view result |
>cgx -v sol.frd |
Max Displacement = 9.52[mm] (3D Result was 9.90[mm] ) |