载荷有冲突,你给的代码不全,不知道具体原因,我转几个例子你比较找找毛病吧
(1)对压电陶瓷加z向电压,产生膨胀、收缩
- /PREP7
- /GRAPHICS,FULL
- !############################## parameter #######################
- *AFUN,DEG
- VOL=0
- VOH=20
- !########################### Material property #################
- !****************** 3---PZT **************************
- et,3,solid5 ! UX,UY,UZ,TEMP,VOLT,MAG
- !mp,nuxy,3,0.27
- mp,dens,3,7.64e3 ! kg/m3
- mp,perx,3,8.93e-9 ! F/m
- mp,pery,3,8.93e-9 ! F/m
- mp,perz,3,6.92e-9 ! F/m
- mp,kxx,3,1
- mp,murx,3,1
- tb,anel,3 ! Anisotropic elastic matrix 72A??? 99.3.4
- tbdata,1,11.565e10,6.489e10,6.229e10,0,0,0
- tbdata,7,11.565e10,6.229e10,0,0,0
- tbdata,12,9.298e10,0,0,0
- tbdata,16,1.786e10,0,0
- tbdata,19,1.786e10,0
- tbdata,21,1.786e10tb,piez,3 ! Piezoelectric matrix 72A??? 99.3.4
- tbdata,1,0,0,-12.3082
- tbdata,4,0,0,-12.3082
- tbdata,7,0,0,20.7574
- tbdata,10,0,0,0
- tbdata,13,0,17.0357,0
- tbdata,16,17.0357,0,0!------------------------ model---------------------------
- NUMSTR,KP,1
- NUMSTR,LINE,1
- NUMSTR,AREA,1block,0,10,0,5,0,1type,3
- mat,3
- esize,0.5
- vmesh,all
- !------------boundary condition-----------------------
- nsel,s,loc,x,0
- d,all,ux
- d,all,uy
- d,all,uz
- ! -------------COUPLE THE TWO FACE OF ELECTRODE -------------
- nsel,s,loc,z,1
- cp,11,volt,all
- *get,ntop,node,,num,minnsel,s,loc,z,0
- cp,21,volt,all
- *get,nbot,node,,num,minALLSfinish
- /eof
- !-------------------------------------------
- /soluantype,staticALLS
- d,nTOP,volt,voh
- d,nBOT,volt,volESEL,S,TYPE,,3
- NSLE
- D,ALL,TEMP
- !D,ALL,MAG
- ALLSnsel,all
- esel,allsavesolvefini/eof
复制代码
|