急问！！！！电压加载问题

cyqwan007

将程序输入以后，得出结果位移都为0，而我已经加载了电压v1,急问高手这是怎么回事？？？
十分感谢！！！！！！！！！！！！！



/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,1.06e11
mp,nuxy,2,0.35
mp,dens,2,8920
esize,1e-3                                         !minimum size for mesh
type,1
mat,3
vmesh,5
vmesh,6
type,2
mat,2
vmesh,1
vmesh,7


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,v1,node                    !define v1         
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,pi,2*asin(1)   
*SET,Fv,74300   
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)
cmsel,s,v1
d,all,volt,volt1   
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

adwph

函数定义有问题吧,如果函数定义正确的话,log文件中应该有个表格文件的.个人认为.

li_1025tian

个人感觉电压加载时有问题

Hallen

我感觉alphad,0.25                          !Defines the mass matrix multiplier for damping       有问题
          betad,0.5                             !Defines the stiffness matrix multiplier for damping