Текст программы

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

ГРАФИЧЕСКАЯ ЧАСТЬ