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)

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