11.12.2012 Views

Nondestructive testing of defects in adhesive joints

Nondestructive testing of defects in adhesive joints

Nondestructive testing of defects in adhesive joints

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

It is an unsymmetrical two po<strong>in</strong>t tensor and is not completely related to material configuration<br />

The force vector dh can be related to its <strong>in</strong>itial configuration <strong>in</strong> a similar manner<br />

dh = F.dH (2.15)<br />

And rearranged as , dh = JF -1 F -T .dA = S.dA (2.16)<br />

Where S is the second piola-kirchh<strong>of</strong>f stress tensor def<strong>in</strong>ed as S=JF -1 F -1 (2.17)<br />

Hyperelasticity is type <strong>of</strong> elasticity where the stress at any po<strong>in</strong>t can be derived from the<br />

deformation gradient and from an energy function. The behaviour <strong>of</strong> the material is said to be path<br />

<strong>in</strong>dependent and P is work conjugate with the rate <strong>of</strong> deformation gradient ₣,a stored elastic potential Ψ<br />

per unit undeformed volume can be established as the work done by the stresses form the <strong>in</strong>itial to the<br />

current position as, Ψ = P : ₣ =S : Ė= (∂Ψ/∂E) Ė (2.18)<br />

Observ<strong>in</strong>g that ½ ċ=Ė is work conjugate to the second PK stress enables a totally lagrangian<br />

constitutive equation to be constructed <strong>in</strong> the same manner as S = ∂Ψ/∂E = 2∂Ψ/∂c (2.19)<br />

Based on the isotropic assumption, the stra<strong>in</strong> energy density function W is expressed as a<br />

function <strong>of</strong> the stra<strong>in</strong> <strong>in</strong>variants .A better approach to model<strong>in</strong>g the response <strong>of</strong> rubbers comes from<br />

assum<strong>in</strong>g existence <strong>of</strong> stra<strong>in</strong> energy which is a function <strong>of</strong> deformation tensor. Stra<strong>in</strong> <strong>in</strong>variants can be<br />

expressed <strong>in</strong> terms <strong>of</strong> the stretch ratios λ1, λ2, and λ3 as follows<br />

Ψ = Ψ (I1, I2, I3)<br />

I1 = λ1 2 + λ2 2 + λ3 2 = tr (c)<br />

I2 = λ1 2 λ2 2 + λ2 2 λ3 2 + λ3 2 λ1 2 = ½ [(tr c) 2 -tr (c 2 )]<br />

I3 = λ1 2 λ2 2 λ3 2 = det c (2.20)<br />

For <strong>in</strong>compressible material I3=1.The stra<strong>in</strong> energy function for Neo-Hookean material is given as<br />

Ψ=C1 (I1-3) (2.21)<br />

In many cases rubber is compressible to some extent and it is important to consider compressibility.<br />

This work deals with the compressible form and <strong>in</strong>troduces a stra<strong>in</strong> energy contribution <strong>in</strong>volv<strong>in</strong>g the<br />

bulk modulus. For compressible form [4], the stra<strong>in</strong> energy should comb<strong>in</strong>e both volumetric and<br />

deviatoric term.<br />

Ψ = Ψv+Ψd<br />

Ψ = ½ λ (ln (detF)) 2 - µ ln (detF)+ ½ µ (tr (c) -3) (2.22)<br />

PK stress and Cauchy stress can be easily be computed from (2.22) and given as<br />

S = λ (ln(detF))c -1 + ½ µ (I –c -1 ) (2.25)<br />

σ= λ/deft (ln (detF)) I+ (µ/detF) (b-I) (2.24)<br />

The components <strong>of</strong> the fourth-order elasticity tensor are obta<strong>in</strong>ed as follows:<br />

Cijkl=λcij -1 ckl -1 + (µ - λ ln (detF) (cik -1 cjl -1 - cil -1 ckj -1 ) i,j,k,l=1,2,3 (2.25)<br />

III. FINITE ELEMENT FORMULATION<br />

The f<strong>in</strong>ite element equations are non-l<strong>in</strong>ear <strong>in</strong> terms <strong>of</strong> displacement because <strong>of</strong> nonl<strong>in</strong>ear<br />

stra<strong>in</strong>-displacement and stress-stra<strong>in</strong> relationships. By the pr<strong>in</strong>ciple <strong>of</strong> virtual work, the weak form <strong>in</strong><br />

the current configuration is given as<br />

∫ ∫ ∫ δe T σ dv = ∫ ∫ δu T q dA + ∫ ∫ ∫ δu T b dv (3.1)<br />

V A V<br />

For problems <strong>in</strong>volv<strong>in</strong>g large displacement [6], consideration should be given to the<br />

configuration. S<strong>in</strong>ce the current configuration is not known, we cannot compute directly from this equ<br />

(3.1). Us<strong>in</strong>g the deformation gradient, the derivatives and <strong>in</strong>tegrals <strong>in</strong> the weak form are transformed to<br />

the <strong>in</strong>itial configuration (Total Langrangian formulation).<br />

∫ ∫ ∫S δe dV o = ∫ ∫ δu T q o dA o + ∫ ∫ ∫ δu T b o dv o<br />

V o A o V o<br />

(3.2)

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

Saved successfully!

Ooh no, something went wrong!