Исследование переходных процессов в механической системе
Курсовой проект - Компьютеры, программирование
Другие курсовые по предмету Компьютеры, программирование
0000,0000,0000,0100,0030,0000,0000,0090,0034,0100,00000,0000,0000,0090,0030,0000,0000,0080,0034,0240,00000,0000,0000,0080,0030,0010,0000,0070,0024,0370,00000,0000,0000,0070,0020,0010,0000,0060,0024,0510,00000,0010,0000,0060,0020,0010,0000,0050,0014,0640,00000,0010,0000,0050,0010,0010,0000,0040,0014,0780,00000,0010,0000,0040,0010,0010,0000,0020,0004,0910,00000,0010,0000,0020,0000,0010,0000,0010,0004,1040,00000,0010,0000,0010,0000,0010,000-0,001-0,0014,1180,00000,0010,000-0,001-0,0010,0010,000-0,002-0,0014,1310,00000,0010,000-0,002-0,0010,0010,000-0,003-0,0024,1450,00000,0010,000-0,003-0,0020,0010,000-0,004-0,0024,1580,00000,0010,000-0,004-0,0020,0010,000-0,005-0,0024,1720,00000,0010,000-0,005-0,0020,0010,000-0,006-0,0034,1850,00000,0010,000-0,006-0,0030,0000,000-0,007-0,0034,1990,00000,0000,000-0,007-0,0030,0000,000-0,007-0,0034,2120,00000,0000,000-0,007-0,0030,0000,000-0,007-0,0034,2260,00000,0000,000-0,007-0,0030,0000,000-0,008-0,0034,2390,00000,0000,000-0,008-0,0030,0000,000-0,008-0,0034,2520,00000,0000,000-0,007-0,0030,0000,000-0,007-0,0034,2660,00000,0000,000-0,007-0,0030,0000,000-0,007-0,0024,2790,00000,0000,000-0,007-0,0020,0000,000-0,006-0,0024,2930,00000,0000,000-0,006-0,0020,0000,000-0,006-0,0024,3060,00000,0000,000-0,006-0,0020,0000,000-0,005-0,0014,3200,00000,0000,000-0,005-0,0010,0000,000-0,004-0,0014,3330,00000,0000,000-0,004-0,0010,0000,000-0,003-0,0014,3470,00000,0000,000-0,003-0,001-0,0010,000-0,0020,0004,3600,0000-0,0010,000-0,0020,000-0,0010,000-0,0010,0004,3740,0000-0,0010,000-0,0010,000-0,0010,0000,0000,0004,3870,0000-0,0010,0000,0000,000-0,0010,0000,0010,0014,4000,0000-0,0010,0000,0010,001-0,0010,0000,0020,0014,4140,0000-0,0010,0000,0020,001-0,0010,0000,0030,0014,4270,0000-0,0010,0000,0030,0010,0000,0000,0030,0024,4410,00000,0000,0000,0030,0020,0000,0000,0040,0024,4540,00000,0000,0000,0040,0020,0000,0000,0050,0024,4680,00000,0000,0000,0040,0020,0000,0000,0050,0024,4810,00000,0000,0000,0050,0020,0000,0000,0050,0024,4950,00000,0000,0000,0050,0020,0000,0000,0050,0024,5080,00000,0000,0000,0050,0020,0000,0000,0050,0024,5220,00000,0000,0000,0050,0020,0000,0000,0050,0024,5350,00000,0000,0000,0050,0020,0000,0000,0050,0024,5490,00000,0000,0000,0050,0020,0000,0000,0050,0024,5620,00000,0000,0000,0050,0020,0000,0000,0040,0014,5750,00000,0000,0000,0040,0010,0000,0000,0040,0014,5890,00000,0000,0000,0040,0010,0000,0000,0030,0014,6020,00000,0000,0000,0030,0010,0000,0000,0030,0014,6160,00000,0000,0000,0030,0010,0000,0000,0020,0004,6290,00000,0000,0000,0020,0000,0000,0000,0010,0004,6430,00000,0000,0000,0010,0000,0000,0000,0010,0004,6560,00000,0000,0000,0010,0000,0000,0000,0000,0004,6700,00000,0000,0000,0000,0000,0000,000-0,001-0,0014,6830,00000,0000,000-0,001-0,0010,0000,000-0,001-0,0014,6970,00000,0000,000-0,001-0,0010,0000,000-0,002-0,0014,7100,00000,0000,000-0,002-0,0010,0000,000-0,003-0,0014,7230,00000,0000,000-0,003-0,0010,0000,000-0,003-0,0014,7370,00000,0000,000-0,003-0,0010,0000,000-0,003-0,0014,7500,00000,0000,000-0,003-0,0010,0000,000-0,004-0,0014,7640,00000,0000,000-0,004-0,0010,0000,000-0,004-0,0014,7770,00000,0000,000-0,004-0,0010,0000,000-0,004-0,0014,7910,00000,0000,000-0,004-0,0010,0000,000-0,004-0,0014,8040,00000,0000,000-0,004-0,0010,0000,000-0,004-0,0014,8180,00000,0000,000-0,004-0,0010,0000,000-0,004-0,0014,8310,00000,0000,000-0,004-0,0010,0000,000-0,003-0,0014,8450,00000,0000,000-0,003-0,0010,0000,000-0,003-0,0014,8580,00000,0000,000-0,003-0,0010,0000,000-0,003-0,0014,8710,00000,0000,000-0,003-0,0010,0000,000-0,002-0,0014,8850,00000,0000,000-0,002-0,0010,0000,000-0,0020,0004,8980,00000,0000,000-0,0020,0000,0000,000-0,0010,0004,9120,00000,0000,000-0,0010,0000,0000,000-0,0010,0004,9250,00000,0000,000-0,0010,0000,0000,0000,0000,0004,9390,00000,0000,0000,0000,0000,0000,0000,0000,0004,9520,00000,0000,0000,0000,0000,0000,0000,0010,0014,9660,00000,0000,0000,0010,0010,0000,0000,0010,0014,9790,00000,0000,0000,0010,0010,0000,0000,0020,0014,9930,00000,0000,0000,0020,0010,0000,0000,0020,0015,0060,00000,0000,0000,0020,0010,0000,0000,0020,0015,0200,00000,0000,0000,0020,0010,0000,0000,0020,0015,0330,00000,0000,0000,0020,0010,0000,0000,0030,0015,0460,00000,0000,0000,0030,0010,0000,0000,0030,0015,0600,00000,0000,0000,0030,0010,0000,0000,0030,0015,0730,00000,0000,0000,0030,0010,0000,0000,0030,0015,0870,00000,0000,0000,0030,0010,0000,0000,0030,0015,1000,00000,0000,0000,0030,0010,0000,0000,0020,0015,1140,00000,0000,0000,0020,0010,0000,0000,0020,0015,1270,00000,0000,0000,0020,0010,0000,0000,0020,0015,1410,00000,0000,0000,0020,0010,0000,0000,0020,0005,1540,00000,0000,0000,0020,0000,0000,0000,0010,0005,1680,00000,0000,0000,0010,0000,0000,0000,0010,0005,1810,00000,0000,0000,0010,0000,0000,0000,0010,0005,1940,00000,0000,0000,0010,0000,0000,0000,0000,0005,2080,00000,0000,0000,0000,0000,0000,0000,0000,000
Как видно из таблицы, полное затухание колебаний происходит на пятой секунде. Это время в 3,8 превышает время нашего наблюдения за системой. Предлагаю убедиться в полном затухании колебаний системы с помощью графиков.
коши excel механический система
График изменения координат центра масс первого тела
График изменения координат центра масс второго тела
График изменения скорости первого тела
График изменения скорости второго тела
6. Обработка исходных данных и получение необходимых значений в PASCAL
.1 Файл ввода
- m1
- m2
.05 - d1
.03 - d2
- начальные условия
- начальные условия
- начальные условия
- начальные условия
- количество шагов
6.2 Программа Метод Эйлера
T2thenT0:=3*T1T0:=3*T2;:=T0/n;:=1;[i]:=0;[i]:=y10;y2[i]:=y20;[i]:=V10;V2[i]:=V20;(f2,ty1y2V1V2);[i+1]:=t[i]+ht;[i+1]:=y1[i]+ht*V1[i];[i+1]:=V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i],t[i]);[i+1]:=y2[i]+ht*V2[i];[i+1]:=V2[i]+ht*fv2(y1[i],y2[i],V1[i],V2[i],t[i]);">program Eiler;odnmas = array [1..201] of real;ht,y10,y20,V10,V20,C1,C2,k,m1,m2,w3,g,d1,d2,T1,T2,T0,di:real;,y1,y2,V1,V2:odnmas;i,n:integer;f1,f2:text;fv1(y1,y2,V1,V2,t:real):real;tT2 then T0:=3*T1T0:=3*T2;:=T0/n;:=1;[i]:=0;[i]:=y10;y2[i]:=y20;[i]:=V10;V2[i]:=V20;(f2, t y1 y2 V1 V2);[i+1]:=t[i]+ht;[i+1]:=y1[i]+ht*V1[i];[i+1]:=V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i],t[i]);[i+1]:=y2[i]+ht*V2[i];[i+1]:=V2[i]+ht*fv2(y1[i],y2[i],V1[i],V2[i],t[i]);
{di:=i/10;(frac(di)=0)or(frac(di)=0.1) then }(f2,t[i]:3:3,y1[i]:8:3,y2[i]:8:3,V1[i]:8:3,V2[i]:8:3);:=i+1;i>n+1;:=y1[1];:=y2[1];:=V1[1];:=V2[1];
i:=2 to n doAbs(y1[i])>maxy1 then:=Abs(y1[i]);Abs(y2[i])>maxy2 then:=Abs(y2[i]);Abs(V1[i])>maxV1 then:=Abs(V1[i]);Abs(V2[i])>maxV2 then:=Abs(V2[i]);;(f2);(f2, maxy1 ,maxy2 ,maxV1 ,maxV2 );(f2,maxy1:8:3,maxy2:8:3,maxV1:8:3,maxV2:8:3);(f1);(f2);.
(f1);(f2);.
6.3 Файл вывода результата
ty1 y2 V1 V2
.000 0.000 0.000 0.000 0.000
.007 0.000 0.000 0.000 0.000
.014 0.000 0.000 0.010 0.000
.020 0.000 0.000 0.029 0.004
.027 0.000 0.000 0.058 0.015
.034 0.001 0.000 0.093 0.035
.041 0.001 0.000 0.136 0.062
.047 0.002 0.001 0.184 0.095
.054 0.003 0.001 0.237 0.131
.061 0.005 0.002 0.292 0.166
.068 0.007 0.003 0.348 0.197
.074 0.009 0.005 0.403 0.222
.081 0.012 0.006 0.456 0.241
.088 0.015 0.008 0.503 0.253
.095 0.019 0.010 0.544 0.259
.101 0.022 0.011 0.577 0.260
.108 0.026 0.013 0.600 0.258
.115 0.030 0.015 0.612 0.251
.122 0.034 0.017 0.612 0.240
.128 0.039 0.018 0.600 0.225
.135 0.043 0.020 0.574 0.205
.142 0.046 0.021 0.536 0.179
.149 0.050 0.022 0.485 0.148
.156 0.053 0.023 0.421 0.111
.162 0.056 0.024 0.346 0.069
.169 0.059 0.025 0.260 0.023
.176 0.060 0.025 0.164 -0.026
.183 0.061 0.024 0.062 -0.076
.189 0.062 0.024 -0.047 -0.126
.196 0.062 0.023 -0.160 -0.175
.203 0.060 0.022 -0.275 -0.223
.210 0.059 0.020 -0.389 -0.267
.216 0.056 0.019 -0.501 -0.308
.223 0.053 0.017 -0.608 -0.345
.230 0.048 0.014 -0.708 -0.377
.237 0.044 0.012 -0.800 -0.403
.243 0.038 0.009 -0.881 -0.424
.250 0.032 0.006 -0.949 -0.439
.257 0.026 0.003 -1.004 -0.446
.264 0.019 0.000 -1.044 -0.446
.270 0.012 -0.003 -1.068 -0.439
.277 0.005 -0.006 -1.076 -0.425
.284 -0.002 -0.009 -1.067 -0.404
.291 -0.010 -0.012 -1.042 -0.376
.298 -0.017 -0.014 -1.001 -0.343
.304 -0.024 -0.016 -0.945 -0.304
.311 -0.030 -0.018 -0.874 -0.261
.318 -0.036 -0.020 -0.791 -0.214
.325 -0.041 -0.022 -0.697 -0.165
.331 -0.046 -0.023 -0.593 -0.114
.338 -0.050