MATLAB в инженерных и научных расчетах


Продолжение таблицы 2.3 - часть 25


a(15,14) = la^ 4*a(13,16); a(15,15) = a(14,14);

a(15,16) = a(13,14); a(16,14) = la^ 4*a(13,15);

a(16,15) = - a(15,14); a(16,16) = a(14,14); b(16,1) = 1; X = a\ b;

X = X \ X(9,1)

Результаты вычислений относительных значений граничных параметров балки при собственных частотах сведены в таблицу 3.2.

Таблица 3.2

№ п/п

Граничные

параметры балки

Относительные значения граничных параметров при частотах

w1=0,4055125

w2=0,7818475

w3=1,1001125

w4=1,4852935

w5=2,4234325

1

0.1677

- 3.5779

0.3375

- 0.4569

- 0.9724

2

- 1.4419

3.5754

- 0.2308

0.3784

- 0.4260

3

- 1.6643

- 1.3379

- 0.0981

1.8113

0.9831

4

- 0.6979

2.6728

- 0.2409

0.5453

- 1.0318

5

0.6626

0.6263

0.5250

0.2441

- 1.6235

6

2.2068

- 1.8784

- 0.0844

0.4282

- 0.4491

7

- 1.1280

- 2.3814

0.4101

- 0.8967

0.0745

8

1.5705

0.6324

- 0.5215

1.7308

- 1.2045

9

1.0000

1.0000

1.0000

1.0000

1.0000

10

- 1.7388

- 1.1200

- 0.4732

0.2324

0.5331

11

- 1.5543

- 0.5464

0.4042

1.2018

- 0.1983

12

0.3185

- 0.4130

- 0.9820

- 1.1665

1.5478

13

1.0055

1.0203

1.0400

1.0723

1.1873

14

0.9850

0.9443

0.8905

0.8024

0.4932

15

- 0.0547

- 0.2021

- 0.3971

- 0.7145

- 1.8141

16

0.0819

0.3020

0.5910

1.0566

2.6174

 

Выражение для прогиба (формы колебания) стержня при собственных колебаниях в соответствии с методом начальных параметров имеет вид

                    (3.17)

Протокол построения форм колебаний балки принимает вид

х1 = 0 : 0.001 : 4.0;  х2 = 0 : 0.001 : 6.0;  х3 = 0 : 0.001 : 3.0; 

х4 = 0 : 0.001 : 1.0;  la = sqrt (0.4055125);

EIv1 = - (X(2,1)*(sinh(la*x1)+sin(la*x1))/(2*la) - …

X(4,1)*(sinh(la*x1)-sin(la*x1))/(2*la^3));




Начало  Назад  Вперед