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

momochong

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

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

我是一个濒

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

Mr_DJ

强烈要求 momochong  做本版版主

momochong

合并起来，命令流文件如下：
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 [sqrt(%r1*%r1*%d*%d-%d1*%d1)-sqrt(%r2*%r2*%d*%d-%d1*%d1)]
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)] [0.815*12.35]
sd 9 sp 0 0 [-sqrt(0.815*0.815*12.35*12.35-6.175*6.175)] [0.88*12.35]
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.0  0.0 1.0 0.0  0.0001 symm ;
plane 2 0.0 0.0 0.0  1.0 0.0 0.0  0.0001 symm ;
c lsdyna keyword
bptol 1 2 -1
bptol 2 3 -1
stp .008
write
当然，还有eos命令没有给出，还有控制接触方式、控制计算时间、输出间隔、计算算法等等参数的命令没有给出。后续一并贴出。

tiny0o0

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

momochong

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

momochong

我来把一些命令添加上。
药型罩紫铜采用采用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

壳体钢材料采用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

炸药材料采用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

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

tiny0o0

加油，继续给补上

mahuaibo

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

momochong

补充上状态方程的定义命令流，如下
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

加油，真的很不错

jyl52lg

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

狼跃冲顶

楼主用的真好

九鼎江湖

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 18:39 发表
parameters
mate0        0            
shapedmate   [%mate0 +1]   
shellmate   [%shapedmate +1]   
explosivemat   [%shellmate +1]   
airmate   [%explosivemate +1]   
c 定义材料参数
lsdymats %s ...

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

haige

强悍  高手 对楼主敬仰

sammy388

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

zyzhx

楼主厉害啊！感谢！