hello, yes I already figured that out (with help of chatGPT)
But why there is no build in function to compute sum?

conta2 = POST_RELEVE_T(ACTION=_F(GROUP_NO=('contact', ),

                                 INTITULE='contact',

                                 NOM_CHAM='CONT_NOEU',

                                 OPERATION=('EXTRACTION', ),

                                 RESULTAT=STA_NL,

                                 TOUT_CMP='OUI'))

Hello, changing the beginning to

#LOADING MODULES
#import sys
#sys.path.append(r'C:\Users\AppData\Local\code_aster\v2021\14.8\lib\aster\Utilitai' )

#from code_aster.MacroCommands.Utils import partition
from code_aster.MacroCommands.Utils.partition import *

#from Utilitai.partition import *
#from partition import *
#from Table import*
import numpy as N
# import Numeric as N
import math

worked for me with 15.4

however only to 5th iteration, then it takes really long to solve.

there is no part 3.9

However if you mean this:

#comment: conta = POST_RELEVE_T(ACTION=_F(GROUP_NO=('contact_tube', ),
#comment:                                 INTITULE='contact',
#comment:                                 NOM_CHAM='CONT_NOEU',
#comment:                                 OPERATION=('EXTRACTION', ),
#comment:                                 RESULTANTE=('RNY', ),
#comment:                                 RESULTAT=STA_NL))

Doesn't work

finer mesh

this is the comm which also calculates total load automatically and works even in graphical mode in salomemeca (The variable is highlighted in red, because of the duplicate name, but it is not a problem )

it is also a little bit different in prepomax

bafybeid6d6irdkqjivza5lnhjiingxegopwc7ttirultkoispvnb3lovbm.ipfs.dweb.link/RQ_1_.zip

sigma_x

The difference with ansys 145 000 MPa may be caused by different averaging or rounding, but from physical point of view it is just a singularity, so no need to worry about it too much.

I guess I have found quite a nice solution how to do that using this post

Easy question: how to set a force on a surface?

The advantage is that Stage_1 and Stage_3 can be in graphical mode.

youtube.com/watch?v=vOfexbjNutE&list=PLBwQ5Llf4Q_Vj8xKcxsbzIfdqrtzcHLbe&index=32

I understand that it is how the numerics work, just wonder why there is no FORCE TOTAL

it could be done just like doing something like this before affe_cara_meca

VOL= POST_ELEM(INST=(0.0, ),
                   INTEGRALE=_F(GROUP_MA=('Fx', ),
                                NOM_CHAM='DEPL',
                                NOM_CMP=('DX', ),
                                TYPE_MAILLE='2D'),
                   MODELE=model,
                   RESULTAT=resnonl)

IMPR_TABLE(NOM_PARA=('VOL', ),
           TABLE=table0,
           UNITE=3)

and then FORCE_FACE / VOL

I can do
VOL= POST_ELEM(INST=(0.0, ),
                   INTEGRALE=_F(GROUP_MA=('Fx', ),
                                NOM_CHAM='DEPL',
                                NOM_CMP=('DX', ),
                                TYPE_MAILLE='2D'),
                   MODELE=model,
                   RESULTAT=resnonl)

and then copy paste, but it is the same as using salome or whatever
or create dummy model and use python but make thing quite complicated for very basic thing.