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


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


Уравнение для граничных параметров балки отличается от уравнения (3.11) только фундаментальными функциями и вектором правой части        (см. (3.16)).

Программа вычисления граничных параметров набирается в отдельном М-файле и принимает вид

a = zeros (16,16);  b =

zeros (16,1); X = zeros

(16,1); la = sqrt (0,4055125);

l1 = 4.0;  l2 =

6.0;  l3 = 3.0;  l4 = 1.0;

a(1,2) = (sinh(la*l1) + sin (la*l1))/(2*la);

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

a(2,2)=(cosh(la*l1)+cos(la*l1))/2; a(2,4) = - (cosh(la*l1) - cos (la*l1))/(2*la^2);

a(2,6) = - 1;  a(3,2) = la^ 4*a(1,4); a(3,4) = a(1,2); a(3,7) = - 1; a(4,1) = - 1;

a(4,2) = la^ 4*a(2,4); a(4,4) = a(2,2); a(5,6) = (sinh(la*l2) + sin (la*l2))/(2*la);

a(5,7)= - (cosh(la*l2) - cos(la*l2))/(2*la^2);

a(5,8) = - (sinh(la*l2) - sin (la*l2))/(2*la^ 3);

a(6,6)=(cosh(la*l2)+cos(la*l2))/2;

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

a(6,10) = -1; a(7,6) =la^ 4*a(5,8); a(7,7) =a(6,6);

a(7,8) = a(5,6); a(7,11) = - 1; a(8,3) = - 1; a(8,6) = la^ 4*a(5,7); a(8,7) = - a(7,6);

a(8,8) = a(6,6); a(9,10) = (sinh(la*l3) + sin (la*l3))/(2*la);

a(9,11)= - (cosh(la*l3) -cos(la*l3))/(2*la^2);

a(9,12) =-(sinh(la*l3)-sin(la*l3))/(2*la^3);

a(10,10)=(cosh(la*l3)+cos(la*l3))/2;

a(10,11) = - a(9,10); a(10,12) = a(9,11); a(10,14) = - 1;  a(11,10) = la^ 4*a(9,12);

a(11,11) =a(10,10); a(11,12)= a(9,10); a(11,15)= -1; a(12,5)= -1; a(12,11) = - a(11,10);

a(12,10)=la^ 4*a(9,11); a(12,12)=a(10,10);

a(13,14)=(sinh(la*l4) + sin (la*l4))/(2*la);

a(13,9) = - 1; a(13,15)= - (cos h(la*l4) -cos(la*l4))/(2*la^2);

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

a(14,14)=(cosh(la*l4)+cos(la*l4))/2; a(14,15) = - a(13,14); a(14,16) =

a(13,15); 




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