Soluzione alle differenze finite problema monoassiale nella tensione con cc cinematiche e statiche

Si consideri un corpo cilindro elastico di lunghezza L = 100 mm, modulo elastico E = 210 GPa e modulo di Poisson n = 0.3.
Si consideri un sistema di riferimento x, y, z con z allineato lungo l'asse del cilindro e origine su una base del cilindro.
Le forze di volume b_z (z) = 1 kN/mm³ z/L variano linearmente con z.
Il cilindro e' in stato monoassiale di tensione e' vincolato all'estremita' z = 0 con spostamento imposti w(0) = 0 e e' libero all'estremita' z = L con forze di superficie imposte f_z(L) = 100 MPa.

Figura 1

Figura 2. Configurazione di riferimento (.), configurazione deformata (+), spostamento.

Figura 3. Tensione normale, deformazione trasversale

Codice Matlab per soluzione mediante metodo delle differenze finite