opti dime 3 elem tet4 dens 2.0 ;
******** cotes
a0 = 10.0 ;
b0 = 15.0 ;
c0 = 25.0 ;
******** géométrie extérieure
p1 = 0.0 0.0 0.0 ;
p2 = a0  0.0 0.0 ;
p3 = a0  b0  0.0 ;
p4 = 0.0 b0  0.0 ;
p5 = a0  0.0 c0  ;
p6 = a0  b0  c0  ;
********
l1 = droi p1 p2 ;
l2 = droi p2 p3 ;
l3 = droi p3 p4 ;
l4 = droi p4 p1 ;
cont1 = l1 et l2 et l3 et l4 ;
l5 = droi p2 p5 ;
l6 = droi p5 p6 ;
l7 = droi p6 p3 ;
cont2 = l5 et l6 et l7 et l2 ;
l8 = droi p6 p4 ;
cont3 = l3 et l7 et l8 ;
l9 = droi p5 p1 ;
cont4 = l1 et l5 et l9 ;
cont5 = l4 et l9 et l6 et l8 ;
******** maillage
mail1 = surf cont1 plan ; trac mail1 ;
mail2 = surf cont2 plan ; trac mail2 ;
mail3 = surf cont3 plan ; trac mail3 ;
mail4 = surf cont4 plan ; trac mail4 ;
mail5 = surf cont5 plan ; trac mail5 ;
mail0 = (mail1 et mail2 et mail3 et mail4 et mail5) volu ;
trac mail0 ;
********
nbno0 = nbno mail0 ; mess 'Nombre de noeuds :' nbno0 ;
********************* modèle (homogène et linéaire sous entendu)
mod0 = mode mail0 mecanique elastique isotrope ;
mat0 = mate mod0 youn 210000.0 nu 0.30 ;
********************* conditions limites
cl0 = bloq depl mail1 ;
******** force répartie constante suivant x sur mail2
frx1 = -10.0 ;
vfr1 = frx1 0.0 0.0 ;
******** pour le calcul
for1 = fsur mass mod0 mail2 vfr1 ;
******** pour le dessin (avec la même amplitude)
amp0 = -3.0e-2 ;
vect1 = vect for1 amp0 fx fy fz vert ;
trac vect1 mail0 titr 'Chargements' ;
********************* matrice de raideur
rig0 = rigi mod0 mat0 ;
******** prise en compte des conditions limites dans la matrice de raideur
rig0 = rig0 et cl0 ;
********************* résolution de [K][q]=[F] avec
******** [K] = rig0 ; [F] = for0 et [q] = res0 contient les ddl (u et v)
******** en chacun des noeuds du maillage
res0 = reso rig0 for1 ;
********************* déformée
def0 = defo res0 mail0 0.0 blanc ;
def1 = defo res0 mail0     rouge ;
trac (def0 et def1) titre 'Déformée amplifiée' ;
def1 = defo res0 mail0 1.0 rouge ;
trac (def0 et def1) titre 'Déformée à l échelle de la structure' ;
******** action aux liaisons
***** reac0 contient les forces nodales pour les noeuds où l'on a imposé
***** les conditions aux limites
reac0 = reac rig0 res0 ;
***** pour le dessin on met la même amplitude
***** que pour les forces extérieurs imposées
vreac0 = vect reac0 amp0 fx fy fz vert ;
trac vreac0 mail0 titr 'Action aux liaisons';
***** on veut avoir la force globale (la resuLTANTE de ces forces nodales)
fg0 = resu reac0 ;
mess 'Force globale :' ; list fg0 ;
*
*********************** calcul des contraintes dans l'élément
****** sig0 contiendra les valeurs du tenseur des contraintes
****** pour chacun des éléments du maillage
sig0 = sigm mod0 mat0 res0 ;
****** et signo0 les valeurs du tenseur des contraintes
****** pour chacun des noeuds du maillage
signo0 = chan chpo mod0 sig0 ;
******** tracé des composantes du tenseur des contraintes
******** moyenné aux noeuds
sxxno0 = exco smxx signo0 ;
trac mail0 sxxno0 titr 'Sigma xx' (cont mail0) ;
sxyno0 = exco smxy signo0 ;
trac mail0 sxyno0 titr 'Sigma xy' (cont mail0) ;
syyno0 = exco smyy signo0 ;
trac mail0 syyno0 titr 'Sigma yy' (cont mail0) ;
syzno0 = exco smyz signo0 ;
trac mail0 syzno0 titr 'Sigma yz' (cont mail0) ;
sxzno0 = exco smxz signo0 ;
trac mail0 sxzno0 titr 'Sigma xz' (cont mail0) ;
szzno0 = exco smzz signo0 ;
trac mail0 szzno0 titr 'Sigma zz' (cont mail0) ;
fin ;