请教双压电驱动悬臂梁振动的仿真！

cyqwan007

小弟是新手，现在做双压电驱动悬臂梁振动的仿真，即一金属悬臂梁，上下两压电片，分别给正弦和余弦信号驱动，由于刚开始学，攒了一段程序，请高手详细指教，不胜感激！
/prep7
thb=1e-3                                          !thickness of beam
thp=0.3e-3                                      !thickness of piezo
wb=20e-3                                         !width of beam
lb=40e-3                                         !length of beam
wp=20e-3                                        !width of piezo
lp=10e-3                                         !length of piezo
/units,si                                             !modelling of structure
block,0,lp,0,wp,0,thp
block,0,lp,0,wb,thp,thb+thp
block,lp,lb,0,wb,thp,thb+thp
block,0,lp,0,wp,thb+thp,thb+thp+thp
vglue,all                              
et,1,solid5,3                                        !define element 1 with material 3 of solid5
mp,dens,3,7600
mp,perx,3,7.124e-9                             !input of permittivity matrix
mp,pery,3,7.124e-9
mp,perz,3,5.841e-9
tb,piez,3                                            !input of piezoelectric matrix [e] or [d]
tbdata,16,10.5
tbdata,14,10.5
tbdata,3,-4.1
tbdata,6,-4.1
tbdata,9,14.1
tb,anel,3                                           !input of elastic coefficient matrix [c]
tbdata,1,13.2e10,7.1e10,7.3e10
tbdata,7,13.2e10,7.3e10
tbdata,12,11.5e10
tbdata,16,3.0e10
tbdata,19,2.6e10
tbdata,21,2.6e10
et,2,solid45                                     !define element 2 with material 2 of solid45
mp,ex,2,1.06e11
mp,nuxy,2,0.35
mp,dens,2,8920
esize,1e-3                           
type,1
mat,3
vmesh,5
vmesh,6
type,2
mat,2
vmesh,1
vmesh,7
nsel,s,loc,z,thp+thb+thp
nsel,r,loc,x,0,lp
cp,1,volt,all                       !add volt 1
nsel,all
nsel,s,loc,z,thp+thb
nsel,r,loc,x,0,lp
cp,2,volt,all
d,all,volt,0                         !to the gnd
nsel,all
nsel,s,loc,z,thp
nsel,r,loc,x,0,lp
cp,2,volt,all
d,all,volt,0                         !to the gnd
nsel,all
nsel,s,loc,z,0
nsel,r,loc,x,0,lp
cp,3,volt,all                        !add volt 3         
nsel,s,loc,x,0
nsel,r,loc,z,thp,thp+thb
d,all,ux,0,,,,uy,uz               !constraint all DOFs at the left
nsel,all
*SET,t,1e-3
*SET,con,1e-7   
*SET,n,60   
*SET,pi,2*asin(1)   
*SET,Fv,74300   
NEQIT,1
finish
/solu   
antyp,trans   
trnopt,full
tintp,,0.25,0.5,0.5   
alphad,0.25           
betad,0.5            
kbc,1                 
*CREATE,LOAD
time,t  
volt1=100*sin(2*pi*Fv*t)
volt2=100*cos(2*pi*Fv*t)
d,1,volt,volt1   
d,3,volt,volt2
alls
t=t+1e-3
OUTPR,BASIC,1
OUTRES,ALL,1
SOLVE   
*END
*DO,I,1,41  
*USE,LOAD   
*enddo
fini

py316

不错。顶

onattacking

这个程序，材料类型和单元类型有问题吧

wg6689

运行出错哈

wangfu2000

运行出错哈