7 LSDYNA 中的完全重启动是什么功能,实际使用过程中,如何实现完全重启动

  • 15068
请先 登录 后评论

4 个回答

王玥琪

LS DYNA 的重启动包含三种类型:

  • 简单重启动
  • 小型重启动
  • 完全重启动

任何类型的重启动都需要首先生成 dump 文件,其为重启动文件,包含计算过程中的所有数据,重启动都是从特定的时间点开始的,重启动的二进制文件包含模型中的所有数据,如应力、应变、变形等,输出重启动文件的命令如下:

*DATABASE_BINARY_D3DUMP
*DATABASE_BINARY_RUNRSF

其中 D3DUMP 为重启动后继续计算,RUNRSF为重启动后覆盖原有计算,前者为默认的重启动文件格式,如果不想在计算结束后生成 D3DUMP 重启动文件,需要在 K 文件中输入如下命令:

d=nodump

RUNRSF重启动保留的文件数由参数 NR 设定,默认情况下为 1,如果 NR>1,将以循环方式写入。

LSDYNA 在进行重启动分析的时候,有两个特征,其一:需要使用相同版本的计算软件开启重启动计算,其二:重启动计算默认使用相同的 CPU 数目

1. 关于简单重启动

简单重启动就是在不改变任何模型和计算文件的情况下,继续展开计算,输入文件没有任何的变化,控制台调用命令如下:

LSDYNA r=d3dump01

2. 小型重启动

新型重启动仅可以修改如下的特定参数,其它参数不可修改:

  • 修改终止时间
  • 修改计算结果输出间隔
  • 修改时间步控制参数
  • 修改加载曲线(加载曲线中的数据点数量不能更改)
  • 添加节点约束
  • 删除 contacts、parts、elements
  • 将 Parts 转为刚性或非刚性

详细内容可以查看用户手册中的 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

使用控制台调用求解器,命令如下:

3. 完全重启动

如果修改 PARTS、加载和接触,则需要使用完全重启动,完全重新启动需要一个二进制转储文件和一个完全重新启动输入文件。在完全重启输入文件中,给出了模型的完整关键字描述:

  • 输入保留节点、元件、零件、材料、触点、负载等原始数据;
  • 上一条中的输入文件中的数据可以修改;
  • 新的输入文件中可以添加和修改 PART、材料和接触;

在计算之前,必须使用命令*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

请先 登录 后评论
Fullmoonlong - 研究生

重启动不显示新加 part,现在还不知道怎么解决呢

请先 登录 后评论