AA2195-AZ31B Friction Stir Welding Temperature Field Numerical Simulation

1. Definition and Fundamental Principles

Friction Stir Welding (FSW) is a solid-state joining process that generates heat through mechanical friction between a rotating tool and the workpiece, avoiding the melting and re-solidification associated with fusion welding. The temperature field numerical simulation of FSW involves the use of finite element analysis (FEA) or computational fluid dynamics (CFD) methods to predict the three-dimensional thermal distribution generated during the welding process. This simulation capability is critical for understanding the thermomechanical behavior of the material, predicting microstructural evolution, and optimizing process parameters prior to physical experimentation.

The specific material system under study—AA2195 (an Al-Cu-Li 2xxx series aluminum alloy) and AZ31B (an Mg-Al-Zn 3xxx series magnesium alloy)—represents a dissimilar lightweight alloy couple of significant interest in aerospace, automotive, and structural applications. The large difference in thermal conductivity (AA2195: ~120 W/m·K; AZ31B: ~72 W/m·K), melting points (AA2195: ~638 °C; AZ31B: ~450 °C), and coefficients of thermal expansion between these materials makes the temperature field highly non-uniform and challenging to predict without rigorous numerical modeling.

The governing equation for the temperature field in FSW is the transient heat conduction equation with a moving heat source:

ρ·Cp·(∂T/∂t + vx·∂T/∂x + vy·∂T/∂y + vz·∂T/∂z) = ∂/∂x(k·∂T/∂x) + ∂/∂y(k·∂T/∂y) + ∂/∂z(k·∂T/∂z) + Q

where ρ is density, Cp is specific heat capacity, v is velocity field, k is thermal conductivity, T is temperature, and Q is the volumetric heat source term representing frictional and plastic deformation heat generation.

2. Category and Business Positioning

This numerical simulation capability belongs to the advanced process engineering and computational design category within Cladding Technology Shanxi Co., Ltd.'s technical portfolio. It serves as an upstream enabling technology that supports all three primary manufacturing routes:

Within the company's qualification and certification framework, this simulation capability positions Cladding Technology Shanxi as a technically advanced provider capable of offering predictive process design services, reducing trial-and-error costs, and accelerating time-to-market for customers in aerospace, energy, and transportation sectors.

3. Technical Purpose and Value

3.1 Process Optimization

The primary purpose of AA2195-AZ31B FSW temperature field numerical simulation is to determine optimal process parameters—including tool rotational speed, traverse speed, tool geometry (shoulder diameter, pin profile), and plunge depth—that achieve sound metallurgical bonding without excessive thermal damage to either parent material. The simulation provides:

3.2 Material Compatibility Assessment

The AA2195-AZ31B system is particularly challenging because:

Numerical simulation enables engineers to map safe processing windows where bonding occurs without exceeding the solidus temperature of AZ31B or inducing excessive intermetallic growth.

3.3 Cost Reduction and Risk Mitigation

Physical FSW trials on AA2195-AZ31B are expensive due to the high cost of Al-Li alloys, the sensitivity of Mg alloys to contamination, and the specialized equipment required. Each physical trial may cost several thousand RMB in material and machine time. Numerical simulation allows pre-screening of parameter combinations, reducing the number of physical trials by 60-80% while maintaining confidence in the final process specification.

4. Key Process and Implementation Points

4.1 Material Property Database Requirements

Property AA2195 (Al-Li Alloy) AZ31B (Mg Alloy) Temperature Range
Density (kg/m³) 2,700 1,810 25-500 °C
Thermal Conductivity (W/m·K) 120 (RT) → 180 (400°C) 72 (RT) → 100 (400°C) 25-500 °C
Specific Heat (J/kg·K) 900 (RT) → 1,100 (400°C) 1,000 (RT) → 1,200 (400°C) 25-500 °C
Young's Modulus (GPa) 73 (RT) → 45 (400°C) 45 (RT) → 25 (400°C) 25-500 °C
Solidus Temperature (°C) ~638 ~450
Yield Strength (MPa) 345 (RT) → 150 (400°C) 175 (RT) → 60 (400°C) 25-500 °C

4.2 Heat Source Modeling

The heat generation in FSW is distributed across two primary sources:

The total heat input rate is estimated as:

Qtotal = η · π · r² · f · μ · σ · v

where η is the efficiency factor (typically 0.8-0.9), r is shoulder radius, f is rotational frequency, μ is friction coefficient, σ is contact pressure, and v is linear velocity.

4.3 Boundary Conditions and Mesh Configuration

Parameter Recommended Value/Setting Rationale
Element type 8-node hexahedral (C3D8RT) Accurate thermal stress coupling; reduced integration for computational efficiency
Mesh density near tool 0.5-1.0 mm element size Capture steep thermal gradients within the thermally affected zone
Mesh density far field 3.0-5.0 mm element size Reduce computational cost while maintaining boundary accuracy
Top surface heat transfer Convective + radiative (h = 5-20 W/m²·K) Accounts for air cooling and oxidation layer radiation
Bottom surface (backing plate) Convective (h = 20-50 W/m²·K) or fixed temperature Represents water-cooled or uninsulated backing plate
Initial temperature 25 °C (ambient) Standard room temperature baseline
Tool temperature Adiabatic or convection-coupled Tool acts as heat source; temperature not independently tracked in thermal-only models

4.4 Recommended FSW Process Parameters for AA2195-AZ31B

Parameter Range (Low) Range (High) Optimal Target Constraint
Tool rotational speed 600 rpm 1,200 rpm 800-1,000 rpm Avoid exceeding AZ31B solidus (~450°C)
Traverse speed 30 mm/min 100 mm/min 50-70 mm/min Balance bonding quality vs. heat input
Speed ratio (ω/v) 100 mm⁻¹ 200 mm⁻¹ 120-160 mm⁻¹ Higher ratio = more heat per unit travel
Shoulder diameter 12 mm 16 mm 14 mm Match total plate thickness (3-5 mm each)
Pin diameter 4 mm 6 mm 5 mm Ensure full penetration without excessive material flow
Pin profile Tapered conical or threaded Enhance material stirring and mixing
Plunge depth Full thickness + 0.5 mm into backing Ensure complete penetration and root bonding

4.5 Simulation Workflow

  1. Geometry Modeling: Create 3D model of AA2195 plate (typically 3-5 mm thick) bonded to AZ31B plate (3-5 mm thick) with FSW tool assembly.
  2. Material Property Assignment: Input temperature-dependent properties for both materials as tabulated above.
  3. Mesh Generation: Generate adaptive mesh with refinement near tool-workpiece interface; use element remeshing or arbitrary Lagrangian-Eulerian (ALE) technique for large deformation.
  4. Heat Source Implementation: Define frictional and plastic deformation heat sources as user-defined subroutines (e.g., DFLUX in ABAQUS).
  5. Boundary Condition Application: Set convective/radiative cooling on exposed surfaces; apply displacement constraints on backing plate.
  6. Solution Procedure: Execute coupled thermal-mechanical analysis or sequential thermal analysis followed by stress evaluation.
  7. Post-Processing: Extract temperature contours, thermal cycles at critical locations, peak temperatures, and cooling rates.
  8. Validation: Compare simulated temperature profiles against thermocouple measurements from physical trials (tolerance: ±15 °C).
  9. Parameter Optimization: Iterate simulation with different parameter sets to identify optimal window.

5. Applicable Standards and Acceptance Criteria

5.1 FSW Process Standards

5.2 Simulation Validation Criteria

Validation Parameter Acceptance Criterion Measurement Method
Peak temperature prediction ±15 °C of thermocouple measurement K-type or N-type thermocouples at defined locations
Thermal cycle shape Correlation coefficient R² > 0.90 Time-temperature history comparison
Thermally affected zone width ±10% of measured value Hardness traverse (HV0.2) across weld cross-section
Cooling rate at weld center ±20% of measured value High-speed thermocouple or thermochromic paint
Weld distortion ±15% of measured value Laser scanning or digital image correlation (DIC)

5.3 Weld Joint Acceptance Standards

6. Common Risks and Controls

6.1 Technical Risks

Risk Description Mitigation Control
AZ31B partial melting Excessive heat input causes local melting of magnesium alloy, leading to porosity and cracking upon solidification Simulate and enforce peak temperature limit of 430°C (20°C below solidus); use lower speed ratios
Intermetallic compound overgrowth Formation of brittle Al-Mg intermetallics at the bonding interface, degrading ductility Limit peak temperature and holding time; simulation predicts time above 350°C threshold
Tool wear and failure High-temperature wear of tungsten carbide tool in contact with both Al and Mg alloys Simulate tool temperature distribution; select tool materials with appropriate hardness retention (e.g., M2 steel for Mg, tungsten carbide for Al)
Simulation-model mismatch Discrepancy between predicted and actual thermal profiles due to inaccurate material properties or boundary conditions Validate with physical trials; update property databases; perform sensitivity analysis on key parameters
Galvanic corrosion in service Electrochemical incompatibility between AA2195 and AZ31B in corrosive environments Specify protective coatings per NACE SP0285; design with electrical isolation; limit joint thickness ratio
Magnesium oxidation and contamination AZ31B is highly reactive with oxygen and moisture; contamination degrades weld quality Specify inert atmosphere (Ar or He) shielding; control ambient humidity; simulation cannot address this—requires procedural controls

6.2 Quality Management Controls

7. Application Scenarios Across Company Technology Routes

7.1 TIG/MIG Weld Overlay Applications

The temperature field simulation methodology developed for AA2195-AZ31B FSW directly transfers to TIG/MIG weld overlay process design for dissimilar aluminum-magnesium cladding systems. Key applications include:

For TIG overlay of magnesium alloys onto aluminum substrates, the simulated temperature profiles guide the selection of backing material (copper backing plate with water cooling) and interpass temperature limits (typically ≤150 °C between passes).

7.2 Hydraulic Explosive Bonding Applications

In hydraulic explosive bonding, the temperature field simulation addresses the thermal effects of high-velocity flyer plate impact. The methodology is adapted as follows:

The temperature field data from simulation guides the selection of explosive charge composition (typically RDX or PETN-based), charge thickness, and stand-off distance to achieve the target interface temperature window of 350-500 °C for AA2195-AZ31B bonding.

7.3 Explosion Welding Applications

For explosion welding of AA2195 onto AZ31B (or vice versa), the temperature field simulation is critical for:

The simulation output directly informs the WPS for explosion welding, including charge parameters, gap distance, and clamping configuration, ensuring compliance with ASME BPV Section II Part D material requirements and NACE MR0175 sour service qualification where applicable.

8. Contribution to Qualification Building, Product Delivery, and Customer Value

8.1 Qualification Building

8.2 Product Delivery

8.3 Customer Value

9. Conclusion

The AA2195-AZ31B Friction Stir Welding Temperature Field Numerical Simulation capability represents a sophisticated computational engineering asset that enhances Cladding Technology Shanxi Co., Ltd.'s ability to deliver high-quality dissimilar metal joining solutions across all three technology routes. By providing predictive, validated thermal process design, this capability reduces development risk, accelerates qualification timelines, and delivers measurable value to customers in aerospace, automotive, energy, and marine industries. The rigorous methodology—combining temperature-dependent material properties, validated heat source models, and systematic post-processing—ensures that simulation outputs are reliable enough to guide production decisions while maintaining compliance with applicable standards including ASTM, ASME, NACE, and ISO frameworks.