;: TLabel; edKsi: TComboBox ;: TLabel; sgValues: TStringGrid ;: TCheckBox; btnCalculate: TButton; chQ: TChart; 1: TLineSeries; chV: TChart ;: TLineSeries; chP: TChart ;: TLineSeries ;: TLineSeries; edTaul: TComboBox ;: TLabel; edMaxQ: TEdit ;: TLabel; edMinQ: TEdit ;: TLabel; edDeltar: TEdit ;: TLabel; edDeltaguar: TEdit ;: TLabel; btnCalculateClick (Sender: TObject); procedure FormCreate (Sender: TObject); private, A1, sigmal, sigmal wave, Taul, Tau, Ksi, X, Z, Zl, Z2, dZ, dZl, DZ2, T, Q, V, Y, Ro, Xstar, Y star ,, MaxQ, MinQ, DeltaY:Extended;, Index: Integer; procedure Calcualte; procedure Setsigmal; procedure SetAl; procedure NextZ; AddPoints; procedure Clear; public {Public declarations} end ;: TfExtreme; implementation uses Math; const C2=-l;
З 1=2;=- 0.5;=1;=0;=0.005;=0.055;=1.5;=l; wave=l;=l;=0.01;=50;
mega=0.01;=l;=- 0.5;=- 1.5;=- 2.5;=0;=0; gty - 0;=i;=l.l;
Х1Ю.055;
{$ R * .dfm}
{Tffixtreme} TfExtreme.FormCreate (Sender: TObject); begin.RowCount:=Ceil (TMax/dT) +1; sgValues.Cells [0,0]: - i A laquo ;; sgValues.Cells [l, 0]:= t laquo ;; sgValues.Cells [2,0] q raquo ;; sgValues.Celk [3,0]: - V laquo ;; sgValues.Cells [4,0]:= p '; end; Tffixtreme.AddPoints; begin. Series [0] .AddXY (T, Q) ;. Series [0]. AddXY (T, V) ;. Series [0] .AddXY (T, Ro) ;. Series [1] .AddXY (Т, -AO) ;. Cells [0, lndex]:=IntToStr (Index); alues. Cells [1, Index]:=FloatToStr (SimpleRoundT o (T, - 3)) ;. Cells [2, Index]:=FloatToStr (SimpleRoundTo (Q, - 3)) ;. Cells [3, Index]:=FloatToStr (SimpleRoundTo (V, - 3)) ;. Cells [4, Index]:=FloatToStr (SimpleRoundTo (Ro, - 3));
(Index) ;; TfExtreme.Calcualte; var, Ro 1: Extended; irstMaxQF ound: Boolean; begin SetAl;:=0;:=X0;:=C2 * X * X + C1;:=Y;:=0;:=Z; irstMaxQF ound:=F alse;:=Q;:=Q;:= 0;:=V 1;:=l;:=0; while T lt; TMax do begin dQ:=(Z-bO * Q)/b 1;:=Q + dQ * dT;:=Q-MaxQ; (Rol gt; -A0) and (Ro lt; -A0) then begin V :=- V ;;:=Ro ;; chbDrift.Checked then case edExtractable.Itemlndex of 0: Y:=C2 * Math.Power (X-Alfal * Sin (Omega * T), 2) + C 1; 1: Y:=C2 * Math.Power (X-Alfal * Sin (Omega * T), 2) + C 1; end:=C2 * X * X + C1 ;;:=T + dT;:=X + V * dT; if MaxQ lt; Q then MaxQ:=Qif not FirstMaxQFound and (MaxQ gt; Q) then begin:=True;:=Q ;; FirstMaxQFound then if MinQ gt; Q then MinQ:=Q; end; end; TfExtreme.NextZ; var: Extended;:=0;:=0; edExtractable.Itemlndex of 0: Power of 1 ::=(YZ)/Taul ;;
::=Zl;:=(Y - 2 * Ksi * Tau * Zl-Z)/(Tau * Tau) ;;
::=Zl;:=Z2;:=(Y- (Tau * Tau + Taul * 2 * Ksi * Taii) * Z2- (Taiil + 2 * Ksi * Tau) * Zl-Z)/(Tau 1 * Tau * T au); end; end;
: Power of 1 ::=Sin (Y-Z) -Z/Taul; end;
::=Zl;:=(Sin (YZ) - 2 * Ksi * Tau * Zl-Z)/(Tau * Tau) ;;
::=Zl;:=Z2;:=(Sin (YZ) - (Tau * Tau + Taul * 2 * Ksi * Tau) * Z2- (Taul + 2 * Ksi * Tau) * Zl- Z)/(T au 1 * T au * T au); end; end; end;:=Z + dZ * dT;:=Zl + dZl * dT ;:=Z2 + dZ2 * dT; end; Tffixtreme.SetA 1; var, Xstroke2m, DeltaGz, DeltaGy, DeltaOmegaG, AO 1, A02, k, h, T auO, BiExtended; Mu: Extended; begin:=Sqrt (1 -Ksi * Ksi)/Ksi; end; DeltaGy 1: Extended; var;:=l/Taul;:=l/Alfa ;; DeltaGy2: Extended; var: Extended; begin Alfa:=Ksi/Tau; if Ksi=l then Result:=l/Alfa else:=2/Alfa * (
/(Mu * Mu + 1) + Exp (ArcTan (Mu)/Mu)/(Sqrt (Mu * Mu + 1) * (Exp (Pi/Mu) - 1))); end; DeltaOmegaGl: Extended; begin Result:=l; end; DeltaOmegaG2: Extended; begin ifKsi=l then Result:=l else:=(Exp (Pi/Mu) + l)/(Exp (Pi/Mu) - 1); end;:=0;:=0;:=() ;; chbDrifi. Checked then begin:=Alfal * Omega; m:=Alfal * Omega * Omega; case Power of 1 ::=0.001;:=DeltaGy 1; DeltaOmegaG:=DeltaOmegaG 1; end;
::=0.001 ;:=DeltaGy2; DeltaOmegaG:=DeltaOmegaG2 end;
::=0.001;:=DeltaGy 1 + DeltaOmegaG 1 * DeltaGy 2; DeltaOmegaG:=DeltaOmegaGl * DeltaOmegaG2; end ;;
АТ 1:=DeltaGy * (gt + gy * (Xl + Xstrokem));:=DeltaGz * DeltaOmegaG * (gt + gy * (X 1 + Xstrokem)); A0:=2 * (A01 + A02);:=- (Math.Power (Xl-Xstrokem, 2) * gyyplus + gtt + 2 * gty * (Xl + Xstrokem)
+ gy * Xstroke2m);:=gt * gt/(2 * k);:=Sqrt (2 * (AO * (Ysmall + 1) + h)/k);
У:=gt/((X 1 -Xstrokem) * Abs (gyyplus)) + (X 1 + Xstrokem) * T auO; DeltaY:=A01 + 0.5 * B * B * Abs (gyyminus); end else:=2 * (V1 * fl * sigma 1 + V1 * f 1 * sigma2 * sigma 1 wave);:=2 * A0 + fl * V 1 * sigma 1; end; TfExtreme.Setsigmal; function sigmal lwave: Extended; begin Result:=- Tau; end; sigmal2wave: Extended; var, Alfa2, Dlcorr: Extended; begincorr:=Sqrt (Ksi * Ksi - 1);:=(- Ksi-Dlcorr)/Tau;:=(- Ksi + D 1 corr)/T au;:=l/Alfal + 1/Alfa2 ;; sigmal l:Extended; begin Result:=Taul; end; sigmal2: Extended; var, Alfa: Extended; ...
