Home GeoForum Forums Geologia & Professione Geotecnica software flac

Salve a tutti!
Sono alle prese con un esercizio che riguarda la filtrazione. ( moti di filtrazione al di sotto di una traversa in c.a). per la risoluzione dell'esercizio non ho nessun problema. L'esercizio però è stato condotto anche con il flac. ( a noi la prof. ci ha fornito anche il listato ) ma per l'esame dovremmo sapere almeno commentare il listato.
Se c'è qualcuno che può aiutarmi, allego il listato.
grazie anticipatamente!

ESERCITAZIONE 4: MOTI DI FILTRAZIONE A REGIME

10.1 LISTATO
STEP 1: DEFINIZIONE GRIGLIA E MODELLO TERRENO
config gwflow
grid 222,40
gen (-40.0,-20.0) (-40.0,0.0) (71.0,0.0) (71.0,-20.0) ratio 1.0,1.0 i=1,223 j=1,41

model elastic
set mech off
set gravity=9.81
model null i 81 142 j 36 40
group 'null' i 81 142 j 36 40
group delete 'null'
model null i 81 j 22 35
group 'null' i 81 j 22 35
group delete 'null'

STEP 2: DEFINIZIONE DEI VINCOLI CINEMATICI E IDRAULICI
fix x i 1
fix x i 223
fix x y j 1
fix saturation
initial saturation 1.0
prop por=0.5 perm=2.0E-10 notnull
apply discharge 0.0 from 1,41 to 1,1
apply discharge 0.0 from 1,1 to 223,1
apply discharge 0.0 from 223,41 to 223,1
apply discharge 0.0 from 81,41 to 143,41
apply pp 100000.0 from 1,41 to 81,41
apply pp 0.0 from 143,41 to 223,41
water bulk=1.0
water density=1000.0
solve

Non è una domanda banale, che si dice sul blog degli ingegneri?

provo a commentare

10.1 LISTATO
STEP 1: DEFINIZIONE GRIGLIA E MODELLO TERRENO
config gwflow (filtrazione)
grid 222,40 (queste sono le celle del tuo dominio 222 direzione x 40 direzione y modello plain strain
gen (-40.0,-20.0) (-40.0,0.0) (71.0,0.0) (71.0,-20.0) ratio 1.0,1.0 i=1,223 j=1,41 (coordinate del dominio, rapporto lunghezza lato celle in questo caso quadrate e numerazione celle secondo i e j sistema coordinate delle celle, risulta evidente che le tue celle sono di 0.5 m di lato)

model elastic (modello costitutivo scelto elastico credo tu abbia impostato almeno un E e un poisson)
set mech off (? qui non so, forse significa che non hai impostato i parametri elastici perché comunque nel listato non li vedo)
set gravity=9.81 (gravita da assegnare alle geometrie)
model null i 81 142 j 36 40 (qui crei la geometria della traversa)
group 'null' i 81 142 j 36 40 (probabilmente disattivando il dominio non so esattamente come lavora Flac)
group delete 'null'
model null i 81 j 22 35 (idem mi pare che venga uno schema a briglia con due taglioni a monte e a valle di lunghezza differente)
group 'null' i 81 j 22 35
group delete 'null'

STEP 2: DEFINIZIONE DEI VINCOLI CINEMATICI E IDRAULICI
fix x i 1 (settaggio vincoli al bordo laterale)
fix x i 223 (altro bordo laterale)
fix x y j 1 (vincolo alla base)
fix saturation
initial saturation 1.0 (dominio saturo?)
prop por=0.5 perm=2.0E-10 notnull (fissi i parametri di porosità, permeabilità)
apply discharge 0.0 from 1,41 to 1,1 (fissi i vincoli di impermeabilità al flusso del modello in questo caso lato)
apply discharge 0.0 from 1,1 to 223,1 (base)
apply discharge 0.0 from 223,41 to 223,1 (altro lato)
apply discharge 0.0 from 81,41 to 143,41 (sotto la traversa)
apply pp 100000.0 from 1,41 to 81,41 (applichi una pressione dell'acqua per generare una differenza di carico idraulico nel tratto di dominio a monte della traversa, di fatto simuli un battente sul paramento che non si vede)
apply pp 0.0 from 143,41 to 223,41 (sul tratto a valle la pressione dell'acqua è ovviamente 0 in superficie)
water bulk=1.0
water density=1000.0
solve

Credo che sia molto importante che la tua insegnante ti faccia fare sforzi per riuscire a comprendere come un software di questo tipo riesca a risolvere un problema geotecnico con una risoluzione di carattere numerico, anzi è fondamentale. Toglimi una curiosità studi geologia o ingegneria? Se a geologia insegnano queste cose non tutto è perduto per la nostra professione.

Un saluto

grazie sei stao super!!!

comunque studio ingegneria civile. si tratta di un esame a scelta, la prof non poteva farci imparare il software, aveva un programma ( molto bello ). ha lasciato però una porta aperta a chi fosse interessato per una tesi con l'utilizzo flac

In bocca al lupo per i tuoi studi.

10.1 LISTATO
STEP 1: DEFINIZIONE GRIGLIA E MODELLO TERRENO
config gwflow (filtrazione) giusto, configuri la griglia e i settaggi del programma per simulare il flusso
grid 222,40 (queste sono le celle del tuo dominio 222 direzione x 40 direzione y modello plain strain
gen (-40.0,-20.0) (-40.0,0.0) (71.0,0.0) (71.0,-20.0) ratio 1.0,1.0 i=1,223 j=1,41 (coordinate del dominio, rapporto lunghezza lato celle in questo caso quadrate e numerazione celle secondo i e j sistema coordinate delle celle, risulta evidente che le tue celle sono di 0.5 m di lato)

model elastic (modello costitutivo scelto elastico credo tu abbia impostato almeno un E e un poisson) assegnati di default, non ha bisogno di specificare in quanto dopo esclude i calcoli meccanici
set mech off (? qui non so, forse significa che non hai impostato i parametri elastici perché comunque nel listato non li vedo) non simula il comportamento meccanico ma solo il moto di filtrazione : mech off
set gravity=9.81 (gravita da assegnare alle geometrie) direi al modello
model null i 81 142 j 36 40 (qui crei la geometria della traversa)
group 'null' i 81 142 j 36 40 (probabilmente disattivando il dominio non so esattamente come lavora Flac)
group delete 'null'
model null i 81 j 22 35 (idem mi pare che venga uno schema a briglia con due taglioni a monte e a valle di lunghezza differente)
group 'null' i 81 j 22 35
group delete 'null'

STEP 2: DEFINIZIONE DEI VINCOLI CINEMATICI E IDRAULICI
fix x i 1 (settaggio vincoli al bordo laterale) fissa in x il modello occero e bloccato rispetto a deformazioni laterali nella colonna i=1 ovvero quella a sinistra
fix x i 223 (altro bordo laterale) fissa in x il modello occero e bloccato rispetto a deformazioni laterali nella colonna i=1 ovvero quella a destra
fix x y j 1 (vincolo alla base) fissa x e y alla base, ovvero nessun grado di libertà
fix saturation fissa la saturazione
initial saturation 1.0 (dominio saturo?) assegna saturazione inziale completa
prop por=0.5 perm=2.0E-10 notnull (fissi i parametri di porosità, permeabilità)
apply discharge 0.0 from 1,41 to 1,1 (fissi i vincoli di impermeabilità al flusso del modello in questo caso lato) assegna una portata pari a 0 sul lato sinistro
apply discharge 0.0 from 1,1 to 223,1 (base)
assegna una portata pari a 0 sulla base
apply discharge 0.0 from 223,41 to 223,1 (altro lato) idem
apply discharge 0.0 from 81,41 to 143,41 (sotto la traversa) idem
apply pp 100000.0 from 1,41 to 81,41 (applichi una pressione dell'acqua per generare una differenza di carico idraulico nel tratto di dominio a monte della traversa, di fatto simuli un battente sul paramento che non si vede)
apply pp 0.0 from 143,41 to 223,41 (sul tratto a valle la pressione dell'acqua è ovviamente 0 in superficie)
water bulk=1.0 modulo di deformazione volumetrica dell'acqua, basso perchè serve a velocizzare il calcolo
water density=1000.0 densità acqua pari a 1000 kg/mc
solve

Anche in geologia viene utilizzato, e insegnato! FLAC è un software un po' macchinoso ma altrettanto geniale, i suoi usi sono i più disparati. Io per esempio l'ho utilizzato molto per per le modellazioni tenso-deformative dei pendii instabili (corso di Stabilità dei Versanti del prof. Salvatore Martino alla Sapienza).