Hello Roozbeh,

Thanks a lot for developing such tools! Previously I tried HYRCAN, and it helped a lot! Recently I have been trying ADONIS. I met a small problem during the analysis.
During one analysis, I used 'solve fos' and it says that the 'fos is less than the lower bound limit 0.5'. Even if I changed the lower bound to 0.001, it still did not work. Then I checked the total displacement and found the displacement only showed at one triangle, as illustrated in the attached image.

I noticed that the geometry was unreasonable at the acute triangle. I had another trial without this acute triangle, and it worked well. However, I still wonder why the displacement only shows in one triangle, without any displacement in other areas of the geometry. Does it mean non-convergence? If yes, is there any notification about non-convergence?

Here is the script. Thank you very much!

Code: [Select]

newmodel()
line('startPoint', 0.000000,-1.000000, 'endPoint', 0.000000,-12.836585)
line('startPoint', 0.000000,-12.836585, 'endPoint', -29.775136,-12.836585)
line('startPoint', -29.775136,-12.836585, 'endPoint', -29.775136,6.434160)
line('startPoint', -29.775136,6.434160, 'endPoint', -29.250000,6.421262)
line('startPoint', -29.250000,6.421262, 'endPoint', -28.740000,6.411262)
line('startPoint', -28.740000,6.411262, 'endPoint', -28.680000,6.441262)
line('startPoint', -28.680000,6.441262, 'endPoint', -28.490000,6.371262)
line('startPoint', -28.490000,6.371262, 'endPoint', -27.170000,6.261262)
line('startPoint', -27.170000,6.261262, 'endPoint', -24.775136,5.986296)
line('startPoint', -24.775136,5.986296, 'endPoint', -24.590117,6.109642)
line('startPoint', -24.590117,6.109642, 'endPoint', -24.395117,5.459642)
line('startPoint', -24.395117,5.459642, 'endPoint', -23.635117,5.459642)
line('startPoint', -23.635117,5.459642, 'endPoint', -23.440117,6.109642)
line('startPoint', -23.440117,6.109642, 'endPoint', -22.940117,6.259642)
line('startPoint', -22.940117,6.259642, 'endPoint', -20.940117,6.310451)
line('startPoint', -20.940117,6.310451, 'endPoint', -20.558233,6.292758)
line('startPoint', -20.558233,6.292758, 'endPoint', -20.182858,6.220315)
line('startPoint', -20.182858,6.220315, 'endPoint', -19.821766,6.094745)
line('startPoint', -19.821766,6.094745, 'endPoint', -19.482436,5.918648)
line('startPoint', -19.482436,5.918648, 'endPoint', -19.171896,5.695670)
line('startPoint', -19.171896,5.695670, 'endPoint', -18.385091,5.040000)
line('startPoint', -18.385091,5.040000, 'endPoint', -17.385091,5.000000)
line('startPoint', -17.385091,5.000000, 'endPoint', -17.024977,4.699905)
line('startPoint', -17.024977,4.699905, 'endPoint', -12.646792,1.051418)
line('startPoint', -12.646792,1.051418, 'endPoint', -11.950000,0.056000)
line('startPoint', -11.950000,0.056000, 'endPoint', -11.450000,0.036000)
line('startPoint', -11.450000,0.036000, 'endPoint', -11.250000,-0.434000)
line('startPoint', -11.250000,-0.434000, 'endPoint', -10.910251,-1.232410)
line('startPoint', -10.910251,-1.232410, 'endPoint', -9.700000,-1.184000)
line('startPoint', -9.700000,-1.184000, 'endPoint', -9.200000,-1.164000)
line('startPoint', -9.200000,-1.164000, 'endPoint', -8.700000,-1.154000)
line('startPoint', -8.700000,-1.154000, 'endPoint', -5.100000,-1.082000)
line('startPoint', -5.100000,-1.082000, 'endPoint', -1.500000,-1.010000)
line('startPoint', -1.500000,-1.010000, 'endPoint', -1.000000,-1.000000)
line('startPoint', -1.000000,-1.000000, 'endPoint', 0.000000,-1.000000)

line('startPoint', -29.775136,2.844160, 'endPoint', -29.250000,2.831262)
line('startPoint', -29.250000,2.831262, 'endPoint', -28.740000,2.821262)
line('startPoint', -28.740000,2.821262, 'endPoint', -28.680000,2.851262)
line('startPoint', -28.680000,2.851262, 'endPoint', -28.490000,2.781262)
line('startPoint', -28.490000,2.781262, 'endPoint', -27.170000,2.671262)
line('startPoint', -27.170000,2.671262, 'endPoint', -24.470000,2.361262)
line('startPoint', -24.470000,2.361262, 'endPoint', -23.270000,2.661262)
line('startPoint', -23.270000,2.661262, 'endPoint', -18.940000,2.771262)
line('startPoint', -18.940000,2.771262, 'endPoint', -15.500000,1.111262)
line('startPoint', -15.500000,1.111262, 'endPoint', -15.340000,1.111262)
line('startPoint', -15.340000,1.111262, 'endPoint', -15.330000,1.111262)
line('startPoint', -15.330000,1.111262, 'endPoint', -15.300000,1.111262)
line('startPoint', -15.300000,1.111262, 'endPoint', -15.260000,1.111262)
line('startPoint', -15.260000,1.111262, 'endPoint', -12.646792,1.051418)

line('startPoint', -29.775136,-0.755840, 'endPoint', -29.250000,-0.768738)
line('startPoint', -29.250000,-0.768738, 'endPoint', -28.740000,-0.778738)
line('startPoint', -28.740000,-0.778738, 'endPoint', -28.680000,-0.748738)
line('startPoint', -28.680000,-0.748738, 'endPoint', -28.490000,-0.818738)
line('startPoint', -28.490000,-0.818738, 'endPoint', -27.170000,-0.928738)
line('startPoint', -27.170000,-0.928738, 'endPoint', -24.470000,-1.238738)
line('startPoint', -24.470000,-1.238738, 'endPoint', -23.270000,-0.938738)
line('startPoint', -23.270000,-0.938738, 'endPoint', -18.940000,-0.828738)
line('startPoint', -18.940000,-0.828738, 'endPoint', -15.500000,-2.488738)
line('startPoint', -15.500000,-2.488738, 'endPoint', -15.340000,-2.488738)
line('startPoint', -15.340000,-2.488738, 'endPoint', -15.330000,-2.488738)
line('startPoint', -15.330000,-2.488738, 'endPoint', -15.300000,-2.488738)
line('startPoint', -15.300000,-2.488738, 'endPoint', -15.260000,-2.488738)
line('startPoint', -15.260000,-2.488738, 'endPoint', -2.160000,-2.788738)
line('startPoint', -2.160000,-2.788738, 'endPoint', -0.830000,-2.938738)
line('startPoint', -0.830000,-2.938738, 'endPoint', 0.000000,-2.984364)

triangle('elemtype', 'T3')
discretize('auto')
triangle('auto')

material('create','Mohr-Coulomb','matid',1, 'matname','Material Layer 1', 'density', 1950,'shear', 2.8e+09,'bulk', 3.9e+09,'coh', 0,'fric', 38,'dil', 0,'tens', 0)
material('create','Mohr-Coulomb','matid',2, 'matname','Material Layer 2', 'density', 1950,'shear', 2.8e+09,'bulk', 3.9e+09,'coh', 5300,'fric', 23,'dil', 0,'tens', 0)
material('create','Mohr-Coulomb','matid',3, 'matname','Material Layer 3', 'density', 1950,'shear', 2.8e+09,'bulk', 3.9e+09,'coh', 7200,'fric', 20,'dil', 0,'tens', 0)

material('assign', 'matid', 1, 'region', -29.075136,3.544160)
material('assign', 'matid', 2, 'region', -29.075136,-0.055840)
material('assign', 'matid', 3, 'region', -29.075136,-12.136585)

watertable('add', 'dens', 1000, 'elev', -7.665096)
applybc('xyfix', 'xlim', -0.100000, 0.100000, 'ylim', -13.836585, 0.000000)
applybc('xyfix', 'xlim', -29.875136, 0.100000, 'ylim', -12.936585, -12.736585)
applybc('xyfix', 'xlim', -29.875136, -29.675136, 'ylim', -13.836585, 7.434160)

set('gravity', 0, 9.8)
solve('fos', 'fosLBLimit', 0.5, 'fosUBLimit', 32.0)


In addition, I also wonder what the difference is between "solve" and "solve fos" except the latter one gives the fos value. Once, I tried to use 'solve' and 'solve fos' on the same model, but the total displacement plots look very different.
Thank you very much once again!

Hi, I'm glad to see you find my programs useful.

regarding your model, when I ran itl with "solve" I realized there is an sharp edge on the geometry which cause the issue (as shown in the attached image). I smoothed the sharp edge and reran the model with "solve fos" and it seems to be working. FOS = 1.014 for the revised geometry. The revised script is attached.

"solve fos" already has "solve" in it so no need to solve the model before fos calculation. Calculation is not triggered when "solve fos" is used because initial fos is out of the specified range. That's why zero displacement is reported. But when you use solve model get cycled without any shear strength reduction consideration hence displacement is nonzero.

thanks for your comments.

-Roozbeh