用某ai写k文件:
写一下lsdyna软件的小球碰撞平板的k文件,小球半径5cm,平板长宽10cm厚2cm,小球速度100m/s。
结果如下~
*KEYWORD
*TITLE
Sphere Impact on Plate (R=5cm, V=100m/s)
*CONTROL_TERMINATION
0.001 ! 仿真时间1ms
*CONTROL_TIMESTEP
0.0, 0.9 ! 时间步长系数0.9
*CONTROL_CONTACT
0.1, 0.1, 2, 2, 0, 1 ! 接触控制参数
*DATABASE_BINARY_D3PLOT
5.0e-5 ! 每0.05ms输出结果
*MAT_PLASTIC_KINEMATIC
1,7850,2.1e11,0.3,4e8,0.0,40.0 ! 小球材料(钢)
*MAT_PLASTIC_KINEMATIC
2,7850,2.1e11,0.3,4e8,0.0,40.0 ! 平板材料
*SECTION_SOLID
1,1 ! 小球单元类型(常应力)
*SECTION_SOLID
2,1 ! 平板单元类型
*PART
1,1,1 ! 小球(材料1,截面1)
*PART
2,2,2 ! 平板(材料2,截面2)
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
1,1,2,0.1 ! 接触定义(摩擦系数0.1)
*BOUNDARY_SPC_SET
1,0,0,0,1,1,1 ! 固定平板底面所有自由度
*INITIAL_VELOCITY_NODE
1,0,0,-100 ! 小球初始速度-Z方向100m/s
*HOURGLASS
1,5 ! 沙漏控制(类型5)
*NODE
!---------------- 小球节点(示例,需细化) ----------------
101, 0.0, 0.0, 7.1 ! 球心(Z=7.1cm)
102, 5.0, 0.0, 7.1 ! 表面节点X+
103, 0.0, 5.0, 7.1 ! 表面节点Y+
... ! (实际需完整球体网格)
!---------------- 平板节点(示例,需细化) ----------------
201, -5.0,-5.0, 0.0 ! 平板角点
202, 5.0,-5.0, 0.0
203, 5.0,5.0, 0.0
204, -5.0,5.0, 0.0
205, -5.0,-5.0, 2.0
... ! (实际需完整长方体网格)
*ELEMENT_SOLID
!---------------- 小球单元(示例) ----------------
1,1,101,102,103,... ! 需完整单元定义
...
!---------------- 平板单元(示例) ----------------
1001,2,201,202,203,...! 需完整单元定义
...
*END