ลูกตุ้มคู่


ลูกตุ้มที่มีลูกตุ้มอีกอันติดอยู่ที่ปลาย
ลูกตุ้มคู่ประกอบด้วยลูกตุ้ม 2 ลูก ที่ติดเข้าด้วยกันปลายทั้ง สองด้าน

ในฟิสิกส์และคณิตศาสตร์ในพื้นที่ของระบบไดนามิกลูกตุ้มคู่หรือที่เรียกอีกอย่างว่าลูกตุ้มสับสนคือลูกตุ้มที่มีลูกตุ้มอีกอันติดอยู่ที่ปลาย โดยก่อให้เกิดระบบฟิสิกส์ ที่เรียบง่าย ซึ่งแสดงพฤติกรรมไดนามิก ที่หลากหลาย พร้อมความไวสูงต่อสภาวะเริ่มต้น [ 1]การเคลื่อนที่ของลูกตุ้มคู่ถูกควบคุมโดยชุดสมการเชิงอนุพันธ์สามัญ ที่เชื่อมโยงกัน และ เป็นแบบสับสน

การวิเคราะห์และการตีความ

ลูกตุ้มคู่สามารถแบ่งออกเป็นหลายรูปแบบ โดยลูกตุ้มทั้งสองอาจมีความยาวและมวลเท่ากันหรือไม่เท่ากัน อาจเป็นลูกตุ้มเดี่ยวหรือลูกตุ้มประกอบ (เรียกอีกอย่างว่าลูกตุ้มประกอบ) และการเคลื่อนที่อาจเป็นสามมิติหรือจำกัดอยู่ในระนาบแนวตั้ง ในการวิเคราะห์ต่อไปนี้ ลูกตุ้มทั้งสองถือเป็นลูกตุ้มประกอบที่เหมือนกันทุกประการ โดยมีความยาวและมวลmและการเคลื่อนที่จำกัดอยู่ในสองมิติ

ลูกตุ้มคู่ผสม
การเคลื่อนที่ของลูกตุ้มประกอบคู่ (จากการบูรณาการเชิงตัวเลขของสมการการเคลื่อนที่)

ในลูกตุ้มประกอบ มวลจะกระจายไปตามความยาว หากมวลของลูกตุ้มคู่กระจายเท่าๆ กัน จุดศูนย์กลางมวลของแขนแต่ละข้างจะอยู่ที่จุดกึ่งกลาง และแขนทั้งสองข้างจะมีโมเมนต์ความเฉื่อยเท่ากับI = 1-12ม. 2เกี่ยวกับจุดนั้น

สะดวกที่จะใช้มุมระหว่างแต่ละแขนงและแนวตั้งเป็นพิกัดทั่วไป ใน การกำหนดค่าของระบบ มุมเหล่านี้แสดงด้วยθ 1และθ 2ตำแหน่งของจุดศูนย์กลางมวลของแกนแต่ละอันอาจเขียนในรูปของพิกัดทั้งสองนี้ หากถือว่าจุดกำเนิดของระบบพิกัดคาร์ทีเซียนอยู่ที่จุดที่ลูกตุ้มลูกแรกแขวนอยู่ จุดศูนย์กลางมวลของลูกตุ้มลูกนี้จะอยู่ที่:

และจุดศูนย์กลางมวลของลูกตุ้มที่ 2 อยู่ที่ นี่เป็นข้อมูลเพียงพอที่จะเขียนลากรองจ์ออกมาได้ เอ็กซ์ 2 - - บาป θ 1 - 1 2 บาป θ 2 - 2 - - คอส θ 1 - 1 2 คอส θ 2 - {\displaystyle {\begin{aligned}x_{2}&=\ell \left(\sin \theta _{1}+{\tfrac {1}{2}}\sin \theta _{2}\right)\\y_{2}&=-\ell \left(\cos \theta _{1}+{\tfrac {1}{2}}\cos \theta _{2}\right)\end{aligned}}}

ลากรองจ์

ลากรองจ์คือ เทอมแรกคือพลังงานจลน์เชิงเส้นของจุดศูนย์กลางมวลของวัตถุ และเทอมที่สองคือ พลังงานจลน์ของ การหมุนรอบจุดศูนย์กลางมวลของแท่งแต่ละแท่ง เทอมสุดท้ายคือพลังงานศักย์ของวัตถุในสนามโน้มถ่วงสม่ำเสมอ สัญกรณ์จุดระบุอนุพันธ์ ของ เวลาของตัวแปรที่เป็นปัญหา - พลังงานจลน์ พลังงานศักย์ - 1 2 ม. - วี 1 2 - วี 2 2 - - 1 2 ฉัน - θ ˙ 1 2 - θ ˙ 2 2 - ม. จี - 1 - 2 - - 1 2 ม. - เอ็กซ์ ˙ 1 2 - ˙ 1 2 - เอ็กซ์ ˙ 2 2 - ˙ 2 2 - - 1 2 ฉัน - θ ˙ 1 2 - θ ˙ 2 2 - ม. จี - 1 - 2 - {\displaystyle {\begin{aligned}L&={\text{พลังงานจลน์}}-{\text{พลังงานศักย์}}\\&={\tfrac {1}{2}}ม.\left(v_{1}^{2}+v_{2}^{2}\right)+{\tfrac {1}{2}}I\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-มก.\left(y_{1}+y_{2}\right)\\&={\tfrac {1}{2}}ม.\left({\dot {x}}_{1}^{2}+{\dot {y}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {y}}_{2}^{2}\right)+{\tfrac {1}{2}}ฉัน\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-มก\left(y_{1}+y_{2}\right)\end{aligned}}}

เนื่องจาก (ดูกฎลูกโซ่และรายการเอกลักษณ์ตรีโกณมิติ ) เอ็กซ์ ˙ 1 - θ ˙ 1 - 1 2 คอส θ 1 - เอ็กซ์ ˙ 1 2 - θ ˙ 1 2 - 1 4 2 คอส 2 θ 1 - ˙ 1 - θ ˙ 1 - 1 2 บาป θ 1 - ˙ 1 2 - θ ˙ 1 2 - 1 4 2 บาป 2 θ 1 - {\displaystyle {\begin{aligned}{\dot {x}}_{1}&={\dot {\theta }}_{1}\left({\tfrac {1}{2}}\ell \cos \theta _{1}\right)&&\rightarrow &{\dot {x}}_{1}^{2}&={\dot {\theta }}_{1}^{2}\left({\tfrac {1}{4}}\ell ^{2}\cos ^{2}\theta _{1}\right)\\[1ex]{\dot {y}}_{1}&={\dot {\theta }}_{1}\left({\tfrac {1}{2}}\ell \sin \theta _{1}\right)&&\rightarrow &{\dot {y}}_{1}^{2}&={\dot {\theta }}_{1}^{2}\left({\tfrac {1}{4}}\ell ^{2}\sin ^{2}\theta _{1}\right)\end{จัดแนว}}}

เอ็กซ์ ˙ 1 2 - ˙ 1 2 - 1 4 θ ˙ 1 2 2 - คอส 2 θ 1 - บาป 2 θ 1 - - 1 4 2 θ ˙ 1 2 - {\displaystyle {\dot {x}}_{1}^{2}+{\dot {y}}_{1}^{2}={\tfrac {1}{4}}{\dot {\theta }}_{1}^{2}\ell ^{2}\left(\cos ^{2}\theta _{1}+\sin ^{2}\theta _{1}\right)={\tfrac {1}{4}}\ell ^{2}{\dot {\theta }}_{1}^{2},}

และ

เอ็กซ์ ˙ 2 - - θ ˙ 1 คอส θ 1 - 1 2 θ ˙ 2 คอส θ 2 - เอ็กซ์ ˙ 2 2 - 2 - θ ˙ 1 2 คอส 2 θ 1 - θ ˙ 1 θ ˙ 2 คอส θ 1 คอส θ 2 - 1 4 θ ˙ 2 2 คอส 2 θ 2 - ˙ 2 - - θ ˙ 1 บาป θ 1 - 1 2 θ ˙ 2 บาป θ 2 - ˙ 2 2 - 2 - θ ˙ 1 2 บาป 2 θ 1 - θ ˙ 1 θ ˙ 2 บาป θ 1 บาป θ 2 - 1 4 θ ˙ 2 2 บาป 2 θ 2 - {\displaystyle {\begin{aligned}{\dot {x}}_{2}&=\ell \left({\dot {\theta }}_{1}\cos \theta _{1}+{\tfrac {1}{2}}{\dot {\theta }}_{2}\cos \theta _{2}\right)&&\rightarrow &{\dot {x}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\cos ^{2}\theta _{1}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \theta _{1}\cos \theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\cos ^{2}\theta _{2}\right)\\{\dot {y}}_{2}&=\ell \left({\dot {\theta }}_{1}\sin \theta _{1}+{\tfrac {1}{2}}{\dot {\theta }}_{2}\sin \theta _{2}\right)&&\rightarrow &{\dot {y}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\sin ^{2}\theta _{1}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\sin \theta _{1}\sin \theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\sin ^{2}\theta _{2}\right)\end{aligned}}} x ˙ 2 2 + y ˙ 2 2 = 2 ( θ ˙ 1 2 cos 2 θ 1 + θ ˙ 1 2 sin 2 θ 1 + 1 4 θ ˙ 2 2 cos 2 θ 2 + 1 4 θ ˙ 2 2 sin 2 θ 2 + θ ˙ 1 θ ˙ 2 cos θ 1 cos θ 2 + θ ˙ 1 θ ˙ 2 sin θ 1 sin θ 2 ) = 2 ( θ ˙ 1 2 + 1 4 θ ˙ 2 2 + θ ˙ 1 θ ˙ 2 cos ( θ 1 θ 2 ) ) , {\displaystyle {\begin{aligned}{\dot {x}}_{2}^{2}+{\dot {y}}_{2}^{2}&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}\cos ^{2}\theta _{1}+{\dot {\theta }}_{1}^{2}\sin ^{2}\theta _{1}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\cos ^{2}\theta _{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}\sin ^{2}\theta _{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \theta _{1}\cos \theta _{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\sin \theta _{1}\sin \theta _{2}\right)\\[1ex]&=\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \left(\theta _{1}-\theta _{2}\right)\right),\end{aligned}}}

การแทนที่พิกัดด้านบนและการจัดเรียงสมการใหม่จะได้ L = 1 2 m 2 ( θ ˙ 1 2 + 1 4 θ ˙ 1 2 + 1 4 θ ˙ 2 2 + θ ˙ 1 θ ˙ 2 cos ( θ 1 θ 2 ) ) + 1 24 m 2 ( θ ˙ 1 2 + θ ˙ 2 2 ) m g ( y 1 + y 2 ) = 1 6 m 2 ( θ ˙ 2 2 + 4 θ ˙ 1 2 + 3 θ ˙ 1 θ ˙ 2 cos ( θ 1 θ 2 ) ) + 1 2 m g ( 3 cos θ 1 + cos θ 2 ) . {\displaystyle {\begin{aligned}L&={\tfrac {1}{2}}m\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{1}^{2}+{\tfrac {1}{4}}{\dot {\theta }}_{2}^{2}+{\dot {\theta }}_{1}{\dot {\theta }}_{2}\cos \left(\theta _{1}-\theta _{2}\right)\right)+{\tfrac {1}{24}}m\ell ^{2}\left({\dot {\theta }}_{1}^{2}+{\dot {\theta }}_{2}^{2}\right)-mg\left(y_{1}+y_{2}\right)\\[1ex]&={\tfrac {1}{6}}m\ell ^{2}\left({\dot {\theta }}_{2}^{2}+4{\dot {\theta }}_{1}^{2}+3{{\dot {\theta }}_{1}}{{\dot {\theta }}_{2}}\cos(\theta _{1}-\theta _{2})\right)+{\tfrac {1}{2}}mg\ell \left(3\cos \theta _{1}+\cos \theta _{2}\right).\end{aligned}}}

สมการออยเลอร์-ลากรองจ์จะให้สมการเชิงอนุพันธ์ลำดับที่สองที่ไม่เป็นเชิงเส้นสองสมการต่อไปนี้ใน[a] : ( θ 1 , θ 2 ) {\displaystyle (\theta _{1},\theta _{2})}

m 2 ( g sin ( θ 1 ) + 2 ( θ ˙ 2 2 sin ( θ 1 θ 2 ) + θ ¨ 2 cos ( θ 1 θ 2 ) ) + 1 θ ¨ 1 ) + m 1 ( g sin ( θ 1 ) + 1 θ ¨ 1 ) = 0 g sin ( θ 2 ) + 1 ( θ ¨ 1 cos ( θ 1 θ 2 ) θ ˙ 1 2 sin ( θ 1 θ 2 ) ) + 2 θ ¨ 2 = 0 {\displaystyle {\begin{aligned}m_{2}\left(g\sin(\theta _{1})+\ell _{2}\left({\dot {\theta }}_{2}^{2}\sin(\theta _{1}-\theta _{2})+{\ddot {\theta }}_{2}\cos(\theta _{1}-\theta _{2})\right)+\ell _{1}{\ddot {\theta }}_{1}\right)+m_{1}\left(g\sin(\theta _{1})+\ell _{1}{\ddot {\theta }}_{1}\right)&=0\\g\sin(\theta _{2})+\ell _{1}\left({\ddot {\theta }}_{1}\cos(\theta _{1}-\theta _{2})-{\dot {\theta }}_{1}^{2}\sin(\theta _{1}-\theta _{2})\right)+\ell _{2}{\ddot {\theta }}_{2}&=0\\\end{aligned}}} ยัง ไม่มี การทราบวิธีแก้ปัญหา รูปแบบปิดสำหรับθ 1และθ 2ในฐานะฟังก์ชันของเวลา ดังนั้นการแก้ระบบจึงสามารถทำได้เพียงเชิงตัวเลข เท่านั้น โดยใช้ วิธี Runge Kuttaหรือ เทคนิคที่คล้ายคลึงกัน

กราฟพาราเมตริกสำหรับวิวัฒนาการของมุมของลูกตุ้มคู่ สังเกตได้ว่ากราฟนี้คล้ายกับการเคลื่อนที่แบบบราวน์

การเคลื่อนไหวที่วุ่นวาย

กราฟแสดงเวลาในการพลิกลูกตุ้มเป็นฟังก์ชันของเงื่อนไขเริ่มต้น
การเปิดรับแสงเป็นเวลานานของลูกตุ้มคู่ที่แสดงการเคลื่อนไหวที่สับสน (ติดตามด้วยLED )

ลูกตุ้มคู่เคลื่อนที่แบบสับสนและแสดงให้เห็นอย่างชัดเจนถึงความสัมพันธ์ที่ละเอียดอ่อนกับสภาวะเริ่มต้นรูปภาพทางขวาแสดงจำนวนเวลาที่ผ่านไปก่อนที่ลูกตุ้มจะพลิกกลับ โดยเป็นฟังก์ชันของตำแหน่งเริ่มต้นเมื่อปล่อยขณะหยุดนิ่ง โดยค่าเริ่มต้นของθ 1จะอยู่ในช่วง ทิศทาง xตั้งแต่ −3.14 ถึง 3.14 ค่าเริ่มต้นθ 2จะอยู่ในช่วง ทิศทาง yตั้งแต่ −3.14 ถึง 3.14 สีของแต่ละพิกเซลระบุว่าลูกตุ้มทั้งสองจะพลิกกลับภายใน:

  • g {\displaystyle {\sqrt {\frac {\ell }{g}}}} (สีดำ)
  • 10 g {\displaystyle 10{\sqrt {\frac {\ell }{g}}}} (สีแดง)
  • 100 g {\displaystyle 100{\sqrt {\frac {\ell }{g}}}} (สีเขียว)
  • 1000 g {\displaystyle 1000{\sqrt {\frac {\ell }{g}}}} (สีฟ้า)หรือ
  • 10000 g {\displaystyle 10000{\sqrt {\frac {\ell }{g}}}} (สีม่วง).
ลูกตุ้มคู่สามอันที่มีเงื่อนไขเริ่มต้นที่เกือบจะเหมือนกันจะแยกออกจากกันเมื่อเวลาผ่านไป แสดงให้เห็นถึงธรรมชาติที่วุ่นวายของระบบ

เงื่อนไขเริ่มต้นที่ไม่นำไปสู่การพลิกภายในจะถูกวางแผนเป็นสีขาว 10000 g {\displaystyle 10000{\sqrt {\frac {\ell }{g}}}}

ขอบเขตของพื้นที่สีขาวตอนกลางถูกกำหนดโดยการอนุรักษ์พลังงานบางส่วนด้วยเส้นโค้งต่อไปนี้: 3 cos θ 1 + cos θ 2 = 2. {\displaystyle 3\cos \theta _{1}+\cos \theta _{2}=2.}

ภายในบริเวณที่กำหนดโดยเส้นโค้งนี้ นั่นก็คือ ถ้า เป็นเช่นนั้น ลูกตุ้มทั้งสองจะพลิกไม่ได้ในเชิงพลังงาน นอกบริเวณนี้ ลูกตุ้มสามารถพลิกได้ แต่การกำหนดว่าจะพลิกเมื่อใดนั้นเป็นเรื่องที่ซับซ้อน พฤติกรรมที่คล้ายกันนี้สังเกตได้จากลูกตุ้มคู่ที่มีมวลจุด สองจุด แทนที่จะเป็นแท่งสองแท่งที่มีมวลกระจาย[2] 3 cos θ 1 + cos θ 2 > 2 , {\displaystyle 3\cos \theta _{1}+\cos \theta _{2}>2,}

การขาดความถี่ในการกระตุ้นตามธรรมชาติทำให้มีการใช้ระบบลูกตุ้มคู่ในการออกแบบเพื่อต้านทานแผ่นดินไหวในอาคาร โดยที่ตัวอาคารนั้นเป็นลูกตุ้มคว่ำหลัก และมีการเชื่อมต่อมวลรองเพื่อทำให้ลูกตุ้มคู่เสร็จสมบูรณ์

ดูเพิ่มเติม

หมายเหตุ

  1. ^ สมการที่ได้มาจาก โค้ด Mathematica ต่อไปนี้ :
    บล็อก[ { m , g , \ [ Theta ], l , L , x , y , v , t }, x = สะสม[{ ตัวห้อย[ l , 1 ], ตัวห้อย[ l , 2 ]} * Sin [{ ตัวห้อย[ \ [ Theta ], 1 ], ตัวห้อย[ \ [ Theta ], 2 ]}]]; y = สะสม[ { ตัวห้อย[ l , 1 ], ตัวห้อย[ l , 2 ]} *- Cos [{ ตัวห้อย[ \ [ Theta ], 1 ], ตัวห้อย[ \ [ Theta ], 2 ]}]]; v = D [{ x , y } /. ตัวห้อย[ \ [ Theta ], i_ ] :> ตัวห้อย[ \ [ Theta ], i ][ t ], t ]; L = บวก@@ ({ ห้อย[ m , 1 ], ห้อย[ m , 2 ]} * ( 1 / 2 แผนที่[ # . # & , ทรานสโพส[ v ]] - ( g y /. ห้อย[ \ [ ธีตา], i_ ] :> ห้อย[ \ [ ธีตา], i ][ t ]))); FullSimplify [ ตาราง[ D [ สร้าง[ ฟังก์ชัน, L /. ห้อย                                                                           [ \ [ Theta ], i ] ' [ t ] -> # ] ' [ ตัวห้อย[ \ [ Theta ], i ] ' [ t ]], t ] == สร้าง[ ฟังก์ชันL /. ตัวห้อย[ \ [ Theta ], i ][ t ] -> # ] ' [ ตัวห้อย[ \ [ Theta ], i ][ t ]] , { i , 2 } ], สมมติฐาน-> { ตัวห้อย[ l , 1 ] > 0 , ตัวห้อย[ l , 2 ] > 0 , ตัวห้อย[ m , 1 ] > 0 , ตัวห้อย[ m , 2 ] > 0 } ] /. h_ [ t ] : > h // คอลัมน์// TeXForm ]                                                    

อ้างอิง

  1. ^ Levien, RB; Tan, SM (1993). "Double Pendulum: An experiment in chaos". American Journal of Physics . 61 (11): 1038. Bibcode :1993AmJPh..61.1038L. doi :10.1119/1.17335.
  2. ^ Alex Small, ตัวอย่างโครงการสุดท้าย: ลายเซ็นแห่งความโกลาหลในลูกตุ้มคู่ (2013) รายงานที่จัดทำขึ้นเพื่อเป็นตัวอย่างสำหรับนักเรียน รวมถึงการหาอนุพันธ์ของสมการการเคลื่อนที่ และการเปรียบเทียบระหว่างลูกตุ้มคู่ที่มีมวล 2 จุดกับลูกตุ้มคู่ที่มีแท่ง 2 แท่ง
  • Meirovitch, Leonard (1986). Elements of Vibration Analysis (พิมพ์ครั้งที่ 2). McGraw-Hill Science/Engineering/Math. ISBN 0-07-041342-8-
  • Eric W. Weisstein, Double pendulum (2005), ScienceWorld (มีรายละเอียดของสมการเชิงซ้อนที่เกี่ยวข้อง)และ "Double Pendulum" โดย Rob Morris, Wolfram Demonstrations Project , 2007 (ภาพเคลื่อนไหวของสมการเหล่านั้น)
  • Peter Lynch , Double Pendulum , (2001). (การจำลองแอพเพล็ต Java)
  • มหาวิทยาลัยนอร์ทเวสเทิร์น, ลูกตุ้มคู่ เก็บถาวร 2007-06-03 ที่ เวย์ แบ็กแมชชีน (การจำลองแอพเพล็ต Java)
  • กลุ่มฟิสิกส์ดาราศาสตร์พลังงานสูงเชิงทฤษฎีที่ UBC ลูกตุ้มคู่ (2005)
  • แอนิเมชั่นและคำอธิบายของลูกตุ้มคู่และลูกตุ้มคู่ทางกายภาพ (แผ่นสี่เหลี่ยมสองแผ่น) โดย Mike Wheatland (มหาวิทยาลัยซิดนีย์)
  • การจำลอง JavaScript ของฟิสิกส์โอเพ่นซอร์สแบบโต้ตอบพร้อมสมการโดยละเอียดของลูกตุ้มคู่
  • การจำลองแบบโต้ตอบของ Javascript ของลูกตุ้มคู่
  • การจำลองฟิสิกส์ลูกตุ้มคู่จาก www.myphysicslab.com โดยใช้โค้ดโอเพ่นซอร์ส JavaScript
  • การจำลองสมการและคำอธิบายของลูกตุ้มของรอทท์
  • วิดีโอเปรียบเทียบลูกตุ้มคู่ที่มีเงื่อนไขเริ่มต้นเดียวกันบนYouTube
  • Double Pendulum Simulator - โปรแกรมจำลองโอเพ่นซอร์สที่เขียนด้วยC++โดยใช้ ชุด เครื่องมือQt
  • เครื่องจำลอง Java ออนไลน์ เก็บถาวร 2022-08-16 ที่ เวย์แบ็กแมชชีนของนิทรรศการจินตนาการ .
Retrieved from "https://en.wikipedia.org/w/index.php?title=Double_pendulum&oldid=1251984686"