论坛 › A18：前后处理软件 › 最近在做的一个数值计算（Dyna-tg）

momochong 发表于 2007-12-4 23:10:47

最近在做的一个数值计算（Dyna-tg）

模型是在tg中建立的。计算效果见图。
命令流见附件

[ 本帖最后由 tiny0o0 于 2008-1-1 16:23 编辑 ]

我是一个濒 发表于 2016-10-21 09:17:40

请问楼主，你的接触是在哪里定义的呢？

Mr_DJ 发表于 2007-12-6 22:52:19

强烈要求 momochong做本版版主

momochong 发表于 2007-12-6 21:55:04

合并起来，命令流文件如下：
parameters
             mate0      0            
             shapedmate   [%mate0 +1]   
             shellmate   [%shapedmate +1]   
             explosivemat   [%shellmate +1]   
             airmate   [%explosivemate +1]   
;
c 定义材料参数
lsdymats %shapedmate 11 brick elfob csb
head shapedchargem(copper solids)
rho 8.93 c copper's density
g0 0.477
sig0 0.12e-02
beta 36.0
n 0.450
gama 0.00
sigm 0.640e-02
b 2.83
bp 2.83
h 0.377e-03
f 0.100e-02
t0 0.179e+04
gam0 2.02
sa 1.50
ispall hydro
pcut -9.00
a 63.5
r 0.00
ivar 0.00
ec0 0.00
ec1 0.00
ec2 0.00
ec3 0.00
ec4 0.00
ec5 0.00
ec6 0.00
ec7 0.00
ec8 0.00
ec9 0.00
;
lsdymats %shellmate 10 brick elfob csb
head shell material(steel solids)
rho 7.86
ispall hydro
g 2.069
sigy 1.72e-02
etan 0.00
pcut -7.030
a1 0.0
a2 0.0
;
lsdymats %explosivemat 8 brick elfob csb
head explosive material(PBX)
rho 1.717
d 7.98e-02
pcj 2.950e-01
beta 0.0
k 0.0
g 0.0
sigy 0.0
;
lsdymats %airmate 15 brick elfob csb
head air material(air null)
rho 1.293e-03
;
parameter
r1 0.88
r2 0.815
d 12.35
d1 [%d/2]
p1
p2 11.43
;
sd 1 cy 0 0 0 0 0 1 6.35
sd 2 cy 0 0 0 0 0 1 6
sd 8 sp 0 0 [-sqrt(0.815*0.815*12.35*12.35-6.175*6.175)]
sd 9 sp 0 0 [-sqrt(0.815*0.815*12.35*12.35-6.175*6.175)]
sd 3 cy 0 0 0 0 0 1 6.7
sd 4 cy 0 0 0 0 0 1 6.175
sd 5 cy 0 0 0 0 0 1 9
c shaped
block 1 11 21;1 11 21;1 7;
-2.2 -2.2 0;
0 2.2 2.2;
0 %p1;
dei 1 2;2 3;;
sfi -1 3;1 -3;;sd 4
sfi ;;-1;sd 8
sfi ;;-2;sd 9
res 1 1 1 3 1 2 i 1
res 3 1 2 3 3 1 j 1
pb 2 2 2 2 2 1 xy -2.7 2.7
mate %shapedmate
endpart
c shell
block 1 3 4 12 20;1 9 17 18 20;1 2 28 30;
-2.2 -2.2 -2.2 -2.2 0;
0 2.2 2.2 2.2 2.2;
-0.175 0 %p2 [%p2+0.35];
dei 1 4;2 5;;
sfi -1 5;1 -5;;sd 3
sfi -2 5;1 -4;;sd 1
sfi -3 5;1 -3;;sd 4
dei 2 5;1 4;2 3;
dei 3 5;1 3;1 2;
res 3 1 1 5 1 4 i 1
res 5 1 4 5 3 1 j 1
pb 4 2 4 4 2 1 xy -2.7 2.7
nri -1;;;
nri ;-5;;
nri 1 3;1 2;-4;
nri 4 5;4 5;-4;
mate %shellmate
endpart
c air and explosive
block 1 8 9 19 29;1 11 21 22 29;1 8 14 34 38;
-2.2 -2.2 -2.2 -2.2 0;
0 2.2 2.2 2.2 2.2;
[-2.0*0.7] 0 %p1 %p2 [%p2+2*0.7];
dei 1 4;2 5;;
sfi -1 5;1 -5;;sd 5
sfi -2 5;1 -4;;sd 1
sfi -3 5;1 -3;;sd 4
sfi 3 5;1 3;-2;sd 8
sfi 3 5;1 3;-3;sd 9
res 3 1 1 5 1 5 i 1
res 5 1 5 5 3 1 j 1
pb 4 2 1 4 2 5 xy -2.7 2.7
dei 3 5;1 3;2 3;
mate %airmate
mti 2 5;1 2;2 4; %explosivemate
mti 4 5;2 4;2 4; %explosivemate
endpart
block 1 11 21;1 11 21;1 31;
-2.2 -2.2 0;0 2.2 2.2 ;-60.0 [-2.0*0.7];
dei 1 2;2 3;;
sfi -1 3;1 -3;;sd 4
res 1 1 1 3 1 2 i 1
res 3 1 2 3 3 1 j 1
nri -1;;;
nri ;-3;;
nri ;;-1;
pb 2 2 2 2 2 1 xy -2.7 2.7
mate %airmate
endpart
lsdyopts
lsdyna keyword
merge
plane 1 0.0 0.0 0.00.0 1.0 0.00.0001 symm ;
plane 2 0.0 0.0 0.01.0 0.0 0.00.0001 symm ;
c lsdyna keyword
bptol 1 2 -1
bptol 2 3 -1
stp .008
write
当然，还有eos命令没有给出，还有控制接触方式、控制计算时间、输出间隔、计算算法等等参数的命令没有给出。后续一并贴出。

tiny0o0 发表于 2007-12-5 21:38:14

哦，这个是做什么分析，看起来有点意思

momochong 发表于 2007-12-6 20:28:27

这个是做爆炸成型弹丸（EFP）战斗部的仿真。
TG对lsdyna非常支持。lsdyna中几乎所有的关键字都支持。
如果真做好了，在TG中建立模型和设置输出lsdyna所要的关键字的命令，输出的lsdyna格式的k文件几乎原封不动（或者稍微做修改即可）就可以进行仿真计算。我正在实现整个命令流。
TG安装目路下的example中有个Brazil文件夹，就是一个对输出结果不需要做任何修改即可在lsdyna求解器中进行计算的典型例子。

momochong 发表于 2007-12-6 21:34:37

我来把一些命令添加上。
药型罩紫铜采用采用11号材料。定义的命令流如下：
lsdymats %shapedmate 11 brick elfob csb
head shapedchargem(copper solids)
rho 8.93 c copper's density
g0 0.477
sig0 0.12e-02
beta 36.0
n 0.450
gama 0.00
sigm 0.640e-02
b 2.83
bp 2.83
h 0.377e-03
f 0.100e-02
t0 0.179e+04
gam0 2.02
sa 1.50
ispall hydro
pcut -9.00
a 63.5
r 0.00
ivar 0.00
ec0 0.00
ec1 0.00
ec2 0.00
ec3 0.00
ec4 0.00
ec5 0.00
ec6 0.00
ec7 0.00
ec8 0.00
ec9 0.00
;

[ 本帖最后由 momochong 于 2007-12-6 21:43 编辑 ]

momochong 发表于 2007-12-6 21:40:44

壳体钢材料采用10号材料。定义的命令流如下：
lsdymats %shellmate 10 brick elfob csb
head shellmate(steel solids)
rho 7.86
ispall hydro
g 2.069
sigy 1.72e-02
etan 0.00
pcut -7.030
a1 0.0
a2 0.0
;
不过TG对10号材料的支持不够完善，在输出的10号材料 定义的关键字中，将关键字*mat_elastic_plastic_hydro改为*mat_elastic_plastic_hydro_spall 即可。

[ 本帖最后由 momochong 于 2007-12-6 21:43 编辑 ]

momochong 发表于 2007-12-6 21:45:27

炸药材料采用8号材料，定义材料的命令如下：
lsdymats %explosivemat 8 brick elfob csb
head explosive material(PBX)
rho 1.717
d 7.98e-02
pcj 2.950e-01
beta 0.0
k 0.0
g 0.0
sigy 0.0
;

momochong 发表于 2007-12-6 21:47:52

空气采用15号材料模型，定义的命令流如下：
lsdymats %airmate 15 brick elfob csb
head air material(air null)
rho 1.293e-03
;

tiny0o0 发表于 2007-12-6 22:20:09

加油，继续给补上

mahuaibo 发表于 2007-12-7 19:49:33

毛毛虫，好厉害啊！

毛毛虫，好厉害啊！
你的QQ多少，有时间多交流啊，我也做了很长时间这方面的仿真了！

momochong 发表于 2007-12-12 22:09:28

补充上状态方程的定义命令流，如下
c 药型罩材料的状态方程定义:采用4号GRUNEISEN状态方程
lsdyeos %shapedmate 4
vci 0.394
s1 1.49
s2 0.0
s3 0.0
gamma 2.02
sa 0.470
e0 0.00
v0 1.0
;
c 壳体材料的状态方程定义:采用4号GRUNEISEN状态方程
lsdyeos %shellmate 4
vci 0.374
s1 1.92
s2 0.0
s3 0.0
gamma 1.69
sa 0.0
e0 0.0
v0 1.0
;
c 炸药材料的状态方程定义:采用2号ＪＷＬ状态方程
lsdyeos %explosivemate 2
a 5.242
b 7.678e-02
r1 4.2
r2 1.10
omega 3.4e-01
e0 8.50e-02
v0 1.0
;
c 空气材料状态方车程定义:采用１号linear polynomial状态方程
lsdyeos %airmate 1
c0 0.0
c1 0.0
c2 0.0
c3 0.0
c4 4.0e-01
c5 4.0e-01
c6 0.0
e0 2.5e-06
v0 1.0
;

XIAOKANG82 发表于 2007-12-29 10:40:41

加油，真的很不错

jyl52lg 发表于 2008-1-8 17:45:28

正在尝试，很有帮助，谢谢

狼跃冲顶 发表于 2008-1-13 21:28:15

楼主用的真好

九鼎江湖 发表于 2008-4-18 18:39:01

parameters
mate0      0            
shapedmate   [%mate0 +1]   
shellmate   [%shapedmate +1]   
explosivemat   [%shellmate +1]   
airmate   [%explosivemate +1]   
c 定义材料参数
lsdymats %shapedmate 11 brick elfob csb
head shapedchargem(copper solids)
rho 8.93 c copper's density
g0 0.477
sig0 0.12e-02
.......
请问：
head后面的语句是什么意思？？？？？
有什么格式、要求需要遵循？？？

mingbai1984 发表于 2008-4-18 21:39:13

原帖由 九鼎江湖 于 2008-4-18 18:39 发表 http://www.simwe.com/forum/images/common/back.gif
parameters
mate0      0            
shapedmate   [%mate0 +1]   
shellmate   [%shapedmate +1]   
explosivemat   [%shellmate +1]   
airmate   [%explosivemate +1]   
c 定义材料参数
lsdymats %s ...

comment，注释，应该是在title里面的内容

haige 发表于 2008-4-23 15:46:42

强悍高手 对楼主敬仰

sammy388 发表于 2009-8-30 17:01:48

感谢牛人楼主，我要努力学习啦！！！

zyzhx 发表于 2016-6-14 07:55:57

楼主厉害啊！感谢！

查看完整版本: 最近在做的一个数值计算（Dyna-tg）