Examples
========

Examples of uses for crystalline materials
-----------------------------------------------------

The following example is a typical use of MONOCRISTAL. It comes from test SSNV171B.

ACIER = DEFI_MATERIAU (

ELAS =_F (E=145200.0, NU=0.3,),

MONO_VISC2 =_F (N=10.0, K=40.0, C=40.0, C=1.0, D=36.68, A=10.0,),

MONO_ISOT2 =_F (R_0=75.5, Q1=9.77, B1=9.77, B1=19.34, H=0.5, Q2=-33.27, B2=5.345,),

MONO_CINE1 =_F (D=36.68,),

)

COMPORT = DEFI_COMPOR (

MONOCRISTAL =_F (MATER = ACIER, ELAS =' ELAS ',

ECOULEMENT =' MONO_VISC2 ',

ECRO_ISOT =' MONO_ISOT2 ',

ECRO_CINE =' MONO_CINE1 ',

FAMI_SYST_GLIS =' OCTAEDRIQUE ',),

)

The following example, implementing POLYCRISTAL, is from the SSNV171B test:

MATPOLY = DEFI_MATERIAU (

ELAS =_F (E=192500.0, NU=0.3,),

MONO_VISC2 =_F (N=10.0, K=40.0, C=40.0, C=6333.0, D=36.68, A=72.21,),

MONO_ISOT2 =_F (R_0=75.5, Q1=9.77, B1=9.77, B1=19.34, H=2.54, Q2=-33.27, B2=5.345,),

MONO_CINE1 =_F (D=36.68,),

)

MONO1 = DEFI_COMPOR (

MONOCRISTAL =_F (MATER = MATPOLY, ELAS =' ELAS ',

ECOULEMENT =' MONO_VISC2 ',

ECRO_ISOT =' MONO_ISOT2 ',

ECRO_CINE =' MONO_CINE1 ',

FAMI_SYST_GLIS =' OCTAEDRIQUE ',),

)

POLY1 = DEFI_COMPOR (

POLYCRISTAL =(

_F (MONOCRISTAL = MONO1, FRAC_VOL =0.025,

ANGL_REP =( -149.67,15.61,154.67,),),

_F (MONOCRISTAL = MONO1, FRAC_VOL =0.025,

ANGL_REP =( -481.7,35.46,188.7),),

),

MU_LOCA = 82500

LOCALISATION =' BETA ', DL=321.5, DA=0.216,

)


Example of use for multifibers
------------------------------------------

The commands below illustrate the use of DEFI_COMPOR for multi-fiber behavior (see for example the SSNL119B test).

GF= DEFI_GEOM_FIBRE (

FIBRE =_F (GROUP_FIBRE =' SACI ', CARA =' DIAMETRE ',

COOR_AXE_POUTRE = (0.,0.,),

VALE =( 0.066, -0.218, 32.E-3,

0.066, -0.218, 32.E-3,

0.066, 0.218, 8.E-3,

0.066, 0.218,8.E-3,),),

),

SECTION =_F (GROUP_FIBRE =' SBET ', MAILLAGE_SECT = MASEC,

TOUT_SECT =' OUI ',

COOR_AXE_POUTRE = (0., 0.,),

),

)

MOPOU = AFFE_MODELE (

MAILLAGE = MAPOU,

AFFE =_F (TOUT =' OUI ', PHENOMENE =' MECANIQUE', MODELISATION =' POU_D_EM ',),

)

# Concrete

EB = 37272.0E+06

BETON = DEFI_MATER_GC (

MAZARS =_F (UNITE_LONGUEUR = 'M', FCJ =40.963E+06, EIJ =EB,

EPSI_C =1.75754E-03, AT=1.0, NU=0.2,),

RHO =2400.0, INFO =2,

)

ACIER = DEFI_MATER_GC (

ACIER =_F (E = 2.0E+11, D_ SIGM_EPSI =3.28E+9, SY=4.E+8,),

RHO =7800.0,

)

MATOR = DEFI_MATERIAU (ELAS =_F (E=2.E11, NU=0.0, RHO =7800.0,))

POUCA = AFFE_CARA_ELEM (

MODELE = MOPOU,

POUTRE =_F (GROUP_MA =' POUTRE ', SECTION =' RECTANGLE',

CARA =( 'HY', 'HZ'), VALE =( 0.2,0.5),

PREC_AIRE =5.0, PREC_INERTIE =10.0,),

ORIENTATION =_F (GROUP_MA =' POUTRE ', CARA =' ANGL_VRIL', VALE =-90.0,),

GEOM_FIBRE =GF,

MULTIFIBRE =_F (GROUP_MA =' POUTRE ', GROUP_FIBRE =(' SBET ',' SACI ')),

)

COMPPMF = DEFI_COMPOR (GEOM_FIBRE =GF, MATER_SECT = MATOR,

MULTIFIBRE =(

_F (GROUP_FIBRE =' SACI ', MATER = ACIER, RELATION =' VMIS_CINE_LINE'),

_F (GROUP_FIBRE =' SBET ', MATER = BETON, RELATION =' MAZARS_UNIL'),

),

)