*INCLUDE,INPUT = B010_2_shaft.inp
*INCLUDE,INPUT = B010_2_snapFit.inp
*Node, NSET=cNode
*ELSET, ELSET=Eall1
EvolumesC3D10
EvolumesC3D10_2
** ■物性値の定義
*MATERIAL, NAME=SolidMaterial
*ELASTIC
2100000, 0.300
*SOLID SECTION,ELSET=EvolumesC3D10,MATERIAL=SolidMaterial
*MATERIAL, NAME=SolidMaterial2
*ELASTIC
23000, 0.370
*SOLID SECTION,ELSET=EvolumesC3D10_2,MATERIAL=SolidMaterial2
*RIGID BODY,ELSET=EvolumesC3D10,REF NODE=22000,ROT NODE=22001
*CONTACT PAIR, INTERACTION=INTFemConstraintContact,TYPE=SURFACE TO SURFACE
face2,face1
*SURFACE INTERACTION, NAME=INTFemConstraintContact
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
 1000000.0
 *FRICTION
0.1,500000.
*BOUNDARY
NodeBC_Group_1,1,3
**NodeBC_Group_2,2,3
**NodeBC_Group_2,1,1,-60
**NodeBC_Group_3,3,3
22000,1,1,-60
22000,2,2
22000,3,3
22001,1,3
*STEP, NLGEOM
*STATIC, SOLVER=SPOOLES,direct
0.05,1.0
** ■出力指定
*NODE FILE, OUTPUT=3d
 U
*EL FILE
 S, E
**出力用の節点セット
*NSET,NSET=NodeBC_Group_4
**出力指定 uは変位
*NODE PRINT,NSET=NodeBC_Group_4
 u
**出力指定　RFは反力　セットに複数節点がある場合にはTOTALSで合計
*NODE PRINT, NSET=NodeBC_Group_4, TOTALS=ONLY
 RF
*END STEP