MATLAB



3.8 - 2


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

([0 1 - 12 12]); grid on

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

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

1-2

= 0 : 0.001 : 1.0; t = 0.5*2.946045; la2 = sqrt (t); q = 10.0;

EIv = - (X(17,1)*(sinh(la2*x)+sin(la2*x))/(2*la2) - X(8,1)*(cosh (

la2*x)-cos(la2*x))/(2*la2^3)-X(9,1)*(sinh(la2*x)-sin(la2*x)) /

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

EIfi

= - (X(17,1)*(cosh(la2*x)+cos(la2*x)) / 2 - X(8,1)*(sinh(la2*x)+

sin(la2*x))/(2*la2) -X(9,1)*(cosh(la2*x) -cos(la2*x))/(2*la2^2)+

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

Q =-X(17,1)*la2^4*(cosh(la2*x)-cos(la2*x))/(2*la2^2)+X(8.1)*la2^4*

(sinh(la2*x)-sin(la2*x))/(2*la2^3)-X(9,1)*(cosh(la2*x)+cos(la2*x)) / 2

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

M=-X(17,1)*la2^4*(sinh(la2*x)-sin(la2*x))/(2*la2^3)+X(8,1)*(cosh(la2*x)+

cos(la2*x)) / 2+X(9,1)*(sinh(la2*x)+sin(la2*x))/(2*la2)-q*(cosh(la1*x)-

cos(la2*x))/(2*la2^2);

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

([0 1 - 0.4 0.4]); grid on

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

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

([0 1 - 28 28]); grid on

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

([0 1 - 12 12]); grid on

2-4

)

= 0 : 0.001 : 0.5; t = 0.5*2.946045; la3 = (2*t^2)^0.25;

EIv = - (X(16,1)*(cosh(la3*x)+cos(la3*x))/2+X(12,1)*(sinh(la3*x)+

sin(la3*x))/(2*la3) - X(13,1)*(cosh(la3*x)-cos(la3*x))/(2*la3^2) -

X(14,1)*(sinh(la3*x)-sin(la3*x))/(2*la3^3));

EIfi

= - (X(16,1)*la3^4*(sinh(la3*x) -sin(la3*x))/(2*la3^3)+X(12,1)*(

cosh(la3*x)+cos(la3*x))/2-X(13,1)*(sinh(la3*x)+sin(la3*x))/(2*la3)-

X(14,1)*(cosh(la3*x)-cos(la3*x))/(2*la3^3));

Q = - X(16,1)*la3^4*(sinh(la3*x)+sin(la3*x))/(2*la3)-X(12,1)*(