4297

CONST g = 9.81; var

Form1: TForm1; {formclik : Boolean;}

a1,b1,l,l1,gamma1,delta,diam,s,mc,m,jc,q,ay,dtime,fi n,fik: extended; sin_beta1,fit,fis,fip,pt,pp,ppr,kr,kkr,kt,tn,time,sqrt1: extended; time_ost,time_shot,sn,sk,s_step,proc_time: extended;

opt_a1,opt_b1,beta_n,opt_gamma1,aa1,aa2: Extended; opt_gamma1_1,opt_a1_1,opt_b1_1,fin1,fik1:Extende d;

Flag2,Flag3,Flag4,Flag5:Boolean;

ii:integer;

aaa1,aaa2:TChartValueList;

bbb1,bbb2:Double; {a,b,c,d: TChartTitle;} implementation

{$R *.DFM}

function ScanStr(S:TEdit):Real; Var i:Integer;

Flag_ok, Flag_ok1: Boolean; begin

Flag_ok:=False;

Flag_ok1:=True;

for i:=1 to Length(S.Text) do Begin

if S.Text[i] in ['0'..'9','E','e','-',','] then Flag_ok:=True else Flag_ok:=False;

if Flag_ok = False then Flag_ok1:=False; End;

if S.Text = '' then Flag_ok1:=False; {if S.Text = '0' then Flag_ok1:=False;} if Flag_ok1 = False then

Begin

MessageDlg('Допустимые значения: цифры, запятая, минус, E (лат.) - показатель степени',

mtInformation,[mbOk],0); S.Text:=''; ScanStr:=0; S.SetFocus;

End;

if Flag_ok1 = True then Begin ScanStr:=StrToFloat(S.Text); End; end;

function myarcsin(x:Extended):extended; Begin if x=1 then myarcsin:=pi/2

else myarcsin:=arctan(x/sqrt(1-x*x)); End; function myarccos(x:Extended):Extended; Begin myarccos:=pi/2-myarcsin(x); End; procedure

TForm1.SeriesChange(number_series:integer;chart_c aption:string);

Var i:integer; begin

for i:=0 to 3 do Chart1.Series[i].Active:=False; Chart1.Series[number_series-1].Active:=True; Chart1.Title.Text.Text:=chart_caption;

end; //----------------------------------

{Начало математических расчетов}; pp:=0; pt:=0; kr:=0; fis:=0; time:=0; sqrt1:=0; sin_beta1:=0; kt:=1; time_shot:=0; s:=pi*diam*diam/4; gamma1:=gamma1*pi/180; delta:=delta*pi/180; fin1:=fin*pi/180; fik1:=fik*pi/180; fit:=fin1; fip:=fin1; Series1.Clear; Series1.Active:=False; Series2.Clear; Series2.Active:=False; Series3.Clear; Series3.Active:=False; Series4.Clear; Series4.Active:=False;

while (fit<=fik1) or (time<time_ost+time_shot) do Begin

sin_beta1:=a1*sin(gamma1+fit)/sqrt(a1*a1+b1*b1- 2*a1*b1*cos(gamma1+fit)); kr:=kkr/(7.28*pp+106);

if time<tn then kt:=time/tn else kt:=1; if Flag5=True then

Begin pt:=(q*kt*dtime-s*b1*sin_beta1*(fit-fip)- ay*dtime*pp)/kr+pp; fis:=(s*pt*b1*sin_beta1-

g*(m*l+mc*l1)*cos(fit+delta))*dtime*dtime/(jc+m*l

*l)

+2*fit-fip; End;

if Flag5=False then Begin

pt:=(g*(m*l+mc*l1)*cos(fit+delta))/(s*b1*sin_beta1 );

fis:=((q*kt-ay*pt)*dtime)/(s*b1*sin_beta1)+fit; End; fip:=fit; fit:=fis; pp:=pt; time:=time+dtime; Series1.AddXY(time,pt,'',ClRed); Series2.AddXY(time,fit*180/pi,'',ClRed); Series3.AddXY(time,(b1*sin_beta1*(fitfip)/dtime),'',ClRed); Series4.AddXY(time,kr*1000000,'',ClRed);

if (fit>=fik1) and (time_shot=0) and (time_ost<>0) then

Begin q:=0; time_shot:=time; fit:=fip; End; End; Chart1.Visible:=True; Chart1.Title.Text.Text:='Зависимость давления (Па) от времени';

Series1.Active:=True;

End;

procedure TForm1.Button2Click(Sender: TObject); {оптимизация положения}

LABEL lab1,lab2,lab3; VAR Flag1:Boolean;

p_m,p_max:Extended; key:Char;

Begin

If Radiobutton2.Checked=True then Flag5:=True; if Radiobutton1.Checked=True then Flag5:=False; Radiobutton1.Visible:=False; Radiobutton2.Visible:=False;

flag1:=False; p_m:=0; p_max:=10e16; flag3:=False; if flag2=false then Begin opt_b1:=s_step; End; if flag2=true then goto lab1; while opt_b1<l do

Begin lab2:

if flag2=true then if flag4=false then if opt_b1>l then opt_b1:=s_step;

{if flag2=true then if flag4=true then if opt_b1<0 then opt_b1:=l;}

{расчет выражения для arccos для бета начальное}

aa1:=4*sn*opt_b1*sin((fik1-fin1)/2); if aa1=0 then goto lab1;

aa1:=(sn*sn-sk*sk+4*opt_b1*opt_b1*(sin((fik1- fin1)/2))*

(sin((fik1-fin1)/2)))/aa1;

if abs(aa1)>1 then goto lab1; {расчет бета начальное}

beta_n:=myarccos(aa1)-0.5*(pi-fik1+fin1); {расчет подкоренного выражения для a1} sqrt1:=sn*sn+opt_b1*opt_b1- 2*sn*opt_b1*cos(beta_n);

if sqrt1<0 then goto lab1; {расчет оптимизационного a1} opt_a1:=sqrt(sqrt1);

{расчет выражения для arccos gamma1} aa2:=(opt_a1*opt_a1+opt_b1*opt_b1sn*sn)/(2*opt_a1*opt_b1);

if abs(aa2)>1 then goto lab1; opt_gamma1:=myarccos(aa2)-fin1;

{Проверка на отсутсвие мертвой точки. Параметр проверки}

if (opt_gamma1+fik1)>=pi then goto lab1; {Конец проверки}

if opt_a1*opt_a1+opt_b1*opt_b1-2*opt_a1*opt_b1* cos(opt_gamma1+fik1)<>sk*sk

then goto lab1; {инициализация переменных} if Flag2=True then Begin

Series1.Clear; Series2.Clear; Series3.Clear; Series4.Clear; End;

pp:=0; pt:=0; kr:=0; fit:=0; fis:=0; time:=0;

p_m:=0; sqrt1:=0; sin_beta1:=0; time_shot:=0; fit:=fin1; fip:=fin1; flag1:=True;

if flag2=true then flag3:=true; {отладка программы}

{MessageDlg(FloatToStr(opt_gamma1*180/pi)+ ' '+ FloatToStr(opt_a1)+ ' '+

FloatToStr(opt_b1)+ ' ', mtInformation,[mbOk],0); {отладка}

{цикл расчета от начального до конечного положения}

while (fit<=fik1) or (time<time_ost+time_shot) do Begin sin_beta1:=sqrt(opt_a1*opt_a1+opt_b1*opt_b1- 2*opt_a1*opt_b1*cos(opt_gamma1+fit)); sin_beta1:=opt_a1*sin(opt_gamma1+fit)/sin_beta1; kr:=(7.28*pp+106);

if kr=0 then goto lab1; kr:=kkr/kr;

if time<tn then kt:=time/tn else kt:=1; if Flag5=True then

Begin pt:=(q*kt*dtime-s*opt_b1*sin_beta1*(fit-fip)- ay*dtime*pp)/kr+pp;

fis:=((s*pt*opt_b1*sin_beta1- g*(m*l+mc*l1)*cos(fit+delta))* (dtime*dtime)/(jc+m*l*l))+2*fit-fip; End; if Flag5=False then

Begin pt:=(g*(m*l+mc*l1)*cos(fit+delta))/(s*b1*sin_beta1 );

fis:=((q*kt-ay*pt)*dtime)/(s*b1*sin_beta1)+fit; End; {Макс. давление за процесс подъема} if pt>=p_m then p_m:=pt;

fip:=fit; fit:=fis; pp:=pt; time:=time+dtime; if Flag2=True then Begin Series1.AddXY(time,pt,'',ClRed); Series2.AddXY(time,fit*180/pi,'',ClRed); Series3.AddXY(time,(b1*sin_beta1*(fitfip)/dtime),'',ClRed); Series4.AddXY(time,kr*1000000,'',ClRed); End;

if (fit>=fik1) and (time_shot=0) and (time_ost<>0) then

Begin q:=0; time_shot:=time; fit:=fip; End; End; if Flag2=True then goto lab1;

{мин. давление за все процессы} if p_m<=p_max then

Begin p_max:=p_m; opt_gamma1_1:=opt_gamma1; opt_a1_1:=opt_a1; opt_b1_1:=opt_b1;

End; {} lab1:

if (flag4=false) and (Flag3=false) then opt_b1:=opt_b1+s_step;

{if (flag4=true) and (Flag3=false) then opt_b1:=opt_b1-s_step;}

if Flag3=True then goto lab3;

if Flag2=true then goto lab2; End; if flag1=False then Begin

MessageDlg('При введенных данных оптимизация невозможна',

mtInformation,[mbOk],0); FormCreate(Form1);goto lab3; End;

Flag2:=True; Flag3:=False; opt_b1:=opt_b1_1; goto lab2; lab3: key:='1'; Flag3:=False; Chart1.Visible:=true; FormKeyPress(Form1,key); End;

procedure TForm1.FormCreate(Sender: TObject); var key: char;

begin with Form1 do Begin caption:='Программа расчета гидроцилндров'; key:='r'; flag2:=False; flag4:=false; Button1.Enabled:=True;

Button1.Visible:=true;

FormKeyPress(Form1,key); Edit22.Enabled:=False; Label22.Enabled:=false; Chart1.Visible:=False;

If Radiobutton2.Checked=True then Flag5:=True; if Radiobutton1.Checked=True then Flag5:=False; End; end;

procedure TForm1.FormKeyPress(Sender: TObject; var Key: Char);

var f: textfile;

f1: file of extended;