MATLAB в инженерных и научных расчетах




Таблица 3.9 - часть 4


(–34.08093)*(x – 4).^3*h2/6 + 10*(x – 4).^4*h2/24 – …

10*(x – 10).^4*h3/24 + 4.36214*(x – 10).^3*h3/6 – …

40*(x – 12).^3*h4/6 + (–51.33059)*(x – 13).^3*h5/6);

EIfi(m,1)= – (–12.79835 – (–1.04938)*x.^2/2 + 20*(x – 2)*h1 + …

(–34.08093)*(x – 4).^2*h2/2 + 10*(x – 4).^3*h2/6 – …

10*(x – 10).^3*h3/6 + 4.36214*(x – 10).^2*h3/2 – …

20*(x – 12).^2*h4 + (–51.33059)*(x – 13).^2*h5/2);

Q(m,1)= (–1.04938) – (–34.08093)*h2 - 10*(x – 4)*h2 + …

10*(x – 10)*h3 - 4.36214*h3 + 40*h4 – …

(-

51.33059)*h5;

M(m,1)= (–1.04938)*x -

20*h1-(34.08093)*(x – 4)*h2 - …

- 5*(x – 4).^2*h2 + 5*(x – 10).^2*h3 - (4.36214)*(x – 10)*h3 + …

40*(x – 12)*h4 - (- 51.33059)*(x – 13)*h5; X(m,1)=x;

x = x + dx; end;

plot(X, EIv); grid on

% plot(X, Q); grid on

% plot(X, EIfi); grid on

% plot(X, M); grid on

[X  EIv  EIfi   Q   M]

На экране компьютера появится окно с первой эпюрой, совпадающей с рис. 3.1, затем, последовательно снимая символ комментария %, строятся другие эпюры, а в окно команд выводится таблица численных значений кинематических и статических параметров балки. Добавим, что совершенно аналогично можно строить эпюры с разрывами 1-го рода и для любых стержней рам. Если вместо оператора plot(X, EIv) использовать оператор stem(X, EIv,¢filled¢), то получится закрашенная эпюра, мало отличающаяся от обычного представления такого рода графиков. Можно сократить число используемых единичных функций до одной. Однако в этом случае единичную функцию необходимо создать в отдельном М-файле, например, такого содержания

function f = H(t)

if  t < = 0.0  f = 0.0;  else f = 1.0; end;

Этот файл должен быть записан в рабочей папке MATLAB (обычно это папка work) под именем Н. Тогда в головной программе построения эпюр можно использовать одну единичную функцию, но от разных аргументов, т.е.

Н(х – 2), Н(х – 4), Н(х – 10) и т.д. Например, выражение для прогиба (3.36) примет вид

EIv(m,1) = - (- 12.79835*x - (- 1.04938)*x.^3/6 + 10*(x – 2).^2*H(x – 2) + …

(-

34.08093)*(x – 4).^3*H(x – 4)/6 + 10*(x – 4).^4*H(x – 4)/24 - …




Содержание  Назад  Вперед