Soluzione alle differenze finite problema monoassiale nella deformazione con cc cinematiche

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 deformazione ed e' vincolato alle estremita' con spostamenti imposti w(0) = 0 e w(L) = w_L = 3 mm.

Figura 1

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

Figura 3. Tensione normale, tensione trasversale

Codice Matlab per soluzione mediante metodo delle differenze finite