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


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


/p>

Рис. 3.2

3.2.1. Определение граничных параметров при статической

           нагрузке

Исходные данные расчета балки:

l1 = 4.0; l2

= 6.0; l3 = 3.0; l4 = 1.0; m = 20.0; q

= 10.0; f1 = - 40.0; f2 = 60.0;          am = 2.0; af1

= 2.0; af2 = 1.0.

Система линейных уравнений для граничных параметров балки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.11)

Программа расчета граничных параметров балки

Протокол с комментариями

% Описание матрицы а(16,16); b(16,1) и Х(16,1);

а = zeros(16,16); ); b = zeros

(16,1); Х = zeros

(16,1);

% Исходные данные расчета балки

l1 = 4.0; l2 = 6.0; l3 = 3.0; l4 = 1.0; m

= 20.0; q = 10.0;

 f1 = - 40.0; f2 = 60.0;  am = 2.0; af1 = 2.0; af2 = 1.0;

% Ввод матриц а(16,16) и b(16,1) с помощью операторов присваивания;

a(1,2) = l1; a(1,4) = -l1^3/6; b(1,1) = -m*(l1-am)^2/2;

a(2,2) = 1.0;  a(2,4) = -l1^2/2;

a(2,6) = -1.0;  b(2,1) = -m*(l1-am);

a(3,4) = l1;  a(3,7) = -1.0;  b(3,1) = m; a(4,1) = -1.0;  a(4,4) = 1.0; 

a(5,6) = l2;  a(5,7) = -l2^2/2;  a(5,8) = -l2^3/6;  b(5,1) = -q*l2 ^4/24;

a(6,6) = 1.0;  a(6,7) = -l2;  a(6,8) = -l2^2/2;  a(6,10) = -1.0; 

b(6,1) = -q*l2 ^3/6;

a(7,7) = 1.0;  a(7,8) = l2;  a(7,11) = -1.0; 

b(7,1) = q*l2^2/2;  a(8,3) = -1.0;  a(8,8) = 1.0;  b(8,1) = q*l2;

a(9,10) = l3;  a(9,11) = -l3^2/2;  a(9,12) = -l3^3/6;  

b(9,1) = -f1*(l3-af1)^3/6;  a(10,10) = 1.0;  a(10,11) = -l3;

a(10,12) = -l1^2/2;  a(10,14) = -1.0;  b(10,1) = -f1*(l3-af1)^2/2;  

a(11,11) = 1.0;  a(11,12) = l3;  a(11,15) = -1.0;  b(11,1) = f1*(l3-af1);  

a(12,5) = -1.0;  a(12,12) = 1.0;  ;  b(12,1) = f1;  a(13,9) = -1.0; 

a(13,14) = l4;  a(13,15) = -l4^2/2;  a(13,16) = -l4^3/6;  

b(13,1) = -f2*(l4-af2)^3/6;  a(14,13) = -1.0;  a(14,14) = 1.0; 

a(14,15) = -l4;   a(14,16) = -l4^2/2;  b(14,1) = -f2*(l4-af2)^2/2;  

a(15,15) = -1.0;   a(15,16) = l4;  b(15,1) = f2*(l4-af2);  a(16,16) = 1.0;  

b(16,1) = f2;

% Решение системы уравнений ах = b и вывод результатов

Х = a\ b

<


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



Книжный магазин