ansys ls/dyna仿真引线键合问题

各位高手，小妹正在做一个有关引线键合仿真建模方面的问题，具体建模如附件里所示，代码如下，但是仿真出来的一直是下面焊盘没有变形，想咨询下各位具体约束和材料接触的定义，代码如下
/filename,3D simulation of wire bonding
/title,3D simulation of wire bonding
/prep7
!define constants
!pad size
padl=120
padh1=1.5
padh2=6.5
padh3=10
!wire&ball&Capillary
wd=20
fa=15
h=24
cd=32
t=100
ca=55
or=10
lh=50
a1=15
bna=15
fab=40
wh=100
h0=60
!others
pi=2*acos(0.0)
!element type
ET,1,shell163
ET,2,SOLID164
r,1,5/6,,2.5e-3
!material model
!pad material
mp,dens,1,2600
mp,ex,1,1.7e10
mp,nuxy,1,0.33
tb,bkin,1,,,,
tbdat,1,3.5e7!屈服应力
tbdat,2,1.35e8!切线模量
!wire&ball material
mp,dens,2,2600
mp,ex,2,7.5e10
mp,nuxy,2,0.42
tb,bkin,2,,,,,
tbdat,1,1.2e8!产生能量
tbdat,2,1.35e8!切线模量
!capillary，rigid劈刀和焊盘单元为163
edmp,rigid,3 !定义材料属性
mp,dens,3,7800
mp,ex,3,2e11
mp,nuxy,3,0.3
!The Pad Modeling and Meshing (Three Layers)3
K,1,-PADL/2,0,0
k,2,-PADL/2,PADH3,0
k,3,-PADL/2,PADH3+PADH2,0
k,4,-PADL/2,PADH3+PADH2+PADH1,0
k,5,PADL/2,PADH3+PADH2+PADH1,0
k,6,PADL/2,PADH3+PADH2,00
k,7,PADL/2,PADH3,0
k,8,PADL/2,0,0
a,3,4,5,6
a,2,3,6,7
a,1,2,7,8
asel,all
!The Wire & Ball Modeling and Meshing
K,9,0,PADH1+PADH2+PADH3+FAB/2,0
K,10,WD/2,PADH1+PADH2+PADH3+FAB/2+((FAB/2)**2-(WD/2)**2)**0.53
k,11,WD/2,PADH1+PADH2+PADH3+FAB/2+((FAB/2)**2-(WD/2)**2)**0.5+WH,0
k,12,FAB/2,PADH1+PADH2+PADH3+FAB/2,0
k,13,0,PADH1+PADH2+PADH3,0
k,14,0,PADH1+PADH2+PADH3+FAB,0
k,15,0,PADH1+PADH2+PADH3+FAB/2+((FAB/2)**2-(WD/2)**2)**0.5+WH,0
l,9,12
l,9,13
l,9,14
l,15,11
l,11,10
l,14,15
larc,10,12,9,FAB/2
larc,12,13,9,FAB/2
larc,14,10,9,FAB/2
al,11,12,18
al,11,13,19,17
al,16,14,15,19
asel,all
!The Capillary Modeling and Meshing
k,16,CD/2,PADH1+PADH2+PADH3+FAB,0
k,17,H/2,PADH1+PADH2+PADH3+FAB+(CD/2-H/2)/tan(pi*CA/180),0
k,18,H/2,PADH1+PADH2+PADH3+FAB+(CD/2-H/2)/tan(pi*CA/180)+LH,0
k,19,H/2+H0*tan(pi*A1/180),PADH1+PADH2+PADH3+FAB+(CD/2-H/2)/tan(pi*CA/180)+LH+H0,0
k,20,T/2-OR+OR*sin(pi*FA/180),PADH1+PADH2+PADH3+FAB+(T/2-OR+OR*sin(pi*FA/180)-CD/2)*tan(pi*FA/180),0
k,21,T/2,PADH1+PADH2+PADH3+FAB+(T/2-CD/2-OR)*tan(pi*FA/180)+OR-OR*sin(pi*BNA/180)
tempx=T/2+((CD/2-H/2)/tan(pi*CA/180)+LH+H0-((T/2-CD/2-OR)*tan(pi*FA/180)+OR-OR*sin(pi*BNA/180)))*tan(pi*BNA/180)
tempy=PADH1+PADH2+PADH3+FAB+(CD/2-H/2)/tan(pi*CA/180)+LH+H0
k,22,tempx,tempy,0
k,23,T/2-OR,PADH1+PADH2+PADH3+FAB+(T/2-CD/2-OR)*tan(pi*FA/180)+OR,0
k,24,T/2-OR,PADH1+PADH2+PADH3+FAB+(T/2-CD/2-OR)*tan(pi*FA/180)+OR,0
l,16,17
l,17,18
l,18,19
l,16,20
l,21,22
l,19,22
larc,21,20,23,OR
al,20,21,22,25,24,26,23
allsel
asel,s,,,1,3,1
vext,all,,,,,-padl/2
allsel
vrotat,4,5,6,7,,,9,14,180,2
!划分网格
vsel,s,,,4,6,1
vSEL,A,,,8,10,1
VATT,2,1,2
vmesh,all
allsel,all
ASEL,S,,,30,36,1
ASEL,A,,,47,53,1
ASEL,A,,,7
AATT,3,1,1
mshape,1,2d
mshkey,0
AMESH,ALL
ALLSEL
asel,s,,,1,3,1
asel,a,,,8,11,1
asel,a,,,13,19,2
asel,a,,,14
asel,a,,,18,20,2
aatt,1,1,1
mshape,1,2d
mshkey,0
amesh,all
allsel,all
nsel,s,loc,y,0
d,all,all
allsel,all
nsel,s,loc,z,0
d,all,uz,0,,,,uy,uz,rotx,roty,rotz
ALLSEL,ALL
edpart,create,1,all
allsel,all
edcgen,asts,2,1,0.13,0.13 !定义接触面
edcgen,asts,1,3,0.13,0.13
ALLSEL,ALL
*dim,time,array,2,1,1,,,
*set,time(1,1,1),0.0
*set,time(2,1,1),1.0
*dim,displacement,array,2,1,1,,,
*set,displacement(1,1,1),0.0
*set,displacement(2,1,1),-130
!use edload to load displacement loading
edload,add,RBUY,0,2,time,displacement,0,,,,,
!deltim,0.151E-03,0.001E-03,0.301E-03,off
allsel,all
time,1!设置求解的结束时间
EDCTS,0,0.9,
EDOPT,ADD,blank,ansys
edrst,100, !定义jobname.rst结果文件的输出步数
edhtime,100, !定义jobname.his结果文件的输出步数
eddump,1 !设置重启动问价d3dump的输出部署
edenergy,1,1,1,1 !设置能量控制选项
edout,all
edwrite,ansys,,,
save
solve