11.02.2013 Views

Computational Methods for Debonding in Composites

Computational Methods for Debonding in Composites

Computational Methods for Debonding in Composites

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

102 C. Schuecker and H.E. Pettermann<br />

5.2 Plasticity Model<br />

The plasticity model proposed here is a phenomenological approach based on the<br />

assumption that plastic de<strong>for</strong>mation occurs <strong>in</strong> the <strong>for</strong>m of shear bands with a specific<br />

orientation. As mentioned previously, experimental results of a study <strong>in</strong>to the<br />

micromechanics of PMC failure mechanisms suggest that the shear bands are precursors<br />

of later cracks [3]. It is there<strong>for</strong>e assumed further that the shear bands<br />

develop <strong>in</strong> planes that have the same orientation as the fracture plane which is predicted<br />

here by Puck’s action plane failure criterion <strong>for</strong> plane stress states (Puck 2D)<br />

[17, 20, 23].<br />

Accord<strong>in</strong>g to Puck, fracture occurs <strong>in</strong> a plane that is parallel to the fiber orientation<br />

and def<strong>in</strong>ed by a fracture plane angle, θfp, as depicted <strong>in</strong> Fig. 5.2, left.<br />

The fracture plane is perpendicular to the lam<strong>in</strong>ate plane (i.e. θfp = 0) <strong>in</strong> the case<br />

of comb<strong>in</strong>ed transverse tensile stresses and <strong>in</strong>-plane shear (mode A) or moderate<br />

transverse compression with <strong>in</strong>-plane shear (mode B). For comb<strong>in</strong>ations of high<br />

transverse compressive stresses with shear, the fracture plane angle is non-zero<br />

(mode C) and can be computed analytically [17,20,23]. For example, uniaxial transverse<br />

compression of epoxy matrix composites typically leads to a predicted fracture<br />

plane angle of approximately θfp = 50 ◦ −56 ◦ , which agrees well with experimental<br />

f<strong>in</strong>d<strong>in</strong>gs [3,17,20] and results of micromechanical models [9]. In light of the considerations<br />

put <strong>for</strong>th above, the plasticity law is <strong>for</strong>mulated with respect to the predicted<br />

fracture plane.<br />

5.2.1 Plastic Stra<strong>in</strong> <strong>for</strong> θθθ fp = 0 (Puck Modes A and B)<br />

In a perpendicular plane, θfp = 0, the only shear stress component under plane<br />

stress conditions is σ12. The relation of shear stress and plastic shear stra<strong>in</strong> <strong>for</strong> <strong>in</strong>-<br />

plane simple shear, γ pl<br />

12 (σ12,σ22 = 0), can be derived from experimental data and<br />

described by any suitable analytical expression. Here, a power law,<br />

1<br />

3<br />

γ12 = γ el pl σ12<br />

12 + γ12 = +<br />

G12<br />

t<br />

l<br />

n<br />

2<br />

θ fp<br />

� �<br />

σ12<br />

n<br />

k<br />

sf<br />

p<br />

l<br />

snn<br />

s ny<br />

sln<br />

t<br />

y<br />

snt<br />

n<br />

(5.1)<br />

Fig. 5.2 Def<strong>in</strong>ition of fracture plane and correspond<strong>in</strong>g coord<strong>in</strong>ate system, l-n-t, with regard to the<br />

ply coord<strong>in</strong>ate system, 1-2-3, by fracture plane angle, θfp (left); tractions on the fracture plane <strong>for</strong><br />

θfp �= 0(right)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!