*INCLUDE,INPUT = T214_2_nodeElem.inp
***要素番号、節点番号の指定
*ELEMENT, TYPE=MASS  , ELSET=EM1
***以下は重力加速度(9810mm/s^2）以下で1000N相当 e-9
*MASS,ELSET=EM1
0.101936799e-3
*****要素番号、節点番号の指定
***ELEMENT, TYPE=MASS  , ELSET=EM2
**10002,14288
**
*****以下は重力加速度(9810mm/s^2）以下で1000N相当 e-9
***MASS,ELSET=EM2
**0.101936799e-3
** Define element set Eall
*ELSET, ELSET=Eall
EfacesB32R,
EvolumesC3D10,
EM1,
***********************************************************
** Materials
** written by write_materials function
** Young's modulus unit is MPa = N/mm2
** FreeCAD material name: CalculiX-Steel
** SolidMaterial
*MATERIAL, NAME=SolidMaterial
*ELASTIC
210000, 0.300
*DENSITY
7.900e-9
*MATERIAL, NAME=SolidMaterial_bar
*ELASTIC
21000000, 0.300
*DENSITY
7.900e-12
***********************************************************
** Sections
** written by write_femelementsets function
**ビーム(縦横サイズ、Vベクトル）
*BEAM SECTION, ELSET=EfacesB32R, MATERIAL=SolidMaterial_bar, SECTION=RECT
 1 , 1
 0.0, 1.0, 0
**
****シェル（板厚）（*解除のこと）
***SHELL SECTION, ELSET=EfacesS6, MATERIAL=SolidMaterial
**10.0
**ソリッド（*解除のこと）
*SOLID SECTION,ELSET=EvolumesC3D10,MATERIAL=SolidMaterial
***********************************************************
** At least one step is needed to run an CalculiX analysis of FreeCAD
** written by write_step_begin function
*STEP
**BUCKLE
**3
*STATIC
***********************************************************
** Displacement constraint applied
** written by write_constraints_displacement function
** FemConstraintDisplacement
*BOUNDARY
NodeBC_Group_1,1
NodeBC_Group_1,2
NodeBC_Group_1,3
NodeBC_Group_1,5
NodeBC_Group_1,4
NodeBC_Group_1,6
***********************************************************
** Node loads Constraints
** written by write_constraints_force function
*CLOAD
** FemConstraintForce
** node load on shape: Edge:Vertex1
NodeBC_Group_2,2,-1000
***DLOAD
**Eall,GRAV,9810,0.0,-1.0,0.0
***********************************************************
** Outputs --> frd file
** written by write_outputs_types function
*NODE FILE, OUTPUT=2d
U
*EL FILE,SECTION FORCES
S, E
**出力指定　RFは反力　セットに複数節点がある場合にはTOTALSで合計
*NODE PRINT, NSET=NodeBC_Group_1, TOTALS=ONLY
 RF
**出力用の要素セット
*ELSET, ELSET=ELGROUP1
EfacesB32R
*EL PRINT,ELSET=ELGROUP1,SECTION FORCES
 S,E
***********************************************************
** written by write_step_end function
*END STEP