ตัวถัง Sears–Haack รูปทรง Sears –Haack เป็นรูปทรงที่มีแรงต้านคลื่น ตามทฤษฎีต่ำที่สุด ในการไหลแบบเหนือเสียง สำหรับวัตถุแข็งทรง กระบอกเรียวที่มีความยาวและปริมาตรที่กำหนด การพิสูจน์ทางคณิตศาสตร์นั้นถือว่าการไหลแบบเหนือเสียงที่มีการรบกวนน้อย (เชิงเส้น) ซึ่งควบคุมโดยสมการ Prandtl–Glauert การพิสูจน์และรูปทรงนี้ได้รับการตีพิมพ์โดยอิสระจากนักวิจัยสองคน ได้แก่Wolfgang Haack ในปี 1941 และต่อมาโดยWilliam Sears ในปี 1947 [ 1 ] [ 2 ] [ 3 ]
ทฤษฎีคาร์มัน-มัวร์ ระบุว่า แรงต้านของคลื่นแปรผันตามกำลังสองของอนุพันธ์อันดับสองของการกระจายพื้นที่ดี คลื่น ~ [ เอส " ( x ) ] 2 {\displaystyle D_{\text{wave}}\sim \left[S''\left(x\right)\right]^{2}} (ดูนิพจน์เต็มด้านล่าง) ดังนั้นสำหรับแรงต้านคลื่นต่ำ จึงจำเป็นต้องมีเอส ( x ) {\displaystyle S(x)} มีลักษณะเรียบเนียน ดังนั้น รูปทรงของ Sears–Haack จึงมีลักษณะแหลมที่ปลายทั้งสองข้าง และค่อยๆ ขยายตัวอย่างราบรื่นไปจนถึงจุดสูงสุด จากนั้นจึงค่อยๆ ลดลงอย่างราบรื่นไปสู่จุดที่สอง
เนื่องจากสูตรเหล่านี้ตั้งอยู่บนสมมติฐานว่าอนุภาคมีรูปร่างเรียว จึงถือว่าL >> 2 R (กล่าวคืออัตราส่วนความละเอียด สูง f = L / 2 R ) แต่ไม่ได้ระบุขีดจำกัดบนที่เฉพาะเจาะจงสำหรับR (หรือขีดจำกัดล่างสำหรับความละเอียด) อาร์ สูงสุด = 4 π วี 3 แอล วี = 3 π 2 16 อาร์ สูงสุด 2 แอล แอล = 16 วี 3 π 2 อาร์ สูงสุด 2 เอฟ := แอล 2 อาร์ สูงสุด = π 8 3 แอล 3 วี = 8 วี 3 π 2 อาร์ สูงสุด 3 เอส สูงสุด = 16 วี 3 π แอล = π อาร์ สูงสุด 2 เอส ( x ) = 128 วี 3 π แอล ( x ( 1 − x ) ) 3 / 2 = 8 π อาร์ สูงสุด 2 ( x ( 1 − x ) ) 3 / 2 ร ( x ) = 8 π 2 วี 3 แอล ( x ( 1 − x ) ) 3 / 4 = 2 2 อาร์ สูงสุด ( x ( 1 − x ) ) 3 / 4 ร ′ ( x ) = 3 อาร์ สูงสุด ( 1 − 2 x ) 2 x ( 1 − x ) 4 , เอส ′ ( x ) = 12 π อาร์ สูงสุด 2 ( 1 − 2 x ) x ( 1 − x ) ร " ( x ) = − 3 อาร์ สูงสุด ( 1 + 4 x ( 1 − x ) ) 4 2 ( x ( 1 − x ) ) 5 / 4 , เอส " ( x ) = − 6 π อาร์ สูงสุด 2 ( 8 ( 1 − x ) x − 1 ) x ( 1 − x ) {\displaystyle {\begin{aligned}R_{\text{max}}&={\frac {4}{\pi }}{\sqrt {\frac {V}{3L}}}\\V&={\frac {3\pi ^{2}}{16}}R_{\text{max}}^{2}L\\L&={\frac {16V}{3\pi ^{2}R_{\text{max}}^{2}}}\\f&:={\frac {L}{2R_{\text{max}}}}={\frac {\pi }{8}}{\sqrt {\frac {3L^{3}}{V}}}&&={\frac {8V}{3\pi ^{2}R_{\text{max}}^{3}}}\\S_{\text{max}}&={\frac {16V}{3\pi L}}&&=\pi R_{\text{max}}^{2}\\S(x)&={\frac {128V}{3\pi L}}{\bigl (}x(1-x){\bigr )}^{3/2}&&=8\pi R_{\text{max}}^{2}{\bigl (}x(1-x){\bigr )}^{3/2}\\r(x)&={\frac {8}{\pi }}{\sqrt {\frac {2V}{3L}}}{\bigl (}x(1-x){\bigr )}^{3/4}&&=2{\sqrt {2}}R_{\text{max}}{\bigl (}x(1-x){\bigr )}^{3/4}\\r'(x)&={\frac {3R_{\text{max}}(1-2x)}{{\sqrt {2}}{\sqrt[{4}]{x(1-x)}}}},&&S'(x)=12\pi R_{\text{max}}^{2}(1-2x){\sqrt {x(1-x)}}\\r''(x)&=-{\frac {3R_{\text{max}}{\bigl (}1+4x(1-x){\bigr )}}{4{\sqrt {2}}{\bigl (}x(1-x){\bigr )}^{5/4}}},&&S''(x)=-{\frac {6\pi R_{\text{max}}^{2}(8(1-x)x-1)}{\sqrt {x(1-x)}}}\\\end{aligned}}}
ที่ไหน:
x คืออัตราส่วนของระยะห่างจากจมูกต่อความยาวลำตัวทั้งหมด (ซึ่งจะมีค่าอยู่ระหว่าง 0 ถึง 1 เสมอ)ร ( x ) {\displaystyle r(x)} คือรัศมีในพื้นที่ร ′ ( x ) {\displaystyle r'(x)} คือความลาดชันในพื้นที่เอส ( x ) {\displaystyle S(x)} คือพื้นที่หน้าตัดเฉพาะที่อาร์ สูงสุด {\displaystyle R_{\text{max}}} คือรัศมีที่มีค่าสูงสุด (เกิดขึ้นที่x = 0.5 ซึ่งเป็นจุดศูนย์กลางของรูปทรง)เอส สูงสุด {\displaystyle S_{\text{max}}} คือพื้นที่หน้าตัดสูงสุดเมื่อx = 0.5 และร ( x ) {\displaystyle r(x)} =อาร์ สูงสุด {\displaystyle R_{\text{max}}} ,V คือปริมาตรL คือความยาวf คือ อัตราส่วน ความละเอียด จากทฤษฎีคาร์มาน–มัวร์ สรุปได้ว่า: ดี คลื่น = − 1 4 π ρ ยู 2 ∫ 0 ℓ ∫ 0 ℓ เอส " ( x 1 ) เอส " ( x 2 ) ln | x 1 − x 2 | ง x 1 ง x 2 , {\displaystyle D_{\text{wave}}=-{\frac {1}{4\pi }}\rho U^{2}\int _{0}^{\ell }\int _{0}^{\ell }S''(x_{1})S''(x_{2})\ln |x_{1}-x_{2}|\mathrm {d} x_{1}\mathrm {d} x_{2},}
หรืออีกทางเลือกหนึ่ง:
ดี คลื่น = − 1 2 π ρ ยู 2 ∫ 0 ℓ เอส " ( x ) ง x ∫ 0 x เอส " ( x 1 ) ln ( x − x 1 ) ง x 1 . {\displaystyle D_{\text{wave}}=-{\frac {1}{2\pi }}\rho U^{2}\int _{0}^{\ell }S''(x)\mathrm {d} x\int _{0}^{x}S''(x_{1})\ln(x-x_{1})\mathrm {d} x_{1}.}
สามารถนำสูตรเหล่านี้มาผสมผสานกันเพื่อให้ได้ผลลัพธ์ดังต่อไปนี้:
ดี คลื่น = 64 วี 2 π แอล 4 ρ ยู 2 = 9 π 3 อาร์ ม เอ x 4 4 แอล 2 ρ ยู 2 , ซี ดี คลื่น = 24 วี แอล 3 = 9 π 2 อาร์ ม เอ x 2 2 แอล 2 , เอ = 64 วี ρ π แอล 4 ρ ข ยู 2 = 12 π อาร์ ม เอ x 2 ρ แอล 3 ρ ข ยู 2 , {\displaystyle {\begin{aligned}D_{\text{wave}}&={\frac {64V^{2}}{\pi L^{4}}}\rho U^{2}&&={\frac {9\pi ^{3}R_{max}^{4}}{4L^{2}}}\rho U^{2},\\C_{D_{\text{wave}}}&={\frac {24V}{L^{3}}}&&={\frac {9\pi ^{2}R_{max}^{2}}{2L^{2}}},\\a&={\frac {64V\rho }{\pi L^{4}\rho _{\text{b}}}}U^{2}&&={\frac {12\pi R_{max}^{2}\rho }{L^{3}\rho _{\text{b}}}}U^{2},\\\end{aligned}}}
ที่ไหน:
เนื่องจากพลังงานจลน์ E ของวัตถุขึ้นอยู่กับปริมาตร ความหนาแน่น และความเร็วเท่านั้น จึงสามารถแสดงได้ดังนี้:
อี = ρ ข วี ยู 2 2 ⇒ ยู 2 = 2 อี ρ ข วี เอ = 128 ρ อี π แอล 4 ρ ข 2 {\displaystyle {\begin{aligned}E&={\frac {\rho _{\text{b}}VU^{2}}{2}}\Rightarrow U^{2}={\frac {2E}{\rho _{\text{b}}V}}\\a&={\frac {128\rho E}{\pi L^{4}\rho _{\text{b}}^{2}}}\\\end{aligned}}}
อนุพันธ์ ตามทฤษฎีของคาร์มัน-มัวร์ แรงต้านของคลื่นมีค่าดังนี้
เอฟ = − ρ คุณ 2 2 π ∫ 0 ล ∫ 0 ล เอส " ( ξ 1 ) เอส " ( ξ 2 ) ln | ξ 2 − ξ 1 | ง ξ 1 ง ξ 2 {\displaystyle F=-{\frac {\rho u^{2}}{2\pi }}\int _{0}^{l}\int _{0}^{l}S''(\xi _{1})S''(\xi _{2})\ln |\xi _{2}-\xi _{1}|d\xi _{1}d\xi _{2}}
ที่ไหนเอส ( x ) {\displaystyle S(x)} คือพื้นที่หน้าตัดของร่างกายที่ตั้งฉากกับแกนของร่างกาย ในที่นี้x = 0 {\displaystyle x=0} แสดงถึงความล้ำหน้าและx = ล {\displaystyle x=l} คือขอบด้านท้าย แม้ว่าทฤษฎีของคาร์มัน-มัวร์จะไม่แยกแยะปลายทั้งสองนี้ เนื่องจากสัมประสิทธิ์แรงต้านไม่ขึ้นอยู่กับทิศทางการเคลื่อนที่ในทฤษฎีเชิงเส้น แทนที่จะเป็นเอส ( x ) {\displaystyle S(x)} เราสามารถกำหนดฟังก์ชันได้เอฟ ( x ) = เอส ′ ( x ) {\displaystyle f(x)=S'(x)} และขยายความต่อไปเป็นชุดๆ
เอฟ ( θ ) = − ล ∑ n = 2 ∞ เอ n บาป n θ x ( θ ) = ล 2 ( 1 − คอส θ ) = ล บาป 2 ( θ 2 ) θ ( x ) = 3 π − 2 อาร์คคอส ( x ล ) = 2 π + 2 อาร์คซิน ( x ล ) {\displaystyle {\begin{aligned}f(\theta )&=-l\sum _{n=2}^{\infty }A_{n}\sin n\theta \\x(\theta )&={\frac {l}{2}}(1-\cos \theta )=l\sin ^{2}\left({\frac {\theta }{2}}\right)\\\theta (x)&=3\pi -2\arccos \left({\sqrt {\frac {x}{l}}}\right)=2\pi +2\arcsin \left({\sqrt {\frac {x}{l}}}\right)\\\end{aligned}}}
ที่ไหน0 ≤ θ ≤ π , 0 ≤ x ≤ ล {\textstyle 0\leq \theta \leq \pi ,\ 0\leq x\leq l} ซีรีส์เรื่องนี้เริ่มต้นจากn = 2 {\displaystyle n=2} เนื่องจากสภาพดังกล่าวเอส ( 0 ) = เอส ( ล ) = 0 {\displaystyle S(0)=S(l)=0} เรามี
เอส ( x ) = ∫ 0 x เอฟ ( x ) ง x , วี = ∫ 0 ล เอส ( x ) ง x = π 16 ล 3 เอ 2 {\displaystyle {\begin{aligned}&S(x)=\int _{0}^{x}f(x)dx,\quad V=\int _{0}^{l}S(x)dx={\frac {\pi }{16}}l^{3}A_{2}\\\end{aligned}}}
โปรดทราบว่าปริมาตรของวัตถุขึ้นอยู่กับค่าสัมประสิทธิ์เท่านั้นเอ 2 {\displaystyle A_{2}} .
ในการคำนวณแรงต้าน เราต้องเขียนสูตรแรงต้านใหม่ก่อน โดยใช้การอินทิเกรตโดยใช้การแยกส่วนหนึ่งครั้ง
เอฟ = พี . วี . ρ ยู 2 2 π ∫ 0 ล ∫ 0 ล เอฟ ( ξ 1 ) เอฟ ′ ( ξ 2 ) ง ξ 1 ง ξ 2 ξ 1 − ξ 2 {\displaystyle F=\mathrm {p.v.} {\frac {\rho U^{2}}{2\pi }}\int _{0}^{l}\int _{0}^{l}f(\xi _{1})f'(\xi _{2}){\frac {d\xi _{1}d\xi _{2}}{\xi _{1}-\xi _{2}}}}
ซึ่งพี . วี . {\displaystyle \mathrm {p.v.} } หมายถึงค่าหลักของโคชี (Cauchy principal value ) ตอนนี้เราสามารถแทนที่การกระจายสำหรับเอฟ {\displaystyle f} และทำการอินทิเกรตนิพจน์โดยใช้เอกลักษณ์สองข้อต่อไปนี้
พี . วี . ∫ 0 π คอส n θ 2 คอส θ 2 − คอส θ 1 ง θ 2 = π บาป n θ 1 บาป θ 1 , ∫ 0 π บาป n θ 1 บาป ม θ 1 ง θ 1 = π 2 { 1 ( ม = n ) , 0 , ( ม ≠ n ) . {\displaystyle \mathrm {p.v.} \int _{0}^{\pi }{\frac {\cos n\theta _{2}}{\cos \theta _{2}-\cos \theta _{1}}}d\theta _{2}={\frac {\pi \sin n\theta _{1}}{\sin \theta _{1}}},\quad \int _{0}^{\pi }\sin n\theta _{1}\sin m\theta _{1}d\theta _{1}={\frac {\pi }{2}}{\begin{cases}1\,\,(m=n),\\0,\,\,(m\neq n).\end{cases}}}
ผลลัพธ์สุดท้ายที่แสดงในรูปของสัมประสิทธิ์แรงต้านC จะได้รับโดย[ 4 ]
ซี ง = 2 เอฟ ρ คุณ 2 ล 2 , ล 2 = พื้นที่ = π 4 ∑ n = 2 ∞ n เอ n 2 {\displaystyle {\begin{aligned}C_{\mathrm {d} }&={\frac {2F}{\rho u^{2}l^{2}}},l^{2}={\text{area}}\\&={\frac {\pi }{4}}\sum _{n=2}^{\infty }nA_{n}^{2}\\\end{aligned}}}
เนื่องจากวี {\displaystyle V} ขึ้นอยู่กับเท่านั้นเอ 2 {\displaystyle A_{2}} ค่าต่ำสุดของเอฟ {\displaystyle F} ถึงเมื่อเอ n = 0 {\displaystyle A_{n}=0} สำหรับn ≥ 3 {\displaystyle n\geq 3} .
ดังนั้น การกำหนดเอ n = 0 {\displaystyle A_{n}=0} สำหรับn ≥ 3 {\displaystyle n\geq 3} เราจึงได้รับ
เอส = ล 3 เอ 2 บาป 3 θ 3 ซี ง = 128 π ( วี ล 3 ) 2 = 9 π 2 ( เอส ม เอ x ล 2 ) 2 ร ( x ) = 8 π ล 2 2 วี 3 [ x ( ล − x ) ] 3 4 = 8 π 2 วี 3 ล [ x ล ( 1 − x ล ) ] 3 4 {\displaystyle {\begin{aligned}S&={\frac {l^{3}A_{2}\sin ^{3}\theta }{3}}\\C_{d}&={\frac {128}{\pi }}\left({\frac {V}{l^{3}}}\right)^{2}={\frac {9\pi }{2}}\left({\frac {S_{\mathrm {max} }}{l^{2}}}\right)^{2}\\r(x)&={\frac {8}{\pi l^{2}}}{\sqrt {\frac {2V}{3}}}[x(l-x)]^{\frac {3}{4}}={\frac {8}{\pi }}{\sqrt {\frac {2V}{3l}}}\left[{\frac {x}{l}}\left(1-{\frac {x}{l}}\right)\right]^{\frac {3}{4}}\\\end{aligned}}}
ที่ไหนร ( x ) {\displaystyle r(x)} รัศมีเป็นฟังก์ชันของx ;x ล {\displaystyle {\frac {x}{l}}} เวอร์ชันนี้แสดงสูตรที่x สัมพันธ์กับl ดังเช่นใน ส่วน สูตรที่มีประโยชน์ ด้านบน
การสรุปโดยทั่วไปโดย อาร์ที โจนส์ การหาค่ารูปร่างของ Sears–Haack นั้นถูกต้องเฉพาะในกรณีที่รูปร่างนั้นเรียวบางเท่านั้น ทฤษฎีนี้ได้รับการขยายไปสู่รูปร่างที่เรียวบางแต่ไม่สมมาตรตามแกนโดยRobert T. Jones ในรายงาน NACA ฉบับที่ 1284 [ 5 ] ในส่วนขยายนี้ พื้นที่เอส ( x ) {\displaystyle S(x)} ถูกกำหนดบนกรวยมัค ซึ่งมีจุดยอดอยู่ที่ตำแหน่งx {\displaystyle x} แทนที่จะอยู่บนx = คงที่ {\displaystyle x={\text{constant}}} ระนาบตามที่ Sears และ Haack สมมติไว้ ดังนั้น ทฤษฎีของ Jones จึงสามารถนำไปใช้กับรูปทรงที่ซับซ้อนกว่า เช่นเครื่องบินความเร็วเหนือเสียง ทั้งลำ ได้
กฎพื้นที่ แนวคิดที่เกี่ยวข้องอย่างผิวเผินคือกฎพื้นที่ของวิทคอมบ์ (Whitcomb area rule ) ซึ่งระบุว่า แรงต้านคลื่นเนื่องจากปริมาตรในกระแสลมความเร็วเหนือเสียงนั้นขึ้นอยู่กับการกระจายตัวของพื้นที่หน้าตัดทั้งหมดเป็นหลัก และเพื่อให้แรงต้านคลื่นต่ำ การกระจายตัวนี้จะต้องราบเรียบ ความเข้าใจผิดที่พบบ่อยคือ รูปทรงของวัตถุแบบ Sears–Haack มีการกระจายพื้นที่ที่เหมาะสมตามกฎพื้นที่ แต่สิ่งนี้ไม่ถูกต้อง สมการ Prandtl–Glauert ซึ่งเป็นจุดเริ่มต้นในการหาที่มาของรูปทรงวัตถุแบบ Sears–Haack นั้นใช้ไม่ได้ในกระแสลมความเร็วเหนือเสียง ซึ่งเป็นบริเวณที่กฎพื้นที่ ใช้ได้
ลิงก์ภายนอก กระสุนปืนไรเฟิล Haack ที่มีแรงต้านน้อยที่สุดถูกเก็บถาวรเมื่อวันที่ 27 กันยายน 2007 บน เว็บไซต์ Wayback Machine ซึ่งเว็บไซต์ดังกล่าว ถูกปิดไปแล้ว – https://web.archive.org/web/20160306044740/http://www.lima-wiederladetechnik.de/englisch/haack_minimum_drag_bullet.htm Geschoßformen kleinsten Wellenwiderstandes โดย W. Haack, Bericht 139 der Lilienthal-Gesellschaft (1941) เครื่องคำนวณรูปร่าง Sears–Haack