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.

184 N. Zobeiry et al.<br />

In compression, on the other hand, this soften<strong>in</strong>g occurs due to matrix crack<strong>in</strong>g<br />

or yield<strong>in</strong>g and fibre <strong>in</strong>stability. In addition, the <strong>in</strong>stability and rotation of the<br />

fibres comb<strong>in</strong>ed with fibre breakage and splitt<strong>in</strong>g contribute to the overall soften<strong>in</strong>g<br />

behaviour under compression. In fact, <strong>in</strong> compression, matrix crack<strong>in</strong>g leads to<br />

the <strong>in</strong>stability of fibres which is not the case <strong>in</strong> tension.<br />

The overall reduction <strong>in</strong> the lam<strong>in</strong>ate tangential stiffness, which leads to the soften<strong>in</strong>g<br />

behaviour, has been modeled us<strong>in</strong>g fuses and spr<strong>in</strong>gs. Fuses fail sequentially<br />

upon application of a remote displacement to the lam<strong>in</strong>ate sub-model. This results<br />

<strong>in</strong> the overall stiffness reduction of the lam<strong>in</strong>ate sub-model and an ensu<strong>in</strong>g soften<strong>in</strong>g<br />

behaviour similar to the composite response under tension.<br />

The rubble sub-model represents the behaviour of the damaged material <strong>in</strong> compression.<br />

Upon reach<strong>in</strong>g the damage <strong>in</strong>itiation stra<strong>in</strong> under compressive loads, the<br />

matrix starts to crack or yield result<strong>in</strong>g <strong>in</strong> the <strong>for</strong>mation of multiple cracked surfaces<br />

or rubble. With further damage growth under compression, the friction between the<br />

newly cracked surfaces results <strong>in</strong> an <strong>in</strong>crease <strong>in</strong> the load carry<strong>in</strong>g capacity of the<br />

material. This damaged material can cont<strong>in</strong>ue to carry compressive loads.<br />

The behaviour of the rubble sub-model can be represented with an analogue<br />

model consist<strong>in</strong>g of gap and spr<strong>in</strong>g elements to model the compressive load carry<strong>in</strong>g<br />

response with an <strong>in</strong>creas<strong>in</strong>g stiffness. The gap element is <strong>in</strong>active (open)<br />

be<strong>for</strong>e damage <strong>in</strong>itiation. After damage <strong>in</strong>itiates, new surfaces <strong>for</strong>m result<strong>in</strong>g <strong>in</strong> an<br />

<strong>in</strong>creased stiffness of the damaged material. This is modeled by sequential closure<br />

of the vary<strong>in</strong>g sized gaps <strong>in</strong> compression. By virtue of the gaps rema<strong>in</strong><strong>in</strong>g open, the<br />

rubble sub-model cannot carry any load <strong>in</strong> tension.<br />

The slider simulates the damage band (k<strong>in</strong>k band) propagation. The yield stra<strong>in</strong><br />

is assumed to co<strong>in</strong>cide with the saturation stra<strong>in</strong> of the rubble sub-model. At the<br />

<strong>in</strong>stant of f<strong>in</strong>al gap closure, the slider activates, <strong>in</strong>dicat<strong>in</strong>g progressive specimen end<br />

shorten<strong>in</strong>g (damage zone broaden<strong>in</strong>g) under a constant applied load. This results<br />

<strong>in</strong> the damage propagation <strong>in</strong>to the undamaged <strong>in</strong>terior material and an <strong>in</strong>creas<strong>in</strong>g<br />

damage height <strong>for</strong> compressive failure.<br />

Together, the Lam<strong>in</strong>ate and Rubble sub-models represent the overall response of<br />

the composite <strong>in</strong> the <strong>for</strong>ce-displacement space dur<strong>in</strong>g progressive damage growth<br />

<strong>in</strong> tension or compression. Figure 9.3 shows a schematic example of the response<br />

of a RVE under controlled displacement load<strong>in</strong>g based on the analogue model. First<br />

displacement is applied <strong>in</strong> the tensile mode until fibre and matrix damage occurs<br />

<strong>in</strong> the RVE. Subsequently the RVE is unloaded and then reloaded <strong>in</strong> the compressive<br />

regime. In this case, the <strong>in</strong>itial reduced modulus (dashed l<strong>in</strong>e) compared<br />

to the undamaged modulus (solid curve) is due to the previous damage <strong>in</strong>curred<br />

<strong>in</strong> tension. After reach<strong>in</strong>g a saturation state of matrix damage <strong>in</strong> the RVE, upon<br />

further application of displacement, the damage broadens <strong>in</strong>to the undamaged material<br />

result<strong>in</strong>g <strong>in</strong> the band broaden<strong>in</strong>g phenomena and the correspond<strong>in</strong>g plateau<br />

stress. The band broaden<strong>in</strong>g then cont<strong>in</strong>ues up to the po<strong>in</strong>t of complete failure (fibre<br />

damage saturation).<br />

The analogue model described above has been used to extend the compression<br />

model<strong>in</strong>g capabilities of a cont<strong>in</strong>uum damage mechanics model (CODAM) previously<br />

developed by Williams et al. [31] <strong>for</strong> gross damage development <strong>in</strong> polymeric

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

Saved successfully!

Ooh no, something went wrong!