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


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


 

 

 

 

 

 

 

1

 

Обозначения переменных, принятых в программе

а – матрица А*(F); п – порядок матрицы А*(F); п1 – число циклов вычисления определителя d = êА*(F) ê; d – величина определителя; f – сжимающая сила F; df – шаг изменения сжимающей силы F; m – счетчик цикла вычисления определителя;

- длины стержней балки; п2 – коэффициент п фундаментальных функций; X, Y – векторы значений сжимающих сил F и определителя d.

Текст программы

п = 16; п1 = 300;  f = 0.01; df = 0.01; X = zeros

(n1,1);  Y = zeros

(n1,1); 

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

for  m = 1 : nп2 = sqrt (f);  a = zeros (n,n); 

a(1,2) = sin (п2*l1) / п2;  a(1,4) = - (п2*l1-sin (п2*l1)) / (п2^3);

a(2,2)=cos (п2*l1); a(2,4) = - (1-cos (п2*l1)) / (п2^2); a(2,6) = - 1; 

a(3,2) = п2*sin(п2*l1);  a(3,7) = - 1; a(4,1) = - 1; a(4,4) = 1;

a(5,6) = sin (п2*l2) / п2; a(5,7)= - (1-cos (п2*l2)) / (п2^2);

a(5,8) = - (п2*l2-sin (п2*l2)) / (п2^ 3); a(6,6)=cos (п2*l2); a(6,7) = - a(5,6);

a(6,8) = a(5,7); a(6,10) = - 1; a(7,6) =

п2*sin(п2*l2); a(7,7) =

a(6,6);

a(7,8) = a(5,6); a(7,11) = - 1; a(8,3) = - 1;  a(8,6) = 1; a(9,10) = (sin (п2*l3) / п2;

a(9,11)= - (1-cos (п2*l3)) / (п2^2); a(9,12) = - (п2*l3-sin (п2*l3)) / (п2^ 3);

a(10,10)=cos (п2*l3); a(10,11) = - a(9,10); a(10,12) = a(9,11); a(10,14) = - 1; 

a(11,10) = п2*sin(п2*l3); a(11,11) = a(10,10); a(11,12) = a(9,10); a(11,15) = - 1;

a(12,5) = - 1; a(12,12) = 1; a(13,9) = - 1; a(13,14) = sin (п2*l4) / п2;

a(13,15)= - (1-cos (la*l4)) / (п2^2); a(13,16) = - (п2*l4-sin (п2*l4)) / (п2^ 3); a(14,13) = - 1;  a(14,14)=cos (п2*l4);  a(14,15) = - a(13,14); a(14,16) = a(13,15);

a(15,14) = п2*sin(п2*l4); a(15,15) = a(14,14); a(15,16) = a(13,14);

a(16,16) = 1; d = det(a); X(m,1) = f; Y(m,1) = d; f = f + df; end;




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



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