-
.è da un po' che non ci smanetto perché non ho tempo però, in linea di principio, adonis funziona come flac. perciò per capire come creare i modelli di calcolo puoi spulciare i manuali di flac.
io per la tesi ho usato plaxis a me sembra molto simile meno userfriendly , ovviamente mi servirebbe un collega geotecnico (anche se non sono a digiunoin generale) per fare due prove di testing..
. -
.
@Rev sto cominciando a vedere di fare un matching tra califfo e adonis, un caso semplicissimo tipo una fondazione superficiale, ovviamente devo comparare la portanza analitica con il campo elastoplastico fem , hai suggerimenti?
proceduralmente con adonis debbo impostare io la litostatica (almeno credo) poi attivare la fondazione (e gravità) ed azzerare gli spost.. poi applico il carico giusto?. -
.
sì, dopo aver creato il modello del suolo, fai girare il programma a vuoto in modo che vengono calcolati gli spostamenti in assenza di carico.
fatto ciò imposti il risultato ottenuto come condizione iniziale e da lì inserisci e calcoli la fondazione.. -
.
Ho fatto una prova con adonis con i dati dell' esempio n.1 di validazione di califfo (fi 25° c'=20Kpa Eed=50Mpa; G shear=20661157 Pa, E bulk =28735632Pa v=0.2 ) con fondazione posta sulla superficie senza approfondimento da cui ho ricavato la q ult (74 Kpa) .
i risultati sono questi :
il cedimento sotto la fondazione è pari a 2.2 mm Adonis , quello con Califfo 2.1 mm ,quindi sono fortunato
; ancora non ho capito bene il significato dei triangolini sullo stato degli elementi ma sembra in condizioni prevalentemente elastico il modello.
Adesso devo capire come modellare una fondazione posta a quota -x dal pc.
qui in link dell'esempio
https://we.tl/t-MRXlyuuJ6E. -
.
la prossima volta è più conveniente che tu posti il testo dello script. così è più semplice selezionarlo e salvarlo in un file di testo con estensione .ajs.
in ogni caso, ora lo posto io:CODICEnewmodel()
set("language","it")
set("unit","stress-pa")
set("useNMD","on")
rect("startPoint",-5,0,"endPoint",5,-5)
delline("id",1)
line("startPoint",-5.000,0.000,"endPoint",4.5,0)
delline("id",5)
line("startPoint",-5.000,0.000,"endPoint",4.5,0)
delline("id",6)
line("startPoint",-5.000,0.000,"endPoint",4.5,0)
delline("id",7)
line("startPoint",-5.000,0.000,"endPoint",-0.5,0)
line("startPoint",-0.500,0.000,"endPoint",0.5,0)
line("startPoint",0.500,0.000,"endPoint",5.000,0.000)
discretize("maxedge",.2)
triangle("maxedge",.2)
material("create","Mohr-Coulomb","matid",1,"matname","prova sand","density",1800,"shear",20661157,"bulk",28735632,"coh",20e03,"fric",25,"dil",0,"tens",0)
material("assign","matid",1,"region",0.899585,-2.60036)
applybc("xyfix","xlim",-5.102,-4.918,"ylim",-5.117,0.153)
applybc("xyfix","xlim",-5.362,5.224,"ylim",-5.316,-4.933)
applybc("xyfix","xlim",4.964,5.454,"ylim",-5.623,0.169)
var grav = 9.81
var dsyy1 = 1800*grav
var isyy1 = 0.0
var k01 = 0.5
initial("syy",0,"yvar",dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("sxx",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("szz",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
set("gravity",0,9.81)
solve()
initial("xydisp",0)
structure("drawbeam","beamid",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
structure("applynode","yforce",-12333,"xlim",-0.550,0.625,"ylim",-0.064,0.159)
solve()
plot("contour","ydisp"). -
.Adesso devo capire come modellare una fondazione posta a quota -x dal pc.
quando definisci il modello (il tuo rettangolo di terreno, per intenderci), devi definire anche il tuo futuro scavo, ma senza scavare, definendo anche al suo interno il terreno in sito.
dopo aver applicato le tensioni litostatiche ed aver annullato i relativi spostamenti, cancellerai quella parte di terreno che avrai già delimitato sin dall'inizio simulando così lo scavo (ma SENZA far di nuovo girare il modello) e definisci lì dentro la tua fondazione.
noi intanto aspettiamo il nuovo script
. -
+1 .
ho ripulito il modello (perché c'erano alcuni comandi poi annullati) e l'ho commentato.
suggerisco a chi vuole approcciarsi a questo programma di copiare ed eseguire nel software una riga per volta. attenzione, perché il primo comando "material" - per la sua lunghezza - è splittato su due righe, che vanno perciò ricopiate assieme.CODICE// script per il calcolo di una fondazione superficiale nastriforme
// su terreno omogeneo
// crea nuovo modello
newmodel()
// impostazioni generali
set("language","it")
set("unit","stress-pa")
set("useNMD","on")
// definizione modello
rect("startPoint",-5,0,"endPoint",5,-5)
// rimozione superficie
delline("id",1)
// creazione nuova superficie suddivisa in 3 tratti
line("startPoint",-5.000,0.000,"endPoint",-0.5,0)
line("startPoint",-0.500,0.000,"endPoint",0.5,0)
line("startPoint",0.500,0.000,"endPoint",5.000,0.000)
// creazione mesh
discretize("maxedge",.2)
triangle("maxedge",.2)
// definizione materiale
material("create","Mohr-Coulomb","matid",1,"matname","prova sand","density",1800,"shear",20661157,"bulk",28735632,"coh",20e03,"fric",25,"dil",0,"tens",0)
material("assign","matid",1,"region",0.899585,-2.60036)
// definizione vincoli di estremità del modello
applybc("xyfix","xlim",-5.102,-4.918,"ylim",-5.117,0.153)
applybc("xyfix","xlim",-5.362,5.224,"ylim",-5.316,-4.933)
applybc("xyfix","xlim",4.964,5.454,"ylim",-5.623,0.169)
// definizione variabili di calcolo
var grav = 9.81
var dsyy1 = 1800*grav
var isyy1 = 0.0
var k01 = 0.5
// condizioni iniziali
initial("syy",0,"yvar",dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("sxx",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("szz",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
// impostazione gravità e risoluzione per tensioni litostatiche
set("gravity",0,9.81)
solve()
// assunzione delle tensioni litostatiche iniziali quali condizioni di partenza
initial("xydisp",0)
// definizione della fondazione e del carico
structure("drawbeam","beamid",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
structure("applynode","yforce",-12333,"xlim",-0.550,0.625,"ylim",-0.064,0.159)
// risoluzione del modello
solve()
// disegno risultati
plot("contour","ydisp")
faccio notare a volsconger un'incongruenza: se si lavora direttamente con lo script, editandolo con un editor di testo, allora ha senso definire le variabili, che poi saranno richiamate in altri punti dello script.
se invece si inseriscono i dati direttamente dall'interfaccia ciò è superfluo.
capisco però che volsconger abbia adattato alcuni degli esempi presenti nei tutorial o sul sito, ma allora, nella riga dove viene impostata la gravità, la corretta linea di codice avrebbe dovuto essere: set("gravity",0,grav).. -
.
noi intanto aspettiamo il nuovo script spero sia un'incitazione non....
cmq ti ringrazio per le tue doverose sottolineature formali reversi, avrei pure apprezzato un commento sulla bontà del risultato, stante la mia autoironia di aver ottenuto per Caso"" o meglio ad casum tale risultato, che cmq è o dovrebbe essere ottenibile col tuo script.
mi sono cmq accorto che c'è un'errore geometrico quindi devo rimetterci mano.. -
....avrei pure apprezzato un commento sulla bontà del risultato, stante la mia autoironia di aver ottenuto per Caso"" o meglio ad casum tale risultato, che cmq è o dovrebbe essere ottenibile col tuo script.
in primo luogo, lo script non è mio ma è il tuo preso dal file che hai condiviso. puoi aprirlo con qualunque editor di testo, compreso notepad.
in secondo luogo, il contour plot degli spostamenti verticali dice poco e niente per dare un giudizio sulla bontà del risultato.
per esempio, ha più senso quello degli spostamenti totali (orizz. + vert.) che dà maggiormente l'idea di come si allontana il terreno da sotto la fondazione. e qualitativamente sembra corretto. il comando è: plot("contour","totdisp")
da vedere anche le deformazioni. interessanti ad esempio quelle verticali: plot("contour","eyy")
non mi convincono però le tensioni perché dovremmo trovare qualcosa di simile alle deformazioni. secondo me c'è qualcosa da definire nei nodi dove finisce la superficie del terreno e inizia la fondazione. in pratica, quello che adesso è un unico nodo dovrebbe potersi sdoppiare, in modo che la fondazione e il terreno possano avere spostamenti diversi. dovresti provare a giocare un po' con le interfacce (puoi fare riferimento al "manuale", cliccando sulla voce "Aiuti" del menu "Aiuti").
e comunque: seguo con interesse il tuo lavoro e lo apprezzo.. -
+1 .
ho apprezzato i tuoi commenti Rev, ed infatti ho adottato uno script con l'interfaccia, purtroppo come funziona non l'ho ben capito l'help è stringatissimo, i parametri jks e jkn forse sono rigidezze , ma adottando l'interfaccia si risolve il problema che hai notato nelle eyy ma gli spostamenti raddoppiano da circa 2 mm a 5 mm ... non capisco il perchè ,
Cmq poi mo sono voluto divertire a fare una push verticale della fondazioni di prova ed effettivamente sembra tornare il calcolo alla terzaghi della portanza della stessa , il modello è senza intefaccia con moltiplicatori e demoltiplicatori della qult che è pari a 74KPa (nel grafico seguente il pallino rosso).
allego lo script della fondazione con liner +
// script per il calcolo di una fondazione superficiale nastriforme
// su terreno omogeneo con interfaccia
// crea nuovo modello
newmodel()
// impostazioni generali
set("language","it")
set("unit","stress-pa")
set("useNMD","on")
// definizione modello
rect("startPoint",-5,0,"endPoint",5,-5)
// rimozione superficie
delline("id",1)
// creazione nuova superficie suddivisa in 3 tratti
line("startPoint",-5.000,0.000,"endPoint",-0.5,0)
line("startPoint",-0.500,0.000,"endPoint",0.5,0)
line("startPoint",0.500,0.000,"endPoint",5.000,0.000)
// creazione mesh
discretize("maxedge",.2)
triangle("maxedge",.2)
// definizione materiale
material("create","Mohr-Coulomb","matid",1,"matname","prova sand","density",1800,"shear",20661157,"bulk",28735632,"coh",20e03,"fric",25,"dil",0,"tens",0)
material("assign","matid",1,"region",0.899585,-2.60036)
// definizione vincoli di estremità del modello
applybc("xyfix","xlim",-5.102,-4.918,"ylim",-5.117,0.153)
applybc("xyfix","xlim",-5.362,5.224,"ylim",-5.316,-4.933)
applybc("xyfix","xlim",4.964,5.454,"ylim",-5.623,0.169)
// definizione variabili di calcolo
var grav = 9.81
var dsyy1 = 1800*grav
var isyy1 = 0.0
var k01 = 0.5
// condizioni iniziali
initial("syy",0,"yvar",dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("sxx",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("szz",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
// impostazione gravità e risoluzione per tensioni litostatiche
set("gravity",0,9.81)
solve()
// assunzione delle tensioni litostatiche iniziali quali condizioni di partenza
initial("xydisp",0)
// definizione della fondazione e del carico
structure("drawbeam","beamid",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
structure("drawliner","beamid",1,"iftype","oneSide","ifid1",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
imaterial("edit","ifid",1,"jkn",1e+8,"jks",1e+8,"friction",16)
structure("applynode","yforce",-12333,"xlim",-0.550,0.625,"ylim",-0.064,0.159)
// risoluzione del modello
solve()
// disegno risultati
plot("contour","totdisp")
e la mia "push" verticale
. -
+1 ....
allego lo script della fondazione con liner +
...
// definizione della fondazione e del carico
structure("drawbeam","beamid",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
structure("drawliner","beamid",1,"iftype","oneSide","ifid1",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
imaterial("edit","ifid",1,"jkn",1e+8,"jks",1e+8,"friction",16)
structure("applynode","yforce",-12333,"xlim",-0.550,0.625,"ylim",-0.064,0.159)
...
c'è una sovrapposizione tra un elemento beam (comando drawbeam) e un elemento con interfaccia (comando drawliner). il secondo comando include il primo (con, in più, l'aggiunta dell'interfaccia). ho quindi riscritto la sezione quotata, portando in cima drawliner al posto di drawbeam:CODICE// definizione della fondazione e del carico
structure("drawliner","beamid",1,"iftype","oneSide","ifid1",1,"xlim",-0.413,0.434,"ylim",-0.064,0.095)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
imaterial("edit","ifid",1,"jkn",1e+8,"jks",1e+7,"friction",16)
structure("applynode","yforce",-12333,"xlim",-0.550,0.625,"ylim",-0.064,0.159)
nota: ho anche cambiato jks perché, facendo riferimento all'esempio del tutorial 06, il jks è di un ordine di grandezza inferiore a jkn.
per il corretto significato di jks e jkn devo andare a guardare sul manuale di flac. l'unità di misura cmq è Pa/m.. -
.
ecco qua il testo del manuale di flac. CITAZIONEThe normal behavior of the liner-zone interface (see Figure 3) is controlled by the normal coupling spring properties of (1) stiffness per unit area, kn, and (2) tensile strength, ft. The shear behavior of the liner-zone interface (see Figure 4) is cohesive and frictional in nature, and is controlled by the shear coupling spring properties of (1) stiffness per unit area, ks, (2) cohesive strength, c, (3) residual cohesive strength, cr, and (4) friction angle, f, and by the interface normal stress, sn. If the liner fails in tension, then the effective cohesion drops from c to cr, and the tensile strength is set to zero. The total relative normal displacement, un, continues to be tracked such that compressive normal stresses will again develop when the gap closes. Note that the coupling spring properties associated with each liner element are averaged at liner nodes.. -
.
grazie per la correzione rev in effetti ora lo spostamento è 2.4 mm non 5 e siamo vicini alla soluzione alla Fox del Vs Califfo. ottimo.
cmq rimane un buco teorico sugli elementi che usa adonis (lo stato degli elementi post analisi con i triangoli e le crocette ad esempio sono plasticizzazioni?),pure l'interfaccia di flac che tu usi funziona differentemente da adonis che è agli elementi finiti a 3 o 6 nodi, flac fa analisi lagrangina del continuo (ad esempio il grigliato è rettangolare se ricordo bene) non vorrei dire castronerie cmq..
adesso cercherò di continuare con gli esempi di validazione di Califfo. tempo permettendo
Edited by volsconger - 22/11/2021, 19:29. -
.
continuo il post anche se forse è meglio che apro in sezione de inutilia...
grazie agli aiutini di reversi e dell'autore del codice , ho provato il confronto tra califfo e adonis con le seguenti ipotesi:
- un materiale tipo MC con 5MPa di modulo di Young 25° e coesione 5kPa
- fondazione 1m x 1m in c.a. spessa 1 m posta a 1 m dal p.c. modellata con elemento beam con interfaccia
- carico da qult alla Terzaghi ottenuto da Califfo pari a 303 kPa
allego lo script di questa prova
newmodel()
set("language","it")
set("unit","stress-pa")
set("useNMD","on")
rect("startPoint",-5,0,"endPoint",5,-5)
line("startPoint",-.500,0.000,"endPoint",-.5,-1)
line("startPoint",-.500,-1,"endPoint",.5,-1)
line("startPoint",.5,-1,"endPoint",.5,0)
discretize("maxedge",.2)
triangle("maxedge",.2)
material("create","Mohr-Coulomb","matid",1,"matname","sand 25° E=5 MPa","density",1800,"shear",1923076,"bulk",4166666,"coh",3000,"fric",25,"dil",0,"tens",0)
material("assign","matid",1,"region",0.899585,-2.60036)
material("assign","matid",1,"region",0,-0.6)
applybc("xyfix","xlim",-5.102,-4.918,"ylim",-5.117,0.153)
applybc("xyfix","xlim",-5.362,5.224,"ylim",-5.316,-4.933)
applybc("xyfix","xlim",4.964,5.454,"ylim",-5.623,0.169)
// impostazione gravità e risoluzione per tensioni litostatiche
var grav = 9.81
var dsyy1 = 1800*grav
var isyy1 = 0.0
var k01 = 0.5
var qu = 302000
// assunzione delle tensioni litostatiche iniziali quali condizioni di partenza
initial("syy",0,"yvar",dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("sxx",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("szz",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
set("gravity",0,9.81)
solve()
// azzero spost litostatici
initial("xydisp",0)
// scavo fondazione
excavate("region",0.216223,-0.516534,"reset","off")
// CC equilibrio scavo
applybc("yfix","xlim",-0.541,-0.444,"ylim",-1.045,0.084)
applybc("yfix","xlim",0.456,0.541,"ylim",-1.045,0.096)
// azzero spost litostatici
initial("xydisp",0)
// definizione della fondazione e del carico
structure("drawliner","beamid",1,"iftype","oneSide","ifid1",1,"xlim",-0.585,0.567,"ylim",-1.061,-0.953)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
imaterial("edit","ifid",1,"jkn",1e+8,"jks",1e+7,"friction",16)
structure("applynode","yforce",-qu/6,"xlim",-0.585,0.567,"ylim",-1.061,-0.953)
solve()
plot("contour","totdisp")
__________________________________________________________
i risultati in termini di abbassamenti non sono in linea con il metodo fox .. per ora ho notato problemi con i vincoli che servono a sostenere lo scavo
. -
.
ho corretto il problema sopra citato dovuto ai fix sbagliati e ho raffinato la soluzione, rimangono spostamenti molto elevati rispetto a quanto ci dicono le previsioni elastiche alla Fox ( non quello degli oroscopi è... )
allego per chi volesse trovare cavolate, lo script
newmodel()
set("language","it")
set("unit","stress-pa")
set("useNMD","on")
rect("startPoint",-5,0,"endPoint",5,-5)
line("startPoint",-.500,0.000,"endPoint",-.5,-1)
line("startPoint",-.500,-1,"endPoint",.5,-1)
line("startPoint",.5,-1,"endPoint",.5,0)
discretize("maxedge",.2)
segment("id",7,"numedge",8)
segment("id",9,"numedge",8)
segment("id",8,"numedge",8)
triangle("maxedge",.2)
material("create","Mohr-Coulomb","matid",1,"matname","sand 25° E=5 MPa","density",1800,"shear",1923076,"bulk",4166666,"coh",0,"fric",25,"dil",0,"tens",0)
material("assign","matid",1,"region",0.899585,-2.60036)
material("assign","matid",1,"region",0,-0.6)
applybc("xyfix","xlim",-5.102,-4.918,"ylim",-5.117,0.153)
applybc("xyfix","xlim",-5.362,5.224,"ylim",-5.316,-4.933)
applybc("xyfix","xlim",4.964,5.454,"ylim",-5.623,0.169)
// impostazione gravità e risoluzione per tensioni litostatiche
var grav = 9.81
var dsyy1 = 1800*grav
var isyy1 = 0.0
var k01 = 0.5
var qu = 302000*1
// assunzione delle tensioni litostatiche iniziali quali condizioni di partenza
initial("syy",0,"yvar",dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("sxx",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
initial("szz",k01*isyy1,"yvar",k01*dsyy1,"xlim",-5,5,"ylim",-5,0)
set("gravity",0,9.81)
solve()
// azzero spost litostatici
initial("xydisp",0)
// scavo fondazione
excavate("region",0.216223,-0.516534,"reset","off")
// CC equilibrio scavo
applybc("xfix","xlim",-0.536,-0.443,"ylim",-0.945,0.034)
applybc("xfix","xlim",0.434,0.562,"ylim",-0.936,0.060)
// azzero spost litostatici
initial("xydisp",0)
// definizione della fondazione e del carico
structure("drawliner","beamid",1,"iftype","oneSide","ifid1",1,"xlim",-0.585,0.567,"ylim",-1.061,-0.953)
structure("material","beamid",1,"area",1,"I",0.0833,"ymod",3e+10)
imaterial("edit","ifid",1,"jkn",1e+8,"jks",1e+7,"friction",16)
structure("applynode","yforce",-qu/9,"xlim",-0.585,0.567,"ylim",-1.061,-0.953)
solve()
plot("contour","totdisp")
----------------------------------
qui delle prove per vedere l'andamento
.
AdonisSoftware geotecnico |
Contacts