Material law in beam nonlinear calculation

Hello,

I’m trying to perform a nonlinear analysis of a corbel, but I can’t figure out the material working law used during the analysis.

+PROG AQUA
head 

ECHO FULL EXTR

NORM DC US NDC AASHTO-2014 COUN US UNIT 0 LANG 1 G 10 

STEE no 1 type s scm0 1.0 scm1 1.0
ssla eps ULTI
SSLA eps   SIG     $ EAE 19.5.1
 10.35   301.17
 1.725   300
  0      0
-1.725  -300
-10.35  -301.17

PROF NO 1 TYPE HEB 300

end


+prog sofimshc
head  Corbel

SYST 3D GDIR ZZ GDIV 10000
CTRL MESH   1
CTRL HMIN    0.2

spt no 1 fix f
spt no 2 z -5.0

sln no 1 npa 1 npe 2 sno 1 grp 1

grp no 1 titl 'Corbel'

end

+prog sofiload 
head

lc no 1 titl 'Point load'
node no 2 p1 100

end

+prog ase
head
lc all
end

+PROG ASE 
HEAD Buckling Eigenvalues for Stresses of LC 1

SYST PLC 1 FACV 0.0 $
                  
EIGE 5 BUCK LC 701 $
               
END

+PROG ASE
HEAD with imperfection

OBLI LC 701 FACV 0.21978

ULTI 50 FAK1 0.1 DFAK 1 FAKE -0.0001 PRO 1 PRIM NO DL NO EPSS 10 $ start at FAC1 0.1 to better recognize predeformation in load-defo plot
SYST PROB TH2 iter 200 tol 0.001 TOL4 0.005
LC  200
 LCC 1
NSTR KMOD S1 KSV UL

END

I wonder why the cross section reaches 345 MPa during the non linear analysis, instead of 300 MPa as specified in AQUA. Besides, isn’t the ultimate strain defined in AQUA taken into account during the ULTI analysis? How could I perform an analysis with different ultimate deformations in tension and compression?

Thanks and regards.

Hello

Thank you for your inquiry.
I looked at your code and am also confused about the results. I also think that the maximum stress should be given by the work law.
In order to solve the problem, I forwarded your problem to the responsible developer.
I will get back to you as soon as we have found a solution.

Best regards
Frederik Höller
Your SOFiSTiK Support Team

Thank you very much!

I’m looking forward to your comments.