- 积分
- 25
- 注册时间
- 2003-3-16
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2012-7-6 22:52:39
|
显示全部楼层
来自 瑞典
K文件如下。和DYNA的格式有些类似,但是比DYNA要简介明了很多,用户需要调整的参数非常少。
如果感兴趣的可以对照 http://manual.impetus-afea.com/solver/3.0b/ 在线用户手册学习K文件。
*UNIT_SYSTEM
SI
*PARAMETER
np = 4000000
sf = 1.0
*INCLUDE
tube.k
*INCLUDE_TRANSFORM
1-coarse-half.k
1.0, 1.0, 1.0, 1000000, 1000000, 0
-0.008, 0.0003, 0.0, -0.008, 0.0, 0.0
1.0, 0.0, 0.0, 0.0, 0.0,-1.0
*TIME
2.0e-3
*OUTPUT
2.0e-5, 1.0e-6
*GPU
0, 1, 1
*CHANGE_P-ORDER
PS, 1, 3
*MAT_RIGID
1, 7800.0
*MAT_JC
2, 8520.0, 1.15e11, 0.31
2.06e8, 5.05e8, 0.42, 0.00, 1.68, 293.0, 1189.0, 1.0
385.0, 0.9, 0, 0, 0, 0, 0, 9.14e8
*MAT_JC
3, 7850.0, 2.1e11, 0.33
1.2e9, 5.0e10, 1.0, 0.0, 1.0, 293.0, 1800.0, 1.0
452.0, 0.9, 0, 0, 0, 0, 0, 1.0e20
*MAT_JC
4, 10660.0, 1.0e9, 0.42
0.24e8, 3.0e8, 1.0, 0.0, 1.0, 293.0, 760.0, 1.0
124.0, 0.9, 0, 0, 0, 0, 0, 1.75e8
*PART
"tube"
1, 1
"jacket"
2, 2, 0, 0, 0, 0.5e-8
"core"
3, 3, 0, 0, 0, 0.5e-8
"lead tip"
4, 4, 0, 0, 0, 0.5e-8
*INITIAL_VELOCITY
PS, 1, 920.0, 0.0, 0.0
*PBLAST
PS, 2, 0, [%np], 0, 0, wet, 0
0, 0, 0, 0, 1, 0, 0, 1.0
1, 2, 0
*GEOMETRY_BOX
1
[-0.1*%sf], [-0.17*%sf], 0.0, [0.6*%sf], [0.17*%sf], [0.17*%sf]
*GEOMETRY_CYLINDER
2
0.0, 0.0, 0.0, [%sf*0.5], 0.0, 0.0, [%sf*0.159]
*CONTACT
PS, 1, PS, 1, 0.1, -1.0e14
1
*SYMMETRY
Z, 0, 0, 0, 1.0e-5
*SET_PART
1
2..4
*SET_PART
2
1..4
*END
|
|