4. Схема алгоритма
Схема головной программы
Из
файла dann.dat
В
файл res.res
В
файл res23.res
Из
файла dann.dat
В
файл res23.res
Схема функции fun
x, f, P, E, J, l: real
5. Таблица идентификаторов
Наименование |
физический смысл |
идентификатор |
Параметр поперечного сечения балки |
a |
a |
Длина балки |
l |
l |
Сосредоточенная нагрузка |
Pk |
P[k] |
Момент инерции сечения |
J |
L |
Модуль упругости |
E |
E |
Число интервалов разбиения |
n |
n |
Прогиб балки |
yi |
y[i] |
6. Текст программы
Program kyr27;
Uses crt;
TYPE vect=array[1..100] of real;
Var P,x,f,y:vect;
i,n,k:integer;
f1,f2:text;
e,a,l,j,dx:real;
function fun(x,f,p,E,J,l:real):real;
begin
fun:=p/(E*J)*exp((3/2)*ln(1+sqr(f)))*(l-x);
end;
begin ClrScr;
assign(f1,'dann.dat');reset(f1);
assign(f2,'res.res');rewrite(f2);
read(f1,a,l,e,n);
writeln(f2,' Kyrsovoi proekt');
writeln(f2,' Issledovanie progiba balki v zavisimosti ot sili');
writeln(f2,' Razrabotal Trashkov A.S.');
writeln(f2,' Isxodnie dannie');
for k:=1 to 4 do read(f1,p[k]);
writeln(f2,'a=',a:3:1,' l=',l:1:0,' P[1]=',P[1]:7:0,' P[2]=',P[2]:7:0,' P[3]=',p[3]:7:0,' P[4]=',P[4]:7:0);
writeln(f2,'E=',E:12:0,' n=',n:2);
writeln(f2,' Opredelenie progiba balki');
dx:=l/n;J:=exp(4*ln(a))/12;
x[1]:=0;y[1]:=0;f[1]:=0;
for k:=1 to 4 do begin
writeln(f2,'Progib balki y(x) pri P=',p[k]:7:0,'H');
writeln(f2,'x y');
writeln(f2,x[1]:5:2,'M ',y[1]:7:5,'M');
for i:=2 to n+1 do begin
x[i]:=(i-1)*dx;
f[i]:=f[i-1]+dx*fun(x[i-1],f[i-1],p[k],E,J,l);
y[i]:=y[i-1]+dx*f[i-1];
writeln(f2,x[i]:5:2,'M ',y[i]:7:5,'M');
end;end;
close(f1);close(f2);
repeat until keypressed
end.
Результаты работы программы
Kyrsovoi proekt
Issledovanie progiba balki v zavisimosti ot sili
Razrabotal Trashkov A.S.
Isxodnie dannie
a=0.8 l=5 P[1]=1500000 P[2]=2000000 P[3]=2500000 P[4]=3000000
E=200000000000 n=20
Opredelenie progiba balki
Progib balki y(x) pri P=1500000H
x y
0.00M 0.00000M
0.25M 0.00000M
0.50M 0.00007M
0.75M 0.00020M
1.00M 0.00040M
1.25M 0.00065M
1.50M 0.00096M
1.75M 0.00132M
2.00M 0.00173M
2.25M 0.00218M
2.50M 0.00268M
2.75M 0.00321M
3.00M 0.00378M
3.25M 0.00437M
3.50M 0.00500M
3.75M 0.00565M
4.00M 0.00632M
4.25M 0.00700M
4.50M 0.00770M
4.75M 0.00841M
5.00M 0.00913M
Progib balki y(x) pri P=2000000H
x y
0.00M 0.00000M
0.25M 0.00000M
0.50M 0.00009M
0.75M 0.00027M
1.00M 0.00053M
1.25M 0.00087M
1.50M 0.00128M
1.75M 0.00176M
2.00M 0.00231M
2.25M 0.00291M
2.50M 0.00357M
2.75M 0.00428M
3.00M 0.00504M
3.25M 0.00583M
3.50M 0.00667M
3.75M 0.00753M
4.00M 0.00842M
4.25M 0.00934M
4.50M 0.01027M
4.75M 0.01122M
5.00M 0.01218M
Progib balki y(x) pri P=2500000H
x y
0.00M 0.00000M
0.25M 0.00000M
0.50M 0.00011M
0.75M 0.00034M
1.00M 0.00066M
1.25M 0.00109M
1.50M 0.00160M
1.75M 0.00220M
2.00M 0.00288M
2.25M 0.00364M
2.50M 0.00446M
2.75M 0.00535M
3.00M 0.00629M
3.25M 0.00729M
3.50M 0.00833M
3.75M 0.00941M
4.00M 0.01053M
4.25M 0.01167M
4.50M 0.01284M
4.75M 0.01402M
5.00M 0.01522M
Progib balki y(x) pri P=3000000H
x y
0.00M 0.00000M
0.25M 0.00000M
0.50M 0.00014M
0.75M 0.00041M
1.00M 0.00080M
1.25M 0.00130M
1.50M 0.00192M
1.75M 0.00264M
2.00M 0.00346M
2.25M 0.00437M
2.50M 0.00536M
2.75M 0.00642M
3.00M 0.00755M
3.25M 0.00875M
3.50M 0.01000M
3.75M 0.01130M
4.00M 0.01263M
4.25M 0.01401M
4.50M 0.01541M
4.75M 0.01683M
5.00M 0.01826M