Users Forum

ADONIS Category => ADONIS Users => Topic started by: R.abed on January 18, 2025, 03:40:08 pm

Title: radial tension should start by zero stress why it doesnt?
Post by: R.abed on January 18, 2025, 03:40:08 pm

Hello Mr. Roozbeh,

i need you help again!
My question ist: why the radial (the tangential also) ist not beginning at zero and has big difference to the analytical solution. We now, at the first it should be zero, but its not. Here you can see the picture.
Here is the script:
set("gravity",0,9.81)
rect("startPoint",-50,-50,"endPoint",50,50)
circle("centerPoint",0,0,"radius",5,"numSeg",20)
discretize("maxedge",1.5)
gmsh("maxedge",1.5,"elemtype","T3","useNMD","on")
applybc("xyfix","xlim",-50,50,"ylim",-50,50)
initial("sxx",-50e6)
initial("syy",-50e6)
initial("szz",-50e6)
material("create","IsoElastic","matid",1,"matname","Material 1","density",2000,"shear",1.92308e+08,"bulk",4.16667e+08)
material("assign","matid",1,"region",36.0031,-3.52607)
excavate("region",0.927914,0.185583)
solve()
plot("contour","sig1")
queryplot("startPoint",5,0,"endPoint",50,0,"numInterval",25,"plName","query1","plID",1)
plot("chart","table","id1",1)
plot("contour","sig3")
queryplot("startPoint",5,0,"endPoint",50,0,"numInterval",25,"plName","query2","plID",2)

Title: Re: radial tension should start by zero stress why it doesnt?
Post by: Roozbeh on January 18, 2025, 08:22:08 pm

You are currently plotting the principal stresses using the provided script, which aligns with the results obtained from other programs.