Scribd
Importer
44 vues2 pages

Code Matlab

Le document présente un script MATLAB utilisant la méthode de Newton-Raphson pour déterminer le diamètre minimal requis d'un élément soumis à un effort axial et un moment de flexion. Les données incluent l'effort axial, le moment de flexion, la limite élastique et le coefficient de sécurité. Le script itère jusqu'à atteindre une tolérance de convergence ou un nombre maximum d'itérations pour calculer le diamètre.

Transféré par

Naminawado Géraud
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
44 vues2 pages

Code Matlab

Le document présente un script MATLAB utilisant la méthode de Newton-Raphson pour déterminer le diamètre minimal requis d'un élément soumis à un effort axial et un moment de flexion. Les données incluent l'effort axial, le moment de flexion, la limite élastique et le coefficient de sécurité. Le script itère jusqu'à atteindre une tolérance de convergence ou un nombre maximum d'itérations pour calculer le diamètre.

Transféré par

Naminawado Géraud
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

clc; clear;

% ==== Données ====


N = 2500; % Effort axial (N)
M = 300; % Moment de flexion (N.m)
Re = 360e6; % Limite élastique (Pa)
gamma_s = 3; % Coefficient de sécurité

% ==== Constantes ====


A = 4 * N; % Numérateur traction
B = 32 * M; % Numérateur flexion
D = pi * (Re / gamma_s); % Terme constant côté droit

% ==== Paramètres Newton-Raphson ====


d = 0.02; % Valeur initiale (m)
tol = 1e-6; % Tolérance de convergence
max_iter = 100; % Nombre max d'itérations

% ==== Méthode de Newton-Raphson ====


for i = 1:max_iter
f = A / d^2 + B / d^3 - D;
df = -2*A / d^3 - 3*B / d^4; % Dérivée de f par
rapport à d

d_new = d - f / df;

if abs(d_new - d) < tol


fprintf('Convergence atteinte après %d
itérations.\n', i);
break;
end

d = d_new;
end

% ==== Résultat ====


fprintf('Diamètre minimal requis : %.6f m\n', d);

Vous aimerez peut-être aussi