26.12.2012 Views

Untitled

Untitled

Untitled

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.

F.E.M. for Prestressed Saint Venant-Kirchhoff Hyperelastic Membranes 129<br />

where λ and µ are known as the Lamé ´ constants. These two constants can be related<br />

to the classical Young modulus E and Poisson’s ratioνas follows:<br />

λ =<br />

νE<br />

(1 + ν)(1 − 2ν)<br />

µ =<br />

E<br />

2(1+ν)<br />

Another important feature which needs to be obtained is the incremental strain<br />

energy accumulated into the structure along the deformation path from the primary<br />

to the secondary states. By performing again a Taylor series expansion truncated<br />

after the second order, the internal strain energy functional per unit of nominally<br />

stressed volume may be developed as:<br />

w ∗ int = wint + ∂wint<br />

(21)<br />

(E<br />

∂EAB<br />

∗ AB −EAB)+ 1 ∂<br />

2<br />

2 wint<br />

(E<br />

∂EAB∂ECD<br />

∗ AB −EAB)(E ∗ CD −ECD) (22)<br />

This above formula can be rewritten as:<br />

where the terms Ω and Υ canbedepictedas:<br />

w ∗ int = wint + Ω + Υ (23)<br />

Ω = SABX pret<br />

X E ∗relat<br />

i,A Xpret<br />

j,B<br />

Υ = 1<br />

pret<br />

CABCDX<br />

i,A 2 Xpret Xj,B E ∗relat<br />

ij X pret<br />

k,C Xpret kl<br />

By substituting (24) and (25) back into (22):<br />

w ∗ int − wint = J[σijE ∗relat<br />

ij<br />

ij<br />

l,D E∗relat<br />

+ 1<br />

2 Cijkl<br />

= JσijE ∗relat<br />

ij<br />

= 1<br />

2 JCijkl C E ∗relat<br />

ij E ∗relat<br />

kl<br />

(24)<br />

(25)<br />

C E ∗relat<br />

ij E ∗relat<br />

kl ]=JΨ (26)<br />

By integrating over the initial undeformed volume V 0 corresponding to the configuration<br />

ℜ0, and by applying the mass conservation principle from this volume V 0<br />

to the prestressed volume V pret corresponding to the configuration ℜt, we obtain<br />

the incremental Helmholtz’s free energy functional, which is given as:<br />

<br />

<br />

<br />

<br />

∆Wint W =<br />

V<br />

(w<br />

0<br />

∗ int−wint)dV =<br />

V<br />

JΨdV =<br />

0<br />

V pret<br />

ΨdV =<br />

V pret<br />

w ∗relat<br />

int dV (27)<br />

Therefore, the internal strain energy functional per unit of volume of the primary<br />

state takes the final form:<br />

w ∗relat<br />

int<br />

= σijE ∗relat<br />

ij<br />

+ 1<br />

2 Cijkl<br />

3 Finite Element Discretization<br />

C E ∗relat<br />

E ∗relat<br />

3.1 From the Strong to the Weak Formulation<br />

The above mentioned primary and secondary states can be understood as an initial<br />

prestressed state ℜ pret and a final in service loading state ℜ due to the consideration<br />

of live and dead load. Henceforth, the coordinates of any body’s particle, in<br />

ij<br />

kl<br />

(28)

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

Saved successfully!

Ooh no, something went wrong!