*** CALCUL NON LINEAIRE PAS A PAS POUTRE ENCASTREE ***
***** DEFINITION DES OPTIONS *****
OPTI ECHO 0 ;
OPTI DIME 2 ELEM TRI3 MODE PLAN CONT ;
N1 N2 = 0.005 0.01 ; COMM DENSITE LOCALE ;
TIMI = 1 ;
*****DEFINITION DE LA GEOMETRIE********
Long = 1.5 ; COMM LONGUEUR POUTRE ;
he = 0.25 ; COMM HAUTEUR ENCASTREMENT ;
hg = 0.15 ; COMM HAUTEUR GOUSSET ;
*****DEFINITION DES POINTS
DENS N1 ;
P1 = 0. 0. ;
P4 = 0. he ;
DENS N2 ;
P2 = long (he – hg) ;
P3 = long he ;
*****DEFINITION DES LIGNES ET CONTOURS
D1 = P1 DROI P2 DROI P3 DROI P4 ;
D2 = P4 DROI P1 ;
CT = D1 ET D2 ;
******DEFINITION DE LA SECTION
ST = SURF CT PLAN ;
TRAC ST ;
*** DECLARATION DES PARAMETRES MATERIAU ***
*** CARACTERISTIQUES DU BETON ***
FT = 3.E6 ;
GF = 100. ;
YUN = 35.E9 ;
SG1 = YUN/10. ;
LCAR1 = 0.02 ; COMM LONGUEUR
CARACTERISTIQUE DES ELEMENTS ;
MODT = MODE ST MECANIQUE ELASTIQUE
ENDOMMAGEMENT MAZARS ;
*** CHAMP ALEATOIRE FONCTION DU
*MODULE D'YOUNG ***
CH1 = ALEA MODT 'BANDES_TOURNANTES' EXPO
'SIGMA' SG1 'MOYENNE' YUN 'LAMBDA' LCAR1 ;
*** INITIALISATION DE LA MATRICE MATERIAU
*POUR LES PARAMETRES BCOM ET BTRA ***
MAT0 = MATE MODT BCOM 1 ;
MAT0 = MATE MODT BTRA 1 ;
*** INTEGRATION SUR LES ELEMENTS FINIS ***
INTG1 = INTG MAT0 MODT ELEM ;
INTG1 = INTG MAT0 MODT ELEM ;
*** CALCUL DES PARAMETRE BETON COMPRESSION BCOM
*BETON TRACTION BTRA ET SEUIL INITIAL EN DEFORMATION
*POUR LA TRACTION KTR ***
BCO1 = INTG1**0.5 ;
BCO2 = BCO1*FT/GF ;
BTR1 = INTG1**0.5 ;
BTR2 = BTR1*FT/GF ;
KTR = FT/YUN ;
MATT = MATE MODT YOUN CH1 NU 0.25 BTRA BTR2 ATRA
0.08 ACOM 1.4 KTR0 1E-4 BCOM 190 BETA 1.06 ;
**** CONDITIONS AUX LIMITES ****
CL1 = BLOQ DEPL ROTA D2 ;
CL2 = BLOQ P3 UY ;
CLT = CL1 ET CL2 ;
**** DEPLACEMENT IMPOSE ****
AMPL = -1*(long/600) ;
DEP1 = DEPI CL2 AMPL ;
VEC1 = ‘VECT’ (DEP1 ‘NOMC’ ‘UY’) DEPL ‘BLEU’ ;
‘TRAC’ VEC1 ST ‘TITR’ ‘DEPLACEMENT IMPOSE’ ;
‘FINS’ ;
*** MATRICE DE RIGIDITE ***
RIG1 = RIGI MODT MATT ;
RIGT = RIG1 ET CLT ;
**** RESOLUTION LINEAIRE ET RESULTATS ****
RES1 = RESOUDRE RIGT DEP1 ;
DXY = EXCO RES1 UX UY ;
DXYABS = ABS DXY ;
DXYMAX = MAXI DXYABS ;
MESS 'DEPLACEMENT EN METRES DEP1 :' DXYMAX ;
*****TRACE DE LA DEFORMEE ********
DEF0 = DEFORM ST RES1 0. BLEU ;
DEF1 = DEFORM ST RES1 TIMI ROUG ;
TRAC (DEF0 ET DEF1) TITR ‘DEFORMEE’ ;
******CHARGEMENT PAR DEPLACEMENT IMPOSE*******
AX1 = PROG 0. TIMI ;
AX2 = PROG 0. TIMI ;
EV1 = EVOL MANU 'temps' AX1 'coefficient multiplicateur' AX2
COUL ORANG ;
@EXCEL1 EV1 TPS ;
*Traçage de EV1
DESS EV1 TITRE 'k=f(t)';
CHA = CHAR DIMP DEP1 EV1 ;
************** TEMPS POUR LES DEPLACEMENTS***************
T0= 0.01 ;
LIS3 = PROG 0. PAS t0 TIMI ;
**** CALCUL PASAPAS ****
TAB1 = TABLE ;
TAB1. ’MOVA' = ‘D’ ;
TAB1. 'CHARGEMENT' = CHA ;
TAB1. ‘MODELE’ = MODT ;
TAB1. 'CARACTERISTIQUES' = MATT ;
TAB1. 'BLOCAGES_MECANIQUES' = CLT ;
TAB1. 'TEMPS_CALCULES' = LIS3 ;
TAB1. 'TEMPS_SAUVES' = LIS3 ;
TAB1. 'ENDOMMAGEMENT' = VRAI ;
*TAB1.'NON_LOCAL'= VRAI ;
TAB1. 'GRANDS_DEPLACEMENTS' = VRAI ;
TAB1. 'MAXITERATION' = 5 ;
*LC = 0.005 ;
TAB1. ’LCAR' = LCAR1 ;
*TAB1. ’CONN' = ‘CONN’ MODT LC
*‘NORMAL’ ‘INTERIEUR’ ;
TAB1. 'LVIAM’ = MOTS EPTI ;
PASAPAS TAB1 ;