MATLAB



3.8 - 4


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

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

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

([0.5 1 - 70 70]); grid on

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

([0.5 1 - 30 30]); grid on

 

1-3

)

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

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

sin(la4*x))/(2*la4) - X(18,1)*(cosh(la4*x)-cos(la4*x))/(2*la4^2) -

X(19,1)*(sinh(la4*x)-sin(la4*x))/(2*la4^3));

EIfi

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

cosh(la4*x)+cos(la4*x)/2-X(18,1)*(sinh(la4*x)+sin(la4*x))/(2*la4)-

X(19,1)*(cosh(la4*x)-cos(la4*x))/(2*la4^2));

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

cosh(la4*x)-cos(la4*x))/(2*la4^2)+X(18,1)*la4^4*(sin(la4*x) -

sinh(la4*x))/(2*la4^3)+X(19,1)*(cosh(la4*x)+cos(la4*x))/2;

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

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

X(19,1)*(sinh(la4*x)+sin(la4*x))/(2*la4);

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

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

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

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

([0 0.25 - 7 7]); grid on

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

([0 0.25 - 18 18]); grid on

)

= 0.25 : 0.001 : 1.0; m = 20.0; t = 0.5*2.946045; la4 = (4*t^2)^0.25;

EIv =-(X(16,1)*(cosh(la4*x)+cos(la4*x))/2+X(17,1)*(sinh(la4*x)+sin(la4*x))

/(2*la4) - X(18,1)*(cosh(la4*x)-cos(la4*x))/(2*la4^2)-X(19,1)*(sinh(la4*x)

-sin(la4*x))/(2*la4^3)+m*(cosh(la4*(x-0.25))-cosh(la4*(x-0.25)))/(2*la4^2));

EIfi

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