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 - Е




—одержание  Ќазад  ¬перед