压电材料系数输入原理与工具

feasky

前面看了那么多人问压电陶瓷的材料系数输入问题，没有一个系统的回答。
这里给出一个完全版的。
它是CSI公司的S. Imaoka写的，很详细，看了都能明白。
如果看不明白，直接用后面的工具，输入IEEE标准的格式就能自动转换成Ansys模式的，方便而且保证正确。：）

附件是转换原理分析

[ 本帖最后由 北极熊甲 于 2007-11-6 13:17 编辑 ]

feasky

这个是转换工具，Excel格式的，方便，一目了然。

feasky

可以结合《耦合场分析指南》一起看，更容易明白了。

lypzsl

好人啊！

yushunlei1982

多谢！！

foxmail3771

好东西，谢谢你！

dqtznn

谢谢

xiaoheshang

谢谢，非常感谢资料分享

powerlong

感谢提供者，切实是好资料。

不过还有几个问题啊：
（1） 关于介电常数
   以上原理与工具中定义为：
  mp,prex,1,reps11
  mp,prex,1,reps11
  mp,prex,1,reps33
可我看了很多资料却这样定义：
  mp,prex,1,reps11
  mp,prex,1,reps33
  mp,prex,1,reps11
为什么有这样的不同？谁对谁错？
（2） 压电应力矩阵问题
    以上原理与工具中定义为：
    tbdata,3,e31
    tbdata,6,e31
   tbdata,9,e33
   tbdata,14,e15
   tbdata,16,e15
   而教材“压电学（孙慷，张福学主编）“却这样定义：
    tbdata,3,e31
    tbdata,6,e31
   tbdata,9,e33
   tbdata,11,e15
   tbdata,14,e15
为什么有这样的不同？谁对谁错？

请教各位！

feasky

powerlong wrote:
感谢提供者，切实是好资料。

不过还有几个问题啊：
（1） 关于介电常数
   以上原理与工具中定义为：
  mp,prex,1,reps11
  mp,prex,1,reps11
  mp,prex,1,reps33
可我看了很多资料却这样定义：
  mp,prex,1,reps11
  mp,prex,1,reps33
  mp,prex,1,reps11
为什么有这样的不同？谁对谁错？

说明前者是以Z轴为极化方向，而后者是Y轴

（2） 压电应力矩阵问题
    以上原理与工具中定义为：
    tbdata,3,e31
    tbdata,6,e31
   tbdata,9,e33
   tbdata,14,e15
   tbdata,16,e15
   而教材“压电学（孙慷，张福学主编）“却这样定义：
    tbdata,3,e31
    tbdata,6,e31
   tbdata,9,e33
   tbdata,11,e15
   tbdata,14,e15
为什么有这样的不同？谁对谁错？

显然后者给出的输入是错的，因为那样输入是IEEE标准，不是ANSYS采纳的
否则要这个工具干什么！需要说明的是IEEE与Ansys采用的e互为转置矩阵。

powerlong

感谢feasky的帮助。
但ANSYS中压电陶瓷参数输入问题学是困扰着我。最近做了几个小计算实验，对比结果更使我迷惑不解。具体计算实验这样：
实验1：在柱坐标，用旋转的方法建立压电陶瓷模型，极化方向定为Z轴；
实验2：在卡式坐标，先做面再旋转面成体而建立压电陶瓷模型，在局部坐标方法中将极化方向定为Z轴；
实验3：在卡式坐标，先做面再旋转面成体而建立压电陶瓷模型，将极化方向定为x轴；
结果是：实验2和实验3结果相同，而与实验1不相同。而我事先的认为是：实验1和实验2结果相同，而与实验3不相同。这是为什么？？困惑啊。实验1和实验2和实验3代码分别如下，请feasky及各位能抽空提示啊。
实验1：
!PZT material is PIC181,solid45 element is used for Horn.  05.5.12
!采用循环方法建立模型，

!05/5/10
/COM, YYP !/filname,Lt-pzt1!jobname
/units,si  ! SI Units
/title,LT-PZT
/NOPR
/PMETH,OFF

/PREP7
               
et,1,solid45 !  horn
et,2,solid5   !PZT
KEYOPT,2,1,3

! taper's material property(steel)
mp,ex,1,1.9e11
mp,ey,1,1.9e11
mp,ez,1,1.9e11
mp,prxy,1,0.3
mp,prxz,1,0.3
mp,pryz,1,0.3
mp,dens,1,7800


!!!!!!!!!!PZT material pameraters  ,polar at Z direction,3-axis.
MP,DENS,2,7850 !PIC-181

! stiffness
c11_E=1.523e11
c12_E=8.909e10
c13_E=8.547e10
c33_E=1.314e11
c44_E=2.830e10
c55_E=2.830e10
c66_E=3.161e10

! rel. Permittivity
eps0=8.854e-12
eps11_S=740
eps33_S=624
mp,perx,2,eps11_S
mp,pery,2,eps33_S    !eps11_S?
mp,perz,2,eps11_S     !eps33_S?

!Piezoelectric "e" matrix
e31=-4.5
e33=14.7
e15=11.0

!Stiffness matrix
TB,ANEL,2
TBDATA,1,c11_E,c12_E,c13_E
TBDATA,7,c11_E,c13_E
TBDATA,12,c33_E
TBDATA,16,c66_E
TBDATA,19,c44_E
TBDATA,21,c55_E

! Piezoelectric Matrix
TB,PIEZ,2
TBDATA,3,e31
TBDATA,6,e31
TBDATA,9,e33
TBDATA,14,e15
TBDATA,16,e15

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!System Modal
L_pzt0=2e-3
N=4
Do_pzt=16e-3
Ro_pzt=Do_pzt/2
Di_pzt=8e-3
Ri_pzt=Di_pzt/2   

L_insteel=6e-3
R_insteel=L_insteel/2  
L_insteel=L_pzt0*N  

L_bsteel=4.8e-3   
R_bsteel=Ro_pzt   

L_fsteel=L_bsteel
R_fsteel=R_bsteel

L_horn=40.49e-3
R_horn=R_fsteel

csys,1  !为了沿Z方向极化

!!!!!!!!!!!!!! Model Setup
!TRANSDUCER

winkel=15

!Back bolt
*do,yi,0,31
CYLIND,R_bsteel,0,0,L_bsteel,yi*winkel,(yi+1)*winkel
VATT,1,,1
VSEL,none
*enddo

!Inner bolt
*do,yi,0,31
CYLIND,R_insteel,0,L_bsteel,L_bsteel+L_insteel,yi*winkel,(yi+1)*winkel
VATT,1,,1
VSEL,none
*enddo

!PZT
! 1
*do,yi,0,31
CYLIND,Ri_pzt,Ro_pzt,L_bsteel,L_bsteel+L_pzt0,yi*winkel,(yi+1)*winkel
VATT,2,,2
VSEL,none
*enddo
! 2
*do,yi,0,31
CYLIND,Ri_pzt,Ro_pzt,L_bsteel,L_bsteel+L_pzt0*2,yi*winkel,(yi+1)*winkel
VATT,2,,2
VSEL,none
*enddo
! 3
*do,yi,0,31
CYLIND,Ri_pzt,Ro_pzt,L_bsteel,L_bsteel+L_pzt0*3,yi*winkel,(yi+1)*winkel
VATT,2,,2
VSEL,none
*enddo
! 4
*do,yi,0,31
CYLIND,Ri_pzt,Ro_pzt,L_bsteel,L_bsteel+L_pzt0*4,yi*winkel,(yi+1)*winkel
VATT,2,,2
VSEL,none
*enddo

!Horn
*do,yi,0,31
CYLIND,R_fsteel,0,L_bsteel+L_insteel,L_bsteel+L_insteel+L_fsteel+L_horn,yi*winkel,(yi+1)*winkel
VATT,1,,1
VSEL,none
*enddo



alls
!NUMMRG,kp
nummrg,all  
numcmp,all

/PSYMB,CS,1
/PSYMB,ESYS,1
ESIZE,,1e-6
MSHKEY,1  
MSHAPE,0,3D
vsel,all
VMESH,ALL

alls
AGLUE,ALL   
LGLUE,ALL   
nummrg,all
numcmp,all

/view,1,1,1
/REPLOt

FINISH  
/SOLU   
allsel,all

csys,1

! Solution
FINISH  
/SOL
ANTYPE,2
MODOPT,LANB,50  
EQSLV,SPAR  
MXPAND,50, , ,1
LUMPM,0
PSTRES,OFF  
MODOPT,LANB,50,10000,100000, ,OFF, ,2   
solve   
finish

/solu
csys,1

!1左
alls
nsel,s,loc,z,L_bsteel
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10

!1右
alls
nsel,s,loc,z,L_bsteel+L_pzt0
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
!2左
alls
nsel,s,loc,z,L_bsteel+L_pzt0
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
!2右
alls
nsel,s,loc,z,L_bsteel+L_pzt0*2
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
!3左
alls
nsel,s,loc,z,L_bsteel+L_pzt0*2
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
!3右
alls
nsel,s,loc,z,L_bsteel+L_pzt0*3
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
!4左
alls
nsel,s,loc,z,L_bsteel+L_pzt0*3
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
!4右
alls
nsel,s,loc,z,L_bsteel+L_pzt0*4
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10

alls
solve
finish
/POST1

实验2：
!PZT material is PIC181,solid45 element is used for Capilarry.  05.3.8   
!05/5/17
/COM,  !/filname,Lt-pzt1!jobname
/units,si ! SI Units
/title,Thermosonic Flip Chip Transducer Design   
/NOPR   
/PMETH,OFF  
/PREP7  

et,3,plane42  !为了旋转成体   
et,1,solid45 !  杆
et,2,solid5   !PZT  
KEYOPT,2,1,3

! taper's material property(steel)  
mp,ex,1,1.9e11  
mp,ey,1,1.9e11  
mp,ez,1,1.9e11  
mp,prxy,1,0.3   
mp,prxz,1,0.3   
mp,pryz,1,0.3   
mp,dens,1,7800  
!!!!!!!!!!PZT material pameraters  ,polar at Z direction,3-axis.
MP,DENS,2,7850 !PIC-181

! stiffness
c11_E=1.523e11
c12_E=8.909e10
c13_E=8.547e10
c33_E=1.314e11
c44_E=2.830e10
c55_E=2.830e10
c66_E=3.161e10

! rel. Permittivity
eps0=8.854e-12
eps11_S=740
eps33_S=624
mp,perx,2,eps11_S
mp,pery,2,eps33_S    !eps11_S?
mp,perz,2,eps11_S     !eps33_S?

!Piezoelectric "e" matrix
e31=-4.5
e33=14.7
e15=11.0

!Stiffness matrix
TB,ANEL,2
TBDATA,1,c11_E,c12_E,c13_E
TBDATA,7,c11_E,c13_E
TBDATA,12,c33_E
TBDATA,16,c66_E
TBDATA,19,c44_E
TBDATA,21,c55_E

! Piezoelectric Matrix
TB,PIEZ,2
TBDATA,3,e31
TBDATA,6,e31
TBDATA,9,e33
TBDATA,14,e15
TBDATA,16,e15

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!System Modal  大小长度参数
L_pzt0=2e-3
N=4
Do_pzt=16e-3
Ro_pzt=Do_pzt/2
Di_pzt=8e-3
Ri_pzt=Di_pzt/2   

L_insteel=6e-3
R_insteel=L_insteel/2  
L_insteel=L_pzt0*N  

L_bsteel=4.8e-3   
R_bsteel=Ro_pzt  

L_fsteel=L_bsteel
R_fsteel=R_bsteel

L_horn=40.49e-3
R_horn=R_fsteel

csys,0  
!!!!!!!!!!!!!! Model Setup  
!TRANSDUCER
k,1,0,0,0   
k,2,L_bsteel,0
k,3,L_bsteel,R_insteel,
k,4,L_bsteel,Ri_pzt,   
k,5,L_bsteel,Ro_pzt,   
!k,6,L_bsteel ,R_bsteel,
k,7,0,R_bsteel,0
k,8,L_bsteel+L_pzt0,Ro_pzt,
k,9,L_bsteel+L_pzt0*2,Ro_pzt,  
k,10,L_bsteel+L_pzt0*3,Ro_pzt,
k,11,L_bsteel+L_pzt0*4,Ro_pzt,
k,12,L_bsteel+L_pzt0*4 ,Ri_pzt,
k,13,L_bsteel+L_pzt0*3,Ri_pzt,
k,14,L_bsteel+L_pzt0*2 ,Ri_pzt,
k,15,L_bsteel+L_pzt0 ,Ri_pzt,  
k,16,L_bsteel+L_pzt0*4,R_insteel,  
k,17,L_bsteel+L_pzt0*4,0,
k,18,L_bsteel+L_pzt0*4+L_horn,0,   
k,19,L_bsteel+L_pzt0*4+L_horn,R_horn,  
!k,20,L_bsteel+L_pzt0*4,R_horn,

a,1,2,3,4,5,7  !Back bolt
a,2,3,16,17  !inner bolt

!pzt
a,4,5,8,15  
a,8,9,14,15
a,9,10,13,14
a,10,11,12,13   

!horn   
a,17,18,19,11   
aplot   
allsel  
aplot   
AGLUE,ALL   
LGLUE,ALL   
nummrg,all  
numcmp,all  
ALLS
CM,BASE,AREA

!   合并线，为面划分做准备
!for backbolt   
lsel,S,line,,4  
lsel,a,line,,3  
lsel,a,line,,2  
lccat,all   !create line26  
!for fontbolt   
lsel,s,line,,20
lsel,a,line,,25
lsel,a,line,,8
lccat,all   !create line27  

!     面划分  
!asel,s,loc,x,0,L_bsteel+L_pzt0*4+L_fsteel  
ALLS
aatt,1,,3   
MSHKEY,1
MSHAPE,0,2D
ESIZE,0.6e-3   
AMESH,base
eplot   
alls
ldelete,26,,,0  
ldelete,27,,,0  


alls
AGLUE,ALL
LGLUE,ALL
nummrg,all
numcmp,all


!   一次旋转成体
allsel  
TYPE,1  
ESIZE,,13              
VROTAT,ALL,,,,,,1,17,360,2 !
ACLEAR,BASE
eplot   
ALLSEL  
EPLOT  

!      修改PZT材料参数

!    pick up  
VSEL,S,VOLU,,3  
VSEL,A,VOLU,,4  
VSEL,A,VOLU,,5  
VSEL,A,VOLU,,6  
VSEL,A,VOLU,,10  
VSEL,a,VOLU,,11  
VSEL,A,VOLU,,12  
VSEL,A,VOLU,,13

eslv,s  
mat,2   
type,2  
emodif,all  
WPSTYL,DEF

/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
ALLSEL,ALL  
VPLOT   
EPLOT  

!电压加载(以Z方向上极化)
local,11,1,,,,90,90,  !YOZ平面（Z轴上）上加电压,采用旋转的圆柱坐标
vsel,s,loc,z,L_bsteel,L_bsteel+L_pzt0*4

!PZT
!1  
!   
alls
nsel,s,loc,z,L_bsteel
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
nsel,all
!   
alls
nsel,s,loc,z,L_bsteel+L_pzt0
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
nsel,all

!2  
alls
nsel,s,loc,z,L_bsteel+L_pzt0
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
nsel,all
!   
alls
nsel,s,loc,z,L_bsteel+L_pzt0*2
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
nsel,all
!3  
alls
nsel,s,loc,z,L_bsteel+L_pzt0*2
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
nsel,all
!   
alls
nsel,s,loc,z,L_bsteel+L_pzt0*3
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
nsel,all

!4  
alls
nsel,s,loc,z,L_bsteel+L_pzt0*3
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,10
nsel,all
!   
alls
nsel,s,loc,z,L_bsteel+L_pzt0*4
nsel,r,loc,x,Ri_pzt,Ro_pzt
D,all,VOLT,0
nsel,all  

alls

!  模态分析
FINISH  
/SOL
ANTYPE,2
MODOPT,LANB,50  
EQSLV,SPAR  
MXPAND,50, , ,1
LUMPM,0
PSTRES,OFF  
MODOPT,LANB,50,10000,100000, ,OFF, ,2   
solve   
finish

/POST1  

实验3：
!PZT material is PIC181,solid45 element is used for Capilarry.  05.3.8   
!XOY平面（X轴上）上加电压
!05/5/17
/COM, !/filname,Lt-pzt1!jobname
/units,si ! SI Units
/title,Thermosonic Flip Chip Transducer Design   
/NOPR   
/PMETH,OFF  
/PREP7  

et,3,plane42  !!为了旋转成体      
et,1,solid45 !  杆
et,2,solid5   !PZT  
KEYOPT,2,1,3

! taper's material property(steel)  
mp,ex,1,1.9e11  
mp,ey,1,1.9e11  
mp,ez,1,1.9e11  
mp,prxy,1,0.3   
mp,prxz,1,0.3   
mp,pryz,1,0.3   
mp,dens,1,7800  
!!!!!!!!!!PZT material pameraters  ,polar at Z direction,3-axis.
MP,DENS,2,7850 !PIC-181

! stiffness
c11_E=1.523e11
c12_E=8.909e10
c13_E=8.547e10
c33_E=1.314e11
c44_E=2.830e10
c55_E=2.830e10
c66_E=3.161e10

! rel. Permittivity
eps0=8.854e-12
eps11_S=740
eps33_S=624
mp,perx,2,eps11_S
mp,pery,2,eps33_S    !eps11_S?
mp,perz,2,eps11_S     !eps33_S?

!Piezoelectric "e" matrix
e31=-4.5
e33=14.7
e15=11.0

!Stiffness matrix
TB,ANEL,2
TBDATA,1,c11_E,c12_E,c13_E
TBDATA,7,c11_E,c13_E
TBDATA,12,c33_E
TBDATA,16,c66_E
TBDATA,19,c44_E
TBDATA,21,c55_E

! Piezoelectric Matrix
TB,PIEZ,2
TBDATA,3,e31
TBDATA,6,e31
TBDATA,9,e33
TBDATA,14,e15
TBDATA,16,e15

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!System Modal  大小长度参数
L_pzt0=2e-3
N=4
Do_pzt=16e-3
Ro_pzt=Do_pzt/2
Di_pzt=8e-3
Ri_pzt=Di_pzt/2   

L_insteel=6e-3
R_insteel=L_insteel/2  
L_insteel=L_pzt0*N  

L_bsteel=4.8e-3   
R_bsteel=Ro_pzt  

L_fsteel=L_bsteel
R_fsteel=R_bsteel

L_horn=40.49e-3
R_horn=R_fsteel

csys,0  
!!!!!!!!!!!!!! Model Setup  
!TRANSDUCER
k,1,0,0,0   
k,2,L_bsteel,0
k,3,L_bsteel,R_insteel,
k,4,L_bsteel,Ri_pzt,   
k,5,L_bsteel,Ro_pzt,   
!k,6,L_bsteel ,R_bsteel,
k,7,0,R_bsteel,0
k,8,L_bsteel+L_pzt0,Ro_pzt,
k,9,L_bsteel+L_pzt0*2,Ro_pzt,  
k,10,L_bsteel+L_pzt0*3,Ro_pzt,
k,11,L_bsteel+L_pzt0*4,Ro_pzt,
k,12,L_bsteel+L_pzt0*4 ,Ri_pzt,
k,13,L_bsteel+L_pzt0*3,Ri_pzt,
k,14,L_bsteel+L_pzt0*2 ,Ri_pzt,
k,15,L_bsteel+L_pzt0 ,Ri_pzt,  
k,16,L_bsteel+L_pzt0*4,R_insteel,  
k,17,L_bsteel+L_pzt0*4,0,
k,18,L_bsteel+L_pzt0*4+L_horn,0,   
k,19,L_bsteel+L_pzt0*4+L_horn,R_horn,  
!k,20,L_bsteel+L_pzt0*4,R_horn,
a,1,2,3,4,5,7  !Back bolt
a,2,3,16,17  !inner bolt
!pzt
a,4,5,8,15  
a,8,9,14,15
a,9,10,13,14
a,10,11,12,13   
!horn   
a,17,18,19,11  
aplot   
allsel  
aplot   
AGLUE,ALL   
LGLUE,ALL   
nummrg,all  
numcmp,all  
ALLS
CM,BASE,AREA
!   合并线，为面划分做准备
!for backbolt   
lsel,s,line,,2  
lsel,a,line,,4  
lsel,a,line,,3  
lccat,all   !create line26  
!for fontbolt   
lsel,s,line,,20
lsel,a,line,,25
lsel,a,line,,8
lccat,all   !create line27  

!     面划分  
!asel,s,loc,x,0,L_bsteel+L_pzt0*4+L_fsteel  
ALLS
aatt,1,,3   
MSHKEY,1
MSHAPE,0,2D
ESIZE,0.6e-3   
AMESH,base
eplot   
alls
ldelete,26,,,0  
ldelete,27,,,0  

!   一次旋转成体
allsel  
TYPE,1  
ESIZE,,13              
VROTAT,all,,,,,,1,17,360,2 !
ACLEAR,BASE
ACLEAR,1
eplot   
ALLSEL  
EPLOT   
!      修改PZT材料参数

!    pick up  
VSEL,S,VOLU,,3  
VSEL,A,VOLU,,4
VSEL,A,VOLU,,5  
VSEL,A,VOLU,,6  
VSEL,A,VOLU,,10  
VSEL,a,VOLU,,11  
VSEL,A,VOLU,,12  
VSEL,A,VOLU,,13

eslv,s  
mat,2   
type,2  
emodif,all  
WPSTYL,DEF
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
ALLSEL,ALL  
VPLOT   
EPLOT   

!电压加载(以x方向上极化)
vsel,s,loc,x,L_bsteel,L_bsteel+L_pzt0*4
!1
!左
alls
vsel,s,volu,,3
vsel,a,volu,,10
aslv,s
asel,r,loc,x,L_bsteel

nsla,s,1
D,all,VOLT,10
nsel,all
!右
alls
vsel,s,volu,,3
vsel,a,volu,,10
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0

nsla,s,1
D,all,VOLT,0
nsel,all

!2左
alls
vsel,s,volu,,4
vsel,a,volu,,11
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0

nsla,s,1
D,all,VOLT,0
nsel,all
!右
alls
vsel,s,volu,,4
vsel,a,volu,,11
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0*2

nsla,s,1
D,all,VOLT,10
nsel,all

!3左
alls
vsel,s,volu,,5
vsel,a,volu,,12
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0*2

nsla,s,1
D,all,VOLT,10
nsel,all
!右
alls
vsel,s,volu,,5
vsel,a,volu,,12
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0*3

nsla,s,1
D,all,VOLT,0
nsel,all

!4左
alls
vsel,s,volu,,6
vsel,a,volu,,13
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0*3

nsla,s,1
D,all,VOLT,0
nsel,all
!右
alls
vsel,s,volu,,6
vsel,a,volu,,13
aslv,s
asel,r,loc,x,L_bsteel+L_pzt0*4

nsla,s,1
D,all,VOLT,10
nsel,all

alls

!  模态分析
FINISH  
/SOL
ANTYPE,2
MODOPT,LANB,50  
EQSLV,SPAR  
MXPAND,50, , ,1
LUMPM,0
PSTRES,OFF  
MODOPT,LANB,50,10000,100000, ,OFF, ,2   
solve   
finish  
/POST1

lgdq006

re

zhangzc80

好文章坚决顶

daven

新来，学习下哈。还好不用积分，要不 又看不了啊

大大象

感谢lz~~~~~~~
谢谢~~~~~~~~~~

tykf

很好，不错，顶，谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢