Formules pour les éléments des matrices orthotropes (orthotropie du matériau)

Convention de signes :

E - Module de Young

v - coefficient de Poisson

G - Module de cisaillement.

Raidisseurs d'un côté dans une direction :

h - épaisseur de la dalle

ha - hauteur des raidisseurs

a - espacement des raidisseurs

a1 - largeur des raidisseurs.

Rigidité - état de membrane :

adoptée : s1 = a*h, s2 = (ha-h)*a1

D_XXXX = E*(s1+s2)/a

D_YYYY = E*h

D_XXYY = v*E*h

D_XYXY = G*h

D_XXXY = 0

D_YYXY = 0

Rigidité de flexion :

adoptée : hz = h + a1/12, xv = 1 - v*v

Dyyyy= E/(12*xv)*a/((a-a1)/(h*h*h)+a1/(hz*hz*hz))

I = moment d’inertie du profilé en T dont les caractéristiques sont : longueur de l’aile a/xv, épaisseur de l’aile h, hauteur totale ha et épaisseur de l’âme a1

K_XXXX = E*I/a

K_YYYY = Dyyyy

K_XXYY = v*Dyyyy

K_XYXY = (1-v)/2*Dyyyy

K_XXXY = 0

K_YYXY = 0

Rigidité de cisaillement :

H_XX = 5/6*G*((a-a1)*h + a1*ha)/a

H_YY = 5/6*G*h

H_XY = 0

 

Raidisseurs de deux côtés dans une direction :

h - épaisseur de la dalle

ha - hauteur des raidisseurs

a - espacement des raidisseurs

a1 - largeur des raidisseurs.

Rigidité - état de membrane :

adoptée : s1 = a*h, s2 = (ha-h)*a1

D_XXXX = E*(s1+s2)/a

D_YYYY = E*h

D_XXYY = v*E*h

D_XYXY = G*h

D_XXXY = 0

D_YYXY = 0

Rigidité de flexion :

adoptée : hz = h + a1/6, xv = 1 - v*v

Dyyyy= E/(12*xv)*a/((a-a1)/(h*h*h)+a1/(hz*hz*hz))

I = moment d’inertie du profilé en croix dont les caractéristiques sont : longueur de l’aile a/xv, épaisseur de l’aile h, hauteur totale ha et épaisseur de l’âme a1

K_XXXX = E*I/a;

K_YYYY = Dyyyy;

K_XXYY = v*Dyyyy;

K_XYXY = (1-v)/2*Dyyyy

K_XXXY = 0;

K_YYXY = 0;

Rigidité de cisaillement :

H_XX = 5/6*G*((a-a1)*h + a1*ha)/a

H_YY = 5/6*G*h

H_XY = 0

 

Raidisseurs de deux côtés dans une direction :

h - épaisseur de la dalle

ha - hauteur des raidisseurs

hb - hauteur des raidisseurs

a - espacement des raidisseurs

a1 - largeur des raidisseurs.

b - espacement des raidisseurs

b1 - largeur des raidisseurs.

Rigidité - état de membrane :

On admet :

s1x = a*h, s2x = (h-ha)*a1

s1y = b*h, s2y = (h-hb)*b1

D_XXXX = E*(s1x+s2x)/a

D_YYYY = E*(s1y+s2y)/b

D_XXYY = v*E

D_XYXY = G*h

D_XXXY = 0

D_YYXY = 0

Rigidité de flexion:

si ha>=hb, to

a2 = a1, a1 = a-a1, h2 = ha, b2 = b1, b1 = b-b2, h1 = hb

dans le cas contraire (hb>ha), alors

a2 = b1, a1 = b-b1, h2 = hb, b2 = a1, b1 = a-b2, h1 = ha

I11 = InertiaTsec(b/xv, b2, hp, h1)

I12 = InertiaTsec(b/xv, b2+a2/6, hp+b2/12, hp+a/12)

I21 = InertiaTsec(a/xv, a2, hp, h2)

I22 = InertiaTsec(a/xv, (a2+a2+b2/6)/2, hp+b2/12, h2)

xv = 1 - v*v

GC1 = G*TorsI(b2,h1)

GC2 = G*TorsI(a2,h2)

D11 = E*a/b/(a1/I11+a2/I12)

D22 = E*b/a/(b1/I21+b2/I22)

K_XXXX = (ha>=hb) D22:D11

K_YYYY = (HA > =hb) D11 :D22

K_XXYY = v*E/12.0*(h3/xv+ h1*h1*h1 *a2*b2/(a*b))

K_XYXY = G*h3/12.0 + (GC1/b+GC2/a)/4.0

K_XXXY = 0.0

K_YYXY = 0.0

xv = 1 - v*v

Rigidité de cisaillement :

H_XX = 5/6*G*((a-a1)*h + a1*ha)/a;

H_YY = 5/6*G*h;

H_XY = 0;

 

Plancher alvéolaire à raidisseurs dans une direction :

h - épaisseur du plancher

h1 - épaisseur de la plaque supérieure

h2 - épaisseur de la plaque inférieure

a - espacement des raidisseurs

a1 - largeur des raidisseurs

Rigidité - état de membrane :

On admet :

s1 = a*h1, s2 = a*h2, s0 = (h-h1-h2)*a1

D_XXXX = E*(s1+s2+s0)/a

D_YYYY = E*(h1+h2)

D_XXYY = v*E*(h1+h2)

D_XYXY = G*(h1+h2)

D_XXXY = 0

D_YYXY = 0

Rigidité de flexion:

On admet :

hz1 = h1 + a1/12, hz2 = h2 + a1/12, xv = 1 - v*v

Dyyyy= a/( (a-a1)/InertiaIsec(1/xv,0,h,h1,h2) + a1/InertiaIsec(1/xv,0,h,hz1,hz2)

Iy = InertiaIsec(a/xv,a1,h,h1,h2)

K_XXXX = E*Iy/a

K_YYYY = Dyyyy

K_XXYY = v*Dyyyy

K_XYXY = (1-v)/2*Dyyyy

K_XXXY = 0

K_YYXY = 0

Rigidité de cisaillement :

H_XX = 5/6*G*((a-a1)*(h1+h2) + a1*ha)/a

H_YY = 5/6*G*(h1+h2)

H_XY = 0

 

Plancher alvéolaire à raidisseurs dans deux directions :

h - épaisseur du plancher

h1 - épaisseur de la plaque supérieure

h2 - épaisseur de la plaque inférieure

a - espacement des raidisseurs

a1 - largeur des raidisseurs.

b - espacement des raidisseurs

b1 - largeur des raidisseurs.

Rigidité - état de membrane :

On admet :

s1 = a*h1, s2 = a*h2, s0x = (h-h1-h2)*a1, s0y = (h-h1-h2)*b1

D_XXXX = E*(s1+s2+s0x)/a

D_YYYY = E*(s1+s2+s0y)/b

D_XXYY = v*E*(h1+h2)

D_XYXY = G*(h1+h2)

D_XXXY = 0.0

D_YYXY = 0.0

Rigidité de flexion:

On admet :

a2 = a1, a1 = a-a2, double b2 = b1, b1 = b-b2, xv = 1 - v*v

GC1 = G*TorsI(a2,h-h1-h2)

GC2 = G*TorsI(b2,h-h1-h2)

Dpl = InertiaIsec(1/xv,0,h,h1,h2)

I11 = InertiaIsec(b/xv, b2, h, h1, h2)

I12 = InertiaIsec(b/xv, b2+a2/6.0, h, h1+a2/12.0, h2+a2/12.0)

I21 = InertiaIsec(a/xv, a2, h, h1, h2)

I22 = InertiaIsec(a/xv, a2+b2/6.0, h, h1+b2/12.0, h2+b2/12.0)

D11 = E*a/b/(a1/I11+a2/I12)

D22 = E*b/a/(b1/I21+b2/I22)

K_XXXX = D22

K_YYYY = D11

K_XXYY = v*E*(Dpl+h3/12.0*a2*b2/(a*b))

K_XYXY = (1.0-v)/2.0*E*Dpl+(GC1/b+GC2/a)/4.0

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = 5.0/6.0*G*((a-a1)*(h1+h2) + a1*h)/a

H_YY = 5.0/6.0*G*((b-b1)*(h1+h2) + b1*h)/b

H_XY = 0.0

 

Radier partiel :

h - hauteur des raidisseurs

a - espacement des raidisseurs

a1 - largeur des raidisseurs.

b - espacement des raidisseurs

b1 - largeur des raidisseurs.

Rigidité - état de membrane :

D_XXXX = E*h*a1/a

D_YYYY = E*h*b1/b

D_XXYY = 0.0

D_XYXY = G*h*(1.-(a-a1)*(b-b1)/(a*b))

D_XXXY = 0.0

D_YYXY = 0.0

Rigidité de flexion:

a2 = a1; a1 = a-a2; b2 = b1; b1 = b-b2

D=E*h3/12, a23=a2*a2*a2, b23=b2*b2*b2

K_XXXX = D * a2/a

K_YYYY = D * b2/b

K_XXYY = v*D * (a2*b2)/(a*b)

K_XYXY = D*(1-v*v*a2*b2/(a*b))/(a*b*(a/a23+b/b23))

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = 5.0/6.0*G*h*a1/a

H_YY = 5.0/6.0*G*h*b1/b

H_XY = 0.0

 

Bandeaux sur bacs acier :

h - épaisseur de la dalle

h1 - hauteur du bac acier

a - espacement des nervures du bac acier

a1 - largeur inférieure des nervures

a2 - largeur supérieure des nervures

Les valeurs des rigidités D, K et H sont calculées de la même façon que pour le type ‘Raidisseur d’un côté dans une direction’ avec les grandeurs suivantes :

a1 = (a1+a2)/2

ha = h+h1.

Orthotropie du matériau :

h - épaisseur de la dalle

n1 - coefficient de rigidité

n2 - coefficient de rigidité

On admet :

E1= E*n1 E2= E*n2

G1= E1*/(2*(1+v)) G2= E2*/(2*(1+v))

G12=sqrt(E1*E2)/(2*(1+v))

H3= (h*h*h)*/(12*(1-v*v))

Rigidité - état de membrane :

D_XXXX = E1 * h/(1-v*v)

D_YYYY = E2 * h/(1-v*v)

D_XXYY = v * h * sqrt ( E1 * E2)/(1- v*v)

D_XYXY = G12 * h

D_XXXY = 0.0

D_YYXY = 0.0

Rigidité de flexion:

K_XXXX = E1 * H3 = K11

K_YYYY = E2 * H3 = K22

K_XXYY = v * sqrt ( K11 * K22)

K_XYXY = (1-v)/2*sqrt(K11*K22) = (h*h*h)/12*G12

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = 5/6 * G1 * h

H_YY = 5/6 * G2 * h

H_XY = 0.0

 

Tôle trapézoïdale Tôle trapézoïdale :

h hauteur de la tôle

t épaisseur de la tôle

a, a1, a2 – dimensions d’un segment de la tôle (espacement, épaisseurs)

On admet :

v2= v*v

B = E*t*t*t /(12*(1-v2))

M = E*t/(1-v2)

b = a-a2

c = (a2-a1)/2

lx = sqrt(h*h+c*c)

ca = c/lx

Ax = a1+b+2*lx

e = h*(lx+a1)/Ax

Ix = e*e*b + 2*lx/3 (h*h+3*e*(e-h))+ (e-h)*(e-h)*a1

As = b+2*ca*ca*lx+a1

Rigidité - état de membrane :

D_XXXX = M*Ax/a

D_YYYY = a/(Ix/B + As/M)

D_XXYY = D_YYYY*v

D_XYXY = G*t*a/Ax

D_XXXY = 0.0

D_YYXY = 0.0

Rigidité de flexion:

K_XXXX = (M*Ix+B*As)/a

K_YYYY = B*a/Ax

K_XYXY = G*t/a*(Ix + t*t/12*(As*(1+Ax/a)-a))

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = G*t*6/5*Ax/kx/a

H_XY = 0.0

kx = Ax/(2*lx)

 

Tôle ondulée Tôle ondulée :

h hauteur de la tôle

t épaisseur de la tôle

a – dimensions d’un segment de la tôle (espacement)

On admet :

L’équation définissant la géométrie de la tôle :

Rigidité - état de membrane :

D_XXXX = M*I1/2p

D_YYYY = 2p M/(I2+ 12*(h/t)*(b/t)*I3)

D_XXYY = D_YYYY*v

D_XYXY = 2p G*t /I1

D_YYXY = 0.0

Rigidité de flexion:

K_XXXX = M* b*b*I3/2?

K_YYYY = 2p B /I1

K_XYXY = G*t*b*b*I3/2p+G*t*t*t/12*(I2/2p *(1+I1/2p)-1)

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = G*t*5/6*I1/2p

H_XY = 0.0

 

Dalle composée avec tôle trapézoïdale

h - épaisseur de la dalle

h1 - hauteur de la tôle trapézoïdale

t - épaisseur de la tôle trapézoïdale

a, a1, a2 - dimensions d’un simple module de la tôle (espacement, largeurs)

Es - module de Young (acier trapézoïdale)

vs - coefficient de Poisson (acier trapézoïdale)

Eb - module de Young (béton ; valeur adoptée sur la base du matériau sélectionné)

vb - coefficient de Poisson (béton ; valeur adoptée sur la base du matériau sélectionné)

Principes :

b = a - a1

e = 0.5*(a-(b+a1))

Remarque : la valeur du paramètre e peut être inférieure à zéro, c’est-à-dire, on admet une dalle / tôle ayant la forme présentée sur la figure ci-dessous

s = sqrt(e*e+h1*h1)

Ab1 = a*h

Ab2 = (a1+e)*h1

As = t*(b+2*s+a1)

Ebx = Eb/(1.0-vb*vb)

Esx = Es/(1.0-vs*vs)

Gb = Eb/(2.0*(1.0+vb))

Gs = Es/(2.0*(1.0+vs))

hx - hauteur équivalente du béton égale :

si h+(a-b)/6.0 < h+h1, to h+(a-b)/12.0

si h+(a-b)/6.0 > h+h1, to h+h1/2

hx3 = hx*hx*hx

y0 = (Ebx*(Ab1*(h1+h/2)+a1* h1* h1/2+2*e*h1* h1/3.0)+Es*t*(b*h1+s*h1))/(Ebx*(Ab1+Ab2)+Es*As)

h3 = h*h*h

d3 = h1*h1* h1

e1 = h1+h/2 -y0

e1 = e1*e1

e2 = y0- h1/2

e2 = e2*e2

e3 = y0-2*h1/3.0

e3 = e3*e3

y1 = Ebx*h*h/(2*(Ebx*h+Esx*t));

hd = h+0.5*h1

hd3 = hd*hd*hd

Rigidité - état de membrane :

D_XXXX = (Ebx*(Ab1+Ab2)+As*Esx)/a

D_YYYY = a/(b/(h*Ebx+t*Esx)+(a-b)/(hx*Ebx))

D_XXYY = D_YYYY*vb

D_XYXY = Gb*a/((c+e)/(h+h1)+(b+e)/h) + Gs*t*a/(As/t)

D_XXXY = 0.0

D_YYXY = 0.0

Rigidité de flexion :

K_XXXX = a/(b/(Ebx*(h3/12.0 + h*(y1-h/2)*(y1-h/2)) + Esx*t*y1*y1) + (a-b)/(Ebx*hx3/12.0))

K_YYYY = a/(b/(Ebx*(h3/12.0 + h*(y1-h/2)*(y1-h/2)) + Esx*t*y1*y1) + (a-b)/(Ebx*hx3/12.0))

K_XYXY = = Gb*a/(b/h3 + (a-b)/hd3)/12.0

K_XXXY = 0.0

K_YYXY = 0.0

Rigidité de cisaillement :

H_XX = Gb*(Ab1+Ab2)/a*5.0/6.0

H_XY = 0.0

 

Plancher alvéolaire

ha = h - h1 - h2 hauteur du raidisseur

Rigidité - état de membrane :

s1 = a*h1, s2 = a*h2, s0 = (h-h1-h2)*a-(3,14*a1*a1/4+c*a1)

D_XXXX = E/xv*(s1+s2+s0)/a

D_YYYY = E/xv *(h1+h2)

D_XXYY = v*E/xv *(h1+h2)

D_XYXY = G*(h1+h2)

D_XXXY = 0

D_YYXY = 0

Rigidité de flexion :

hz1 = h1 + a1/12, hz2 = h2 + a1/12, xv = 1 - v*v, r = a1/2

Dyyyy= E*a/((a-a1)/(InertiaIsec(1/xv,0,h,h1,h2) + a1/InertiaIsec(1/xv,0,h,hz1,hz2))

Iy = InertiaIsec(a/xv,a1,h,h1,h2) + 4*(0,00755*r*r*r*r + 0,2146*r*r*(c/2+0,7766*r)*(c/2+0,7766*r))

Si la rigidité de la flexion dans la direction Y a été ignorée, alors Dyyyy = 1,0*e-4*Dyyyy

Sinon, Dyyyy = Dyyyy

K_XXXX = E*Iy/a

K_YYYY = Dyyyy

K_XXYY = v*Dyyyy

K_XYXY = (1-v)/2*Dyyyy

K_XXXY = 0

K_YYXY = 0

Rigidité de cisaillement :

H_YY = 5/6*G*(h1+h2)

H_XY = 0