Structural Analysis Workflow

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 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 (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:

1. Do stress concentrations align with geometric features?
2. Are boundary conditions realistic (no rigid body motion)?
3. Is the deformation pattern physically reasonable?
4. Are peak stresses at expected locations?
5. 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 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 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:

1. Preliminary Design: Use hand calculations and simplified FEA
2. Design Refinement: Detailed FEA with accurate geometry and loading
3. Design Verification: Final analysis with all load cases and safety factors
4. Fatigue Assessment: If cyclic loading is significant
5. 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.

Source & license

This open-source skill is cataloged on AgentStack and links to its original source — we do not rehost the code.

• Author: Soljourner
• Source: Soljourner/claude-engineering-skills
• License: MIT

Install and usage instructions live in the source repository linked above.