压电瞬态分析加载问题！

cyqwan007

我做的项目，悬臂梁上下各有一个压电片。在/solu中，先给上面压电片电信号，从而使梁发生弯曲；当弯曲最大时候，再给下面的压电片电信号。我想请教的是，在/solu中，如何加载（即一组循环，在每一步中，先激发上压电片，然后获取某点位移信息，和前一步位移信息比较，一旦比前面小，则激发下面的压电片）下面是我初步写的程序，请大家指教修改，十分感谢！
/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                              !cementation of each parts together
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,2.1e11
mp,nuxy,2,0.3
mp,dens,2,7800
esize,1e-3                            !minimum size for mesh
type,1
mat,3
vmesh,1
vmesh,6
type,2
mat,2
vmesh,5
vmesh,7
nsel,s,loc,z,thp+thb+thp
nsel,r,loc,x,0,lp
cm,v1,node              !add v1
nsel,all
nsel,s,loc,z,thp+thb
nsel,r,loc,x,0,lp
cm,gnd,node           !to the gnd
nsel,all
nsel,s,loc,z,thp
nsel,r,loc,x,0,lp
cm,gnd,node              !to the gnd
nsel,all
nsel,s,loc,z,0
nsel,r,loc,x,0,lp
cm,v2,node              !define v2         
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
cmsel,s,gnd
d,all,volt,0
*SET,t,1e-3
*SET,con,1e-7   
*SET,n,60   
*SET,pi,2*asin(1)   
*SET,Fv,500   
!NEQIT,1
!f,1322,fz,-0.1
finish
/solu   
antyp,trans   
trnopt,full
tintp,,0.25,0.5,0.5    !Defines transient integration parameters
alphad,0.25            !Defines the mass matrix multiplier for damping
betad,0.5              !Defines the stiffness matrix multiplier for damping
kbc,1                  !Specifies stepped or ramped loading within a load step
*CREATE,LOAD
time,t  
!volt1=100*sin(2*pi*Fv*t)
volt2=100*sin(2*pi*Fv*t)
!cmsel,s,v1
!d,all,volt,volt1   
cmsel,s,v2
d,all,volt,volt2
alls
t=t+1e-3
OUTPR,BASIC,1
OUTRES,ALL,1
SOLVE   
*END
*DO,I,1,41  
*USE,LOAD   
*enddo
fini
/post26
NUMVAR,20   
nsol,2,1322,u,x,u1322_x
nsol,3,1322,u,y,u1322_y
nsol,4,1322,u,z,u1322_z

daicy

回复 1 # cyqwan007 的帖子

你做出来没有？316185329，能否交流一下？

xuliujuan1987

请问用ansys如何在构件端部施加加速度载荷或者是施加冲击力，非常感谢！