LS-DYNA边界条件-粘弹性边界

这是ansys粘弹性边界条件,如何转化为lsdyna中二维模型粘弹性边界条件?
!以下建立底边界法向和切向弹簧阻尼单元
nsel,s,loc,y,0.
*get,np,node,,count !得到选中的结点数，存入np
*get,npmax,node,,num,maxd !得到已经定义的最大结点数，存入npmax
*do,ip,1,np
npnum=node((ip-1)*dxyz-L/2.,0.,0.)
x=nx(npnum)
y=ny(npnum)
z=nz(npnum)
npmax=npmax+1
n,npmax,x.,y-dxyz/2,z !定义底边界法向结点以便与边界点形成法向单元
type,3
real,3
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
npmax=npmax+1
n,npmax,x-dxyz/2.,y,z !定义底边界切向结点以便与边界点形成切向单元
type,2
real,2
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
*enddo
!以下建立左边界法向和切向弹簧阻尼单元
nsel,s,loc,x,-L/2
*get,np,node,,count !得到选中的结点数，存入np
*get,npmax,node,,num,maxd !得到已经定义的最大结点数，存入npmax
*do,ip,2,np !侧边界最下面一个点按底边界上处理
npnum=node(-L/2,(ip-1)*dxyz,0.)
x=nx(npnum)
y=ny(npnum)
z=nz(npnum)
npmax=npmax+1
n,npmax,x-dxyz/2.,y,z !定义左边界法向结点以便与边界点形成法向单元
type,3
real,3
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
npmax=npmax+1
n,npmax,x,y-dxyz/2.,z !定义左边界切向结点以便与边界点形成切向单元
type,2
real,2
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
*enddo
!以下建立右边界法向和切向弹簧阻尼单元
nsel,s,loc,x,L/2
*get,np,node,,count !得到选中的结点数，存入np
*get,npmax,node,,num,maxd !得到已经定义的最大结点数，存入npmax
*do,ip,2,np !侧边界最下面一个点按底边界上处理
npnum=node(L/2,(ip-1)*dxyz,0.)
x=nx(npnum)
y=ny(npnum)
z=nz(npnum)
npmax=npmax+1
n,npmax,x+dxyz/2.,y,z !定义右边界法向结点以便与边界点形成法向单元
type,3
real,3
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
npmax=npmax+1
n,npmax,x,y-dxyz/2.,z !定义右边界切向结点以便与边界点形成切向单元
type,2
real,2
e,npnum,npmax
d,npmax,all,0. !约束新生成的点
*enddo
allsel,all
/pnum,type,1
/number,1
eplot
finish
/solu
ANTYPE,trans
!*
TRNOPT,FULL
LUMPM,0
btime=0.02
etime=16.00
dtime=0.02
*DO,itime,btime,etime,dtime
TIME,itime
nsel,s,loc,y,H !选中需要加荷载的点
nsel,r,loc,x,-L/4,L/4
*if,itime,lt,1.,then
f,all,fy,1*itime
*elseif,itime,ge,1.0,and,itime,le,2.0
f,all,fy,1*(2-itime)
*else
f,all,fy,0.0
*endif
allsel,all
SOLVE
*ENDDO