Bonjour,
Dans le cadre de calculs sur un plancher en struture sandwich (matériaux composites), je cherche à calculer une déflection (w) en fonction de la longueur du plancher (x) et de sa largeur (y). Je bute simplement sur lécriture de l'algorithme et ce serait vraiment très gentil de votre part si vous pouviez m'apporter un peu d'aide.
La formule mathématique simplifiée est:
w= somme (pour m=1 jusqu'à 1000) somme (pour n=1 jusqu'à 1000) sin(m*pi*x)*sin(n*pi*y).
Je souhaiterais calculer w en fonction de x et de y. Ci dessous, j'ai copié-collé ce que j'ai écrit. J'ai tracé en 3D, w(x,y) et cela me donne n'importe quoi: je n'ai pas une belle symétrie et la déflextion maximale ne se trouve pas au centre!
Merci beaucoup.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26 clear all close all clc a=2060 b=1400 x=0:a/10:a y=0:b/10:b w=zeros(length(x),length(y)) for i=1:length(x) for j=1:length(y) for m=1:50 for n=1:50 w_temp=sin(m*pi*x(i))*sin(n*pi*y(j)); w(i,j)=w(i,j)+w_temp end end end end %Plot the deflection (w) versus (x,y) [x,y] = meshgrid(x,y); surf(x,y,-w)
Beauvaim.
Partager