22.11.2014 Views

Download full text - ELSA - Europa

Download full text - ELSA - Europa

Download full text - ELSA - Europa

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.

Geometric non-linearities (3)<br />

Set up following incrementally objective scheme to update<br />

the Cauchy stress (2D case for simplicity) in three phases:<br />

• Let α be the angle of rotation over ∆ t and let θ = α/2<br />

n+<br />

1/2<br />

tan θ = ( ∆t/ 2) ⋅W<br />

12<br />

1. Apply first half of the rotation increment:<br />

n* n T<br />

⎡ cosθ<br />

sinθ⎤<br />

σ = Rσ<br />

R with R= ⎢<br />

−sinθ<br />

cosθ<br />

⎥<br />

⎣<br />

⎦<br />

2. Apply the constitutive equation:<br />

σ<br />

= σ + C⋅∆t⋅D<br />

( n + 1)* n* n + 1/2<br />

3. Apply second half of the rotation increment:<br />

σ<br />

= Rσ<br />

R<br />

n+ 1 ( n+<br />

1)* T<br />

27<br />

Geometric non-linearities (4)<br />

For structural elements (bars, beams, shells) use co-rotational<br />

formulation:<br />

• The stress is measured in a reference frame that rotates with the element<br />

• This greatly simplifies the stress increment procedure: the stress may be<br />

incremented directly by applying the constitutive law.<br />

Example (bar element). In the longitudinal direction:<br />

∆L<br />

∆ ε = → ∆σ<br />

L<br />

∆L L dL L<br />

∑∆ ε = ∑ ln<br />

L<br />

∫ =<br />

L0<br />

L L<br />

A small-strain formulation would be:<br />

∆L L−L<br />

∆ ε = → ∑ ∆ ε =<br />

L<br />

L<br />

0 0<br />

0<br />

0<br />

y<br />

L<br />

x<br />

(Natural or logarithmic strain)<br />

(Engineering strain)<br />

∆L<br />

28<br />

σε ,<br />

14

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

Saved successfully!

Ooh no, something went wrong!