solve this quastion with this data and formula i send ypu
وزن هواپیما ( W = 4,600 ) پوند
نسبت میانگین نیروی پیشران به وزن ( \left(\frac{T}{W}\right){\text{ave}} = 0.14 )
سرعت باد ( V_w = 0 ) (بدون باد)
شیب باند ( \phi = 0^\circ ) (بدون شیب)
انحراف فلپ ( \delta_f = 15^\circ )
مساحت بال ( S = 175 \text{ft}^2 )
طول بال ( b = 35 \text{ft} )
نسبت Aspect ( A = 7 )
ضریب اصطکاک زمین ( \mu_g = 0.03 )
ضریب بالابری حداکثر ( C{L_{\text{max}}} = 1.69 )
فاکتور اوزولد ( e = 0.80 )
شتاب گرانشی ( g = 32.174 \text{ft/s}^2 )
چگالی هوا در شرایط استاندارد سطح دریا ( \rho = 0.0023769 \text{slugs/ft}^3 )
:نیروی پیشران هنگام برخاست
نیروی پیشران در (V=0)،(TV=0=2,000 پوند)
نیروی پیشران در (65=V_R)،(TV=VR=1,200 پوند)
A 7
agave 8.7883
agV0 12.1445
agVR 6.1160
b 35
CD0 0.0620
CD_muCL 0.0572
CD_out_of_ge 0.1012
CDg 0.0839
CL_max 1.6900
CL_out_of_ge 0.8300
CLg 0.8900
CLOF 1.2779
dCDi -0.0173
dCL 0.1200
e 0.8000
g 32.2000
h 3.6000
h_screen 50
hTR 47.4389
k 0.7236
mu_g 0.0300
phi 0
qLOF 20.5697
qR 18.8199
rho 0.0024
RTR 5.7240e+03
S 175
SCL 19.7690
sigma 0.4800
SNGR 900.9520
SR 125.8400
STO 1.7820e+03
STR 735.4132
T_D_VLOF 591
t_rotate 1
theta_CL 0.1288
tNGR 14.3190
tTCL 0.1471
tTO 21.0560
tTR 5.5899
TV_V0 2000
TV_VR 1200
V50 137.2800
VLOF 131.5600
VLOF_VS_ratio 1.1500
VR 125.8400
VS 114.4000
Vw 0
W 4600
"The ground speed \left(V_G\right), which is the time rate of change of ground distance (s), is illustrated in Figure 9.6. It is defined as follows:
VG=ds dt=V-Vw#(9.1)
 to provide a reasonably accurate representation of jet engine thrust at low speeds:
FNFN0=1-k1V+k2V2#(9.2)
From the definitions of the lift and drag coefficients, given by Equations 7.1 and 7.5, respectively, the lift and drag forces can be expressed as follows:
L=qSCL=12ρSCLV2#(9.3)D=qSCDg=12ρSCDgV2#(9.4)
With the slats and flaps set for takeoff and the landing gear extended, the zero-lift drag coefficient can be determined from the following expression:
CD0=CD0clean +ΔCD0slats +ΔCD0flaps +ΔCD0gear #(9.5)
C_{D_g}=C_{D_0}+\lambda KC_L^2
\mu=\frac{F_{fr}}{R}
The total friction force that acts on the airplane can be expressed as follows:
Ffr=μRMW+μRNW#(9.8)
Ffr=μWcos⁡γG-L#(9.9)
0.0.8. Net Acceleration Force
For the sake of brevity, the net acceleration force acting on the airplane in the direction of motion is herein denoted as \sum F, where
 F=FN-D-Ffr-Wsin⁡γG#(9.10)
or {\ }\sum F=F_N-D-\mu\left(W\cos\funcapply\gamma_G-L\right)-W\sin\funcapply\gamma_G{\ } (from Equation 9.9)
As the runway gradient is a very small quantity, \ ^6 the small angle approximations \cos\funcapply\gamma_G\cong1 and \sin\funcapply\gamma_G\cong\gamma_G (for \gamma_G measured in radians) can be used without loss of accuracy.
Thus \sum F=F_N-D-\mu(W-L)-W\gamma_G
Substituting Equations 9.3 and 9.4 into Equation 9.12 results in the required expression for the net acceleration force:
 F=FN-12ρV2SCDg-μCL-μ+γGW#(9.13)
 The airplane's acceleration (with respect to the ground) can be expressed as follows:
a=ddtds dt=ddtV-Vw=dV dt  if dVw/dt=0#(9.14)
a=\frac{dV}{\mathrm{\ }dt}=\frac{dV}{\mathrm{\ }ds}\frac{ds}{\mathrm{\ }dt}=\frac{dV}{\mathrm{\ }ds}\left(V-V_w\right)
The ground distance measured from the start of the takeoff to the point of rotation \left(s_R\right) can be determined from the following expression:
sR=start end   ds=VwVR  V-Vwa dV#(9.16)
At any instant during the ground run, the airplane's acceleration can be obtained from Newton's second law (see Section 2.4.5):
a= Fm=g FW#(9.17)
s_R=\int_{V_w}^{V_R}\hairsp\hairsp\frac{W\left(V-V_w\right)}{g\sum\ F}\mathrm{\ }dV=\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV}{g\sum\ F}\mathrm{\ }dV-\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV_w}{g\sum\ F}\mathrm{\ }dV
The mean thrust FN is calculated at a reference speed, herein denoted as V^\ast, which is determined as follows:
V*2=12VR-Vw2#(9.19) or  V*=12VR-Vw=0.707VR-Vw#(9.20)
An expression for the net accelerating force ( \sum F ) can be obtained by substituting Equation 9.6 into Equation 9.13, that is,
 F=FN-12ρV2SCD0+λKCL2-μCL-μ+γGW#(9.21)
The required distance \left(s_R\right) is obtained by a summation of the two integrals that are given in Equation 9.18. For the sake of convenience, the two integral expressions will be evaluated separately, where
1  =VwVR  WVg F dV  and  2  =VwVR  WVwg F dV#(9.22)
For the first integral, the following variable substitution is made:
f=V2#(9.23)
thus df=2V\mathrm{\ }dV and dV=\frac{1}{2V}\mathrm{\ }df
hence \int_1\hairsp=\int_{V_w}^{V_R}\hairsp\frac{WV}{g\sum F}\mathrm{\ }dV=\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{WV}{g\sum F}\left(\frac{1}{2V}\right)df=\frac{1}{2g}\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{W}{\sum F}\mathrm{\ }df
Now substitute \sum F from Equation 9.21 into Equation 9.24:
1  =12gfwfR  1FNW-μ-γG-ρS2WCD0+λKCL2-μCLf df#(9.25)
To simplify the analysis, arbitrary constants A_2,B_2, and C_2 are now introduced, where
A_2=FNW-μ-γG  and  B2=ρS2WCD0+λKCL2-μCL
and {\ }C_2=\sqrt{\frac{A_2}{B_2}}
Thus \int_1\hairsp=\frac{1}{2g}\int_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{1}{A_2-B_2f}\mathrm{\ }df=\frac{-1}{2gB_2}\left[\ln\funcapply\left(A_2-B_2f\right)\right]_{\sqrt{f_w}}^{\sqrt{f_R}}
1  =-12gB2ln⁡A2-B2V2VwVR=12gB2ln⁡A2-B2Vw2A2-B2VR2#(9.27)
The second integral is now considered:
\begin{matrix}\ &\ \int_{2}\hairsp\hairsp=\frac{V_w}{g}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{A_2-B_2V^2}\mathrm{\ } dV=\frac{V_W}{gB_2}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{C_2^2-V^2}\mathrm{\ } dV\\\end{matrix}
(\int_{2}\ =\frac{V_w}{2gB_2C_2}\left[\ln{\left(\frac{C_2+V}{C_2-V}\right)}\right]_{V_w}^{V_R}=\frac{V_w}{2gB_2C_2}\ln{\left\{\frac{\left(C_2+V_R\right)\left(C_2-V_w\right)}{\left(C_2-V_R\right)\left(C_2+V_w\right)}\right\}})
The two results are now added to give the distance s_R :
sR=12gB2ln⁡A2-B2Vw2A2-B2VR2-VwC2ln⁡C2+VRC2-VwC2-VRC2+Vw#(9.29)
where A_2,B_2, and C_2 are defined in Equation 9.26 .
uniform acceleration:
sR=1aVwVR  V-VwdV=1aV22-VVwVwVR=12aVR-Vw2#(9.30)
the following integral expression is obtained:
sR=VwVR  V-Vwa dV=WgVwVR  V-Vw F dV#(9.31)
To integrate Equation 9.31 numerically, it is possible to divide the takeoff distance s_R into n segments (or velocity intervals) and to use the trapezoidal rule for incremental values of V :
sR=W2gi=0n  (V)i+1-Vw Fi+1+(V)i-Vw FiΔV#(9.32)
The time that the airplane takes to accelerate to the rotation speed \(\left(V_R\right)\) may be determined using similar approaches to those presented for the distance calculation (Sections 9.4.1 to 9.4.4). The time-velocity relationship is given by
\[
\mathrm{d}t=\frac{1}{a}\mathrm{d}V\emsp\left(\mathrm{\ based\ on\ Equation\ }9.15\right)
Thus\ t=\int_{\mathrm{start\ }}^{\mathrm{end\ }}\mathrm{d}t=\int_{V_w}^{V_R}\frac{1}{a}\mathrm{d}V)
t_R=\frac{1}{\bar{a}}\int_{V_w}^{V_R}\mathrm{d}V=\frac{V_R-V_w}{\bar{a}}
s_{RL}=\Delta t_{RL}\left(\frac{V_R+V_{LOF}}{2}-V_w\right)
s_g=s_R+s_{RL}
U\cong\int\left(F_N-D\right)\mathrm{d} s=\left(F_N-D\right)_{\mathrm{avg\ }}s_{L35}
U=\frac{W}{2g}\left(V_{35}^2-V_{LOF}^2\right)+h_{sc}W
By combining Equations 9.37 and 9.38, the following expression for the distance from liftoff to the 35 ft point is obtained (for zero wind):
s_{L35}=\left(\frac{W}{\left(F_N-D\right)_{avg}}\right)\left(\frac{V_{35}^2-V_{LOF}^2}{2g}+h_{sc}\right)
An alternative-and relatively straightforward-method to estimate the distance from liftoff to the 35 ft point is based on the product of the average time and the average ground speed for this segment of the takeoff:
s_{L35}=\Delta t_{L35}\left(\frac{V_{LOF}+V_{35}}{2}-V_W\right)
s35=V3522a#(9.41)
where s_{35} is the distance from brake release (V=0) to the 35 ft screen height (typical units: m,ft ); and
V_{35} is the TAS at 35 ft (equal to the ground speed as V_w=0 ) (typical units: m/s,ft/s ).
As the greatest part of the takeoff distance will be the ground run, Equation 9.21 can be used to derive a relationship for the mean acceleration. Assuming zero slope (i.e., \gamma_G=0 ), it follows that
a= Fm=gWFN-12ρV2SCD-μRCL-μRW#(9.42)
a=gFNW(1-k)
s_{35}=\frac{V_{35}^2W}{2g(1-k)F_{N_{0.7}}}
where 𝐹𝑁0.7 is the thrust evaluated at 0.7𝑉35 (typical units: N,lb ).
The speed term in Equation 9.44 is now eliminated using the lift coefficient at the screen height 𝐶𝐿35, that is,
$$
𝐶𝐿35=2𝐿𝜎𝜌0𝑉352𝑆=2𝑛𝑧𝑊𝜎𝜌0𝑉352𝑆#(9.45)

تعیین مسافت (s_R) با فرض ثابت بودن (F_N)

این بخش روش ریاضی برای برآورد مسافت (s_R) برای یک هواپیما که در باند خشک یا خیس با شیب (\gamma_G)و باد مخالف (V_w) برخاست می‌کند، ارائه می‌دهد. برای سادگی، فرض می‌شود که کشش موتور \( F_N \) یک مقدار میانگین ثابت است. این مقدار میانگین در سرعت مرجع خاص (V^\ast) تعیین می‌شود، نه در شروع یا پایان حرکت روی زمین، که اطمینان از برآورد دقیق مسافت برخاستن را تضمین می‌کند.

سرعت مرجع (V^\ast):
سرعت مرجع به صورت زیر محاسبه می‌شود:
V^{\ast2}=\frac{1}{2}\left(V_R-V_w\right)^2\emsp\mathrm{یا V*=12VR-Vw=0.707VR-Vw
که در آن (V_R) سرعت چرخش و (V_w) سرعت باد است.

نیروی شتاب خالص (\sum F):
با جایگزینی مقادیر از معادلات قبلی، نیروی شتاب خالص به صورت زیر است:
\sum F=F_N-\frac{1}{2}\rho V^2S\left(C_{D_0}+\lambda K C_L^2-\mu C_L\right)-\left(\mu+\gamma_G\right)W
مسافت برخاستن (s_R) با جمع دو انتگرال محاسبه می‌شود:
\int_{1}\hairsp\hairsp=\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV}{g\sum\ F}\mathrm{\ }dV{\ }\mathrm{\ and\ }{\ }\int_{2}\hairsp\hairsp=\int_{V_w}^{V_R}\hairsp\hairsp\frac{WV_w}{g\sum\ F}\mathrm{\ }dV
جایگزینی انتگرال اول:
f=V2#
1  =12gfwfR  1FNW-μ-γG-ρS2WCD0+λKCL2-μCLf df#
A_2=FNW-μ-γG  and  B2=ρS2WCD0+λKCL2-μCL
{\ }C_2=\sqrt{\frac{A_2}{B_2}}\bigm{\int\ }_1\hairsp=\frac{1}{2g}{\int\ }_{\sqrt{f_w}}^{\sqrt{f_R}}\hairsp\frac{1}{A_2-B_2f}\mathrm{\ }df=\frac{-1}{2gB_2}\left[ln{\left(A_2-B_2f\right)}\right]_{\sqrt{f_w}}^{\sqrt{f_R}}
1  =-12gB2lnA2-B2V2VwVR=12gB2ln⁡A2-B2Vw2A2-B2VR2#
\begin{matrix}\ &\ \int_{2}\hairsp\hairsp=\frac{V_w}{g}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{A_2-B_2V^2}\mathrm{\ } dV=\frac{V_W}{gB_2}\int_{V_w}^{V_R}\hairsp\hairsp\frac{1}{C_2^2-V^2}\mathrm{\ } dV\\\end{matrix}
(\int_{2}\ =\frac{V_w}{2gB_2C_2}\left[\ln{\left(\frac{C_2+V}{C_2-V}\right)}\right]_{V_w}^{V_R}=\frac{V_w}{2gB_2C_2}\ln{\left\{\frac{\left(C_2+V_R\right)\left(C_2-V_w\right)}{\left(C_2-V_R\right)\left(C_2+V_w\right)}\right\}})

#### محاسبه انتگرال دوم:
برای انتگرال دوم:
\int_{2}\ =\frac{V_w}{2gB_2C_2}\ln{\left\{\frac{\left(C_2+V_R\right)\left(C_2-V_w\right)}{\left(C_2-V_R\right)\left(C_2+V_w\right)}\right\}}
که در آن (C_2) به صورت (\sqrt{\frac{A_2}{B_2}}) تعریف شده است.
#### معادله نهایی برای (s_R):
مسافت کل (s_R) مجموع دو انتگرال است:
s_R=\frac{1}{2gB_2}\left[ln{\left(\frac{A_2-B_2V_w^2}{A_2-B_2V_R^2}\right)}-\frac{V_w}{C_2}ln{\left\{\frac{\left(C_2+V_R\right)\left(C_2-V_w\right)}{\left(C_2-V_R\right)\left(C_2+V_w\right)}\right\}}\right]
این روش یک راه حل تحلیلی برای مسافت برخاستن تحت فرض ثابت بودن کشش ارائه می‌دهد.
تعیین مسافت (s_R) با استفاده از شتاب میانگین
 
یک روش ساده‌تر و جایگزین برای برآورد (s_R) بر اساس استفاده از شتاب میانگین \(\ \ \bar{a}\ ) است. معادله برای مسافت روی زمین با شتاب یکنواخت به صورت زیر است:
sR=1aVwVR  V-VwdV=1aV22-VVwVwVR=12aVR-Vw2#
این نشان می‌دهد که باد مخالف مسافت برخاستن را کاهش می‌دهد، در حالی که باد دمی آن را افزایش می‌دهد.
مراحل محاسبه:
1. سرعت مرجع : تعیین (V^\ast) با استفاده از معادله {\ }V^\ast=\frac{1}{\sqrt2}\left(V_R-V_w\right)=0.707\left(V_R-V_w\right).
2. کشش: محاسبه (F_N) در (V^\ast).
3. شتاب: محاسبه (\bar{a}) با استفاده از معادلات \ a=\frac{\sum\ F}{m}=\frac{g\sum\ F}{W}و \sum\ F=F_N-\frac{1}{2}\rho V^2S\left(C_{D_g}-\mu C_L\right)-\left(\mu+\gamma_G\right)W.
4. مسافت برخاستن: استفاده از معادله s_R=1aVwVR  V-VwdV=1aV22-VVwVwVR=12aVR-Vw2 برای یافتن (s_R).
این روش، اگرچه کمتر دقیق از روش‌های عددی است، تخمینی معقول از مسافت برخاستن ارائه می‌دهد و تأثیر باد و شتاب بر عملکرد برخاستن را نشان می‌دهد.
ارزیابی عددی مسافت (s_R)
روش دقیق‌تری برای ارزیابی، ادغام عددی است، به ویژه زمانی که داده‌های واقعی کشش موتور در دسترس باشند. معادله حاکم بر مسافت برخاستن به صورت زیر است:
s_R=\int_{V_w}^{V_R}\frac{V-V_w}{a}dV=\frac{W}{g}\int_{V_w}^{V_R}\frac{V-V_w}{\sum F}dV
که در آن (\sum F) نیروی شتاب خالص است که با معادله 9.13 داده شده است.

مثال از ادغام عددی
مسافت برخاستن (s_R) می‌تواند با تقسیم محدوده سرعت به (\ n\ ) بخش و استفاده از قاعده ذوزنقه‌ای برای ادغام عددی محاسبه شود:
s_R=\frac{W}{2g}\sum_{i=0}^{n}\left[\left(V_{i+1}-V_w\right)F_{i+1}+\left(V_i-V_w\right)F_i\right]\Delta V
برآورد زمان رسیدن به (V_R)
زمان مورد نیاز برای افزایش سرعت هواپیما به سرعت چرخش (V_R) می‌تواند با استفاده از روش‌های مشابه به آنچه در محاسبه مسافت استفاده شده است، تعیین شود. رابطه بین زمان و سرعت به صورت زیر است:
\mathrm{d}t=\frac{1}{a}\thinsp\mathrm{d}V\emsp

بنابراین،
t=\int_{\mathrm{start}}^{\mathrm{end}}\mathrm{d}t=\int_{V_w}^{V_R}\frac{1}{a}\thinsp\mathrm{d}V
توسعه یک راه‌حل تحلیلی برای انتگرال در معادله  پیچیده است. با این حال، یک روش ساده‌تر که فرض می‌کند شتاب میانگین ثابت است، اغلب نتایج قابل قبولی را ارائه می‌دهد. با شتاب ثابت، زمان برای حرکت زمینی (t_R) به صورت زیر است:
t_R=\frac{1}{\bar{a}}\int_{V_w}^{V_R}\mathrm{d}V=\frac{V_R-V_w}{\bar{a}}
برآورد مسافت چرخش
مسافت از چرخش تا برخاستن \left((s_{RL})\right) به دلیل تغییرات در (C_L) به طور دقیق برآورد کردن دشوار است. با این حال، این مسافت نسبت به (s_R) نسبتاً کوچک است. به طور معمول، داده‌های آزمایش پروازی برای تحلیل‌های عملکردی استفاده می‌شود. تخمین مسافت چرخش می‌تواند با ضرب زمان چرخش معمولی در سرعت متوسط زمینی به دست آید:

s_{RL}=\Delta t_{RL}\left(\frac{V_R+V_{LOF}}{2}-V_w\right)
تخمین سرعت برخاستن (V_{LOF})می‌تواند با توجه به اینکه مقررات سازگاری FAR/CS 25.105 نیاز به (V_{LOF}) به میزان حداقل 1.1 برابر سرعت حداقل جداسازی (V_{MU}) می‌شود، به دست آید.

پس از دانستن مسافت چرخش، کل حرکت زمینی \left((s_g)\right) می‌تواند با جمع کردن دو بخش زمینی به دست آید:
s_g=s_R+s_{RL}
روش انرژی برای برآورد مسافت از برخاستن تا 35 فوت
مسافت برخاستن تا پاک کردن صفحه نمایش 35 فوت می‌تواند با در نظر گرفتن کار انجام شده برای تسریع هواپیما از سرعت برخاستن تا سرعت در نقطه 35 فوت و در عین حال بالا بردن آن 35 فوت از سطح باند تخمین زده شود. این تقریب فرض می‌کند که شیب باند و باد صفر است.
کار انجام شده به طور تقریبی به صورت زیر داده می‌شود:
U\approx\int\left(F_N-D\right)\thinsp\mathrm{d}s=\left(F_N-D\right)_{\mathrm{avg}}\thinsp s_{L35}

کار انجام شده می‌تواند به مجموع تغییرات در انرژی جنبشی و پتانسیل معادل شود:
U=\frac{W}{2g}\left(V_{35}^2-V_{LOF}^2\right)+h_{sc}W
که در آن:
- (h_{sc})ارتفاع صفحه نمایش (یعنی 35 فوت یا 10.7 متر).
ترکیب معادلات منجر به عبارت زیر برای مسافت از برخاستن تا نقطه 35 فوت (با فرض باد صفر) می‌شود:
s_{L35}=\left(\frac{W}{\left(F_N-D\right)_{\mathrm{avg}}}\right)\left(\frac{V_{35}^2-V_{LOF}^2}{2g}+h_{sc}\right)
روش مبتنی بر زمان برای برآورد مسافت از برخاستن تا 35 فوت
روش جایگزین و نسبتاً ساده برای برآورد مسافت از برخاستن تا نقطه 35 فوت شامل حاصل‌ضرب زمان میانگین و سرعت زمینی میانگین برای این بخش از برخاستن است:
s_{L35}=\Delta t_{L35}\left(\frac{V_{LOF}+V_{35}}{2}-V_W\right)
که در آن:
- (\Delta t_{L35}) زمان از برخاستن تا پاک کردن صفحه نمایش 35 فوت (واحد: ثانیه).
"
try to use all formula i send you and expand them to solve and find all parametrs can you find