Modélisation et simulation par éléments finis : cas d'un tablier de pont.

Chapitre 3

19 format(A,x,I1,2x,A,x,$)

do j=1,3

write(20,21)c_noeud(j,i),"|"

21 format(F8.4,x,A,$)

enddo

write(20,22);22 format(2/:)

enddo

!!!!!!!!!!!! écriture la matrice Ke

write(20,23)"Ke =",E/vol,"*"

23 format(A,x,F18.4,x,A)

do i=1,12

write(20,11)"|"

11 format(x,A,x,$)

do j=1,12

write(20,12)ke(i,j)

12 format(F8.4,x,"|",$)

enddo

write(20,13);13 format(/:)

enddo

close(20)

call system('start c:\matrice_K_B\metrice_K.txt')

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

end program

!

!procédure pour la table de localisation pour les matrices Bi

!cette sous routine permet d'insérer la matrice identité à la

!place de Bi

!

subroutine table_Bi(col,tab)

integer i,j

real::col

real, dimension(1:3,1:12)::tab

tab=0;j=3*col-2

do i=1,3

tab(i,j)=1;j=j+1

enddo

return

end

!procédure pour calculer les termes ji la matrice jacobienne inverse

!

! procédure pour calculer le déterminant de la matrice

! de la matrice jacobienne de transformation

!warning: cette procédure est pour ce cas spécificique

!il ne peut être utilisé pour calculer le déterminant d'1e autre matrice

function det_jab(x_1,y_1,z_1,x_2,y_2,z_2,x_3,y_3,z_3,x_4,y_4,z_4)

real:: x_1,y_1,z_1,x_2,y_2,z_2,x_3,y_3,z_3,x_4,y_4,z_4,det_jab

det_jab=((x_2-x_1)*(y_3-y_1)*(z_4-z_1)+(y_2-y_1)*(z_3-z_1)*(x_4-x_1)+&

&(z_2-z_1)*(x_3-x_1)*(y_4-y_1)-(x_2-x_1)*(z_3-z_1)*(y_4-y_1)-(y_2-&

&y_1)*(x_3-x_1)*(z_4-z_1)-(z_2-z_1)*(y_3-y_1)*(x_4-x_1))

return

end function det_jab

!

subroutine produit(mat1,mat2,n,c,n1,c1,res)

integer::i,j,k,n,c,n1,c1

real, dimension(1:n,1:c1)::res

real, dimension(1:n,1:c)::mat1

real, dimension(1:n1,1:c1)::mat2

do i=1,n

do j=1,c1

res(i,j)=0.

do k=1,c

res(i,j)=res(i,j)+mat1(i,k)*mat2(k,j)

enddo

enddo

Chapitre 3

enddo

end subroutine produit

3.3.5.3. Exemple de l'exécution du programme Input :

,

Noeu n°	x	y	z
1	0.75	0.33	0.1
2	1.01	0.1	0.05
3	0.02	1.11	0.12
4	0.105	0	1.2

Tableau 3.2 : Valeurs d'essai pour le test du programme matrice_K

Figure 3.7 : Capture d'écran des données en Input du programme matrice_K

88 /176

Output :

Chapitre 3

89 /176

Figure 3.8 : Capture d'écran des données en Output du programme matrice_K

90 /176

Chapitre 4

Chapitre 4 : Modélisation et Simulation

numérique d'un tablier de pont

Sommaire

4.1.	Matériels employés pour la simulation	91
4.2.	Présentation générale de l'ouvrage	92
4.3.	Caractéristiques du Tablier	94
4.4.	Définition des charges et actions appliquées à la structure	99
4.5.	Définition du flux de travail	109
4.6.	Définition de la structure	110
4.7.	Construction du modèle éléments finis (EF)	116
4.8.	Introduction des conditions de fixations (Appuis)	118
4.9.	Choix des normes et règlements à utiliser	119
4.10.	Définitions des charges	120
4.11.	Lancement des calculs de la structure	128
4.12.	Résultats de calcul	129
4.13.	Définition des combinaisons d'actions	134
4.14.	Exploitation des résultats	136

91 /176