| name | structural-analysis-workflow |
| description | FEA workflow for pump casings and impellers under fluid loads |
| category | thinking |
| domain | structural |
| complexity | advanced |
| dependencies | [] |
Structural Analysis Workflow for Pump Components
This skill provides a systematic approach to performing Finite Element Analysis (FEA) on pump components, specifically casings and impellers under fluid loads.
FEA Workflow
1. Geometry and Material Definition
Geometry Preparation:
- Import or create CAD geometry of pump casing or impeller
- Simplify geometry by removing non-structural features (small fillets, bolt holes, cosmetic features)
- Identify critical regions: high-stress areas, pressure boundaries, mounting interfaces
- Define coordinate system aligned with pump axis (Z-axis typically along shaft)
Material Definition:
- Assign material properties to all components
- Define isotropic or anisotropic properties as needed
- Input temperature-dependent properties if thermal effects exist
- Consider material orientation for cast or forged components
Critical Parameters:
- Young's Modulus (E)
- Poisson's Ratio (ν)
- Yield Strength (σ_y)
- Ultimate Tensile Strength (σ_u)
- Density (ρ)
- Thermal expansion coefficient (α)
2. Load Application
Pressure Loads:
- Apply internal pressure to casing wetted surfaces
- Use static pressure at design point or maximum operating pressure
- For impellers, apply pressure distribution from CFD analysis
- Consider pressure pulsations for fatigue analysis
Centrifugal Forces (Impellers):
- Apply rotational velocity to impeller body
- Centrifugal force: F_c = mω²r
- Ensure correct direction and magnitude
- Consider overspeed conditions (typically 120% of rated speed)
Hydraulic Forces:
- Radial thrust on impeller
- Axial thrust from pressure difference
- Tangential forces from fluid momentum change
- Unbalanced forces at off-design conditions
Thermal Loads:
- Temperature distribution from thermal analysis
- Differential expansion between components
- Thermal gradients during startup/shutdown
3. Boundary Conditions (Constraints)
Casing Constraints:
- Fixed supports at mounting feet
- Symmetry planes if applicable
- Pipe connections (nozzles) - typically allow thermal expansion
- Anchor points: fully constrained in all directions
Impeller Constraints:
- Hub constrained to shaft (cylindrical constraint)
- Contact between impeller and shaft (may use bonded or frictional contact)
- Symmetry conditions if analyzing only a sector
Connection Types:
- Bonded: No relative motion (welded joints)
- Frictional: Coulomb friction model (shrink fits)
- No separation: Contact that can slide but not separate
- Frictionless: Free sliding contact
4. Mesh Generation
Element Types:
- Tetrahedral (Tet10): General purpose, automatic meshing
- Hexahedral (Hex20): Higher accuracy, structured regions
- Wedge (Prism): Transition between tet and hex meshes
- Shell elements: Thin-walled sections of casings
Mesh Quality Criteria:
- Element aspect ratio < 5:1 (ideally < 3:1)
- Skewness < 0.85
- Jacobian ratio > 0.6
- Minimum of 3 elements through thickness
- Finer mesh at stress concentrations (fillets, holes, welds)
Mesh Refinement:
- Use sphere of influence around critical features
- Create local mesh controls at:
- Fillet radii
- Nozzle-to-shell junctions
- Impeller eye (inlet)
- Impeller blade-to-hub junction
- Wear ring gaps
- Perform mesh convergence study (verify results converge as mesh density increases)
Recommended Element Sizes:
- Impeller blade thickness: minimum 4 elements
- Casing wall thickness: minimum 3 elements
- Fillet radius: minimum 6 elements around arc
- Global element size: 5-10% of characteristic dimension
5. Solver Selection
Linear Static Analysis:
- Use when: Deformations are small, material is elastic, no time-dependent effects
- Suitable for: Initial design verification, factor of safety calculations
- Fast computation, good for parametric studies
Nonlinear Analysis:
- Use when: Large deformations, material plasticity, contact nonlinearity
- Required for: Plastic collapse analysis, contact pressure distribution
- Slower computation, requires incremental loading
Modal Analysis:
- Determine natural frequencies and mode shapes
- Avoid resonance with pump operating speed and harmonics
- Check: 1st critical frequency > 1.5 × operating frequency
Harmonic/Transient Analysis:
- Pressure pulsations at blade passing frequency
- Startup and shutdown transients
- Water hammer events
Fatigue Analysis:
- High-cycle fatigue (> 10^4 cycles): Impellers, shafts
- Low-cycle fatigue: Thermal cycling, startup/shutdown
- Use S-N curves or strain-life methods
Solver Settings:
- Convergence criteria: typically 0.1% force/moment convergence
- Large deflection: ON for impellers with thin blades
- Contact settings: program-controlled for initial runs
- Stabilization: may be needed for very thin features
6. Results Interpretation
Primary Results:
Von Mises Stress (σ_vm):
- Equivalent stress for ductile materials
- Compare against material yield strength
- σ_vm = √[(σ₁-σ₂)² + (σ₂-σ₃)² + (σ₃-σ₁)²] / √2
- Most commonly used failure criterion
Principal Stresses (σ₁, σ₂, σ₃):
- Maximum and minimum normal stresses
- Critical for brittle materials
- Check maximum principal stress against ultimate strength
Deformation:
- Total deformation: overall displacement magnitude
- Directional deformation: check clearances (impeller-to-casing)
- Critical gaps: wear ring clearances, axial thrust bearing clearances
Safety Factor:
- Factor of Safety (FOS) = σ_allowable / σ_actual
- Color-coded: FOS < 1 (red/failure), FOS > 1 (safe)
Secondary Results:
Shear Stress:
- Critical at welds, keyways, shaft connections
- Maximum shear stress theory (Tresca criterion)
Strain:
- Elastic strain vs. plastic strain
- Strain energy density
- Equivalent plastic strain for ductile failure
Contact Pressure:
- Shrink fit interfaces
- Bearing surfaces
- Sealing faces
Critical Review Checklist:
- Do stress concentrations align with geometric features?
- Are boundary conditions realistic (no rigid body motion)?
- Is the deformation pattern physically reasonable?
- Are peak stresses at expected locations?
- Does the solution satisfy equilibrium (sum of reactions = applied loads)?
7. Safety Factor Evaluation
Design Codes and Standards:
- ASME Section VIII Division 1: Pressure vessels (casings)
- ASME Section VIII Division 2: Higher pressure, FEA-based design
- API 610: Centrifugal pumps for petroleum industry
- ISO 13709: Centrifugal pumps for petroleum, petrochemical and natural gas industries
Required Safety Factors:
Static Loads:
- Against yield: FOS ≥ 1.5 (general machinery)
- Against yield: FOS ≥ 2.0 (pressure vessels per ASME)
- Against ultimate: FOS ≥ 2.5 (static loads)
- Against ultimate: FOS ≥ 4.0 (dynamic/impact loads)
Fatigue Loads:
- Infinite life design: Stress < endurance limit
- Finite life: Use cumulative damage theory (Miner's rule)
- FOS on fatigue life: typically 2.0 or higher
Safety Factor Calculation:
For ductile materials:
FOS_yield = σ_y / σ_vm
FOS_ultimate = σ_u / σ_vm
Design FOS = min(FOS_yield, FOS_ultimate)
For brittle materials:
FOS = σ_u / σ₁ (based on maximum principal stress)
Material-Specific Adjustments:
- Cast iron: Reduce allowable stress by 25%
- Welded joints: Apply joint efficiency factor (0.7-1.0)
- High temperature: Apply creep reduction factors
- Corrosive environment: Add corrosion allowance
Applications to Pumps
Pump Casing Pressure Containment
Analysis Objectives:
- Verify casing can withstand maximum allowable working pressure (MAWP)
- Check stress at nozzle-to-shell junctions
- Evaluate flange sealing capability
- Assess casing distortion effects on wear ring clearances
Critical Regions:
- Volute cutwater (point of maximum pressure gradient)
- Nozzle reinforcement pads
- Casing split lines (horizontal or vertical split)
- Mounting foot attachments
Load Cases:
- Maximum operating pressure
- Hydrostatic test pressure (typically 1.5 × design pressure)
- Pressure + thermal expansion
- Pressure + piping loads (nozzle forces and moments)
Acceptance Criteria:
- Von Mises stress < allowable stress
- Flange face distortion < gasket capability
- No yielding at bolt holes
- Maximum deformation < 0.5% of casing diameter
Common Failure Modes:
- Casing rupture at cutwater
- Flange leakage due to uneven loading
- Cracking at weld toes
- Bolt failure under pressure cycling
Impeller Centrifugal Stress
Analysis Objectives:
- Calculate maximum stress in blades, shrouds, and hub
- Evaluate stress at blade-to-shroud junctions
- Check deflection to ensure clearance maintenance
- Assess fatigue life at blade passing frequency
Critical Regions:
- Blade trailing edge (maximum bending moment)
- Blade-to-hub fillet (stress concentration)
- Blade-to-shroud junction (closed impellers)
- Hub bore (stress concentration from keyway)
Load Cases:
- Centrifugal force at rated speed
- Centrifugal force at overspeed (120% of rated)
- Combined centrifugal + hydraulic pressure
- Thermal loads from fluid temperature
Stress Components:
Centrifugal Stress (Hoop Stress):
σ_hoop = ρω²r² (for thin disk)
where:
ρ = material density (kg/m³)
ω = angular velocity (rad/s)
r = radius (m)
Blade Bending Stress:
- From pressure difference across blade
- Maximum at blade root (hub junction)
- Increases with blade length and pressure differential
Acceptance Criteria:
- Von Mises stress at rated speed < 0.5 × σ_y
- Von Mises stress at overspeed < 0.67 × σ_y
- Blade tip deflection < 50% of clearance gap
- No natural frequencies within ±20% of operating speed range
Design Optimizations:
- Increase hub diameter to reduce stress
- Taper blade thickness from hub to shroud
- Add fillets at blade-to-hub junction (r ≥ 3mm typical)
- Use materials with higher strength-to-density ratio
Shaft Deflection
Analysis Objectives:
- Calculate shaft deflection under radial and axial loads
- Verify bearing alignment is maintained
- Ensure coupling alignment tolerances are met
- Check shaft critical speeds vs. operating speed
Critical Regions:
- Midspan between bearings (maximum deflection)
- Impeller location (affects wear ring clearance)
- Coupling location (misalignment causes vibration)
- Bearing journals (stress concentration at shoulders)
Load Cases:
- Radial hydraulic thrust (maximum at shutoff)
- Axial hydraulic thrust
- Impeller weight + fluid weight
- Thermal expansion
Hydraulic Radial Thrust:
F_radial = K_r × ρ × g × H × D₂ × b₂
where:
K_r = radial thrust coefficient (0.2-0.4 for single volute)
H = head (m)
D₂ = impeller diameter (m)
b₂ = impeller outlet width (m)
Acceptance Criteria:
- Maximum shaft deflection < 0.0005 × shaft span
- Deflection at impeller < 25% of wear ring clearance
- Shaft stress < 0.3 × σ_y (for fatigue resistance)
- 1st critical speed > 1.5 × maximum operating speed
Deflection Calculation:
- Use beam theory for preliminary estimates
- FEA for complex geometries and load distributions
- Include gyroscopic effects for high-speed pumps
Fatigue Analysis
Analysis Objectives:
- Predict fatigue life in cycles or years
- Identify locations prone to crack initiation
- Evaluate stress concentration factors
- Design for infinite life or safe finite life
Fatigue-Critical Locations:
- Impeller blade trailing edges (pressure pulsations)
- Blade-to-hub fillets (stress concentration)
- Shaft keyway (stress concentration factor ~3)
- Shaft shoulder at bearing locations
- Weld toes on casings
Load Cycles:
- Startup/shutdown cycles: low-cycle fatigue
- Blade passing frequency: high-cycle fatigue
- Pressure pulsations: high-cycle fatigue
- Rotor imbalance: high-cycle fatigue
Fatigue Analysis Methods:
S-N Curve Method (High-Cycle Fatigue):
- Applicable for N > 10⁴ cycles
- Use material S-N curves (stress vs. cycles to failure)
- Apply mean stress correction (Goodman or Soderberg)
- Calculate cumulative damage using Miner's rule
Strain-Life Method (Low-Cycle Fatigue):
- Applicable for N < 10⁴ cycles
- Uses plastic strain range
- Coffin-Manson equation
- Required for startup/shutdown analysis
Stress Concentration Factors:
- Sharp corners: K_t = 2-3
- Keyways: K_t = 2-3
- Fillets: K_t = 1.5-2.5 (radius dependent)
- Threads: K_t = 2-4
Acceptance Criteria:
- Infinite life: alternating stress < endurance limit
- Finite life: calculated life > 2 × required life
- Factor of safety on stress amplitude ≥ 1.5
- Safety factor on life ≥ 2.0
Mean Stress Effects:
Goodman correction:
σ_a / S_e + σ_m / σ_u = 1 / FOS
where:
σ_a = alternating stress amplitude
σ_m = mean stress
S_e = endurance limit
σ_u = ultimate tensile strength
Material Selection Criteria
Pump Casing Materials
Cast Iron (ASTM A48, A278):
- Applications: Low-pressure water service, non-corrosive fluids
- Properties: σ_u = 150-400 MPa, low cost, good castability
- Limitations: Brittle, poor fatigue resistance, limited to <250°C
- Use for: Municipal water, HVAC, pressures < 20 bar
Ductile Iron (ASTM A536):
- Applications: General industrial service, moderate pressures
- Properties: σ_y = 275-550 MPa, σ_u = 400-800 MPa, good machinability
- Grades: 65-45-12, 80-55-06 (σ_u-σ_y-elongation)
- Use for: Process water, oil transfer, pressures < 40 bar
Carbon Steel (ASTM A216 WCB):
- Applications: High-pressure, high-temperature service
- Properties: σ_y = 250 MPa, σ_u = 485 MPa, good weldability
- Temperature range: -29°C to 400°C
- Use for: Boiler feed, power plants, pressures < 100 bar
Stainless Steel 316 (ASTM A743 CF-8M):
- Applications: Corrosive fluids, seawater, chemicals
- Properties: σ_y = 205 MPa, σ_u = 485 MPa, excellent corrosion resistance
- Temperature range: -196°C to 400°C
- Use for: Chemical processing, marine, food & beverage
Bronze (ASTM B584):
- Applications: Seawater service, small pumps
- Properties: Good corrosion resistance, moderate strength
- Use for: Marine, desalination, pump internals
Impeller Materials
Ductile Iron (ASTM A536):
- Best strength-to-cost ratio
- Density: 7,100 kg/m³
- σ_y = 275 MPa minimum
- Use for: General service, water
316 Stainless Steel:
- Corrosion resistance
- Density: 8,000 kg/m³
- σ_y = 275 MPa (cast), 290 MPa (investment cast)
- Use for: Chemical, food, pharmaceutical
Duplex Stainless (CD4MCu):
- High strength and corrosion resistance
- σ_y = 450 MPa
- Expensive but long-lasting
- Use for: Seawater, aggressive chemicals
Nickel-Aluminum Bronze (NAB):
- Excellent cavitation resistance
- σ_y = 240 MPa
- Best for seawater
- Use for: Marine, desalination
Titanium (Ti-6Al-4V):
- Highest strength-to-weight ratio
- Density: 4,430 kg/m³
- σ_y = 880 MPa
- Very expensive
- Use for: High-speed, aerospace
Shaft Materials
Carbon Steel (AISI 1045, 4140):
- General purpose
- σ_y = 400-650 MPa
- Good machinability
- Requires corrosion protection
Stainless Steel 416, 17-4PH:
- Corrosive environments
- σ_y = 520-1,170 MPa (17-4PH heat treated)
- Good fatigue resistance
- 17-4PH for high loads
Alloy Steel (AISI 4340):
- High-power applications
- σ_y = 860-1,380 MPa (heat treated)
- Excellent fatigue properties
- Requires protection coating
Material Selection Process
Step 1: Service Conditions
- Fluid type (corrosive, abrasive, clean)
- Temperature range
- Pressure level
- Environmental exposure
Step 2: Performance Requirements
- Required strength (stress levels)
- Fatigue resistance (cyclic loading)
- Wear resistance (sliding contact)
- Thermal expansion compatibility
Step 3: Manufacturing Considerations
- Casting vs. machining from bar stock
- Weldability requirements
- Heat treatment feasibility
- Inspection requirements (NDT)
Step 4: Cost-Benefit Analysis
- Material cost
- Processing cost
- Expected lifetime
- Maintenance costs
Step 5: Code Compliance
- ASME allowable stresses
- Temperature limits
- Impact testing requirements (low temperature)
- Corrosion allowance
Stress Concentrations
Theoretical Stress Concentration Factor (K_t)
The stress concentration factor relates the peak local stress to the nominal stress:
σ_max = K_t × σ_nominal
Common Stress Concentrations in Pumps
Fillets:
K_t = 1 + 2√(d/r)
where:
d = smaller diameter
r = fillet radius
Design rule: r ≥ 0.1 × d for K_t ≤ 2.0
Shoulder (step change in diameter):
- r/d = 0.05: K_t ≈ 2.5
- r/d = 0.10: K_t ≈ 2.0
- r/d = 0.20: K_t ≈ 1.6
Circular hole in plate:
- Single hole: K_t = 3.0
- Multiple holes: interaction increases K_t
Keyway on shaft:
- Profile keyway: K_t = 2.0-2.5
- Sled runner keyway: K_t = 2.5-3.0
- End-milled keyway: K_t = 3.0-3.5
Threads:
- Coarse thread: K_t = 2.2-3.0
- Fine thread: K_t = 2.8-4.0
- Stress relief groove reduces K_t
Welds:
- Toe of weld: K_t = 2.0-3.0
- Root of weld: K_t = 2.5-4.0
- Grinding flush reduces K_t by 30-50%
Stress Concentration Mitigation Strategies
Increase Fillet Radius:
- Double the radius → reduce K_t by 20-30%
- Optimize using FEA parametric study
- Blended radii better than constant radius
Add Relief Grooves:
- Undercut at shoulder changes
- Reduces stress at critical section
- Common at bearing fits
Eliminate Sharp Corners:
- All internal corners should have radius
- Minimum radius = 2 × wall thickness
- Check manufacturability
Relocate Transitions:
- Move stress concentrations away from high nominal stress regions
- Place diameter changes away from high bending moment locations
Use Material Locally:
- Increase local thickness at stress concentration
- Add boss or reinforcement pad
- Ensure smooth transition
Surface Treatment:
- Shot peening introduces compressive residual stress
- Improves fatigue life by 2-4x
- Effective for shafts, blades
Fatigue Notch Factor (K_f)
For fatigue analysis, the effective stress concentration is reduced:
K_f = 1 + q(K_t - 1)
where:
q = notch sensitivity (0 to 1)
K_t = theoretical stress concentration factor
Notch sensitivity depends on:
- Material ductility (lower for ductile materials)
- Notch radius (lower for larger radii)
- Material grain size
Typical values:
- Steel, r = 0.5 mm: q ≈ 0.9
- Steel, r = 2.0 mm: q ≈ 0.95
- Aluminum: q ≈ 0.7-0.8
Verification with FEA
Local Mesh Refinement:
- Mesh size at notch < r/4 (fillet radius)
- Use at least 6-8 elements around fillet arc
- Verify stress convergence with mesh refinement
Linearization:
- Extract stress through section thickness
- Separate membrane and bending components
- Use for pressure vessel code compliance
Submodeling:
- Global model with coarse mesh
- Local model with fine mesh around stress concentration
- Apply displacement boundary conditions from global to local
Usage Guidelines
This structural analysis workflow should be applied iteratively:
- Preliminary Design: Use hand calculations and simplified FEA
- Design Refinement: Detailed FEA with accurate geometry and loading
- Design Verification: Final analysis with all load cases and safety factors
- Fatigue Assessment: If cyclic loading is significant
- Documentation: Report stress levels, safety factors, and compliance with codes
Always validate FEA results against:
- Theoretical calculations where possible
- Previous designs with field experience
- Industry standards and guidelines
- Physical testing when available
Remember: FEA is a tool to aid engineering judgment, not replace it. Always perform sanity checks on results and ensure boundary conditions represent actual installation and operating conditions.
