LS DYNA 的重启动包含三种类型:
任何类型的重启动都需要首先生成 dump 文件,其为重启动文件,包含计算过程中的所有数据,重启动都是从特定的时间点开始的,重启动的二进制文件包含模型中的所有数据,如应力、应变、变形等,输出重启动文件的命令如下:
*DATABASE_BINARY_D3DUMP
*DATABASE_BINARY_RUNRSF
其中 D3DUMP 为重启动后继续计算,RUNRSF为重启动后覆盖原有计算,前者为默认的重启动文件格式,如果不想在计算结束后生成 D3DUMP 重启动文件,需要在 K 文件中输入如下命令:
d=nodump
RUNRSF重启动保留的文件数由参数 NR 设定,默认情况下为 1,如果 NR>1,将以循环方式写入。
LSDYNA 在进行重启动分析的时候,有两个特征,其一:需要使用相同版本的计算软件开启重启动计算,其二:重启动计算默认使用相同的 CPU 数目。
简单重启动就是在不改变任何模型和计算文件的情况下,继续展开计算,输入文件没有任何的变化,控制台调用命令如下:
LSDYNA r=d3dump01
新型重启动仅可以修改如下的特定参数,其它参数不可修改:
详细内容可以查看用户手册中的 RESTART 部分。例如小型重启动中删除部分 PART,可以如下修改 K 文件:
*DELETE_PART
1
*DELETE_CONTACT
1
重启小型重启动需要一个二进制文件 dump 和 K 文件,K 文件如下所示:
*keyword
*CONTROL_TERMINATION
15e-03
*DATABASE_BINARY_d3plot
1e-5
*DELETE_PART
4,5
*DELETE_CONTACT
3
*END
使用控制台调用求解器,命令如下:
如果修改 PARTS、加载和接触,则需要使用完全重启动,完全重新启动需要一个二进制转储文件和一个完全重新启动输入文件。在完全重启输入文件中,给出了模型的完整关键字描述:
在计算之前,必须使用命令*STRESS_INITIALIZATION 初始化新输入的 PARTS,该命令会初始化应力分布等参数。但完全重启动不能改变单元间的连接和网格的拓扑,*NODE 节点坐标保持不变,*DELETE 命令只用于小型重启动中,在完全重启动中不能使用,在完全重启动中,则需要删除 PART 和 ELEMENT 下的数据;改变初始速度不能使用 *INITIAL_VELOCITY,而应该使用 *CHANGE_VELOCITY_OPTION,保留原有的 CONTACTS 时影该保留 ID。
注意,在使用 MPP 版本的 LS DYNA 和独立版本的 LS DYNA 版本启动重启动时,控制台而关键字有所不同,其中 MPP 版本的关键字为:
n=d3full01
而独立版本的 LS DYNA 版本为:
r=d3dump01
在小型重启动中,MPP 和独立版本的 LS DYNA 重启动的调用没有区别。
下面是基于 ANSYS LS-DYNA 9.0 的完全重启动模型。注意不要任何突然的位移变化都会导致加速度无穷大。
fini
/clear
/title, ANSYS LS-DYNA 9.0 Full Restart Example
/filnam,file_1
/plopts,info,1
/pnum,node,1
/pnum,elem,1
/view,,1,2,3
/prep7
et,1,SHELL163,10
r,1,,5,0.10
edint,5,4
mp,ex,1,30.0e6
mp,nuxy,1,0.30
mp,dens,1,0.00074
tb,biso,1
tbdata,1,60.0e3
tbdata,2,30.0e4
edmp,hgls,1,5
mat,1
type,1
real,1
n,1, 0.0, 0.0, 0.0
n,2, 1.0, 0.0, 0.0
ngen,4,2,1,2,1,0.0,1.0,0.0
e,1,2,4,3
e,3,4,6,5
e,5,6,8,7
eplot
edpart,create
eddamp,all,,0.02
eddamp,1,,1.0e-4
fini
/solu
nsel,s,loc,y,0.0
cm,nbase,node
esln
cm,ebase,elem
nsel,s,loc,y,3.0
cm,ntip,node
esln
cm,etip,elem
esel,all
nsel,all
eplot
d,nbase,ux,0.0,,,,uy,uz,rotx,roty,rotz
*dim,etime,,3
*dim,zdisp,,3
etime(1)=0.0,0.05,0.101
zdisp(1)=0.0,0.50,0.50
edload,add,uz,,ntip,etime(1),zdisp(1)
time,0.100
edrst,20
edhtime,200
edhist,ntip
edhist,nbase
edhist,ebase
edenergy,1,1,1,1
edout,glstat
edout,spcforc
edopt,add,,both
save
solve
save
fini
/post1
file,file_1,rst
set,last
/edge,,1
/dscale,,1
plnsol,u,z,2
/wait,2
fini
/post26
numvar,30
file,file_1,his
nsol,2,8,u,z,n8uz
prvar,2
plvar,2
/wait,2
fini
/solu
edstart,3,,,d3dump01 ! full restart - Jobname changed automatically ……
edis ! stress initialization for all parts ……
*dim,etime2,,3
*dim,zdisp2,,3
etime2(1)=0.10,0.15,0.201
zdisp2(1)=0.50,1.00,0.75
edload,add,uz,,ntip,etime2(1),zdisp2(1)
time,0.200
edrst,40
edhtime,400
edhist,ntip
edhist,nbase
edhist,ebase
edenergy,1,1,1,1
edout,glstat
edout,spcforc
edopt,add,,both
save
solve
save
fini
/post1
file,file_2,rst
set,last
/edge,,1
/dscale,,1
plnsol,u,z,2
/wait,2
fini
/post26
numvar,30
file,file_2,his
nsol,2,8,u,z,n8uz
prvar,2
plvar,2
/eof