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
