!#!KAPITEL MATERIAUX ET GEOMETRIES
+PROG AQUA urs:1
HEAD Définition des caractéristiques des matériaux et de la section des éléments de fondations
ECHO FULL VAL FULL
CTRL WARN 12343 ; CTRL WARN 12303
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DU SYSTÈME D'UNITÉ !!!!
$=========================================================================================================================================================$
UNIT 5
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES NORMES DE CALCUL !!!!
$=========================================================================================================================================================$
NORM 'NF' 'en199x-200x' COUN 33 UNIT 5
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES VARIABLES !!!!
$=========================================================================================================================================================$
// Définition des variables pour béton //
// Valeur caractéristique de la résistance à la compression fck
STO#fck 18.3299999237061
// Valeur caractéristique de la résistance à la traction fctm
STO#fctm_etoile 2.0855929851532
// Valeur caractéristique de la résistance fctk,0.05
STO#fctk_005_etoile 1.45991504192352
// Valeur du module d'élasticité Ecm
STO#Ecm_etoile 29414.26171875
// Définition des variables pour les aciers //
STO#fe 263.32
STO#classe 'Classe A'
// Définition des variables pour la section transversale des éléments de fondation //
// Diamètre des éléments de fondation s'ils ont une section circulaire (pieux ou micropieux)
STO#Diametre 0.3
// Largeur des éléments de fondation s'ils ont une section rectangulaire (barrettes)
STO#TB_Lbarrettes 0.5
// Hauteur des éléments de fondation s'ils ont une section rectangulaire (barrettes)
STO#TB_Hbarrettes 0.5
// Section A (m²) dans le cas où les éléments de fondation sont de type micropieux
STO#Amicropieux 50
// Inertie I (m²) dans le cas où les éléments de fondation sont de type micropieux
STO#Imicropieux 50
// Variable de type BOOLEAN (true ou false) qui reflète si les éléments de fondation sont de type barrettes ou non
STO#OB_Barrettes 'Vrai'
// Variable de type BOOLEAN (true ou false) qui reflète si les éléments de fondation sont de type pieux ou non
STO#OB_Pieux 'Faux'
// Définition des variables pour les lois de comportement décrivant l'interaction sol-structure. //
// Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction x à court terme
STO#valueskfcourtX 43643.3586194397, 87286.7172388794, 87286.7172388794, 87286.7172388794, 87286.7172388794, 87286.7172388794, 60849.8716821696, 263306.82765317, 559318.257858671, 886554.251045694 $$
443277.125522847
$Fin variable
// Palier de pression r1 dans la direction x
STO#valuespalierr1X 381.25, 762.5, 762.5, 762.5, 762.5, 762.5, 576.25, 2220, 3606.25, 4712.5 $$
2356.25
$Fin variable
//Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction x à long terme
STO#valueskflongX 21821.6793097199, 43643.3586194397, 43643.3586194397, 43643.3586194397, 43643.3586194397, 43643.3586194397, 30424.9358410848, 131653.413826585, 279659.128929335, 443277.125522847 $$
221638.562761423
$Fin variable
// Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction x sous combinaisons sismique
STO#valueskfsismX 130930.075858319, 261860.151716638, 261860.151716638, 261860.151716638, 261860.151716638, 261860.151716638, 182549.615046509, 789920.48295951, 1677954.77357601, 2659662.75313708 $$
1329831.37656854
$Fin variable
// Palier de pression r2 dans la direction x
STO#valuespalierr2X 787.5, 1575, 1575, 1575, 1575, 1575, 913.75, 3641.25, 4062.5, 5000 $$
2500
$Fin variable
// !!!!!!!! Direction YY !!!!!!!! //
// Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction y à court terme
STO#valueskfcourtY 43643.3586194397, 87286.7172388794, 87286.7172388794, 87286.7172388794, 87286.7172388794, 87286.7172388794, 60849.8716821696, 263306.82765317, 559318.257858671, 886554.251045694 $$
443277.125522847
$Fin variable
// Palier de pression r1 dans la direction y
STO#valuespalierr1Y 381.25, 762.5, 762.5, 762.5, 762.5, 762.5, 576.25, 2220, 3606.25, 4712.5 $$
2356.25
$Fin variable
// Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction y à long terme
STO#valueskflongY 21821.6793097199, 43643.3586194397, 43643.3586194397, 43643.3586194397, 43643.3586194397, 43643.3586194397, 30424.9358410848, 131653.413826585, 279659.128929335, 443277.125522847 $$
221638.562761423
$Fin variable
// Module de réaction de mobilisation de la réaction frontale et tangentielle dans la direction y sous combinaisons sismique
STO#valueskfsismY 130930.075858319, 261860.151716638, 261860.151716638, 261860.151716638, 261860.151716638, 261860.151716638, 182549.615046509, 789920.48295951, 1677954.77357601, 2659662.75313708 $$
1329831.37656854
$Fin variable
// Palier de pression r2 dans la direction y
STO#valuespalierr2Y 787.5, 1575, 1575, 1575, 1575, 1575, 913.75, 3641.25, 4062.5, 5000 $$
2500
$Fin variable
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES MATÉRIAUX !!!!
$=========================================================================================================================================================$
// Définition du matériau béton pour les éléments de fondation //
CONC 1 C '20' FC #fck FCT #fctm_etoile FCTK #fctk_005_etoile EC #Ecm_etoile MUE 0.200000 GAM 25 ALFA 1.00000E-05 FCR 28 ECR 1.05*#Ecm_etoile $$
FCTD #fctk_005_etoile/1.5 GMOD #Ecm_etoile/2.4 KMOD #Ecm_etoile/1.80 RHO 2400 TITL "=Béton élément de fondation"
SSLA EPS SERV SIG 1.500000[-] TYPE LIM
STO#EPS_SERV 0.00, -0.19, -0.39, -0.58, -0.77, -0.96, -1.16, -1.35, -1.54, -1.74, -1.93, -2.10, -2.28, -2.45, -2.63, -2.80, -2.98, -3.15, -3.33, -3.50
$Fin variable
STO#SIG_SERV 0.00, -5.55, -10.32, -14.37, -17.75, -20.51, -22.69, -24.33, -25.46, -26.12, -26.33, -26.16, -25.67, -24.88, -23.79, -22.43, -20.82, -18.95, -16.85, -14.52
$Fin variable
loop#i SIG_SERV
SSLA EPS #EPS_SERV(#i) SIG #SIG_SERV(#i) TYPE 'POL'
endloop
SSLA EPS ULTI SIG 1.500000[-] TYPE LIM
STO#EPS_ULTI 0.00, -0.19, -0.39, -0.58, -0.77, -0.96, -1.16, -1.35, -1.54, -1.74, -1.93, -2.10, -2.28, -2.45, -2.63, -2.80, -2.98, -3.15, -3.33, -3.50
$Fin variable
STO#SIG_ULTI 0.00, -2.58, -4.79, -6.67, -8.24, -9.52, -10.53, -11.29, -11.82, -12.12, -12.22, -12.14, -11.91, -11.55, -11.04, -10.41, -9.66, -8.79, -7.82, -6.74
$Fin variable
loop#i SIG_ULTI
SSLA EPS #EPS_ULTI(#i) SIG #SIG_ULTI(#i) TYPE 'POL'
endloop
SSLA EPS CALC SIG -1.500000[-] TYPE LIM
STO#EPS_CALC 0.00, -0.20, -0.40, -0.60, -0.80, -1.00, -1.20, -1.40, -1.60, -1.80, -2.00, -3.50
$Fin variable
STO#SIG_CALC 0.00, -3.48, -6.60, -9.35, -11.73, -13.75, -15.40, -16.68, -17.60, -18.15, -18.33, -18.33
$Fin variable
loop#i SIG_CALC
SSLA EPS #EPS_CALC(#i) SIG #SIG_CALC(#i) TYPE 'POL'
endloop
HMAT 1 TYPE FOUR NSP 0.030000 A 1
// Définition du matériau pour la semelle de liaison //
CONC 2 C '100' FC 100 EC 1e10 GAM 0 TYPR B FCR 33 FCTD 1.200000 TITL "=C 100/115 (EN 1992)_Semelle"
HMAT 2 TYPE FOUR NSP 0.030000 A 1
// Définition du matériau béton pour les éléments de fondation //
STEE 3 B #classe FY #fe FT #fe+40 FP #fe ES 200000 MUE 0.300000 GAM 78.50000 ALFA 1.20000E-05 EPST 50 FYC #fe FTC #fe+40 TMAX 32 GMOD 76923 KMOD 166667 TITL "=B 500 B (EN 1992) OMRI"
SSLA EPS SERV SIG 1.150000[-] TYPE EXT
SSLA EPS ULTI SIG -1.150000[-] TYPE EXT
SSLA EPS CALC SIG -1.150000[-] TYPE EXT
HMAT 3 TYPE FOUR NSP 0 A 1
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DE LA SECTION TRANSVERSALE DES ÉLÉMENTS DE FONDATION !!!!
$=========================================================================================================================================================$
if #OB_Barrettes == 'Vrai'
SREC 1 H #TB_Lbarrettes*1000 B #TB_Hbarrettes*1000 SO 90 SU 90 SS 90 MNO 1 MRF 3 RTYP 'ASYM' ASO 12.56000 ASU 12.56000 DASO 20 DASU 20 $$
DASS 20 A 150 INCL 90 REF C IT 100[o/o] AY 100[o/o] AZ 100[o/o] BCYZ '0' SPT 0 TITL "B/H = #TB_Lbarrettes*1000 / #TB_Hbarrettes*1000 mm"
elseIF #OB_Pieux == 'Vrai'
SCIT 1 D #Diametre*1000 SA 90 MNO 1 MRF 3 RTYP BEND ASA 18.84000[cm2] DAS 32 A 20[-] ASL 0 IT 100[o/o] AY 100[o/o] AZ 100[o/o] TITL "D #Diametre mm"
else
$SCIT 1 D #Diametre*1000 SA 90 MNO 1 MRF 3 RTYP BEND ASA 18.84000[cm2] DAS 32 A 20[-] ASL 0 IT 100[o/o] AY 100[o/o] AZ 100[o/o] TITL "D #Diametre mm"
SVAL NO 1 MNO 3 A #Amicropieux Iy #Imicropieux Iz #Imicropieux
endif
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES LOI DE COMPORTEMENT !!!!
$=========================================================================================================================================================$
// !!!!!!!! Direction XX !!!!!!!! //
//Définition des loi de comportement à court terme _ direction XX
LOOP#i valueskfcourtX
let#k #i+1
if #valuespalierr1X(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure à court terme n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskfcourtX(#i)/#valuespalierr1X(#i)))[mm] F #valuespalierr1X(#i)[kN] S POL
SFLA U (1000/(#valueskfcourtX(#i)/#valuespalierr1X(#i))+25)[mm] F #valuespalierr1X(#i)[kN] S POL
endif
ENDLOOP
//Définition des loi de comportement à long terme _ direction XX
LOOP#i valueskfcourtX
let#k #i+101
if #valuespalierr1X(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure à long terme n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskflongX(#i)/#valuespalierr1X(#i)))[mm] F #valuespalierr1X(#i)[kN] S POL
SFLA U (1000/(#valueskflongX(#i)/#valuespalierr1X(#i))+25)[mm] F #valuespalierr1X(#i)[kN] S POL
endif
ENDLOOP
//Définition des loi de comportement sous combinaisons sismique _ direction XX
LOOP#i valueskfcourtX
let#k #i+301
if #valuespalierr1X(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure sous situation sismique n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskfsismX(#i)/#valuespalierr1X(#i)))[mm] F #valuespalierr1X(#i)[kN] S POL
SFLA U (1000*(2*(#valuespalierr2X(#i)-#valuespalierr1X(#i))/#valueskfsismX(#i)+#valuespalierr1X(#i)/#valueskfsismX(#i)))[mm] F #valuespalierr2X(#i)[kN] S POL
SFLA U (1000*(2*(#valuespalierr2X(#i)-#valuespalierr1X(#i))/#valueskfsismX(#i)+#valuespalierr1X(#i)/#valueskfsismX(#i))+25)[mm] F #valuespalierr2X(#i)[kN] S POL
endif
ENDLOOP
// !!!!!!!! Direction YY !!!!!!!! //
//Définition des loi de comportement à court terme _ direction YY
LOOP#i valueskfcourtY
let#k #i+401
if #valuespalierr1Y(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure à court terme n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskfcourtY(#i)/#valuespalierr1Y(#i)))[mm] F #valuespalierr1Y(#i)[kN] S POL
SFLA U (1000/(#valueskfcourtY(#i)/#valuespalierr1Y(#i))+25)[mm] F #valuespalierr1Y(#i)[kN] S POL
endif
ENDLOOP
//Définition des loi de comportement à long terme _ direction YY
LOOP#i valueskfcourtY
let#k #i+501
if #valuespalierr1Y(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure à long terme n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskflongY(#i)/#valuespalierr1Y(#i)))[mm] F #valuespalierr1Y(#i)[kN] S POL
SFLA U (1000/(#valueskflongY(#i)/#valuespalierr1Y(#i))+25)[mm] F #valuespalierr1Y(#i)[kN] S POL
endif
ENDLOOP
//Définition des loi de comportement sous combinaisons sismique _ direction YY
LOOP#i valueskfcourtY
let#k #i+601
if #valuespalierr1X(#i) <> 0
SMAT #k LTYP STD MTYP PLAS ALPH 0 TITL "Loi d'interaction sol structure sous situation sismique n°#k "
SFLA U 0[mm] F 0[kN] S POL TYPE PT
SFLA U (1000/(#valueskfsismY(#i)/#valuespalierr1Y(#i)))[mm] F #valuespalierr1X(#i)[kN] S POL
SFLA U (1000*(2*(#valuespalierr2Y(#i)-#valuespalierr1Y(#i))/#valueskfsismY(#i)+#valuespalierr1Y(#i)/#valueskfsismY(#i)))[mm] F #valuespalierr2Y(#i)[kN] S POL
SFLA U (1000*(2*(#valuespalierr2Y(#i)-#valuespalierr1Y(#i))/#valueskfsismY(#i)+#valuespalierr1Y(#i)/#valueskfsismY(#i))+25)[mm] F #valuespalierr2Y(#i)[kN] S POL
endif
ENDLOOP
END
+PROG SOFIMSHA urs:2
HEAD Définition de la géométrie des éléments de fondation
ECHO FULL VAL FULL
CTRL WARN 10752
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DU SYSTÈME D'UNITÉ !!!!
$=========================================================================================================================================================$
UNIT 5
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES PARAMETRES DU SYSTEME !!!!
$=========================================================================================================================================================$
SYST SPAC GDIV 10000 GDIR NEGZ
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES VARIABLES !!!!
$=========================================================================================================================================================$
// Caractéristiques géométriques du systéme de fondation //
// Niveau supérieur des éléments de fondation
STO#Zinf_semelle 100
// Longueur des éléments de fondation
STO#Lfondation 25
// Nombre des éléments de fondation
STO#Nbr_fondation 6
// Nombre des points de déscritisation
STO#Nbr_descri 10
// Longueur des segments de barre suite à la discrétisation
STO#Linf_noeuds #Lfondation/#Nbr_descri
// Cordonnée x des différents éléments de fondation
STO#x -1.5, -1.5, -1.5, 1.5, 1.5, 1.5
$Fin variable
// Cordonnée y des différents éléments de fondation
STO#y -3, 0, 3, -3, 0, 3
$Fin variable
// Inclinaison des différents éléments de fondation autour de l'axe x
STO#Ix 0, 0, 0, 0, 0, 0
$Fin variable
// Inclinaison des différents éléments de fondation autour de l'axe y
STO#Iy 0, 0, 0, 0, 0, 0
$Fin variable
// Epaisseur de la semelle
STO#E_semelle 2
STO#semellex -2.5, 2.5, 2.5, -2.5
$Fin variable
STO#semelley -2.5, -2.5, 2.5, 2.5
$Fin variable
// Rigidité verticale (axiale) (terrain seul: frottement+pointe) modélisée par un ressort à la base des pieux en kN/m
STO#Raxiale 4542518.5
//Rigidité en rotation autour de l'axe X sous la pointe d'un pieu en kN.m/rd, à calculer pour des barrettes
STO#Rrotx 8594.38337430909
//Rigidité en rotation autour de l'axe Y sous la pointe d'un pieu en kN.m/rd, à calculer pour des barrettes
STO#Rroty 8594.38337430909
// Caractéristiques géométriques de la semelle de liaison //
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES NOEUDS !!!!
$=========================================================================================================================================================$
// Noeuds pour éléments de fondation //
Loop#i #Nbr_fondation
Loop#j #Nbr_descri+1
NODE no #i*100+(#j+1) X #x(#i)+#Linf_noeuds*sin(#Ix(#i))*#j Y #y(#i)+#Linf_noeuds*sin(#Iy(#i))*#j Z #Zinf_semelle-(#Linf_noeuds/((1+tan(#Ix(#i))**2+tan(#Iy(#i))**2)**0.5))*#j $ ((1+tan(#Ix(#i)*tan(#Ix(#i)+tan(#Iy(#i)*tan(#Iy(#i))**0.5))*#j
$ NODE (#i*200)+(#j+1) X #x(#i)+#Linf_noeuds*sin(#Ix(#i)) Y #y(#i)+#Linf_noeuds*sin(#Iy(#i)) Z #Zinf_semelle+#Linf_noeuds*(cos(#Ix(#i))+cos(#Ix(#i)))
Endloop
Endloop
// Noeuds pour semelle //
Loop#i semellex
NODE no 10000+#i X #semellex(#i) Y #semelley(#i) Z #Zinf_semelle+#E_semelle/2
Endloop
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES LIAISONS RIGIDES !!!!
$=========================================================================================================================================================$
Loop#i #Nbr_fondation
GRP 999
Node 70001 x 0 y 0 Z #Zinf_semelle+#E_semelle/2
NODE NO #i*100+1 FIX KF NR1 70001
Endloop
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES RESSORTS D'INTERACTION SOL STRUCTURE!!!!
$=========================================================================================================================================================$
//Ressort pour la loi d'interaction sol structure à courte durée
Loop#j #Nbr_fondation
GRP #j+1
LOOP#i #Nbr_descri+1
if #valuespalierr1X(#i) <> 0
let#k1 #i+1
let#k2 #i+401
SPRI (#i+1) #j*100+(#i+1) DX 1 0 0 MNO #k1 AR 1
SPRI (#i+101) #j*100+(#i+1) DX 0 1 0 MNO #k2 AR 1
ENDIF
ENDLOOP
ENDLOOP
$//Ressort pour la loi d'interaction sol structure à longue durée
$Loop#j #Nbr_fondation
$GRP #j+100
$ LOOP#i #Nbr_descri+1
$ if #valuespalierr1X(#i) <> 0
$ let#k1 #i+101
$ let#k2 #i+501
$ SPRI (#i+1) #j*100+(#i+1) DX 1 0 0 MNO #k1 AR 1
$ SPRI (#i+101) #j*100+(#i+1) DX 0 1 0 MNO #k2 AR 1
$ ENDIF
$ ENDLOOP
$ENDLOOP
$$//Ressort pour la loi d'interaction sol structure sous combinaison sismique
$Loop#j #Nbr_fondation
$GRP #j+200
$ LOOP#i #Nbr_descri+1
$ if #valuespalierr1X(#i) <> 0
$ let#k1 #i+301
$ let#k2 #i+601
$ SPRI (#i+1) #j*100+(#i+1) DX 1 0 0 MNO #k1 AR 1
$ SPRI (#i+101) #j*100+(#i+1) DX 0 1 0 MNO #k2 AR 1
$ ENDIF
$ ENDLOOP
$ENDLOOP
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES CONDITIONS D'APPUI A LA BASE DES ELEMENTS DE FONDATION!!!!
$=========================================================================================================================================================$
Loop#j #Nbr_fondation
GRP #j+1
SPRI 201 #j*100+(#Nbr_descri+1) DX 0 0 1 CP #Raxiale
SPRI 202 #j*100+(#Nbr_descri+1) DX 1 0 0 CM #Rrotx
SPRI 203 #j*100+(#Nbr_descri+1) DX 0 1 0 CM #Rroty
ENDLOOP
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES BARRES !!!!
$=========================================================================================================================================================$
loop#k1 #Nbr_fondation
loop#k2 #Nbr_descri
GRP #k1+1 $REF #i
BEAM (#k1*100)+(#k2+1) #k1*100+(#k2+1) #k1*100+(#k2+2) NCS 1 NP -1 $DRX -0.364851 -0.931066 0.0
endloop
endloop
END
!#!KAPITEL INTRODUCTION DES ACTIONS ET CAS DE CHARGES ELEMENTAIRES
+PROG SOFILOAD urs:3
HEAD Définition des différentes cas de charge pour l'ensemble des combinaisons de calcul
ECHO OPT ACT VAL FULL
ECHO OPT LOAD VAL FULL
$=========================================================================================================================================================$
!*! !!!! DÉFINITION DES VARIABLES !!!!
$=========================================================================================================================================================$
//Combinaisons ELU_Fondamental //
STO#FX_ELU_FOND_LD 120,120,-13.1,13.1
$Fin variable
STO#FY_ELU_FOND_LD 0,0,-1.86,-1.86
$Fin variable
STO#FZ_ELU_FOND_LD 0,-0,-486.11,-486.11
$Fin variable
STO#MX_ELU_FOND_LD 0,-0,1.94,1.94
$Fin variable
STO#MY_ELU_FOND_LD 0,0,-14.32,14.32
$Fin variable
STO#MZ_ELU_FOND_LD 0,0,0.03,-0.03
$Fin variable
STO#FX_ELU_FOND_CD -13.1,13.1,-13.1,13.1
$Fin variable
STO#FY_ELU_FOND_CD 1.86,1.86,-1.86,-1.86
$Fin variable
STO#FZ_ELU_FOND_CD -486.11,-486.11,-486.11,-486.11
$Fin variable
STO#MX_ELU_FOND_CD -1.94,-1.94,1.94,1.94
$Fin variable
STO#MY_ELU_FOND_CD -14.32,14.32,-14.32,14.32
$Fin variable
STO#MZ_ELU_FOND_CD -0.03,0.03,0.03,-0.03
$Fin variable
LOOP#i FX_ELU_FOND_LD
LET#K1 #i+1
LET#K #i+1
LC #K1 TYPE NONE TITL 'Comb.ELU fond_LD n°#K'
POIN NODE '70001' WIDE 0 TYPE PX #FX_ELU_FOND_LD(#i)
POIN NODE '70001' WIDE 0 TYPE PY #FY_ELU_FOND_LD(#i)
POIN NODE '70001' WIDE 0 TYPE PZ #FZ_ELU_FOND_LD(#i)
$POIN NODE '70001' WIDE 0 TYPE MX #MX_ELU_FOND_LD(#i)
$POIN NODE '70001' WIDE 0 TYPE MY #MY_ELU_FOND_LD(#i)
$POIN NODE '70001' WIDE 0 TYPE MZ #MZ_ELU_FOND_LD(#i)
ENDLOOP
END
!#!KAPITEL ANALYSIS DES CAS DE CHARGES
-prog ASE urs:6
echo load,forc,line,lsum yes
CTRL WARN 860 ; ctrl solv 4
lc all
end
+PROG ASE URS:4
CTRL WARN 860 ; ctrl solv 4
SYST PROB nonl; LC 1
END
-PROG ASE URS:8
CTRL WARN 860 ; ctrl solv 4
SYST PROB nonl; LC 10 ; LCC 1 fact 10
END
$-PROG SOFIMSHC urs:341.4
$HEAD Géométrie de la semelle de liaison
$unit 5
$syst SPAC gdiv 10000 gdir NEGZ
$ctrl TOPO 0
$ctrl TOLG val 0.05[m]
$ctrl MESH 1
$ctrl EDRL 1
$ctrl HMIN val 0.30[m]
$ctrl FINE val 0.30
$ctrl PROG val 1.5
$ctrl EFAC val 1.4
$$sar 1 t (#E_semelle*1000) mno 2 mctl auto nra 3 nx 0 0 1 grp 8000 h1 0.20 h2 0.20 $ thickness must be set explicitly to 0.0!
$$LOOP#i 4
$$ IF #i< 4
$$ sarb out ;slnb x1 #semellex(0) y1 #semelley(0) z1 #Zinf_semelle+#E_semelle/2 x2 #semellex(#i+1) Y2 #semelley(#i+1) Z2 #Zinf_semelle+#E_semelle/2
$$ ELSE
$$ sarb out ;slnb x1 #semellex(#i) y1 #semelley(#i) z1 #Zinf_semelle+#E_semelle/2 x2 #semellex(0) Y2 #semelley(0) Z2 #Zinf_semelle+#E_semelle/2
$$ ENDIF
$$ LET#i #i+1
$$ENDLOOP
$sar 1 t (#E_semelle*1000) mno 2 mctl auto nra 3 nx 0 0 1 grp 8000 h1 0.20 h2 0.20 $ thickness must be set explicitly to 0.0!
$sarb out ; slnb x1 0 0 10.75 x2 5 0 10.75
$sarb out ; slnb x1 5 0 10.75 x2 5 5 10.75
$sarb out ; slnb x1 5 5 10.75 x2 0 5 10.75
$sarb out ; slnb x1 0 5 10.75 x2 0 0 10.75
$end
!#!CHAPTER TEMPLATE-Principle
+prog template urs:14
head TEMPLATE-Principle
DEL#REPERAGE ; DEL#B ; DEL#A ; DEL#Nr_LC ; DEL#Ux_LC ; DEL#Uy_LC ; DEL#Uz_LC ; DEL#Urx_LC ; DEL#Ury_LC ; DEL#Urz_LC ; DEL#Px_LC ; DEL#PY_LC ; DEL#PZ_LC ; DEL#Mx_LC ; DEL#My_LC ; DEL#Mz_LC
@key mat_stee 3 $ Steel Material Nr. 2
let#fy @fy $ define variable #fy and plot variable value
txb fy = #(#fy,7.2) kN/m**2
let#cdb_len 0 $ actual length of record
let#cdb_ier 1 $ error flag
let#node 5
let#lc 1
nr ux uy uz urx ury urz Px Py A REPERAGE
[-] [mm] [mm] [mm] [rad] [rad] [rad] [kN] [kN]
@key N_DISP #lc $#node $ record with KWH/KWL 2 c'est le LC
loop#i $ loop
let#nr @(nr)
LET#ux @1*1000 //displacement X
LET#uy @2*1000 //displacement Y
LET#uz @3*1000 //displacement R
LET#urx @4*1000 //rotation X
LET#ury @5*1000 //rotation Y
LET#urz @6*1000 //rotation Z
let#px @8 $ ou @PX $ read value Px and save it to variable #px nodal support X
let#py @9 $ ou @PY $ read value PY and save it to variable #py nodal support Y
let#pz @10 $ ou @PZ $ read value PY and save it to variable #py nodal support Z
let#mx @11 $ support moment X
let#my @12 $ support moment y
let#mz @13 $ support moment z
STO#REPERAGE(#i) #nr
STO#A(#i) #py
STO#Nr_LC(#i) #nr
STO#Ux_LC(#i) #ux
STO#Uy_LC(#i) #uy
STO#Uz_LC(#i) #uz
STO#Urx_LC(#i) #urx
STO#Ury_LC(#i) #ury
STO#Urz_LC(#i) #urz
STO#Px_LC(#i) #px
STO#Py_LC(#i) #py
STO#Pz_LC(#i) #pz
STO#Mx_LC(#i) #mx
STO#My_LC(#i) #my
STO#Mz_LC(#i) #mz
if #cdb_ier<1 $ check
$endif
txb #(#nr,5.0) #(#ux,7.2) #(#uy,7.2) #(#uz,7.2) #(#urx,7.2) #(#ury,7.2) #(#urz,7.2) #(#px,7.2) #(#py,7.2) #(#A(#i), 7.3) #(#REPERAGE(#i), 7.3)
$txb PX = #(#px,5.3) kN PY = #(#py,5.3) kN $ text at the beginning of the output with formated printout
$txb PX = #(#px,5.3) kN $ text at the beginning of the output with formated printout
$txb A(#i) = #(#A(#i), 7.3)
$txb REPERAGE(#i) = #(#REPERAGE(#i), 7.3)
$txb B(#i) = #(#B(#i), 7.3)
$txb B(#i) = #(#B(#i), 7.3)
$txb B(#i) = #(#B(#i), 7.3)
endif
endloop #cdb_ier<2 $ Ende-Bedingung
end
+prog template urs:12
head Récupération des résultats dans les springs
$ N_Disp récupération des résultats des noeuds
let#cdb_len 0 $ actual length of record
let#cdb_ier 1 $ error flag
let#lc 1
nr P v
[-] [kN] [mm]
@key SPRI_RES #lc
loop#i
let#nr @(nr)
let#p @2 $ ou @P
let#v @8*1000 $ ou @v
if #cdb_ier<1 $ check
txb #(#nr,5.0) #(#p,7.2) #(#v,7.2)
endif
endloop #cdb_ier<2 $ Ende-Bedingung
end
-prog template urs:5
head filer only one node
let#lc 2
let#node 5
let#cdb_ier 0
nr urx ury urz urb
[-] [rad] [rad] [rad] [rad]
@key n_disp #lc $#node
loop
let#nr @(nr)
let#urx @(urx)
let#ury @(ury)
let#urz @(urz)
let#urb @(urb)
if #cdb_ier<2
txb #(#nr,15.0) #(#urx,15.8) #(#ury,15.8) #(#urz,15.8) #(#urb,15.8)
endif
endloop #cdb_ier<2
end
+prog template urs:26
head @NAME ; selector KWL not necessary
@key BEAM
let#nr @NR
txb selection ok with BEAM and without KWL , 1.Beam: #(#nr,5.0)
end
!#!KAPITEL INTRODUCTION DES ACTIONS ET CAS DE CHARGES ELEMENTAIRES COURTE DUREE
+PROG SOFILOAD urs:7
HEAD Définition des différentes cas de charge pour l'ensemble des combinaisons de calcul
ECHO OPT ACT VAL FULL
ECHO OPT LOAD VAL FULL
ctrl
$LOOP#i FX_ELU_FOND_LD
$LET#K1 #i+1
$LET#K #i+1
LC 1000 TYPE NONE TITL 'Essai'
LOOP#i REPERAGE
if (#REPERAGE(#i)<>0)&(#REPERAGE(#i)<>70001)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE PXX -#Px_Lc(#i)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE pyy -#Py_LC(#i)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE pzz -#Pz_LC(#i)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE mxx -#Mx_LC(#i)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE myy -#My_LC(#i)
$POIN NODE #REPERAGE(#i) WIDE 0 TYPE mzz -#Mz_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE Wxx -#Ux_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE wyy -#Uy_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE wzz -#Uz_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE dxx -#Urx_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE dyy -#Ury_LC(#i)
POIN NODE #REPERAGE(#i) WIDE 0 TYPE dzz -#Urz_LC(#i)
endif
ENDLOOP
LC 1001
copy 1 fact 1000000000000
copy 1000 fact 1
END
+PROG ASE URS:9
CTRL WARN 860 ; CTRL WARN 174; CTRL WARN 191; ctrl solv 4; SYST PROB nonl ; GRP 999 val yes; grp (1 6 1) facs 1;LC 1000
END
+PROG ASE URS:10
CTRL WARN 860 ; CTRL WARN 174; CTRL WARN 191;ctrl solv 3 ; SYST PROB nonl ; LC 1001
END
+PROG ASE URS:11
CTRL WARN 860 ; CTRL WARN 174; CTRL WARN 191;ctrl solv 4; SYST PROB nonl ; GRP 999 val yes; grp (1 6 1) facs 1;LC 2
END