TESI DI DOTTORATO Modellazione e analisi non lineare - LabMec ...

Capitolo 6. Metodo di analisi
START
INPUT
Dati geometrici e meccanici
Condizioni iniziali: {u1}, {ε1},{f1}
LSW = 3
NSTEP = 0
LOADS
Condizione di carico iniziale: {p1}
ASSEM
Assemblaggio della matrice di
rigidezza elastica [KE]
SOLVE
Fattorizzazione della matrice [KE]
([H] = [KE] -1 )
Stampa dei parametri dello stato
tensionale e deformativo corrente
(NSTEP, {u1}, {p1}, {s1}, {ε1},{f1})
LSW = LSW - 4
{u0} = {u1}
{p0} = {p1}
LOADS
Incremento di carico {Δp}
({p1} = {p1} + {Δp})
Processo di
carico
completato?
SI
Stampa dei parametri dello stato
tensionale e deformativo finale
(NSTEP, {u1}, {p1}, {s1}, {ε1},{f1},
{εmax}, {εmin},{fmax} ,{fmin})
STOP
NO
A
154
A
NSTEP = NSTEP + 1
Aggiornamento
fattorizzazione della
matrice di rigidezza?
ASSEM
Assemblaggio della matrice di
rigidezza tangente [KT]
SOLVE
Fattorizzazione della matrice [KT]
([H] = [KT] -1 )
STRUCT
Calcolo reazione strutturale
{s1 (k) } = {s1(u1 (k) )}
> 4
LSW
< 4
{r (k) } = {s1 (k) } - {p1}
SOLVE
Aggiornamento dello spostamento
{u1 (k+1) } = {u1 (k) } + [H] {r (k) }
||{r (k) }|| < ftol ?
LSW = 4
Figura 6.20. Diagramma di flusso del programma di calcolo automatico.
SI
SI
NO
NO