opti dime 3 elem tet4 dens 5.0 ;
******** cotes
a0 = 10.0 ;
b0 = 15.0 ;
c0 = 45.0 ;
e0 = 2.0 ;
******** géométrie extérieure
p1 = 0.0         0.0         0.0 ;
p2 = a0          e0          0.0 ;
p3 = (2.0*a0)    0.0         0.0 ;
p4 = (2.0*a0+e0) b0          0.0 ;
p5 = (2.0*a0)    (2.0*b0)    0.0 ;
p6 = a0          (2.0*b0-e0) 0.0 ;
p7 = 0.0         (2.0*b0)    0.0 ;
p8 = ((-1.0)*e0) b0          0.0 ;
q1 = 0.0         0.0         c0 ;
q2 = a0          ((-1.0)*e0) c0 ;
q3 = (2.0*a0)    0.0         c0 ;
q4 = (2.0*a0-e0) b0          c0 ;
q5 = (2.0*a0)    (2.0*b0)    c0 ;
q6 = a0          (2.0*b0+e0) c0 ;
q7 = 0.0         (2.0*b0)    c0 ;
q8 = ((-1.0)*e0) b0          c0 ;
********
l1 = cer3 p1 p2 p3 ;
l2 = cer3 p3 p4 p5 ;
l3 = cer3 p5 p6 p7 ;
l4 = cer3 p7 p8 p1 ;
m1 = cer3 q1 q2 q3 ;
m2 = cer3 q3 q4 q5 ;
m3 = cer3 q5 q6 q7 ;
m4 = cer3 q7 q8 q1 ;
n1 = droi p3 q3 ;
n2 = droi p5 q5 ;
n3 = droi p7 q7 ;
n4 = droi p1 q1 ;
********
mail1 = dall l1 l2 l3 l4 quel ; trac mail1 ;
mail2 = dall m1 m2 m3 m4 quel ; trac mail2 ;
mail3 = dall l2 n1 m2 n2 quel ; trac mail3 ;
mail4 = dall l1 n4 m1 n1 quel ; trac mail4 ;
mail5 = dall l4 n3 m4 n4 quel ; trac mail5 ;
mail6 = dall l3 n2 m3 n3 quel ; trac mail6 ;
********
mail0 = (mail1 et mail2 et mail3 et mail4 et mail5 et mail6) volu ;
trac mail0 ;
trac mail0 (cont 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
cl1 = bloq uz mail1 ;
cl2 = bloq depl p1 ;
cl0 = cl1 et cl2 ;
******** force répartie constante suivant z sur mail2
fr1 = -1000.0 ;
vfr1 = 0.0 0.0 fr1 ;
******** pour le calcul
for1 = fsur mass mod0 mail2 vfr1 ;
******** pour le dessin (avec la même amplitude)
amp0 = -3.0e-3 ;
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 ;