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
