Hello, I'm trying to model a culvert. The idea is to check the Axial, Moment and Shear in the structure at different backfill stages.
The script below is in the final condition but there are already some lines for the phases.
Even in the final stage only I'm struggling. Also I tried to excavate all the material and to check the arch only since I know the results I should get. But still no convercenge (or no correct effort laws)
MODEL IN THE FINAL CONDITION
line("startPoint",13.25,0,"endPoint",30,0)
line("startPoint",30,0,"endPoint",30,37)
line("startPoint",30,37,"endPoint",-30,37)
line("startPoint",-30,37,"endPoint",-30,0)
line("startPoint",-30,0,"endPoint",-13.25,0)
line("startPoint",-13.25,0,"endPoint",-11.925,1.68743549710113)
line("startPoint",-11.925,1.68743549710113,"endPoint",-10.6,2.99839811368318)
line("startPoint",-10.6,2.99839811368318,"endPoint",-9.275,4.0465209209307)
line("startPoint",-9.275,4.0465209209307,"endPoint",-7.95,4.89125454233673)
line("startPoint",-7.95,4.89125454233673,"endPoint",-6.625,5.56792546719726)
line("startPoint",-6.625,5.56792546719726,"endPoint",-5.3,6.09908589562254)
line("startPoint",-5.3,6.09908589562254,"endPoint",-3.975,6.49966388605202)
line("startPoint",-3.975,6.49966388605202,"endPoint",-2.65,6.77958212850985)
line("startPoint",-2.65,6.77958212850985,"endPoint",-1.325,6.94518046749565)
line("startPoint",-1.325,6.94518046749565,"endPoint",0,7)
line("startPoint",0,7,"endPoint",1.325,6.94518046749565)
line("startPoint",1.325,6.94518046749565,"endPoint",2.65,6.77958212850985)
line("startPoint",2.65,6.77958212850985,"endPoint",3.975,6.49966388605202)
line("startPoint",3.975,6.49966388605202,"endPoint",5.3,6.09908589562254)
line("startPoint",5.3,6.09908589562254,"endPoint",6.625,5.56792546719726)
line("startPoint",6.625,5.56792546719726,"endPoint",7.95,4.89125454233673)
line("startPoint",7.95,4.89125454233673,"endPoint",9.275,4.0465209209307)
line("startPoint",9.275,4.0465209209307,"endPoint",10.6,2.99839811368318)
line("startPoint",10.6,2.99839811368318,"endPoint",11.925,1.68743549710113)
line("startPoint",11.925,1.68743549710113,"endPoint",13.25,0)
line("startPoint",-30,1,"endPoint",-12.4647848007961,1);line("startPoint",30,1,"endPoint",12.4647848007961,1)
line("startPoint",-30,2,"endPoint",-11.6090886527941,2);line("startPoint",30,2,"endPoint",11.6090886527941,2)
line("startPoint",-30,3,"endPoint",-10.597974951642,3);line("startPoint",30,3,"endPoint",10.597974951642,3)
line("startPoint",-30,4,"endPoint",-9.33381011252398,4);line("startPoint",30,4,"endPoint",9.33381011252398,4)
line("startPoint",-30,5,"endPoint",-7.73706380589127,5);line("startPoint",30,5,"endPoint",7.73706380589127,5)
line("startPoint",-30,6,"endPoint",-5.54717355562253,6);line("startPoint",30,6,"endPoint",5.54717355562253,6)
line("startPoint",-30,6.8,"endPoint",-2.4866307301744,6.
;line("startPoint",30,6.8,"endPoint",2.4866307301744,6.
set("unit","stress-kpa")
discretize("size","auto")
gmsh("size","auto","elemtype","T3","useNMD","off")
material("create","Mohr-Coulomb","matid",1,"matname","Material 1","density",2500,"shear",2.8e+09,"bulk",3.9e+09,"coh",3.45e+06,"fric",30,"tens",1e+10)
material("assign","matid",1,"region",-20.6954,0.646731)
material("assign","matid",1,"region",-20.7672,1.65276)
material("assign","matid",1,"region",-20.7672,2.37135)
material("assign","matid",1,"region",-20.8391,3.44923)
material("assign","matid",1,"region",-20.7672,4.52712)
material("assign","matid",1,"region",-20.6954,5.31757)
material("assign","matid",1,"region",-20.6954,6.39545)
material("assign","matid",1,"region",-20.0487,8.98237)
material("assign","matid",1,"region",9.77282,6.53917)
material("assign","matid",1,"region",10.4914,5.605)
material("assign","matid",1,"region",11.7849,4.3834)
material("assign","matid",1,"region",13.5095,3.30551)
material("assign","matid",1,"region",14.8748,2.37135)
material("assign","matid",1,"region",15.8808,1.72462)
material("assign","matid",1,"region",16.8869,0.431154)
structure("drawbeam","beamid",1,"xlim",-13.502,13.604,"ylim",-2.030,7.716)
structure("material","beamid",1,"area",0.5,"I",.010417,"ymod",2e+11,"dens",2500)
applybc("xfix","xlim",29.361,35.190,"ylim",-1.871,38.860)
applybc("xfix","xlim",-34.686,-29.505,"ylim",-1.439,40.443)
applybc("xyfix","xlim",-30.590,-13.763,"ylim",-2.441,0.260)
applybc("xyfix","xlim",14.126,31.576,"ylim",-3.012,0.467)
structure("applynode","pin","xlim",-13.508,-12.609,"ylim",-0.591,0.205)
structure("applynode","pin","xlim",12.994,13.559,"ylim",-0.436,0.200)
set("gravity",0,9.81)
solve()
MODEL ARCH ONLY
line("startPoint",13.25,0,"endPoint",30,0)
line("startPoint",30,0,"endPoint",30,37)
line("startPoint",30,37,"endPoint",-30,37)
line("startPoint",-30,37,"endPoint",-30,0)
line("startPoint",-30,0,"endPoint",-13.25,0)
line("startPoint",-13.25,0,"endPoint",-11.925,1.68743549710113)
line("startPoint",-11.925,1.68743549710113,"endPoint",-10.6,2.99839811368318)
line("startPoint",-10.6,2.99839811368318,"endPoint",-9.275,4.0465209209307)
line("startPoint",-9.275,4.0465209209307,"endPoint",-7.95,4.89125454233673)
line("startPoint",-7.95,4.89125454233673,"endPoint",-6.625,5.56792546719726)
line("startPoint",-6.625,5.56792546719726,"endPoint",-5.3,6.09908589562254)
line("startPoint",-5.3,6.09908589562254,"endPoint",-3.975,6.49966388605202)
line("startPoint",-3.975,6.49966388605202,"endPoint",-2.65,6.77958212850985)
line("startPoint",-2.65,6.77958212850985,"endPoint",-1.325,6.94518046749565)
line("startPoint",-1.325,6.94518046749565,"endPoint",0,7)
line("startPoint",0,7,"endPoint",1.325,6.94518046749565)
line("startPoint",1.325,6.94518046749565,"endPoint",2.65,6.77958212850985)
line("startPoint",2.65,6.77958212850985,"endPoint",3.975,6.49966388605202)
line("startPoint",3.975,6.49966388605202,"endPoint",5.3,6.09908589562254)
line("startPoint",5.3,6.09908589562254,"endPoint",6.625,5.56792546719726)
line("startPoint",6.625,5.56792546719726,"endPoint",7.95,4.89125454233673)
line("startPoint",7.95,4.89125454233673,"endPoint",9.275,4.0465209209307)
line("startPoint",9.275,4.0465209209307,"endPoint",10.6,2.99839811368318)
line("startPoint",10.6,2.99839811368318,"endPoint",11.925,1.68743549710113)
line("startPoint",11.925,1.68743549710113,"endPoint",13.25,0)
line("startPoint",-30,1,"endPoint",-12.4647848007961,1);line("startPoint",30,1,"endPoint",12.4647848007961,1)
line("startPoint",-30,2,"endPoint",-11.6090886527941,2);line("startPoint",30,2,"endPoint",11.6090886527941,2)
line("startPoint",-30,3,"endPoint",-10.597974951642,3);line("startPoint",30,3,"endPoint",10.597974951642,3)
line("startPoint",-30,4,"endPoint",-9.33381011252398,4);line("startPoint",30,4,"endPoint",9.33381011252398,4)
line("startPoint",-30,5,"endPoint",-7.73706380589127,5);line("startPoint",30,5,"endPoint",7.73706380589127,5)
line("startPoint",-30,6,"endPoint",-5.54717355562253,6);line("startPoint",30,6,"endPoint",5.54717355562253,6)
line("startPoint",-30,6.8,"endPoint",-2.4866307301744,6.
;line("startPoint",30,6.8,"endPoint",2.4866307301744,6.
set("unit","stress-kpa")
discretize("size","auto")
gmsh("size","auto","elemtype","T3","useNMD","off")
material("create","Mohr-Coulomb","matid",1,"matname","Material 1","density",2500,"shear",2.8e+09,"bulk",3.9e+09,"coh",3.45e+06,"fric",30,"tens",1e+10)
material("assign","matid",1,"region",-20.6954,0.646731)
material("assign","matid",1,"region",-20.7672,1.65276)
material("assign","matid",1,"region",-20.7672,2.37135)
material("assign","matid",1,"region",-20.8391,3.44923)
material("assign","matid",1,"region",-20.7672,4.52712)
material("assign","matid",1,"region",-20.6954,5.31757)
material("assign","matid",1,"region",-20.6954,6.39545)
material("assign","matid",1,"region",-20.0487,8.98237)
material("assign","matid",1,"region",9.77282,6.53917)
material("assign","matid",1,"region",10.4914,5.605)
material("assign","matid",1,"region",11.7849,4.3834)
material("assign","matid",1,"region",13.5095,3.30551)
material("assign","matid",1,"region",14.8748,2.37135)
material("assign","matid",1,"region",15.8808,1.72462)
material("assign","matid",1,"region",16.8869,0.431154)
structure("drawbeam","beamid",1,"xlim",-13.502,13.604,"ylim",-2.030,7.716)
structure("material","beamid",1,"area",0.5,"I",.010417,"ymod",2e+11,"dens",2500)
applybc("xfix","xlim",29.361,35.190,"ylim",-1.871,38.860)
applybc("xfix","xlim",-34.686,-29.505,"ylim",-1.439,40.443)
applybc("xyfix","xlim",-30.590,-13.763,"ylim",-2.441,0.260)
applybc("xyfix","xlim",14.126,31.576,"ylim",-3.012,0.467)
structure("applynode","pin","xlim",-13.508,-12.609,"ylim",-0.591,0.205)
structure("applynode","pin","xlim",12.994,13.559,"ylim",-0.436,0.200)
set("gravity",0,9.81)
excavate("region",-24.1273,0.709627,"reset","off")
excavate("region",-24.1273,1.52063,"reset","off")
excavate("region",-24.0259,2.93988,"reset","off")
excavate("region",-24.2287,3.64951,"reset","off")
excavate("region",-24.2287,4.46051,"reset","off")
excavate("region",-24.2287,5.57564,"reset","off")
excavate("region",-24.2287,6.38664,"reset","off")
excavate("region",-22.6067,9.63065,"reset","off")
excavate("region",19.2613,6.38664,"reset","off")
excavate("region",19.2613,5.77839,"reset","off")
excavate("region",19.3627,4.46051,"reset","off")
excavate("region",19.3627,3.75088,"reset","off")
excavate("region",19.5654,2.83851,"reset","off")
excavate("region",19.5654,1.622,"reset","off")
excavate("region",19.5654,0.506876,"reset","off")
Culvert(Read 15 times)