Gopalkrishna Advanced Rural Research Foundation (GARRF)
Advanced Virtual EMI & SHM Research Laboratory
Step 1: Interfacial Adhesive Shear-Lag & Viscoelastic Boundary Dissipation
1. Specimen & Transducer Configuration
Specimen Geometry (Aspect Ratio Constrained)
🧪 Interfacial Adhesive Shear-Lag Parameters
• PZT-5H Technical Specifications
Dimensions: 10 × 10 × 0.3 mmd31 = -274 × 10-12 m/V, d33 = 593 × 10-12 m/V
Y11E = 62 GPa, Y33E = 50 GPa, ρ = 7750 kg/m³
ε33T / ε0 = 3400, Loss factor η = 0.015
• Crack Geometric Progression
| State | Length (a) | Depth (d) | Width | ΔK Loss |
|---|---|---|---|---|
| Pristine | 0.0 mm | 0.0 mm (0%) | 0.0 mm | 0.0 % |
| Crack 1 | 2.0 mm | 10% t_s | 0.1 mm | 4.0 % |
| Crack 2 | 5.0 mm | 30% t_s | 0.2 mm | 10.0 % |
| Crack 3 | 10.0 mm | 50% t_s | 0.4 mm | 18.0 % |
| Crack 4 | 18.0 mm | 80% t_s | 0.8 mm | 30.0 % |
2. Spectral Measurement Sequence (Active When Locked)
Screen 1: Specimen Support, Bond Layer & Dynamic Crack Kinematics
Setup Unlocked • Configure parameters & click 'Lock Parameters'
Screen 2: Real Conductance Spectrum G(ω) [10 kHz – 1 MHz]
Pristine
Crack 1
Crack 2
Crack 3
Crack 4
Screen 3: Imaginary Susceptance Spectrum B(ω) [10 kHz – 1 MHz]
Pristine
Crack 1
Crack 2
Crack 3
Crack 4
1. Interfacial Adhesive Shear-Lag & Viscoelastic Boundary Mechanics
In physical structural health monitoring installations, the piezoceramic transducer is bonded to the host structure through an adhesive layer (typically cyanoacrylate or high-shear epoxy) of finite thickness tb. The adhesive acts as a viscoelastic boundary interface that transmits dynamic stress into the substrate through interfacial shear stress τ(x):
τ(x) = ( G̅b / tb ) • [ usubstrate(x) − uPZT(x) ]
Complex Shear Modulus: G̅b = Gb ( 1 + j ηb )
Complex Shear-Lag Parameter: Γ = √( G̅b / ( tb • Y̅11E • h ) )
Here, Gb is the storage shear modulus (0.5 to 3.0 GPa), ηb is the viscoelastic loss factor (0.02 to 0.30), Y̅11E = Y11E(1 + jη) is the complex Young's modulus of PZT, and h is transducer thickness.
Modified Effective Structural Impedance Z̅s,eff
Due to interfacial shear lag, the mechanical drive-point impedance Zs presented by the host substrate is filtered before reaching the electro-active ceramic core. The effective impedance Z̅s,eff coupled to the PZT element is governed by:
Z̅s,eff = ( G̅b / ( Γ • tb ) ) • tanh( Γ • l ) • [ Zs / ( Zs + Zbond ) ]
Bond Transmission Efficiency Factor: αbond = tanh( Γ • l ) / ( Γ • l )
- Thick or Soft Adhesive ($t_b \uparrow$, $G_b \downarrow$): A thick adhesive layer increases mechanical compliance between PZT and substrate. This attenuates high-frequency shear stress transmission, reducing resonance peak heights and broadening the half-power bandwidth (reducing effective system $Q$-factor).
- Viscoelastic Dissipation ($\eta_b \uparrow$): Higher adhesive loss factors introduce damping into the boundary impedance, elevating the baseline conductance floor without altering structural modal frequencies.
2. Structural Geometry & Analytical Modal Dynamics
The natural frequencies and mechanical drive-point impedances of the host structure are calculated strictly from analytical elastodynamics based on the selectable dimensions:
- Length ($L$): 50 mm to 1000 mm (metals), 50 mm to 500 mm (concrete).
- Width ($W$): Constrained between 1/5 L and 1/2 L (0.20 L ≤ W ≤ 0.50 L), ensuring valid beam and plate aspect ratios.
- Thickness ($t_s$): 2 mm to 6 mm for 1D/2D thin formulations ($t_s \ll L, W$); 10% L to 100% L for 3D volumetric solids.
1D Euler-Bernoulli Flexural Base:
fbase,1D = ( βn2 / ( 2π L2 ) ) • √( E • I / ( ρ • A ) ) = ( βn2 / ( 2π L2 ) ) • √( E / ρ ) • ( ts / √12 )
2D Kirchhoff-Love Thin Plate Rigidity:
D = ( E • ts3 ) / [ 12 • ( 1 − ν2 ) ],
fbase,2D = ( π / ( 2 L2 ) ) • √( D / ( ρ • ts ) )
3D Volumetric Continuum Base:
cs = √( E / [ 2 • ρ • ( 1 + ν ) ] ),
fbase,3D = cs / ( 2 • √( L2 + W2 + ts2 ) )
3. Exact Quantitative PZT Placement Coordinates
Transducer locations are defined using Cartesian coordinates (x, y, z) rather than qualitative descriptors:
| Configuration | Location Designation | Exact Mathematical Coordinates (x, y, z) | Modal Sensitivity & Strain Energy Coupling |
|---|---|---|---|
| 1D Beam | Near Fixed Clamped Root | x = 0.15 L, y = 0.50 W, z = ts | Maximum dynamic bending curvature; maximum modal sensitivity. |
| 1D Beam | Mid-Span Center | x = 0.50 L, y = 0.50 W, z = ts | Maximum displacement for odd modes; zero curvature for even modes. |
| 1D Beam | Near Free Tip End | x = 0.85 L, y = 0.50 W, z = ts | Zero bending moment boundary; minimal dynamic strain transmission. |
| 2D Plate | Plate Geometric Center | x = 0.50 L, y = 0.50 W, z = ts | Symmetric biaxial radial expansion coupled with (1,1) plate modes. |
| 2D Plate | Plate Edge Midpoint | x = 0.35 L, y = 0.20 W, z = ts | Mixed boundary reflection; asymmetric modal spectrum. |
| 2D Plate | Plate Corner Region | x = 0.20 L, y = 0.20 W, z = ts | Near boundary constraint; localized high-frequency edge diffraction. |
| 3D Volumetric | Top Surface Center | x = 0.50 L, y = 0.50 W, z = ts | Surface normal dilatation coupled with planar Rayleigh wave modes. |
| 3D Volumetric | Top Surface Offset Edge | x = 0.25 L, y = 0.25 W, z = ts | Corner and quarter-span edge reflections; asymmetric 3D wave scattering. |
| 3D Volumetric | Centrally Embedded Aggregate | x = 0.50 L, y = 0.50 W, z = 0.50 ts | Full 6-face triaxial solid confinement (T1, T2, T3 ≠ 0). |
4. Comparative Admittance Formulations: 1D, 2D, and 3D
A. 1D Formulation [Liang, Sun, & Rogers, 1993, 1994]
Assumes pure uniaxial extension along length l (T2 = T3 = 0):
Y(ω) = jω ( w l / h ) • [ ε̅33T − d312 Y̅11E + d312 Y̅11E • ( Za,1D / ( Z̅s,eff + Za,1D ) ) • ( tan(κl) / (κl) ) ]
B. 2D Formulation [Zhou et al., 1996; Bhalla & Soh, 2004]
Biaxial plane stress (T3 = 0) with Poisson's ratio ν cross-coupling and adhesive shear-lag:
Y(ω) = 2 jω ( w l / h ) • [ ε̅33T − ( 2 d312 Y̅E / ( 1 − ν ) ) + ( 2 d312 Y̅E / ( 1 − ν ) ) • ( Za,2D / ( Z̅s,eff + Za,2D ) ) • ( tan(κl) / (κl) ) ]
C. 3D Directional Sum Impedance [Annamdas & Soh, 2007, 2008]
Captures full 3D stress-strain tensors across all six faces without planar assumptions:
Y(ω) = jω ∫∫A D • n dA = G(ω) + j B(ω)
Total Boundary Impedance: Zs,3D = Zs,xx + Zs,yy + Zs,zz + Zshear
Directional Wavenumbers: κx = ω √( ρ / C11E ),
κy = ω √( ρ / C22E ),
κz = ω √( ρ / C33E )
High-Frequency Thickness Mode Proof: When impedance sweeps reach 600–900 kHz, physical measurements display a pronounced thickness dilatation resonance peak. Because 1D and 2D models enforce T3 = 0, their wave equations lack the ∂2uz/∂z2 eigenvalue, rendering them mathematically incapable of predicting this mode. The 3D Annamdas-Soh formulation incorporates κz = ω √(ρ / C33E), predicting the thickness peak within 2% to 4% of measured experimental data.
5. Peer-Reviewed Scholarly References
- Liang, C., Sun, F. P., & Rogers, C. A. (1993). "Coupled Electro-Mechanical Analysis of Adaptive Material Systems Determination of the Actuator Power Consumption and System Energy Transfer." AIAA Conference Proceedings, pp. 3573–3582.
- Liang, C., Sun, F. P., & Rogers, C. A. (1994). "Coupled Electro-Mechanical Analysis of Adaptive Material Systems—Determination of the Actuator Power Consumption and System Energy Transfer." Journal of Intelligent Material Systems and Structures, 5(1), pp. 12–20.
- Zhou, S. W., Liang, C., & Rogers, C. A. (1996). "An Dynamic Model of Piezoelectric Actuator-Driven Thin Plates." Journal of Intelligent Material Systems and Structures, 7(4), pp. 434–443.
- Bhalla, S., & Soh, C. K. (2004). "Electromechanical Impedance Modeling for Adhesively Bonded Piezo-Transducers." Journal of Aerospace Engineering (ASCE), 17(4), pp. 154–165.
- Annamdas, V. G. M., & Soh, C. K. (2007). "Three-dimensional electromechanical impedance model for piezoceramic transducer–structure interaction, Part I: Formulation of directional sum impedance." Smart Materials and Structures, 16(5), pp. 1801–1813.
- Annamdas, V. G. M., & Soh, C. K. (2007). "Three-dimensional electromechanical impedance model for piezoceramic transducer–structure interaction, Part II: Numerical implementation and validation." Smart Materials and Structures, 16(5), pp. 1814–1821.
- Annamdas, V. G. M., & Soh, C. K. (2008). "Three-dimensional electromechanical impedance model for multiple piezoceramic transducers-structure interaction." Smart Materials and Structures, 17(5), p. 055009.