Cross sectional analysis Aqua/Aqb

Hello Guys

I am trying to do an elastic/plastic analysis of a given cross section with aqua and aqb. i am just trying to get the correct stress - state for given acting forces in aqb.

however, as far as i understand it - the full plastic strain state is not calculated correctly in my file. you can see it as the areas above and below the neutral axis are not equal but the cross section is nearly in full plastic mode. hence there has to be an acting normal force which is not in equlibrium with the forces given by me.

i am suspicious that sofistik is calculating something in regards to “effective cross sectional depth” with steel profiles. however i am not sure if thats the problem and if it is like that how to solve it.

maybe someone can help me.
best regards

code below:

+prog aqua urs:1
KOPF U-Profil Test
ECHO VOLL EXTR

STEU STYP FEMX !SAVE THE MESH
STEU HDIV 0.5[-] !REFINE MESH BY A FACTOR
NORM DC EN 1993-2005
UNIT TYPE 7

STO#fy 636
STO#eniro 145000
STAH NR 1 ART B FY #fy FT #fy ES #eniro EPSY #fy/#eniro*1000 EPST 1000 SCM 1.0

$ define U-Profil
QNR 1 MNR 1 LTAU 100
QPOL ART U
QP NR 1 Y 0.0000 Z -0.0000
QP NR 2 Y 4.3000 Z -0.0000
QP NR 3 Y 4.3000 Z -15.7000
QP NR 4 Y 4.9699 Z -18.2000
QP NR 5 Y 6.8000 Z -20.0301
QP NR 6 Y 9.3000 Z -20.7000
QP NR 7 Y 15.0000 Z -20.7000
QP NR 8 Y 20.7000 Z -20.7000
QP NR 9 Y 23.2000 Z -20.0301
QP NR 10 Y 25.0301 Z -18.2000
QP NR 11 Y 25.7000 Z -15.6135
QP NR 12 Y 25.7000 Z -0.0000
QP NR 13 Y 30.0000 Z -0.0000
QP NR 14 Y 30.0000 Z -15.6135
QP NR 15 Y 28.7540 Z -20.3500
QP NR 16 Y 25.3500 Z -23.7540
QP NR 17 Y 20.7000 Z -25.0000
QP NR 18 Y 15.0000 Z -25.0000
QP NR 19 Y 9.3000 Z -25.0000
QP NR 20 Y 4.6500 Z -23.7540
QP NR 21 Y 1.2460 Z -20.3500
QP NR 22 Y 0.0000 Z -15.7359
ENDE

PROG AQB urs:2
KOPF PLASTISCHE SG VERTEILUNG
ECHO VOLL EXTR
give acting forces Plastic Moment from ReportBrowser
S NQ 1 NR 01 X 0 N 0.00 MY 1.165
S NQ 1 NR 01 X 1 N 0.00 MY -1.165

DEHN KMOD Gebr KSV SL

$ now the problem occurs - see result viewer - calculated stresses are not in equilibrium
ENDE

Above code was edited 17.06.2020 as a copy&paste mistake occured.

However i think my initial guess was right, as a change of NORM to
"norm din 18800 "
solved the “problem” in this case

best regards