{"id":238721,"date":"2022-10-26T16:06:54","date_gmt":"2022-10-26T16:06:54","guid":{"rendered":"\/forum\/forums\/topic\/how-to-get-the-linearity-in-frequency-response-of-piezoelectric-transducer\/"},"modified":"2022-10-26T16:06:54","modified_gmt":"2022-10-26T16:06:54","slug":"how-to-get-the-linearity-in-frequency-response-of-piezoelectric-transducer","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/how-to-get-the-linearity-in-frequency-response-of-piezoelectric-transducer\/","title":{"rendered":"How to get the linearity in frequency response of piezoelectric transducer?"},"content":{"rendered":"

I am trying to model a piezoelectric sensor in ANSYS and I would like to get the frequency response of the sensor to calculate its usable frequency range and linearity. I first tried using harmonic analysis through the mode superposition method. I get a frequency response that includes the natural frequency and a linear range can be obtained through this plot. But I realized, this method is not accurate because the resistant element (CIRCU94) is ignored in the mode superposition method (please correct me if I am wrong). By changing the resistance value, the amplitude does not change. We need to have resistance to model charge dissipation, otherwise this open circuit voltage that does not properly change with excitation frequency won’t give me an accurate frequency range and linearity. Please see the following plot:<\/p>\n

<\/p>\n

 <\/p>\n

Then I used full harmonic analysis which I could see the effect of resistance in my voltage. But my expecation was that at least for low frequencies in the frequency response I see nonlinear behavior, and once the frequency increases, a more linear behavior in the voltage should be observed like this:<\/p>\n

\"Piezoelectric<\/p>\n

 <\/p>\n

But  what I get from the full harmonic analysis is this for low frequencies up to 500 Hz:<\/p>\n

<\/p>\n

 <\/p>\n

I am sure over frequencies like 20 Hz there must be almost a linear voltage amplitude. <\/p>\n

I attached my code for mode superposition and full harmonic. Can anyone check where I am making a mistake? or this is correct? The bottom of the PZT sensor is fixed, and excitation is from the top. The middle disc is PZT and I selected the top node at the end of the code to easily get results for this node. <\/p>\n

 <\/p>\n

I hope someone can help me soon because I need to report this simulation very soon. Thank you in advance, everyone! <\/p>\n

Mode superposition:<\/strong><\/p>\n

finish<\/span><\/div>\n
\/clear<\/span><\/div>\n
 <\/div>\n
\/PREP7  <\/span><\/div>\n
et,1,solid226,0 <\/span><\/div>\n
keyopt,1,1,1001<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
\/com Stiffness<\/span><\/div>\n
TB, ANEL, 1 , 1 , 0<\/span><\/div>\n
TBDATA, 1, 1.327E+5 , 0.866E+4 , 0.856E+4<\/span><\/div>\n
TBDATA, 7, 1.327E+5 , 0.856E+4<\/span><\/div>\n
TBDATA, 12, 1.192E+5<\/span><\/div>\n
TBDATA, 16, 2.12E+4<\/span><\/div>\n
TBDATA, 19, 2.12E+4<\/span><\/div>\n
TBDATA, 21, 2.99E+4<\/span><\/div>\n
 <\/div>\n
\/com Piezo matrix<\/span><\/div>\n
TB,PIEZ,1,,,1<\/span><\/div>\n
TBDATA, 3, -1.86e-10<\/span><\/div>\n
TBDATA, 6, -1.86e-10<\/span><\/div>\n
TBDATA, 9, 4.00e-10<\/span><\/div>\n
TBDATA, 11, 6.17e-10<\/span><\/div>\n
TBDATA, 13, 6.17e-10<\/span><\/div>\n
 <\/div>\n
\/com Permittivity<\/span><\/div>\n
tb,DPER,,,,1 ! Relative permittivity at const <\/span><\/div>\n
tbdata,1,1852,1852,1751<\/span><\/div>\n
TB,DPER,1,,,1852,1852,1751,1<\/span><\/div>\n
EMUNIT, EPZRO, 8.85E-15<\/span><\/div>\n
MP, DENS, 1 , 7.8e-9<\/span><\/div>\n
 <\/div>\n
 <\/div>\n
!==============================================<\/span><\/div>\n
et,2,solid185<\/span><\/div>\n
MP,EX,2,350e3 <\/span><\/div>\n
MP,PRXY,2,0.22  <\/span><\/div>\n
MP,DENS,2,3.9e-9 <\/span><\/div>\n
MP,BETD,2,1e-4<\/span><\/div>\n
!==============================================<\/span><\/div>\n
CYL4,0,0,2.5, , , ,1.5  <\/span><\/div>\n
BLC4,0,-1.25,7.1,2.5,1.5<\/span><\/div>\n
VSBV,       2,       1  <\/span><\/div>\n
CYL4,0,0,2.5, , , ,1.5  <\/span><\/div>\n
vglue,all<\/span><\/div>\n
k,21,0,0,1.5<\/span><\/div>\n
csys,4<\/span><\/div>\n
wpcsys,-1,0<\/span><\/div>\n
kwpave,21<\/span><\/div>\n
CYL4,0,0,1.5, , , ,2  <\/span><\/div>\n
k,,0,0,2<\/span><\/div>\n
kwpave,26<\/span><\/div>\n
CYL4,0,0,1.5, , , ,1.5<\/span><\/div>\n
csys,4<\/span><\/div>\n
wpcsys,-1,0<\/span><\/div>\n
save<\/span><\/div>\n
!==============================================<\/span><\/div>\n
Type,2<\/span><\/div>\n
Mat,2<\/span><\/div>\n
esize,0.3<\/span><\/div>\n
SMRT,1<\/span><\/div>\n
vsweep,2<\/span><\/div>\n
vsweep,4<\/span><\/div>\n
vsweep,3<\/span><\/div>\n
 <\/div>\n
Type,1<\/span><\/div>\n
Mat,1<\/span><\/div>\n
vsweep,1<\/span><\/div>\n
save<\/span><\/div>\n
!================================================<\/span><\/div>\n
DA,7,all,0<\/span><\/div>\n
DA,18,all,0<\/span><\/div>\n
 <\/div>\n
lsel,s,,,39,42<\/span><\/div>\n
DL,all,,ux,0<\/span><\/div>\n
DL,all,,uy,0<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
 <\/div>\n
asel,s,,,1<\/span><\/div>\n
lsla,s<\/span><\/div>\n
ksll,s<\/span><\/div>\n
nsla,s<\/span><\/div>\n
nsll,a<\/span><\/div>\n
nslk,a<\/span><\/div>\n
!D,all,volt,0<\/span><\/div>\n
cp,next,volt,all<\/span><\/div>\n
*get,n_ground,node,0,num,min<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
asel,s,,,2<\/span><\/div>\n
lsla,s<\/span><\/div>\n
ksll,s<\/span><\/div>\n
nsla,s<\/span><\/div>\n
nsll,a<\/span><\/div>\n
nslk,a<\/span><\/div>\n
cp,next,volt,all<\/span><\/div>\n
*get,n_supply,node,0,num,min<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
d,n_ground,volt,0<\/span><\/div>\n
!======= CONTACT ==========<\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – START <\/span><\/div>\n
CM,_NODECM,NODE <\/span><\/div>\n
CM,_ELEMCM,ELEM <\/span><\/div>\n
CM,_KPCM,KP <\/span><\/div>\n
CM,_LINECM,LINE <\/span><\/div>\n
CM,_AREACM,AREA <\/span><\/div>\n
CM,_VOLUCM,VOLU <\/span><\/div>\n
\/GSAV,cwz,gsav,,temp<\/span><\/div>\n
MP,MU,1,<\/span><\/div>\n
MAT,1   <\/span><\/div>\n
MP,EMIS,1,7.88860905221e-31 <\/span><\/div>\n
R,3 <\/span><\/div>\n
REAL,3  <\/span><\/div>\n
ET,3,170<\/span><\/div>\n
ET,4,174<\/span><\/div>\n
R,3,,,1.0,0.1,0,<\/span><\/div>\n
RMORE,,,1.0E20,0.0,1.0, <\/span><\/div>\n
RMORE,0.0,0,1.0,,1.0,0.5<\/span><\/div>\n
RMORE,0,1.0,1.0,0.0,,1.0<\/span><\/div>\n
RMORE,,,,,,1.0  <\/span><\/div>\n
KEYOPT,4,4,0<\/span><\/div>\n
KEYOPT,4,5,0<\/span><\/div>\n
KEYOPT,4,7,0<\/span><\/div>\n
KEYOPT,4,8,0<\/span><\/div>\n
KEYOPT,4,9,0<\/span><\/div>\n
KEYOPT,4,10,0   <\/span><\/div>\n
KEYOPT,4,11,0   <\/span><\/div>\n
KEYOPT,4,12,5   <\/span><\/div>\n
KEYOPT,4,14,0   <\/span><\/div>\n
KEYOPT,4,18,0   <\/span><\/div>\n
KEYOPT,4,2,0<\/span><\/div>\n
KEYOPT,3,5,0<\/span><\/div>\n
KEYOPT,4,1,5<\/span><\/div>\n
! Generate the target surface   <\/span><\/div>\n
ASEL,S,,,17  <\/span><\/div>\n
CM,_TARGET,AREA <\/span><\/div>\n
TYPE,3  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
ESLL,U  <\/span><\/div>\n
ESEL,U,ENAME,,188,189   <\/span><\/div>\n
NSLE,A,CT2  <\/span><\/div>\n
ESURF   <\/span><\/div>\n
CMSEL,S,_ELEMCM <\/span><\/div>\n
! Generate the contact surface  <\/span><\/div>\n
ASEL,S,,,1  <\/span><\/div>\n
CM,_CONTACT,AREA<\/span><\/div>\n
TYPE,4  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8\/2008) <\/span><\/div>\n
ESURF   <\/span><\/div>\n
ALLSEL  <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,3  <\/span><\/div>\n
ESEL,A,TYPE,,4  <\/span><\/div>\n
ESEL,R,REAL,,3  <\/span><\/div>\n
\/PSYMB,ESYS,1   <\/span><\/div>\n
\/PNUM,TYPE,1<\/span><\/div>\n
\/NUM,1  <\/span><\/div>\n
EPLOT   <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,3  <\/span><\/div>\n
ESEL,A,TYPE,,4  <\/span><\/div>\n
ESEL,R,REAL,,3  <\/span><\/div>\n
CMSEL,A,_NODECM <\/span><\/div>\n
CMDEL,_NODECM   <\/span><\/div>\n
CMSEL,A,_ELEMCM <\/span><\/div>\n
CMDEL,_ELEMCM   <\/span><\/div>\n
CMSEL,S,_KPCM   <\/span><\/div>\n
CMDEL,_KPCM <\/span><\/div>\n
CMSEL,S,_LINECM <\/span><\/div>\n
CMDEL,_LINECM   <\/span><\/div>\n
CMSEL,S,_AREACM <\/span><\/div>\n
CMDEL,_AREACM   <\/span><\/div>\n
CMSEL,S,_VOLUCM <\/span><\/div>\n
CMDEL,_VOLUCM   <\/span><\/div>\n
\/GRES,cwz,gsav  <\/span><\/div>\n
CMDEL,_TARGET   <\/span><\/div>\n
CMDEL,_CONTACT  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – END   <\/span><\/div>\n
APLOT   <\/span><\/div>\n
!*  <\/span><\/div>\n
!*  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – START <\/span><\/div>\n
CM,_NODECM,NODE <\/span><\/div>\n
CM,_ELEMCM,ELEM <\/span><\/div>\n
CM,_KPCM,KP <\/span><\/div>\n
CM,_LINECM,LINE <\/span><\/div>\n
CM,_AREACM,AREA <\/span><\/div>\n
CM,_VOLUCM,VOLU <\/span><\/div>\n
\/GSAV,cwz,gsav,,temp<\/span><\/div>\n
MP,MU,1,0   <\/span><\/div>\n
MAT,1   <\/span><\/div>\n
MP,EMIS,1,7.88860905221e-31 <\/span><\/div>\n
R,4 <\/span><\/div>\n
REAL,4  <\/span><\/div>\n
ET,5,170<\/span><\/div>\n
ET,6,174<\/span><\/div>\n
R,4,,,1.0,0.1,0,<\/span><\/div>\n
RMORE,,,1.0E20,0.0,1.0, <\/span><\/div>\n
RMORE,0.0,0,1.0,,1.0,0.5<\/span><\/div>\n
RMORE,0,1.0,1.0,0.0,,1.0<\/span><\/div>\n
RMORE,,,,,,1.0  <\/span><\/div>\n
KEYOPT,6,4,0<\/span><\/div>\n
KEYOPT,6,5,0<\/span><\/div>\n
KEYOPT,6,7,0<\/span><\/div>\n
KEYOPT,6,8,0<\/span><\/div>\n
KEYOPT,6,9,0<\/span><\/div>\n
KEYOPT,6,10,0   <\/span><\/div>\n
KEYOPT,6,11,0   <\/span><\/div>\n
KEYOPT,6,12,5   <\/span><\/div>\n
KEYOPT,6,14,0   <\/span><\/div>\n
KEYOPT,6,18,0   <\/span><\/div>\n
KEYOPT,6,2,0<\/span><\/div>\n
KEYOPT,5,5,0<\/span><\/div>\n
KEYOPT,6,1,5<\/span><\/div>\n
! Generate the target surface   <\/span><\/div>\n
ASEL,S,,,2  <\/span><\/div>\n
CM,_TARGET,AREA <\/span><\/div>\n
TYPE,5  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
ESLL,U  <\/span><\/div>\n
ESEL,U,ENAME,,188,189   <\/span><\/div>\n
NSLE,A,CT2  <\/span><\/div>\n
ESURF   <\/span><\/div>\n
CMSEL,S,_ELEMCM <\/span><\/div>\n
! Generate the contact surface  <\/span><\/div>\n
ASEL,S,,,11<\/span><\/div>\n
CM,_CONTACT,AREA<\/span><\/div>\n
TYPE,6  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8\/2008) <\/span><\/div>\n
ESURF   <\/span><\/div>\n
ALLSEL  <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,5  <\/span><\/div>\n
ESEL,A,TYPE,,6  <\/span><\/div>\n
ESEL,R,REAL,,4  <\/span><\/div>\n
\/PSYMB,ESYS,1   <\/span><\/div>\n
\/PNUM,TYPE,1<\/span><\/div>\n
\/NUM,1  <\/span><\/div>\n
EPLOT   <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,5  <\/span><\/div>\n
ESEL,A,TYPE,,6  <\/span><\/div>\n
ESEL,R,REAL,,4  <\/span><\/div>\n
CMSEL,A,_NODECM <\/span><\/div>\n
CMDEL,_NODECM   <\/span><\/div>\n
CMSEL,A,_ELEMCM <\/span><\/div>\n
CMDEL,_ELEMCM   <\/span><\/div>\n
CMSEL,S,_KPCM   <\/span><\/div>\n
CMDEL,_KPCM <\/span><\/div>\n
CMSEL,S,_LINECM <\/span><\/div>\n
CMDEL,_LINECM   <\/span><\/div>\n
CMSEL,S,_AREACM <\/span><\/div>\n
CMDEL,_AREACM   <\/span><\/div>\n
CMSEL,S,_VOLUCM <\/span><\/div>\n
CMDEL,_VOLUCM   <\/span><\/div>\n
\/GRES,cwz,gsav  <\/span><\/div>\n
CMDEL,_TARGET   <\/span><\/div>\n
CMDEL,_CONTACT  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – END   <\/span><\/div>\n
 <\/div>\n
!============================================================<\/span><\/div>\n
\/prep7<\/span><\/div>\n
ET,7,CIRCU94<\/span><\/div>\n
R,5,1, , , , ,<\/span><\/div>\n
RMORE, , , , , ,<\/span><\/div>\n
RMORE, , , ,<\/span><\/div>\n
type,7<\/span><\/div>\n
real,5<\/span><\/div>\n
E,n_ground,n_supply <\/span><\/div>\n
!============================================================<\/span><\/div>\n
\/solu<\/span><\/div>\n
ANTYPE,2 <\/span><\/div>\n
MODOPT,LANB,50 <\/span><\/div>\n
EQSLV,SPAR <\/span><\/div>\n
MODCONT,ON <\/span><\/div>\n
MXPAND,50, , ,0,,YES  <\/span><\/div>\n
SFA,12,1,PRES,1\/7.069<\/span><\/div>\n
LUMPM,0 <\/span><\/div>\n
PSTRES,0 <\/span><\/div>\n
MODOPT,LANB,50,0,4e10, ,OFF <\/span><\/div>\n
\/solu<\/span><\/div>\n
SOLVE <\/span><\/div>\n
finish <\/span><\/div>\n
!============================================================<\/span><\/div>\n
 <\/div>\n
!===========================================================<\/span><\/div>\n
\/solu<\/span><\/div>\n
ANTYPE,3 <\/span><\/div>\n
HROPT,MSUP, , ,0<\/span><\/div>\n
HROUT,ON<\/span><\/div>\n
LUMPM,0  <\/span><\/div>\n
HROPT,MSUP,all, ,0<\/span><\/div>\n
HROUT,ON,OFF,0 <\/span><\/div>\n
LVSCAL,1,1 <\/span><\/div>\n
 <\/div>\n
 <\/div>\n
HARFRQ,0,4e5,<\/span><\/div>\n
NSUBST,4e2,<\/span><\/div>\n
KBC,1    <\/span><\/div>\n
!BETAD,1e-6, <\/span><\/div>\n
DMPRAT,0.03 <\/span><\/div>\n
MDAMP,1, , , , , , ,  <\/span><\/div>\n
SOLVE  <\/span><\/div>\n
 <\/div>\n
save,pzt_harmonic,db<\/span><\/div>\n
finish<\/span><\/div>\n
save<\/span><\/div>\n
\/SOLU<\/span><\/div>\n
EXPASS,1 ! Harmonic Analysis Expansion Pass<\/span><\/div>\n
NUMEXP,ALL,0,4e5,0 ! Frequency range<\/span><\/div>\n
HREXP,ALL, <\/span><\/div>\n
SOLVE  <\/span><\/div>\n
 <\/div>\n
nsel,s,,,n_supply<\/span><\/div>\n

 <\/p>\n

Full harmonic analysis: <\/strong><\/p>\n

finish<\/span><\/div>\n
\/clear<\/span><\/div>\n
 <\/div>\n
\/PREP7  <\/span><\/div>\n
et,1,solid226,0 <\/span><\/div>\n
keyopt,1,1,1001<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
\/com Stiffness<\/span><\/div>\n
TB, ANEL, 1 , 1 , 0<\/span><\/div>\n
TBDATA, 1, 1.327E+5 , 0.866E+4 , 0.856E+4<\/span><\/div>\n
TBDATA, 7, 1.327E+5 , 0.856E+4<\/span><\/div>\n
TBDATA, 12, 1.192E+5<\/span><\/div>\n
TBDATA, 16, 2.12E+4<\/span><\/div>\n
TBDATA, 19, 2.12E+4<\/span><\/div>\n
TBDATA, 21, 2.99E+4<\/span><\/div>\n
 <\/div>\n
\/com Piezo matrix<\/span><\/div>\n
TB,PIEZ,1,,,1<\/span><\/div>\n
TBDATA, 3, -1.86e-10<\/span><\/div>\n
TBDATA, 6, -1.86e-10<\/span><\/div>\n
TBDATA, 9, 4.00e-10<\/span><\/div>\n
TBDATA, 11, 6.17e-10<\/span><\/div>\n
TBDATA, 13, 6.17e-10<\/span><\/div>\n
 <\/div>\n
\/com Permittivity<\/span><\/div>\n
tb,DPER,,,,1 ! Relative permittivity at const <\/span><\/div>\n
tbdata,1,1852,1852,1751<\/span><\/div>\n
TB,DPER,1,,,1852,1852,1751,1<\/span><\/div>\n
EMUNIT, EPZRO, 8.85E-15<\/span><\/div>\n
MP, DENS, 1 , 7.8e-9<\/span><\/div>\n
 <\/div>\n
 <\/div>\n
!==============================================<\/span><\/div>\n
et,2,solid185<\/span><\/div>\n
MP,EX,2,350e3 <\/span><\/div>\n
MP,PRXY,2,0.22  <\/span><\/div>\n
MP,DENS,2,3.9e-9 <\/span><\/div>\n
MP,BETD,2,1e-4<\/span><\/div>\n
!==============================================<\/span><\/div>\n
CYL4,0,0,2.5, , , ,1.5  <\/span><\/div>\n
BLC4,0,-1.25,7.1,2.5,1.5<\/span><\/div>\n
VSBV,       2,       1  <\/span><\/div>\n
CYL4,0,0,2.5, , , ,1.5  <\/span><\/div>\n
vglue,all<\/span><\/div>\n
k,21,0,0,1.5<\/span><\/div>\n
csys,4<\/span><\/div>\n
wpcsys,-1,0<\/span><\/div>\n
kwpave,21<\/span><\/div>\n
CYL4,0,0,1.5, , , ,2  <\/span><\/div>\n
k,,0,0,2<\/span><\/div>\n
kwpave,26<\/span><\/div>\n
CYL4,0,0,1.5, , , ,1.5<\/span><\/div>\n
csys,4<\/span><\/div>\n
wpcsys,-1,0<\/span><\/div>\n
save<\/span><\/div>\n
!==============================================<\/span><\/div>\n
Type,2<\/span><\/div>\n
Mat,2<\/span><\/div>\n
esize,0.3<\/span><\/div>\n
SMRT,1<\/span><\/div>\n
vsweep,2<\/span><\/div>\n
vsweep,4<\/span><\/div>\n
vsweep,3<\/span><\/div>\n
 <\/div>\n
Type,1<\/span><\/div>\n
Mat,1<\/span><\/div>\n
vsweep,1<\/span><\/div>\n
save<\/span><\/div>\n
!================================================<\/span><\/div>\n
DA,7,all,0<\/span><\/div>\n
DA,18,all,0<\/span><\/div>\n
 <\/div>\n
lsel,s,,,39,42<\/span><\/div>\n
DL,all,,ux,0<\/span><\/div>\n
DL,all,,uy,0<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
 <\/div>\n
asel,s,,,1<\/span><\/div>\n
lsla,s<\/span><\/div>\n
ksll,s<\/span><\/div>\n
nsla,s<\/span><\/div>\n
nsll,a<\/span><\/div>\n
nslk,a<\/span><\/div>\n
!D,all,volt,0<\/span><\/div>\n
cp,next,volt,all<\/span><\/div>\n
*get,n_ground,node,0,num,min<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
asel,s,,,2<\/span><\/div>\n
lsla,s<\/span><\/div>\n
ksll,s<\/span><\/div>\n
nsla,s<\/span><\/div>\n
nsll,a<\/span><\/div>\n
nslk,a<\/span><\/div>\n
cp,next,volt,all<\/span><\/div>\n
*get,n_supply,node,0,num,min<\/span><\/div>\n
allsel<\/span><\/div>\n
 <\/div>\n
d,n_ground,volt,0<\/span><\/div>\n
!======= CONTACT ==========<\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – START <\/span><\/div>\n
CM,_NODECM,NODE <\/span><\/div>\n
CM,_ELEMCM,ELEM <\/span><\/div>\n
CM,_KPCM,KP <\/span><\/div>\n
CM,_LINECM,LINE <\/span><\/div>\n
CM,_AREACM,AREA <\/span><\/div>\n
CM,_VOLUCM,VOLU <\/span><\/div>\n
\/GSAV,cwz,gsav,,temp<\/span><\/div>\n
MP,MU,1,<\/span><\/div>\n
MAT,1   <\/span><\/div>\n
MP,EMIS,1,7.88860905221e-31 <\/span><\/div>\n
R,3 <\/span><\/div>\n
REAL,3  <\/span><\/div>\n
ET,3,170<\/span><\/div>\n
ET,4,174<\/span><\/div>\n
R,3,,,1.0,0.1,0,<\/span><\/div>\n
RMORE,,,1.0E20,0.0,1.0, <\/span><\/div>\n
RMORE,0.0,0,1.0,,1.0,0.5<\/span><\/div>\n
RMORE,0,1.0,1.0,0.0,,1.0<\/span><\/div>\n
RMORE,,,,,,1.0  <\/span><\/div>\n
KEYOPT,4,4,0<\/span><\/div>\n
KEYOPT,4,5,0<\/span><\/div>\n
KEYOPT,4,7,0<\/span><\/div>\n
KEYOPT,4,8,0<\/span><\/div>\n
KEYOPT,4,9,0<\/span><\/div>\n
KEYOPT,4,10,0   <\/span><\/div>\n
KEYOPT,4,11,0   <\/span><\/div>\n
KEYOPT,4,12,5   <\/span><\/div>\n
KEYOPT,4,14,0   <\/span><\/div>\n
KEYOPT,4,18,0   <\/span><\/div>\n
KEYOPT,4,2,0<\/span><\/div>\n
KEYOPT,3,5,0<\/span><\/div>\n
KEYOPT,4,1,5<\/span><\/div>\n
! Generate the target surface   <\/span><\/div>\n
ASEL,S,,,17  <\/span><\/div>\n
CM,_TARGET,AREA <\/span><\/div>\n
TYPE,3  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
ESLL,U  <\/span><\/div>\n
ESEL,U,ENAME,,188,189   <\/span><\/div>\n
NSLE,A,CT2  <\/span><\/div>\n
ESURF   <\/span><\/div>\n
CMSEL,S,_ELEMCM <\/span><\/div>\n
! Generate the contact surface  <\/span><\/div>\n
ASEL,S,,,1  <\/span><\/div>\n
CM,_CONTACT,AREA<\/span><\/div>\n
TYPE,4  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8\/2008) <\/span><\/div>\n
ESURF   <\/span><\/div>\n
ALLSEL  <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,3  <\/span><\/div>\n
ESEL,A,TYPE,,4  <\/span><\/div>\n
ESEL,R,REAL,,3  <\/span><\/div>\n
\/PSYMB,ESYS,1   <\/span><\/div>\n
\/PNUM,TYPE,1<\/span><\/div>\n
\/NUM,1  <\/span><\/div>\n
EPLOT   <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,3  <\/span><\/div>\n
ESEL,A,TYPE,,4  <\/span><\/div>\n
ESEL,R,REAL,,3  <\/span><\/div>\n
CMSEL,A,_NODECM <\/span><\/div>\n
CMDEL,_NODECM   <\/span><\/div>\n
CMSEL,A,_ELEMCM <\/span><\/div>\n
CMDEL,_ELEMCM   <\/span><\/div>\n
CMSEL,S,_KPCM   <\/span><\/div>\n
CMDEL,_KPCM <\/span><\/div>\n
CMSEL,S,_LINECM <\/span><\/div>\n
CMDEL,_LINECM   <\/span><\/div>\n
CMSEL,S,_AREACM <\/span><\/div>\n
CMDEL,_AREACM   <\/span><\/div>\n
CMSEL,S,_VOLUCM <\/span><\/div>\n
CMDEL,_VOLUCM   <\/span><\/div>\n
\/GRES,cwz,gsav  <\/span><\/div>\n
CMDEL,_TARGET   <\/span><\/div>\n
CMDEL,_CONTACT  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – END   <\/span><\/div>\n
APLOT   <\/span><\/div>\n
!*  <\/span><\/div>\n
!*  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – START <\/span><\/div>\n
CM,_NODECM,NODE <\/span><\/div>\n
CM,_ELEMCM,ELEM <\/span><\/div>\n
CM,_KPCM,KP <\/span><\/div>\n
CM,_LINECM,LINE <\/span><\/div>\n
CM,_AREACM,AREA <\/span><\/div>\n
CM,_VOLUCM,VOLU <\/span><\/div>\n
\/GSAV,cwz,gsav,,temp<\/span><\/div>\n
MP,MU,1,0   <\/span><\/div>\n
MAT,1   <\/span><\/div>\n
MP,EMIS,1,7.88860905221e-31 <\/span><\/div>\n
R,4 <\/span><\/div>\n
REAL,4  <\/span><\/div>\n
ET,5,170<\/span><\/div>\n
ET,6,174<\/span><\/div>\n
R,4,,,1.0,0.1,0,<\/span><\/div>\n
RMORE,,,1.0E20,0.0,1.0, <\/span><\/div>\n
RMORE,0.0,0,1.0,,1.0,0.5<\/span><\/div>\n
RMORE,0,1.0,1.0,0.0,,1.0<\/span><\/div>\n
RMORE,,,,,,1.0  <\/span><\/div>\n
KEYOPT,6,4,0<\/span><\/div>\n
KEYOPT,6,5,0<\/span><\/div>\n
KEYOPT,6,7,0<\/span><\/div>\n
KEYOPT,6,8,0<\/span><\/div>\n
KEYOPT,6,9,0<\/span><\/div>\n
KEYOPT,6,10,0   <\/span><\/div>\n
KEYOPT,6,11,0   <\/span><\/div>\n
KEYOPT,6,12,5   <\/span><\/div>\n
KEYOPT,6,14,0   <\/span><\/div>\n
KEYOPT,6,18,0   <\/span><\/div>\n
KEYOPT,6,2,0<\/span><\/div>\n
KEYOPT,5,5,0<\/span><\/div>\n
KEYOPT,6,1,5<\/span><\/div>\n
! Generate the target surface   <\/span><\/div>\n
ASEL,S,,,2  <\/span><\/div>\n
CM,_TARGET,AREA <\/span><\/div>\n
TYPE,5  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
ESLL,U  <\/span><\/div>\n
ESEL,U,ENAME,,188,189   <\/span><\/div>\n
NSLE,A,CT2  <\/span><\/div>\n
ESURF   <\/span><\/div>\n
CMSEL,S,_ELEMCM <\/span><\/div>\n
! Generate the contact surface  <\/span><\/div>\n
ASEL,S,,,11<\/span><\/div>\n
CM,_CONTACT,AREA<\/span><\/div>\n
TYPE,6  <\/span><\/div>\n
NSLA,S,1<\/span><\/div>\n
ESLN,S,0<\/span><\/div>\n
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8\/2008) <\/span><\/div>\n
ESURF   <\/span><\/div>\n
ALLSEL  <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,5  <\/span><\/div>\n
ESEL,A,TYPE,,6  <\/span><\/div>\n
ESEL,R,REAL,,4  <\/span><\/div>\n
\/PSYMB,ESYS,1   <\/span><\/div>\n
\/PNUM,TYPE,1<\/span><\/div>\n
\/NUM,1  <\/span><\/div>\n
EPLOT   <\/span><\/div>\n
ESEL,ALL<\/span><\/div>\n
ESEL,S,TYPE,,5  <\/span><\/div>\n
ESEL,A,TYPE,,6  <\/span><\/div>\n
ESEL,R,REAL,,4  <\/span><\/div>\n
CMSEL,A,_NODECM <\/span><\/div>\n
CMDEL,_NODECM   <\/span><\/div>\n
CMSEL,A,_ELEMCM <\/span><\/div>\n
CMDEL,_ELEMCM   <\/span><\/div>\n
CMSEL,S,_KPCM   <\/span><\/div>\n
CMDEL,_KPCM <\/span><\/div>\n
CMSEL,S,_LINECM <\/span><\/div>\n
CMDEL,_LINECM   <\/span><\/div>\n
CMSEL,S,_AREACM <\/span><\/div>\n
CMDEL,_AREACM   <\/span><\/div>\n
CMSEL,S,_VOLUCM <\/span><\/div>\n
CMDEL,_VOLUCM   <\/span><\/div>\n
\/GRES,cwz,gsav  <\/span><\/div>\n
CMDEL,_TARGET   <\/span><\/div>\n
CMDEL,_CONTACT  <\/span><\/div>\n
\/COM, CONTACT PAIR CREATION – END   <\/span><\/div>\n
 <\/div>\n
!============================================================<\/span><\/div>\n
\/prep7<\/span><\/div>\n
ET,7,CIRCU94<\/span><\/div>\n
R,5,1e6, , , , ,<\/span><\/div>\n
RMORE, , , , , ,<\/span><\/div>\n
RMORE, , , ,<\/span><\/div>\n
type,7<\/span><\/div>\n
real,5<\/span><\/div>\n
E,n_ground,n_supply <\/span><\/div>\n
!============================================================<\/span><\/div>\n
\/solu<\/span><\/div>\n
ANTYPE,2 <\/span><\/div>\n
MODOPT,LANB,50 <\/span><\/div>\n
EQSLV,SPAR <\/span><\/div>\n
MODCONT,ON <\/span><\/div>\n
MXPAND,50, , ,0,,YES  <\/span><\/div>\n
!!!!!!!!!!!!!!!!!!!!!!!!!!!!SFA,12,1,PRES,1\/7.069<\/span><\/div>\n
LUMPM,0 <\/span><\/div>\n
PSTRES,0 <\/span><\/div>\n
MODOPT,LANB,50,0,4e10, ,OFF <\/span><\/div>\n
\/solu<\/span><\/div>\n
SOLVE <\/span><\/div>\n
finish <\/span><\/div>\n
!============================================================<\/span><\/div>\n
\/solu<\/span><\/div>\n
asel,s,,,12<\/span><\/div>\n
nsla,s<\/span><\/div>\n
nsel,r,loc,x,-0.25,0.25<\/span><\/div>\n
nsel,r,loc,y,-0.25,0.25<\/span><\/div>\n
*get,n_force,node,0,num,min<\/span><\/div>\n
allsel<\/span><\/div>\n
f,n_force,fz,-1<\/span><\/div>\n
!===========================================================<\/span><\/div>\n
\/SOL <\/span><\/div>\n
ANTYPE,3<\/span><\/div>\n
HROPT,FULL  <\/span><\/div>\n
HROUT,ON<\/span><\/div>\n
LUMPM,0  <\/span><\/div>\n
EQSLV, ,0,  <\/span><\/div>\n
PSTRES,0<\/span><\/div>\n
HARFRQ,0,500,  <\/span><\/div>\n
NSUBST,100,  <\/span><\/div>\n
KBC,0    <\/span><\/div>\n
DMPRAT,0.1<\/span><\/div>\n
MDAMP,1, , , , , , ,  <\/span><\/div>\n
solve<\/span><\/div>\n
 <\/div>\n
vsel,s,,,1<\/span><\/div>\n
eslv,s<\/span><\/div>\n
nsle,s<\/span><\/div>\n

 <\/p>\n","protected":false},"template":"","class_list":["post-238721","topic","type-topic","status-closed","hentry","topic-tag-Modal_Analysis-1","topic-tag-frequency-response","topic-tag-pzt-1"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["180812","50680"],"_bbp_author_ip":["23.56.168.190"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_btv_view_count":["1279"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["238721"],"_bbp_forum_id":["27793"],"_bbp_engagement":["50680","180812"],"_bbp_voice_count":["2"],"_bbp_reply_count":["1"],"_bbp_last_reply_id":["240147"],"_bbp_last_active_id":["240147"],"_bbp_last_active_time":["2022-11-08 16:34:21"]},"test":"ali_safian"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/238721","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/238721\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=238721"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}