MATLAB


2.3 - 30


X(4,1)*(cosh(la*x)-

cos(la*x))/(2*la^2))+

m*(sinh(la*(x-am))+sin(la*(x-am))) / (2*la));

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*(sinh(la*(x-am))-

sin(la*(x-am)))*

la / 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))- m*(cosh(la*(x-am)+

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

subplot

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

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

subplot

(2,2,2), plot (x, EIfi); axis ([2 4 - 40 40]); grid on

subplot

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

([2 4 - 2 2]); grid on

subplot

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

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

1-2

= 0 : 0.0001 : 6.0; q = 10.0; la

= sqrt (0.1*0.315365);

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

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

X(8,1)*(sinh(la*x)-sin(la*x)

) / (2*la^3)) + q* (cos(la*x) + cos(la*x) - 2) / (2*la^4));

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

sin(la*x)) / (2* la) - X(8,1)* (cosh(la*x) - cos(la*x)) / (2*la^2) +

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

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

la^4 *(sinh(la*x)-sin(la*x))/(2*la^3) + X(8,1)* (cosh(la*x) + cos(la*x)

) / 2 - q*(sinh(la*x)+sin(la*x))/(2*la);

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

(cosh(la*x)+cos(la*x)) / 2 + X(8,1)*(sinh(la*x)+sin(la*x))/(2*la) -

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

subplot

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

([0 6 - 110 110]); grid on

subplot

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

subplot

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

([0 6 - 40 40]); grid on

subplot

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

([0 6 - 40 40]); grid on

2-3

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

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