- 积分
- 0
- 注册时间
- 2006-12-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
小弟是新手,现在做双压电驱动悬臂梁振动的仿真,即一金属悬臂梁,上下两压电片,分别给正弦和余弦信号驱动,由于刚开始学,攒了一段程序,请高手详细指教,不胜感激!
/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 |
|