#1 2022-11-01 15:54:42

bay bach

Member

Registered: 2022-11-01

Posts: 9

about RC

hi guys ! Is here anyone experienced by Reinforced Concrete.

mecour

Member

From: Ostrava (Czech)

Registered: 2011-04-04

Posts: 199

Re: about RC

Hello,

you need to be more specific. I use CA possibilities for RC concrete and also prestressed concrete.

mecour

bay bach

Member

Registered: 2022-11-01

Posts: 9

Re: about RC

Thank u for replying. Can CA modelize steel and concrete in 3D? I need that example of code.

mecour

Member

From: Ostrava (Czech)

Registered: 2011-04-04

Posts: 199

Re: about RC

Hello,

of course you can simulate in 3D.  You can use many nonlinear material models and also interface model JOINT_BA for steel-concrete connection.

https: // code-aster.org/V2/UPLOAD/DOC/Formations/02-concretece.pdf

mecour

Member

From: Ostrava (Czech)

Registered: 2011-04-04

Posts: 199

Re: about RC

Hello,

try to start with simplier model. In 3D you need to use 3D_INTERFACE not PLAN_INTERFACE.
JOINT_BA is possible to use only in 2D models, in 3D use CZM_LAB_MIX (maybe I caused this misunderstand, I am sorry).
In your model are collision of groups, group t13 is also in group t16.
Wthat is the material with elastic modulus over 600GPa?

bay bach

Member

Registered: 2022-11-01

Posts: 9

Re: about RC

Hello again
Trying to use 3D_INTERFACE and CZM_LAB_MIX . But something wrong in my code. Result graphic is not good.
Can u help me please.

Bay Bach

DEBUT(identifier='0:1',
      )

MAIL = LIRE_MAILLAGE(identifier='1:1',
                     FORMAT='MED',
                     UNITE=2)

MODE = AFFE_MODELE(identifier='2:1',
                   AFFE=(_F(GROUP_MA=('bond', ),
                            MODELISATION=('3D_INTERFACE', ),
                            PHENOMENE='MECANIQUE'),
                         _F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI')),
                   MAILLAGE=MAIL)

ST = DEFI_MATERIAU(identifier='3:1',
                   ECRO_LINE=_F(D_SIGM_EPSI=30000.0,
                                SY=3000.0),
                   ELAS=_F(E=210000.0,
                           NU=0.3))

mazars = DEFI_MATERIAU(identifier='4:1',
                       ELAS=_F(E=26823.0,
                               NU=0.2),
                       MAZARS=_F(AC=0.79095138,
                                 AT=0.9,
                                 BC=2339.0,
                                 BT=18247.0,
                                 EPSD0=5.48027e-05,
                                 EPSI_LIM=0.0035,
                                 K=0.7,
                                 SIGM_LIM=8.7))

bond = DEFI_MATERIAU(identifier='5:1',
                     CZM_LAB_MIX=_F(ALPHA=0.5,
                                    BETA=1.0,
                                    CINEMATIQUE='UNILATER',
                                    GLIS_C=0.08,
                                    PENA_LAGR=100.0,
                                    SIGM_C=10.0),
                     ELAS=_F(E=26823.0,
                             NU=0.2))

MATE = AFFE_MATERIAU(identifier='6:1',
                     AFFE=(_F(GROUP_MA=('conc', ),
                              MATER=(mazars, )),
                           _F(GROUP_MA=('rebar', ),
                              MATER=(ST, )),
                           _F(GROUP_MA=('bond', ),
                              MATER=(bond, ))),
                     MAILLAGE=MAIL)

ramp = DEFI_FONCTION(identifier='7:1',
                     NOM_PARA='INST',
                     PROL_DROITE='LINEAIRE',
                     PROL_GAUCHE='LINEAIRE',
                     VALE=(0.0, 0.0, 1.0, 0.1, 2.0, 0.2, 3.0, 0.3, 4.0, 0.4, 5.0, 0.5, 6.0, 0.6, 7.0, 0.7, 8.0, 0.8, 9.0, 0.9, 10.0, 1.0))

listr = DEFI_LIST_REEL(identifier='8:1',
                       DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=1.0,
                                     PAS=0.1))

times = DEFI_LIST_INST(identifier='9:1',
                       DEFI_LIST=_F(LIST_INST=listr),
                       ECHEC=_F(ACTION='DECOUPE',
                                EVENEMENT='ERREUR',
                                SUBD_METHODE='MANUEL',
                                SUBD_NIVEAU=3,
                                SUBD_PAS=5,
                                SUBD_PAS_MINI=1e-10),
                       METHODE='MANUEL')

fix = AFFE_CHAR_MECA(identifier='10:1',
                     DDL_IMPO=_F(DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,
                                 GROUP_MA=('fix', )),
                     MODELE=MODE)

load = AFFE_CHAR_MECA(identifier='11:1',
                      DDL_IMPO=_F(DZ=0.01,
                                  GROUP_MA=('load', )),
                      MODELE=MODE)

resnonl = STAT_NON_LINE(identifier='12:1',
                        CHAM_MATER=MATE,
                        COMPORTEMENT=(_F(GROUP_MA=('bond', ),
                                         RELATION='CZM_LAB_MIX'),
                                      _F(GROUP_MA=('conc', ),
                                         RELATION='MAZARS'),
                                      _F(GROUP_MA=('rebar', ),
                                         RELATION='VMIS_ISOT_LINE'),
                                      _F(RELATION='ELAS')),
                        CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                       RESI_GLOB_MAXI=0.001,
                                       RESI_GLOB_RELA=1e-05),
                        EXCIT=(_F(CHARGE=fix),
                               _F(CHARGE=load,
                                  FONC_MULT=ramp)),
                        INCREMENT=_F(LIST_INST=times),
                        METHODE='NEWTON',
                        MODELE=MODE,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  REAC_INCR=1,
                                  REAC_ITER=1))

resnonl = CALC_CHAMP(identifier='13:1',
                     reuse=resnonl,
                     CHAM_MATER=MATE,
                     CONTRAINTE=('SIEF_ELNO', 'SIEF_NOEU'),
                     DEFORMATION=('EPSI_ELNO', 'EPSI_NOEU'),
                     FORCE=('REAC_NODA', 'FORC_NODA'),
                     MODELE=MODE,
                     RESULTAT=resnonl)

IMPR_RESU(identifier='14:1',
          FORMAT='MED',
          RESU=_F(NOM_CHAM=('DEPL', 'SIEF_NOEU', 'EPSI_NOEU', 'REAC_NODA', 'FORC_NODA'),
                  RESULTAT=resnonl),
          UNITE=80)

FIN()

bay bach

Member

Registered: 2022-11-01

Posts: 9

Re: about RC

Med file is here.
h ttps://drive.google.com/file/d/1Xt35_nFcxRK1eJC0SsPYfjtQBeQgCrkT
I need a bond-slip graphic between steel and concrete. The analysis runs. But something wrong with code. Thank you

mecour

Member

From: Ostrava (Czech)

Registered: 2011-04-04

Posts: 199

Re: about RC

Hello,

please study the documentation a little (r3.06.13, u3.13.14, ). You can't use cohesive zone model with volume element. Interface elements must be zero thickness.  So in your study the interface can be only face between rebar and concrete created by degenerated zero thickness elements.
You can create zero thickness elements this way:
-- create group of faces on between rebar and concrete
-- than use in mesh module Modification->Transformation->Duplicate Nodes or Elements, this operation create zero thickness volume element for CZM model.

mecour

bay bach

Member

Registered: 2022-11-01

Posts: 9

Re: about RC

Thank you very much.
Can you share some simple examples please.