orthotropic-plate

name	orthotropic-plate
description	Orthotropic plate theory — governing equation (Huber differential equation), bending stiffness (D_x, D_y, D_xy, D_1), Navier/Lévy solutions, composite laminates as equivalent orthotropic plates, grillage analogy, stiffened deck plates (bridges/ships), interaction formulas for biaxial loading, vibration of orthotropic plates, and AASHTO/AISC bridge deck design.
metadata	{"priority":7,"promptSignals":{"phrases":["orthotropic plate","Huber plate equation","orthotropic deck","anisotropic plate","composite plate theory","stiffened plate"],"minScore":3}}

Orthotropic Plate Theory — Complete Skill

Fundamentals

Definition and Assumptions

Orthotropic plate: Different bending and in-plane stiffness in two perpendicular directions (x and y) Special case of anisotropic material; principal axes aligned with plate axes Kirchhoff thin-plate assumptions apply: small deflection; linear elastic; plane sections remain plane

Comparison: Isotropic plate: E_x = E_y; ν_xy = ν_yx; one D Orthotropic plate: D_x ≠ D_y; D_1 (coupling); D_xy (twisting)

Sources of orthotropy: Material orthotropy: wood (grain direction), composites (0/90 laminate), rolled steel (weak minor) Structural orthotropy: stiffened plates (ribs in one direction only); corrugated decks; grillage of beams

Bending Stiffness Components

Plate stiffness parameters: D_x = E_x × h³ / (12(1 - ν_xy × ν_yx)) [bending stiffness in x-direction; h = plate thickness] D_y = E_y × h³ / (12(1 - ν_xy × ν_yx)) [bending stiffness in y-direction] D_1 = ν_yx × D_x = ν_xy × D_y [coupling stiffness; D₁ = ν D for isotropic] D_xy = G_xy × h³ / 12 [torsional/twisting stiffness; G_xy = shear modulus in xy-plane] Effective twisting stiffness: D_t = D_1 + 2 D_xy [appears in governing equation]

Reciprocal relation: ν_xy / E_x = ν_yx / E_y [Maxwell reciprocal; ν_xy × E_y = ν_yx × E_x]

For isotropic plate: D_x = D_y = D = Eh³/(12(1-ν²)); D_1 = νD; D_xy = (1-ν)D/2; D_t = D

Governing Equation

Huber Differential Equation

Governing equation for orthotropic plate (Huber, 1923): D_x × ∂⁴w/∂x⁴ + 2D_t × ∂⁴w/∂x²∂y² + D_y × ∂⁴w/∂y⁴ = q(x,y)

Where: w = transverse deflection [m] q = distributed transverse load [Pa = N/m²] D_t = D_1 + 2D_xy [effective twisting rigidity]

Isotropic reduction: D_x = D_y = D; D_t = D → D × ∇⁴w = q (biharmonic; standard isotropic)

Moments from curvature: M_x = -D_x × ∂²w/∂x² - D_1 × ∂²w/∂y² [moment per unit width in x-direction] M_y = -D_y × ∂²w/∂y² - D_1 × ∂²w/∂x² [moment per unit width in y-direction] M_xy = -2D_xy × ∂²w/∂x∂y [twisting moment]

Shear forces: Q_x = -∂/∂x(D_x × ∂²w/∂x² + D_t × ∂²w/∂y²) Q_y = -∂/∂y(D_y × ∂²w/∂y² + D_t × ∂²w/∂x²)

Analytical Solutions

Navier Double-Fourier Series (Simply Supported All Edges)

Boundary conditions: w = 0 and M = 0 on all edges (simply supported; SSSS) Plate dimensions: a × b (x: 0 to a; y: 0 to b)

Deflection series: w(x,y) = Σ_m Σ_n W_mn × sin(mπx/a) × sin(nπy/b)

Load series: q(x,y) = Σ_m Σ_n q_mn × sin(mπx/a) × sin(nπy/b) q_mn = (4/ab) × ∫∫ q(x,y) sin(mπx/a) sin(nπy/b) dx dy

Amplitude: W_mn = q_mn / [D_x(mπ/a)⁴ + 2D_t(mπ/a)²(nπ/b)² + D_y(nπ/b)⁴]

Uniform load q₀: q_mn = 16q₀/(π² × m × n) [odd m, n only; even terms zero] W = 16q₀/π⁶ × Σ_m Σ_n (1/mn) / [D_x(m/a)⁴ + 2D_t(m²n²)/(a²b²) + D_y(n/b)⁴] × sin(mπx/a)sin(nπy/b)

Maximum deflection (center of SSSS plate, uniform load): w_max ≈ q₀ × a⁴ / (D_eff × π⁴ × K) [K from Navier series; converges rapidly] For isotropic: w_max = q₀a⁴/(75.7D) for a = b (square plate)

Lévy Solution (Three Sides SS; One Side Arbitrary)

For plate with x-edges (x=0, x=a) simply supported; y-edges arbitrary: w(x,y) = Σ_m W_m(y) × sin(mπx/a)

W_m(y) satisfies 4th-order ODE in y: D_y × d⁴W_m/dy⁴ - 2D_t × (mπ/a)² × d²W_m/dy² + D_x × (mπ/a)⁴ × W_m = q_m

General solution: W_m = (particular) + A₁ cosh(α_m y) + A₂ sinh(α_m y) + A₃ cosh(β_m y) + A₄ sinh(β_m y) α_m, β_m from characteristic equation of ODE (two distinct real roots for orthotropic plates)

Apply remaining 4 BCs (y=0 and y=b edges; each contributes 2 conditions → 4 equations for A₁–A₄) Solve 4×4 system → closed-form

Equivalent Isotropic Plate (Huber Approximation)

For some orthotropic plates (wood, composites): Replace with equivalent isotropic using geometric mean: D_eff = √(D_x × D_y) [Huber approximation for D when D_t ≈ √(D_x × D_y)] Use standard isotropic plate tables with D_eff and equivalent aspect ratio

When valid: not valid for strongly orthotropic plates; use full Huber equation for D_x/D_y > 4

Structural Orthotropy — Stiffened Plates

Equivalent Orthotropic Stiffnesses

Stiffened plate (ribs in x-direction only; uniform spacing s_x): D_x = (E × I_x) / s_x [I_x = moment of inertia of rib+plate strip per unit width] D_y = E × t³ / (12(1-ν²)) [plate only; no ribs in y-direction; t = plate thickness] D_xy = G × (I_T_rib/s_x + t³/6) [twisting; I_T_rib = torsional inertia of rib cross-section] D_1 ≈ E × ν × t³ / (12(1-ν²)) [same as isotropic for coupling; rib contribution small]

Rib cross-section I_x (T-section): I_x = I_plate_strip + I_rib + A_rib × y_centroid² [parallel axis theorem; composite section] A_rib = h_rib × t_rib; I_rib = t_rib × h_rib³ / 12 + A_rib × (y_centroid - h_rib/2)²

Typical bridge orthotropic deck: 12 mm plate + 300 mm ribs at 300 mm spacing (trapezoidal or U-ribs) D_x/D_y ≈ 10–30 (highly orthotropic; ribs dominate x-direction)

Grillage Analogy

Replace orthotropic plate with grillage of beam members: Longitudinal beams (in x): EI_x per unit width = D_x; GJ_x per unit width = D_xy Transverse beams (in y): EI_y per unit width = D_y; GJ_y per unit width = D_xy Cross-rigidity: coupling D_1 handled approximately by Poisson ratio of each beam set

Grillage mesh: one beam per rib spacing; transverse beams at 1–2× rib spacing Used in bridge deck analysis (AASHTO); results compared to plate theory

Composite Laminates as Orthotropic Plates

Classical Laminate Theory (CLT)

Bending stiffness matrix [D]: D_ij = Σ_k (Q̄_ij)k × (z_k³ - z(k-1)³) / 3 [sum over plies k; z_k = distance from mid-plane to ply top] Q̄_ij = transformed reduced stiffness (function of E₁, E₂, G₁₂, ν₁₂, and ply angle θ)

Orthotropic laminate (0°/90° symmetric, e.g., [0/90]_s): D_11 corresponds to D_x; D_22 to D_y; D_66 to D_xy; D_12 to D_1 D_16 = D_26 = 0 for orthotropic symmetric laminate (no bending-twisting coupling)

Apparent engineering constants (for Huber equation): D_x ≈ D_11; D_y ≈ D_22; D_t = D_12 + 2D_66; D_xy = D_66

Quasi-isotropic laminate ([0/±45/90]_s): D_11 ≈ D_22 ≈ D_66 + D_12 → nearly isotropic; small D_16, D_26 terms → can use isotropic plate formulas

Vibration of Orthotropic Plates

Natural Frequencies (Free Vibration)

Frequency equation (SSSS boundary condition): ω_mn² = π⁴ × [D_x(m/a)⁴ + 2D_t(mn)²/(a²b²) + D_y(n/b)⁴] / (ρ_plate × h)

Where: ρ_plate = plate density [kg/m³]; h = plate thickness [m]; m, n = mode numbers

Fundamental frequency (m=n=1): f₁₁ = π²/(2π) × √[D_x/a⁴ + 2D_t/(a²b²) + D_y/b⁴] / (ρh) [Hz]

Isotropic check (D_x = D_y = D; D_t = D): f₁₁ = π/(2) × (1/a² + 1/b²) × √(D/(ρh)) [matches Leissa 1969 Table 4.1 for SSSS]

Buckling (in-plane compression): N_x_cr = π² × [D_x m²/a² + 2D_t n²/b² + D_y n⁴a²/(m²b⁴)] (n=1 governs for long plates) For uniaxial compression (N_x only; N_y = 0): N_x_cr = π² × min_m [D_x m²/a² + 2D_t/b² + D_y a²/(m²b⁴)] [per unit width, N/m]

Bridge Orthotropic Deck Design (AASHTO)

AASHTO LRFD Bridge Design

Orthotropic steel deck (OSD) design: AASHTO LRFD BDS Article 9.8: specific requirements for orthotropic deck Deck plate minimum: 14 mm (0.5625 in) per AASHTO 9.8.2.5 Rib thickness: ≥ 6 mm; rib height: 200–350 mm typical (closed U-ribs preferred) Rib spacing: 300–600 mm (typical 300 mm for truck wheel load distribution)

Fatigue design: Critical details: rib-to-deck weld (Category C or C'); rib-to-floor-beam connection AASHTO fatigue categories (Table 6.6.1.2.5-3): Cat. C: ΔF_n = 10 MPa (infinite life); Cat. C': ΔF_n = 7 MPa (infinite life) Rib-to-deck penetration weld: Cat. C' (typical); requires full penetration or enhanced detail

IMO and AISC application for ship decks: IACS CSR: structural analysis of container ship deck grillages AISC Design Guide 9: stiffened plate design for buildings (storage racks, bridge piers)

Standards and References

Standard	Scope
AASHTO LRFD Bridge Design Spec. Art. 9.8	Orthotropic steel deck bridge
Timoshenko & Woinowsky-Krieger "Theory of Plates and Shells"	Classical plate theory reference
Leissa (1969) NASA SP-160	Vibration of plates (tabulated frequencies)
Huber (1923) Czasopismo Techniczne	Original orthotropic plate equation
ISO 19902	Fixed steel offshore structures (plate analysis)
AISC Design Guide 9	Torsional analysis (includes plate components)

Output

Provide: plate geometry (a [m] × b [m] × h [mm]), source of orthotropy (material / structural / composite layup), bending stiffness values (D_x, D_y, D_t, D_xy [N·m]), stiffness ratio D_x/D_y (degree of orthotropy), loading (q₀ [kPa], concentrated P [kN], boundary conditions), solution method (Navier/Lévy/FEA), maximum deflection w_max [mm] and location, moment distribution (M_x_max [kN·m/m] and M_y_max [kN·m/m] at critical section), natural frequency f₁₁ [Hz] (fundamental) and mode shape, buckling load N_x_cr [kN/m] under in-plane compression (if applicable), equivalent orthotropic parameters (if composite: from CLT D-matrix), bridge application (rib geometry; fatigue category; AASHTO check), and applicable reference (Timoshenko & Woinowsky-Krieger, AASHTO LRFD 9.8, Leissa 1969).