Bobby Nordlund | > iak Taye Newman 00% =r ae Captain O'Reilly ¢ eee ‘ + ob 21 Apri 2010 aie Rocket Competitive Success Criteria (CSC) Analysis pe Costly components Lo ae. ‘component Rate Quantity, Cost Engine Adaptor System | $5000 2.00 $10,000.00 Body Tube | $20000/in 10.50" {$105,000.00 Balsa Wood '$100/in’ 10.03in* $1,003.00 B6-0 & A8-3 $1,000,000 1.00 ‘Actual payload: 32.48 Cost-to-payload ratio: $34,444.54/¢ ‘Two tools were used to predict the rockets performance. Firstly, the commercial RockSim software (by ‘Apogee) was ands a highly-accurate rocketry modeling program that does just about anything rocketry. It was used to model each of the four configurations in order to determine which design held a KanKowed competitive advantage in the competition. Ukimately, our modeling predicted the B-A design would have acost-to-payload ratio more than 12% less than any other design. rahe MATLAB proved a useful, albeit cumbersome, tool for predicting our rocket’s performance. Unfortunately, it was hamstrung by our measured thrust data which lacked the accuracy of Estes’ Boy es published thrust profiles. The -12% error in Impulse seen in the Thrust Lab was carried forward and reflected clmost exactly in the discrepancy between RockSim’s predictions and MATLAB's predictions of our rocket’s performance: the difference between 305' and 260'.~ Th's 4 O6™ Coinc: derdal foes easly The margin we built into our design was swallowed by Skt winds out gf the SE. Our over-stable rocket experienced significant weather-cocking that tured precious vertiatglstance inte horizontal distance, resulting ina fight just beneath the altitude predicted by RockSim: 292 Xe ellens MATLAB Inputs 8M pays Temperature: 293 K 9 Pressure: 101325 Pa Payload-less initial mass: 0.1151 kg. Booster mass: 0.0810 kg Payload: .0324 kg Drag coefficient: 0.80 Cross-sectional area: 7.45*10* m? ‘5 “_g ge PreaghMATLAB Outputs © Takeoff mass: 147.5 g © Atend of boost. (© Velocity: 21.7603 m/s © Height: 9.6785 m ‘+ Atend of upper stage... (© Velocity: 36.0384 m/s © Height: 24.8395 m © Atend of coast. (© Velocity: 0 m/s (stopped) © Height: 79.2591 m (260.0363 ft) Flight Profle dct Tine Toe oct ii) 2 wane = Tvs TneBobby Nordlund Taylor Newman Captain O'Reilly A201-5001 21 April 2010 Rocket Project: MATLAB CodeA204 Rocket Project EA204 Rocket Project Contents i Nordlund and Newman 21 March 2010 This script seeks to model the performance of a two-stage rocket fitted with an Estes BO and A8 serial-stage propulsion system. The values of input parameters are freely adjusted, although the thrust data is the direct result of experimentation and is less-readily modified. Thrust, mass, weight, and drag (all functions of time) are accounted for in the script. Assumptions are limited to the use of published burn- time values. format short, clear all; cle Enter Parameters temperature=293; $Temperature (K) pressure=101325; SPressure (Pa) m0=.1151; $Payload-iess initial mass (kg) mb=.0410; Booster mass (kg) payload=.0324; sPayload mass. (kg) Cod-.8; iDrag’ coefficient (approximate) As. i5e-4; SEffective cross-sectional area (m*2) Preliminary Calculations Baire2672, 83/ (arK) Rye rho= (pressure) /(R_air*temperature) ;_tkg/n” mem Ospayload; Yactual initial mass S-pisplay Kelavent Pareneters disp(cat (2, "Payload mass: *,num2str (payload*1000) ," 9") disp (cat (2, "Booster mass: ';mumastr(m 8"1000) , disp (cat (2; ‘Takeoff mass: ',num2ste (mat i ay) a) Booster Stage Calculations ‘Assume three forces: Thrust (T_B6), Weight (W), Drag (D) Note: Worthless data in 1001b6-4a, 2001a8-3a, 2001a8-3b and data from 5001b6-4a (C465:C545) appears to be from’a different size engine. L)sxlsread(‘thrustdata!, '2001b6-4a" , 'C225:C305'); 2)=xisread(*thrustdata* , '3001b6-6a", 'C685:C765"); 3)sxlsread(*thrustdata’, '3001b6-6b" , '315:C395"); y ) 4)=xIsread(*thrustdata' , '5001b6-6a", 'C553:C633")7 S)=xlsread(*thrustdata', *6001b6-4a", 'C605:C6a5" 2,6) =axlsread(*thrustdata’ , '6001b6~4b", 'C463:C543"); ‘erence: [Link] commonwealth .net /récket store/estes/[Link] Ns impulse, -8 8 burn time, 6.24 g propellant [Link] [Link]A204 Rocket Project I_Béstrapz(T_B6,1)/100; Ns 17Bé=mean (I B6, 2); SAverage total impulse T3p_B6=1 B67 (.00624*9.81); SAverage specific impulse n_dot=7 B6/(Isp_B6*9.61); tMass Flow Rate (Sellers 14-7) t=.01; Ttime interval © Mass Lost (trapezoidal method) m_lostat®(~.5* (m dot (1)4m_dot (Length (m_dot) )) +cumsum (m_dot (1: length (m_dot))) molost(1)=0; tinitial mass lost. due to thrust miiost (2) tt (m dot (1) +m dot (2))/27 Strapezoidal integration of m dot from 0 to t msm_O-m. lost, Wen.*9.81; Weight v boost (1)=0; initial velocity “boost (1)=(T BG (1)-W(1))./m(1); Yacceleration at t=0 (D=0) cost (2)=(TB6 (2) -W(2)) :/m(2); Aacceleration at t-0 (assume D=0) foe ke2: length (T_B6) =1 ¥ boost (k) =v boost (k-1) +t poin D(k)=.S*rhoty boost (k)~2*C d*a; tbrag boost (ik) = (TB6 (k) “W(k) -Dik)) «/m(k)} Sredefine acceleration with Drag acboost (+1) =(T_B6 (k+l) ~W (k#1)-D(k))./m(k+1) ; Stemporary D assumption 1 boost (k-1) +a_boost (k))/2; Strapezoidal integration point by end v_bew boost (length(v_boost)); Svelocity at burnout hoboost=t* (—. 5+ (v boost (1) +v boost (length (v boost) ))+cunsum(v_boost (1:Jength(v_boost))))7 hoboost (1)=0; singeial height hoboost (2) =t* (v boost (1) + boost (2))/2; Strapezoidal integration of v from 0 to t hob=h boost (length (h_boost})? disp(*At end of boost disp(cat(2,' Velocity is ‘,num2ste(v b)," m/s')) disp(cat(2;" Height is *,ndmauee (h by," m')) At end of boost. Velocity 18 2127603 m/s Height is 9.6785 m Upper Stage Calculations Assume three forces: Thrust (T_B6), Weight (W), Drag (D) Reference: Estes A8-3 (2.5 Ns impulse, .5 s burn time, 3.12 g llant) Note: Worthless data in 1001b6- 4a, 2001a8-3a, 2001a8-3b and data from 5001b6-4a (C465:C545) appears to be from a different size engine. sA)mxisread(' thrustdata', '1001a8-3a!, '¢305, '2)axisread(*thrustdata', *2001b6-4b", "C121 (3)=xisread *3001a8-3a", "C260 ‘4)sxisread *5001a8-3a", 1C550 {S)exisread *5001a8-3b' , "C470. 16)=xlsread "e001a8-3a", "C510. FTAG(2) 7) aularead(*thrustdata’y "6001a6-3b"» "C475:C525") + ‘Refetence: [Link] commonwealth: not/rocketstore/eates/[Link] § Estes Ae-3 (2.5 No inpuise, «5 8 burn ime, 3.12 q propellant) EAB=trape (78,1) /1007"aNs TAB-neat (7 8,2); daverage total impulse Tep_AB=1 G7 (.50312"9,81)7 SAverage specific impulse iStet NO/(1sp_A8*9,01); wave Flow Rate (Sellers 14-7) 01; Stine interval Viigo Lose (trapezoidal myLiod) m_Lostet*(~.5* (im dot (1) rm dot (length (m dot) )) +cunsum(m_dot (1:ength (m_dot))))7 milost(1)-0; tinfeial mass lost. due to thrust matost (2)=t* in dot (1)+m_dot (2))/27 Strapezoidal integration of mdot from 0 to t men O-m b-m Losey fens? 5-81) SWeight Upper (1)=v'b; Sinitial velocity Du}=.S*rhoty upper (1) "2c dea; Atnitial drag a_upper (1)= (FAB) -W(2) -D12))'./m(1)7 Qacceleration, D-D(v_b) ‘acupper (2)= (T7A8 (2) -W (2) -D(1}) ./mi2); tacceleration, assume D is still D(v_b) Tor? he siengen (rae) -1 v_upper (k) =v upper (k-1) +t* (a_upper(k-1) +a_upper (k)}/27 trapezoidal integration point by point D(k)=.S*zhoty_upper (k)*2*C_d*A; $rag upper (k) = (TAB (K) -W(k) -Dk)) -/m(k)? Sredefine acceleration with Drag acupper (ke1)=(T_AS (kel) -W(k+1) -D(x)) «/miks1) 7, temporary D assumption v_uev_upper (length (v_upper)); tvelocity at burnout sppSe=h_b+t* (=. 5* (¥_upper (1) +v_upper (length (v_upper) )) tcumsum(v_upper (1: length (v_upper))))+ [Link] [Link]A204 Rocket Project b_upper (1)=h by initial height houpper (2) =h-b+t* (v_upper (1) #v_upper (2)) /2, housh upper (Tength (f_upper) ) disp(TAt end of upper stage, ..') disp(cat (2,' velocity is ‘;numastr(v_u),* m/s") disp(cat(2;' Height 1s ',namaste(h_uJ,'’m')) Strapezoidal integration of v from 0 to t At end of upper stage. Velocity 15.36.0384 m/s Reight is 24,8395 m Coast Stage Calculations mem(length(m)); tmass at burnout Wom+9.81; aWeight’ Byers av ur Mnitial, velocity DU)=.5*rho*v coast (1)*24C_d*A; SInitial drag; a_coast (1)=(-W-D(1)) /m; saéceleration, D=D(v 0) accoast (2)=(-W-D(1))/m; tacceleration, assum® D is still D(v_u) j,hiterating factor while v_coast (length(v_coast))>=0 ‘y_cBast (k)=v_coast (k~1) +t* (a_coast (k-1)+a_coast (k))/2; Strapezoidal integration point by point D(k)=.5*rhoty coast (k) ~2*C_d*A; Drag a_coast (k)=(-W-D(k))/mj_8¥edefine acceleration with Drag cast (kt1)=(-W-D(k)}/mz Stemporary 0 assumption end. vle=v_coast (Length (y coast)) Treoasech ust® ("599 const (1) +v_coast (Length (v_coast))) +cunsum(v_coast ( foeoast(iJ=n us jinstral height hrcoast (2) Hust! (y-coast (1) sv coast (2))/2; ttrapezoidal integration of v from 0 to t fresh coast (Zength (eeoast)) 7 Ire fgeteh_c*106/(2.54°12)7 theight in feet Gisp(cat end of coaste--*) Gisplcat gr, Velocity ‘2s num2atecy.e)," m/e, (otogped)’)) displcat (2,' Height Le *,ndmastethieTs*’m or * numistr(h « length (v_coast)))) ty,' £0")) At end of coast... Velocity is 0.035923 m/s (stopped) Height is 79.2591 m or 260.0365 ft Concatenated Profiles Thrust vs. Time ‘Tecat (1,mean(T_B6,2),mean(T_A8,2))"; figure (i) sthrust’ Pict hold on Lot (0: ,01:1.31,7,'-k') S$Plot of averages legend (‘average’ thrust", 3) plot (0:.01:.8,7_86,".")’AB6-0 Plots plot (.8:.01:11377 Ae, ".*) SA8-3 Plots title(*Thrust vs.—Tine*), ylabel ("Thrust (M)'),xlabel (Timo (3) get (figure (1), "cotor', ‘whlee") SPeLight brorite ; velocity=cat (2,v boost, v_upper,v_coast) height=cat (2,h’boost,h’upper,h coast); figure (a) | subplot (1,2,1) Velocity Profile plot (0;.04; (length (velocity) ~1) /100, velocity, '-k") axis({0°5 0 max(velocity) +5}) title('velocity vs. Time’), ylabel (‘Velocity (m/s)'),xlabel('Time (s)') subplot (1,2,2) sHeight Profile plot (0: 01; (length (height) ~1) /100,height, ‘axis ({0'5 0 max (height) +10) ) title(*Height vs. Tine") ,ylabel (‘Height (m)'),xlabel (‘Time (s) set (figure (2), color’, !white!) suptitie(*Flignt Profile") eat) [Link]A204 Rocket Project Thnust vs. Time Thuust (N) —— Average thrust 02 O4 O86 08 1 12 14 Time (8) Flight Profile Velocity vs. Time Height vs. Time locity (mvs) Bom 8 RS Height (m) sos 8 8 irae ape iatiia si) spina ta aiG Time (s) Time (8) [Link] [Link])Bobby Nordlund Taylor Newman Captain O'Reilly A201-5001 21 April 2010 Rocket Project: RockSim Model, Predictions, and Parts Listpha Newmar Nordlund) a Saraeeaereeeees ameter 0.5780 Ih, Span ameter: 3970 In ted stage mass 145.63 ny Margin: 3.71 Overstable pes Alpha (Newman/Nordlund Daeter 0.9760 I, Span dlametr: 35760 In 39 9 Selected Stage mse 104 9.0992 Iny Margin: 4.49 OverstableModified Estes Alpha (Newman/Nordlund) - Simulation summary Sim. Ress Er Alt. Ft, [Link]/s Acc. Ft./s/s Time App. Vel. deploy ft/s Alt. deploy Ft. CSimulation index 495 301.28,(zy uno) sssee" * (els) £0-2985¢8' ‘epi0u Jo WeUOH (wD) ESSH5"T '(H) BSSHSTO‘O “UOREAS Jo srPeY 6 1407 :S5eH "O00 °99 © “UL. 0000'0 914 Inds "Ut 0900°7T 21a PEPE g:adeus ‘snupered T « ada 8ua/AineAIod :1eHeVeW ‘OASeId “Ul ZT - Ul ZT - Sansa BINY>eIEG (cyu06) gece * (@ wi) £0-279840°2 :en40U Jo waLIOW (UD) ££9062'0 * (W) €2906200°0 :[Link] Jo sypeY 6 SZE'E :S8eH "UL OSZT'O 9D « uno [Link] 0 ase BA WOsy "UI 000S'0 :uOReDOT "Ut O0Sz"9 :4AT “Ul 29EC'O :aI “UL 0096°0 :aOBUH BULERRD soded :JeueneW 'SS-07-UD - SZTOE - S0ISy Buy! BujLEWUOD (eyuo6) eoscoe * (ew) £0-27992072 :ens0u Jo weWO (uD) €£9062'0 *(W) €/906L00'0 voReAAS jo srpey 6 Sze“C:SSEW "UE OSZTO :9D © UNO J0.0W Jo UO kw “UT 0005"0 :UONeDeT “UT 0062 “URN “UI ZSEL' :aI “UT 0096°0 140.8UH BuLERURD « soded 1eHeNeW ‘55-07-49 - SZTOE - Sorsg Buys GujLeqU2D (296) cease * (26x) 60-22£9976 ‘ens8u Jo wowOH (uD) 21262990 (Ww) 21262900'0 ;vOReNA Jo snpeY 6 Spz"0 'SeN '“U SBGO'D 9D « {UNO Ja}OW Jo seq 2X LOLS “UT 0005'Z {LORIN “Ut O/T {UBT UL ¥SES'OT : 7T# AIH “UT 0OT' :ACSPOR AUR « sodeg :Jeu0IeW ‘07-19 P14 “Gua - Z-Z9TOE - saIS3 >POIg BU6U (0 :°N9 * “ut 0000"0 :Nx wisypeW *(ZyW96) HO06T'L * (sub) c0-240061' sensu Jo woUdH (uD) £2021'z () £20710"0 ‘vOReAKS 0 SKIpeY 6 66ST :SSEW "UL OSLET “9D © ‘21 Apo Jo 3569 a4 HOH “UT 000" UOC) “UE 0OSL" [UAT “UL ODT" ZAI“ OSEL'O :00 © sodeg yeuoreW ‘OZ-L8 - SBOE 183 - 50953 IUNOW 10; (0 Fen * “Ur 0000"0 :Nx wuSIPOW *(Z 96) BOLB6T * (@ iby) $0-280486"F 40440 WoUON (uD) BH45U'S *() ByEELSO‘O “UOREAS Jo SFIPEY B T00'9 :S5EW "UL OSLA'E ‘SD © "ul 00S¢'c UR1 "Ut 00S6°0 ‘aT * "Ul 0526°0 :00 + sodeg s1euoveW ‘0S-19 - 980€ 18a - saxsa 0qm Apos, (2106) 0 (evitBx) 0 zeps0u 0 yWaWOy (uD) 0 (ww) 9 :uOReAB 0 snpeY 6 OOL'ST 858M “UI 0000.0 5 « op HeUaeW “aBWINIY - ZAY OM/ASTAIV - SULPe}V04 JOWALY (226) 0* (2 vw) 0 zen.0U 0 WOLOW (w2) 0 (W) 0 :uOReAAB Jo srupey B 000'0S :358W * "UI 0000'0 :59 « ‘esjeg :1eueveW 'BET-ONM - O806T - e9Body peoed 22ND "Ul 2002's °NX WHS>POY ' (v6) 6>EO'HY * (Zw) 90-26460F'y :esoU jo weWON (uD) TZHEY'E * (WH) TzHHYED'O :UOREAAB J Shipey 6 C1" :SSEW "Ul LT6E'Z $99 "UL 0926°0 :UB1 “UT 006° [00 954 APO "UT 0946'0 £210 “UL OOSC2 "UR "BNO PHOS fades 35ON « sje ‘reuDyeW “Wwoysn9 - U0 950N syed seureysnsLaunch lug Estes - "1/8" per #0: 0.1720 In, 1D: 0.1560 In, Len: 2.0000 In, Loc: 0.3750 In CG: 1.0000 In, Mas: 0.162 g Radu of gyration: 0.0147552 (m) , 1.47552 (cm) Moment of inertia: 3.53335e-08 (kgm2) , (0.353335 (gcm*2) "=, Material: Fin set - Custom, Material: Balsa * OG: 2.0298 In. , Mass: 1.551 g Radius of gyration: 0.0252479 (m) , 2.52479 (cm) Moment of inertia: 9.8882e-07 (kam*2), ‘9.8882 (gem™2) , RockSim XN: 9.9263 In. , CNa: 18.9332 Booster parts Body tube Estes - EST 3086 - BT-50, Material: Paper 00: 0.9760 In, ID: 0,9500 In. , Len: 2.7500 In. ‘CG: 1.3750 In.» Mass: 1.987 g Radius of gyration: 0.0219652 (m) , 2.19652 (cm) Moment of inertia: 9.588316-07 kgm™2) , 9.58831 (gem42) , RockSim XN: 0,0000 In. ,CNa: 0 Motor Mount Estes - EST 3085 - BT-20, Material: Paper * 00: 0.7360 In, , 1D: 0.7100 In. Len: 2.7500 In. Location: 0.0000 In, From the front of Body tube © CG: 1.3750 In. « Mass: 1.492.g Radius of gyration: 0.0212077 (m) , 2.12077 (cm) Moment of inertia: 6.710722-07 (kgm2) , 6.71072 (gam*2), RockSim XN: 0.0000 In. , CHa: Engine hook Estes - Engine hook - Standard size, Material: * 0G: 0.0000 In. , Mas: 1.400 g Racus of gyration: 0 (rm) , 0 (cm) Moment of inertia: © (kam/2) , 0 (gcm”2) Centering ring Estes - 30125 - CR-20-55, Material: Paper ‘Centering ringOD: 0,9600 In, ID: 0.7362 In, Len: 0.2500 In. Location: 0.7500 In. From the base of Motor Mount "CG: 0.1250 In, , Mass: 3.325 q Radius of gyration: 0.00790673 (m) , 0.790673 (cm) Moment of inertia: 2.07867e-07 (kam“2) , 2.07967 (gem*2) ‘Tube Coupler Estes - 30125 - CR-20-55, Material: Paper Centering ringOD: 0.9600 n, 1D: 0.7362 In, Len: 1.0000 In. Location: -0.5000 In, From the front of Motor Mount * CG: 0.5000 In, , Mass: 13.300 g Radius of gyration: 00106347 (rm) , 1.05317 (cm) Moment of inertia: 1.50334e-06 (karn2), 15.0834 (gom*2) Fin set - Custom, Material: Balsa * CG: 1.4694 In. , Mass: 1.299 g Radius of gyration: 0.017447 (m), 1.73447 (cm) Moment of inertia: 3.90662¢-07(hgm*2) 3.90662 (gem™2) ,RockSim XN: 11.6108 In , CNa: 18.3099Bobby Nordlund Taylor Newman Captain O'Reilly A201-5001 21 April 2010 Rocket Wind Tunnel Lab Data Reduction a. Calculated air density: 1.1804 kg/m? b. Rocket's cross-sectional area: 4.82710 m? [see MATLAB script} d. Calculated drag coefficient: 1.2384 Drag va. Velocity Squared "Zoo —so0 tooo 200 1400 ve00 Tao 7000 —Z200 2400700 Volocty Squared (nis?) €. Calculated r value: 0.99845, f. [see MATLAB script) ‘Notes: The obviously incorrect value calculated for the drag coefficient demands a reassessment of the Rocket's actual cross-sectional area. With a cross-sectional area of 7.45x10* m’, the rocket yields a drag coefficient equal to 0.80237—a much more realistic result. \ rt Win Arh ogy Exell + : cise up wth tas?A204 Drag Lab EA204 Drag Lab Contents ‘Lab Data ‘Data Reduction ‘Drag Calculations Plot (for drag coefficient) Nordlund and Newman 21 March 2010 ‘This script seeks to determine drag coefficient from collected drag data. format short, clear all; ele Lab Data v=60, 70,80, 90,100,110); aWind Tunnel Velocity (mph) vev.*5286/3600#12*5.54/100;_ 4(m/a} drag (1,:)=[.05, .08, -1,.13,.16,.19]; $Axial force while streamlined (1b) rag (2, :)=[ 106, :08, :10,.13,.16, .19]7 saxial force 'w/ misaligned fins (1b) Data Reduction ‘Air Density Calculations (rho) P_ambient=29.8; Sinlig Plambient=P_ambient/29,92*1.0132e5; %Pa ‘Tambient=536.87. 88 ‘Tlambient-T_ambient*20a/519; x Raire287; 13/ (Kg"K) rho=(P_ambient)/(R_air*T ambient); kg/m”3 disp(cat (2, "Air Derisity:~", numzst® (rho), " kg/m*3")) Air Density: 1.1804 kg/m*3 Drag Calculations Cross-sectional area d=.976" or 2.48cm r=,0123952; radius (m) a=pitr*2; tcross-sectional area (m*2) disp (cat (2, ‘Cross-sectional area: ",num2str(a),' m°2")) % Calculate v2 dragemean (drag, 1) .*4.448; vaverage drag (to account for slop) (i) coettpolytit (¥."3,drag;1)7 Neoetfsesente of polynomial Jesum ((polyval (coeff, v.42) “drag) .*2) 7 musmean (drag) 5 Sesum ( (drag-mu) .*2) 7 z2e1-J/S; disp(cat (2, 'r°2 value (for D vs. v-2 curve fit) Calculate Drag Coefficient ¢ decoeft (1) /(a*rho/2); disp(eat (2, "Coefficient of drag: ‘,numastr(c a), M(enis clearly denotes troubled data) )) *-mumastr (x2) )) sectional area: 0.00040260 m2 2 value (for D'vs. V*2 curve fit): 0.99845 efficient of drag: 1.2384 (this clearly denotes troubled data) file /1N:/Astro\ntm/EA204_ [Link][4/29/2010 [Link])A704 Drag Lab Assume effective area nearer 7.45e-4 m*2 better=7.45e-4; more likely cross-sectional area (m~2) d_betterscoeff(1)/(a better*rho/2) ; disp (cat (2, ‘Likely coefficient of drag: ‘snumaste (C_d better), ... (assumes a=",num2str(a better), m-2 with £n3)7)) Likely coefficient of drag: 0.80237 (assumes a=0.000745 m*2 with fins) Plot (for drag coefficient) figure (a) Plot (v:~2, drag, '*k", (min (v) :max (v)) .*2,polyval (coeff, (min(v) :max(v)) .*2,'~k")) title(‘prag vs. Velocity Squared") ylabel ("Drag (i) ") xlabel ("Velocity Squared (m*2/s*2)") legend ("Data", ‘Linear Interpotation® ,1) text (2000, .50, 'm=c dta*zho/2") text (2000, .45, cat (2, 'C_d=',numgatr (Cd) )) reat (2, 'r™25" ynumastr (2))) 1); "color", ‘white’) Drag vs. Velocity Squared + Data Linear Interpolation + A boo eto 000 a0 Tao 1600 1800 m0 2200 Dan] aN Velocity Squared (m/s?) file//N:/Astrotfhtml/EAZ04_Draglab.htmil4/29/2010 [Link]}Bobby Nordlund Taylor Newman Captain O'Reilly EA201-5001 21 April 2010 Rocket Thrust Lab Data Reduction a 1804 kg/m? b. Thrust data taken every 0.01s for the published duration of burn. [see MATLAB script} ‘c._ Integration done in MATLAB. [see MATLAB script) Estes B6 Results: ‘Average total impulse: 4.3665 Ns Average specific impulse: 71.3318 s ‘Average thrust: 5.4582 N ‘Average maximum thrust: 10.8837 N Average standard deviation: 0.5635 N Estes 06 Tet Pi Te) Estes 86 Error Analy: “Published Value Measured Value |___% Error ‘Maximum Thrust (N) 12.10 10.88 7-10.05 Average Thrust (N) 6.25 5.46 Impulse (Ns) 5.00 437Estes AB Results: Published Value Measured Value ‘* Average total impulse: 2.1898 Ns ‘Average specific impulse: 35.7734 s + Average thrust: 4.3797 N '* Average maximum thrust: 12.0857 N ‘+ Average standard deviation: 0.72021 N aoa 0 Tee) Estes AB Error Analysis. ‘Maximum Thrust (N) 11.80 Average Thrust (N) 5.00 Impulse (Ns) 2.50 3 Error 2.42 “12.41 2.41 Conclusion: With motors in the A and B range, even the smallest standard deviation will result in significant variability in motor performance. Adding margin between the performance requirements and predicted performance is the best way to mitigate the demonstrated effects of motor inaccuracy. y “sg Nott values were used instead. ropellant mass comparison i Pee rele impossible due to unmeasured ejection charge mass. PublishedEA204 Thrust Lab 1of3 + Prcliminary Cteulatons + Estes A Data Reduction Nordlund and Newman 21 March 2010 ‘This script seeks to aggregate data fom test cell experiments and publish averaged results ‘This script will also display errr analysis results to compare published and measured values, format short, clear ally ele ‘Thrust Data ‘Note: Worthless data in 100106-4a, 2001a8-3, 2001a8-3b and data from S001b6-4a, (C465:C545) appears to be from a different size engine, Estes B6 Aperdsread(*chrusedaca’, “2001b6—4a", 1622510308") {2yexdsread(*thrustdata’ !001b6~6a", *C6BSiC765") 7 1) 3)eelaread(*thrustdata’, *3001b6-€b', "C318:6395) ‘A)sxlsread(*thrustdata’, "5001b6~€a",°C3532C633") 7 ‘Syexdsread(*chruseaaea’, "6001b6 60308854); 6 {:, )nalaread("ehrustdata', "6001b6—48","C463s0543") 5 Estes Re fr 1)=xLszead (*ehrustdata ‘2o01a8~3a', 1030820355") fad‘ ehruatdata',*3001a8-34"s"e260:0310"97 yad(*thrustdata', *S001a8-34", "€550:0600") + fad" thrustdata’,*5001a8-38',*C420:€520") yead|"thraatdata’, €004a8-38", *CS10;0560"93 ‘Dyeadsread("thrastdata’,*6001a8~3b", *C475:525°) 7 Preliminary Calculations Pamosont=29. ambient? fambient~536-87 aR ‘Tanbient=T anbiont*200/S197 4K 177437 kgK) Fhow (P_anbient) /(Rair*T ambient); thg/n°3 Gtsplcae (2, "Air Seaeity:*youmbatx (eho) k9/0°2")) sat Estes B6 Data Reduction Reference: 5.0 Ns impulse, 8 s bum time, 6.24 g propellant dtepl86 weoute:) rape (F_B6,1)/100; ans 0 total impulse Average total inpulse:"ymum2ats(1_B6)," Ns")) 46/(.00624"9.81); saverage apecttic inp *sgunastr (Eep_B6)," ©*)) thrust: *;mumaste(T (0,2),2); "tae Slopleaeray average macinia thrusts fomumdeee(@ BE max)" MD) TBE stdengan(eed(? B6,0,2)) 1 Aispieat (2s Average standard deviation: *,mumaste(3_B6_sta)," N")) Areay Propellant mass comparison impossible due to unmeasured ejection Error=(‘category’, "Published Value", ‘Measured Value’, "¥ Exror')7 [Link] hil 4/29/2010 13:12EA204 Thrust Lab 20f3 Beror (2s Error (2! ADT max. H's °T ave. (M's Iapulse (Hs) "2 2)=(12.176-25/5.0)7 Error (2:4,3)=(7_B6 maxi?_B6 average; 1 86); pereent=100*((1F_B6 maxi _BG_average/I_B6]~(12.1/6.25;5.0}) «/{12-1;6.25/5.01); Error (2:4,4)={percent (1) ;percant (2) spezeant (3) 2 atsp('a6 trror analysis") tepterzer) Wthrost Plot figured) old on plot (0: .01:.8,nean(7_B6,2),"-k*,0:-01:-8,medtan(1_B6,2),":4") Yogend\ oan thruct';"moaien erst") plot (0!.019.8,7.B6,"<") Eitle (‘Estes BE Thrust Profile"), xlabel ("Time (s)"),ylabel (Thrust (4) *) set (figured), ‘color", “wnite') Average total impulse: 4.3665 Ns Average specific inpuise! 71.3318 # Avorage thrust: 5.4562 N Average maximin thrast: 10.8637 8 eaory” ‘Published Value’ ‘Measured Value' "8 Err trae. wt 2.1000) { 10.8837) {-10.0522) saves asf 612800) | Si4se2] (71216698), Meputae (We)! 3s. 3665) (12/6604), Estes B& Thst Profle 001 02 03 04 05 06 07 08 Tie (9) Estes AS Data Reduction Reference: 2.5 Ns impulse, .5 s bum time, 3.12 g propellant tap (cas Resuieas') TAbserape(f AU, 1)/100) ms TAtenean(Z_AB,2); tAverage total impulse Gispreat (2,7 Average total inpsleet-mum2ate(Z_A8),* Ne") Tap ABsI_AG/(.00624"9.81); bAverage specific. tnpal aisplcatt2, "Average specific impulses." numdetx (Zap_A@)," ©") |_A8/.33 Sinpulse cividea oy burn tine ‘Avezaye thrust: ";nuadate(T AG average), N')) TAB saxenoan (max(T AG, (]/1),2))"ehvorage Aaximum thrust pleat (2, Average masini thrust: senumaste(2.A® max)" N*)) TMM staemean(sta(? A8,0,2))7 dispfeat 2," Average’ standard deviation: ‘,nua2ete(RAG_std)," N*)) error Array § Note: propelant mass comparison impossible due to unneasured ejection § charge’ mass file:/N:/Astro\hrnV/EA204_ThrustLab bul 4/29/2010 13:12EA204 Thrust Lab [Link] Thrust ab bil 30f3 Error={"Catogory', ‘Published Value', ‘Measured Valve's "8 Exxox") Error (2:4,2)=(°T tas fo. iw)" "Impulse (is); Error (2:4,2)=(11. | Error (2:4, 3)=(1 AB maxyt AB avorage/T_A®); percent~100*(((7_A8_nax/7_A8_average;I_A0}~ (12. Ercor (2:4, 4)={pexcent (1) :perGent (2) spezeent (3) GispU’Ab Error Analyese") isp (Error) throat, Plot figure (2) hold on lot (0: 01; .5,mean(t_A8,2),'-X',0:.01:.5,median(T_) Tegend(?mean thrust"; median thrust") LOE (0: 015.8) T.A8,".") Set (figure 2), "Color", white") title ("estes AS Thrust Profiie'),zlabel (Tine (2)),ylabel(' Thrust (M)") 5.052.519 ./112.85.0;2.8)) 5 29, Average total impulse: 2.1898 ms Average thrust: 4.3797 Average maximum thrust: 12.0857 6 Average standard deviation: 0.72021 8 KCategory’ “Published Value’ “Measured Value * Error" Tmax. iW)" 11-8000) 12.0857) (3.4216) rave. mM)" 31 43797) [-12/4060) Mmpule imey* 2.8000) 251680) {-12 4060), Estes Ab Tavs Profle “ —Wean trot 2 Median thrust 10 8 € % 6 4 “9 G05 OT 018 02 0m 03 om 04 0% OS Time (8) 4/29/2010 13:12