Zaco intop1(p)/intop1(pabe.iomega*(w+eps))/(rho0*c0) "比声阻抗"
pext_1 at3_spatial(0,0,-1[m],pabe.p_t,'minc') "1 m 处的声压"
prms sqrt(0.5*pext_1*conj(pext_1)) "1 m 处的均方根压力"
TVR 20*log10(prms/Vrms/1[uPa/V]) "发射电压响应 (TVR)"
pext_Zeval at3_spatial(0,0,Zeval,pabe.p_t,'minc') "Zeval 处的声压"
Ifront 0.5*pext_Zeval*conj(pext_Zeval)/(rho0*c0) "Zeval 处的轴上强度"
Ptot -intop1(pabe.I_bndx*pabe.nx+pabe.I_bndy*pabe.ny+pabe.I_bndz*pabe.nz) "总辐射功率"
Iave Ptot/(4*pi*Zeval^2) "Zeval 处的单极源平均强度"
Di Ifront/Iave "强度方向性"
DI 10*log10(Di) "Tonpilz 型换能器的指向性指数"
k0 2*pi*freq/c0 "波数"
pzt_stress F_pre/intop2(1) "PZT 处的名义压应力"
a sqrt(intop3(1)/pi) "等效带凸缘活塞的半径"
DI_fl_pist 10*log10((k0*a)^2/(1-2*besselj(1,2*k0*a)/(2*k0*a))) "带凸缘活塞的指向性指数"