MATLAB



3.8


1

= (1,1) = - 27.6953

2

= (2,1) = 1.2487 2

3

= (3,1) = - 68.9032

4

= (4,1) = 10.1048

5

= (5,1) = 15.2920

6

= (6,1) = - 12.7248

7

= (7,1) = - 12.8813

8

= (8,1) = - 11.6889

9

= (9,1) = 25.1927

10

= (10,1) = 0.5814

11

= (11,1) = - 25.7170

12

= (12,1) = -

1.3367 2

13

= (13,1) = 8.6147

14

= (14,1) = - 0.5814

15

= (15,1) = 15.2920

16

= (16,1) = 3.7219 3

17

= (17,1) = -

2.0502 2

18

= (18,1) = 16.4936

19

= (19,1) = 0.5814

20

= (20,1) = - 25.7170

 

 

 

 

3.3.6. -

EIJ(x), EIj(x), Q(x), (x) (3.21), (3.22). EIJ(x), EIj(x), Q(x) (x)

0-1

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

EIv = - (X(2,1)*(sinh(la1*x)+sin(la1*x))/(2*la1) - X(4,1)*(sinh

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

EIfi

= - (X(2,1)*(cosh(la1*x)+cos(la1*x)) / 2 - X(4,1)* (cosh(la1*x) -

cos(la1*x))/(2*la1^2)+q*(sinh(la1*x)-sin(la1*x))/(2*la1^3));

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

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

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

X(4,1)*(sinh(la1*x)+sin(la1*x))/(2*la1)-q*(cosh(la1*x)-cos(la1*x))/(2*la1^2);

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

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

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