Moisture Simulation Model for Polymer Composites

Moisture Simulation Model for Polymer Composites

OBJECTIVES

• Develop 3D finite element model for moisture absorption of fiber-reinforced polymer composites.

• Perform case studies for glass fiber-reinforced polymer composites over long time periods.

• Validate the simulation results with the analytical and experimental results.

FINITE ELEMENT MODELING & SIMULATION

FINITE ELEMENT RESULTS

CONCLUSIONS AND FUTURE WORK• Simulation results match well with 1D analytical results,

which validates the finite element model.

• Results are presented for both 2D and 3D models.

• Study polymer composites moisture evolution under complex environments.

• Evaluate the effects of temperature and humidity on polymer composites moisture behavior subjected to long time exposure.

ACKNOWLEDGEMENTS

Student(s): Zhen Huo Mechanical & Aerospace Engineering Department

Faculty Advisor(s): Dr. K. Chandrashekhara, Mech. & Aero. Engineering Dept. Dr. Thomas Schuman, Chemistry Dept.

This project is funded by Bell Helicopter Inc. and Intelligent Systems Center.

GOVERNING EQUATIONS

X=0 X=1.7703 mm

B.C.=SATURATION

Face B.C. is

applied

Diffusivity (D11) 4.44158E-5 mm2/h

Saturation concentration 4170 ppm

Half thickness 1.7703 mm

ABAQUS 2D half modeling ABAQUS 3D half modeling

where c is moisture concentration, t is time span, n is outward normal to S. J is concentration flux and is given by:

where s is solubility, ∅ is the normalized concentration, which equals c/s, ks is “soret effect” factor, providing temperature gradient driven diffusion; kp is the pressure stress factor,

providing pressure gradient driven diffusion; θz is absolute zero temperature.

ABAQUS 3D half modeling moisture content after 168 hours

D11= 4.44158E-5 mm2/h; D22= D11 / 2; D33 = D11 / 4

ABAQUS 2D half modeling moisture content after 168 hours

INTRODUCTION

• Composite materials degrade due to environmental exposure such as moisture diffusion, thermal spike, UV, and thermal oxidation.

• Moisture absorption of composites decreases the strength and stiffness, degrade the fiber/matrix interface strength, swell and plasticize the resin lowering its glass transition temperature.

• The relative degree of degradation is related to chemistry of the reinforcement and matrix as well as the exposure time.

• The determination of moisture concentration profile through the thickness of the composite structure requires the knowledge of 3D moisture diffusion coefficients.

FINITE ELEMENT RESULTS (CONT. )

Convergence study with half modeling after 168 soaking

hours

Aspect ratio effects study with half modeling using

50 nodes mesh

VARIANCE STUDY

2D half modeling moisture content results comparison

after 168 hours

3D half modeling moisture content results comparison

after 168 hours

VALIDATION

0 0.02 0.04 0.06 0.080

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

X position to the midpoint of whole length (inch)

Mo

istu

re C

on

ten

t (%

)

10 nodes mesh30 nodes mesh50 nodes mesh

0 0.02 0.04 0.06 0.080

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

X position to the midpoint of whole length (inch)

Mo

istu

re C

on

ten

t (%

)

ASPECT RATIO 2.8ASPECT RATIO 10ASPECT RATIO 20

0 0.02 0.04 0.06 0.080

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

X position to the midpoint of whole length (inch)

Mo

istu

re C

on

ten

t (%

)

1D analytical solution2D finite element half model

0 0.01 0.02 0.03 0.04 0.05 0.06 0.070

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

X position to the midpoint of whole length (inch)

Mo

istu

re C

on

ten

t (%

)

1D analytical solution3D finite element half model