MATLAB



3.8 - 3


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

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

M = - X(16,1)*la3^4*(cosh(la3*x)-cos(la3*x))/(2*la3^2)-X(12,1)*la3^4*

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

X(14,1)*(sinh(la3*x)+sin(la3*x))/(2*la3);

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

([0 0.5 - 4 4]); grid on

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

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

([0 0.5 - 10 10]); grid on

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

([0 0.5 - 10 10]); grid on

)

= 0.5 : 0.001 : 1.0; t = 0.5*2.946045; la3 = (2*t^2)^0.25; f = 60.0;

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)+f*(sinh(la3*(x-0.5)) -

sinh(la3*(x-0.5)))/(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^2)+

f *(cosh(la3*(x-0.5))-cos(la3*(x-0.5)))/(2*la3^2));

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

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

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

f *(cosh(la3*(x-0.5))+cos(la3*(x-0.5)))/2;

M = - X(16,1)*la3^4*(cosh(la3*x)-cos(la3*x))/(2*la3^2)-X(12,1)*la3^4* (

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

X(14,1)*(sinh(la3*x)+sin(la3*x))/(2*la3) - f *(sinh(la3*(x-0.5))+

sinh(la3*(x-0.5)))/(2*la3);

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