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


Продолжение таблицы 2.3 - часть 29


/p>

 

 

 

3.2.6. Построение эпюр напряженно-деформированного состояния

                    балки при вынужденных колебаниях

 

Эпюры напряженно-деформированного состояния можно построить для каждого стержня в отдельности, используя соотношения метода начальных параметров.

 

 

 

 

 

 

 

 

(3.21)

 

где выражения от внешней нагрузки имеют вид

 

 

 

 

 

 

(3.22)

 

 

 

 

 

 

 

(3.22)

 

 

- единичная функция Хевисайда.

 

M, F, q – сосредоточенные момент, сила и распределенная нагрузка.

аМ, аF, аH, аK – координаты внешней нагрузки.

Протоколы построения эпюр EIJ(x), EIj(x), M(x) и Q(x) примут вид

стержень 0-1

х = 0 : 0.0001 : 2.0;  la = sqrt (0.1*0.315365);

EIv = - (X(2,1)*(sinh(la*x)+sin(la*x))/(2*la) - X(4,1)* …

(sinh(la*x)-sin(la*x))/(2*la^3));

EIfi = - (X(2,1)*(cosh(la*x)+cos(la*x)) / 2 - X(4,1)* …

(cosh(la*x)-cos(la*x))/(2*la^2));

Q = - X(2,1)*la^4*(cosh(la*x)-cos(la*x)) / (2*la^2) + …

X(4,1)* (cosh(la*x)+cos(la*x)) / 2;

M = - X(2,1)*la^4*(sinh(la*x)-sin(la*x))/(2*la^3) + …

X(4,1)*(sinh(la*x)+sin(la*x))/(2*la);

subplot

(2,2,1), plot (x, EIv); axis

([0 2 - 30  30]); grid on

subplot

(2,2,2), plot (x, EIfi); axis ([0 2 - 30  30]); grid on

subplot

(2,2,3), plot (x, Q); axis

([0 2 - 10  10]); grid on

subplot

(2,2,4), plot (x, M); axis

([0 2 - 3    3 ]); grid on

Перед выполнением этого протокола необходимо поместить в окно команд вектор граничных параметров Х (см. п. 3.2.5), значения нагрузки и их координаты.

х = 2.0 : 0.0001 : 4.0; m

= 20.0; am = 2.0; la = sqrt (0.1*0.315365);

EIv = - (X(2,1)*(sinh(la*x)+sin(la*x))/(2*la) - X(4,1)* …

(sinh(la*x)-sin(la*x))/(2*la^3)+m*(cosh(la*(x-am)) - …

cos (la*(x-am))) / (2*la^2));

EIfi = - (X(2,1)*(cosh(la*x)+cos(la*x)) / 2 -




Начало  Назад  Вперед