高 藝,郭志榮,侯 娟,史 芹
(石河子大學(xué) 理學(xué)院物理系,新疆 石河子 832003)
國際青年物理學(xué)家競標(biāo)賽(IYPT)第35屆(2022年)賽題[1]第七題“Three-Sided Dice:To land a coin on its side is often associated with the idea of a rare occurrence.What should be the physical and geometrical characteristics of a cylindrical dice so that it has the same probability to land on its side and one of its faces?”(三面骰子:在拋硬幣時,硬幣能夠豎著立起來的情況是罕見的。如果一個圓柱狀骰子側(cè)面立地和兩底面立地的概率相等,那么這個骰子應(yīng)當(dāng)具備什么樣的物理和幾何特征?)。
骰子是具有某種“對稱性”的物體,有質(zhì)量分布和幾何形狀兩個要素。最簡單的是兩要素非耦合的骰子,即質(zhì)量分布均勻的規(guī)則體。另外,投擲過程中受作用(如落地撞擊等)而變形也會影響某種概率。本文不討論變形情況,即質(zhì)量分布均勻剛體圓柱狀骰子。
滿足賽題要求的圓柱狀骰子的幾何特征主要是高H和底面半徑R的關(guān)系決定。聯(lián)想到投擲球狀物體時,球面上任意點立地的概率相等,設(shè)想在柱體外接一球面,如圖1所示,側(cè)面和兩底映射到球面上面積相等時,可確定骰子的幾何特征。不難得到
圖1
(1)
這與文獻[2]用計算機模擬得到結(jié)果有甚微相差,這個差異是由棱先著地,然后翻向底面或側(cè)面的概率不同引起的。
接下來圓柱狀骰子棱上某點觸地會導(dǎo)致概率不均,這也是下文運動特征討論的重點。
棱觸地后,骰子最終哪面立地,是由觸地初態(tài)和環(huán)境(重力、阻尼)作用下的動力學(xué)過程決定,討論簡單,下面在不考慮阻力的情況下,分析柱體狀骰子的穩(wěn)定轉(zhuǎn)動,然后由轉(zhuǎn)速與章動角的函數(shù)單調(diào)性定性地給出因為運動而導(dǎo)致的概率差異。
2.1.1 動力學(xué)方程
如圖2所示,建立固定坐標(biāo)系O-XYZ和本體坐標(biāo)系o-xyz。觸地點為O點,O與o重合,旋轉(zhuǎn)軸作為oz。不計阻力,并且節(jié)線與X軸重合時計時,柱體狀骰子的運動學(xué)方程為[3]。
圖2 繞對角線自旋
(2)
2.1.2 自旋與進動
采用前文的幾何特征,設(shè)柱體繞x,y,z軸的轉(zhuǎn)動慣量分別是Ix,Iy,Iz。柱體平動動能T1和轉(zhuǎn)動動能T2分別為
(3)
(4)
取XOY面為零勢面,系勢能為
V=mgHcosθ.
(5)
則圓柱狀骰子運動的拉格朗日方程[4]為
(6)
其中
f(θ)=(mH2+Iy-Iz)cosθ.
(7)
拉氏函數(shù)中不顯含φ和ψ,由φ和ψ循環(huán)坐標(biāo)可以得到兩個初值決定的積分常數(shù)如下。
(8)
(9)
式(7)可求得進動角速度
(10)
2.1.3 章動
設(shè)體系總能量E(初始條件確定)。定義新的能量和有效勢能[5]并對能量式簡化。
(11)
其中
(12)
(13)
由(11)式知V′≤E′,如圖3所示,可得θ1<θ<θ2,這里
圖3 章動角范圍
(14)
2.1.4 概率分析
(15)
本文討論了國際青年物理學(xué)家競標(biāo)賽(IYPT)第35屆(2022年)賽題的第七題,圓柱形骰子的幾何參數(shù)和運動特征。在均質(zhì)剛體圓柱形骰子的假設(shè)下,本文認為骰子立地概率主要由幾何特征決定,而運動導(dǎo)致概率變化主要是柱體棱邊著地的時候才出現(xiàn),因此在不考慮阻力的情況下,就柱棱觸地時能穩(wěn)定不倒的情況進行了轉(zhuǎn)動運動的分析。結(jié)果顯示,棱點觸地后且要滿足很嚴(yán)格的運動(轉(zhuǎn)動)狀態(tài)才能保持穩(wěn)定不倒。在對穩(wěn)定轉(zhuǎn)動時由角速度與章動角的函數(shù)關(guān)系的單調(diào)性分析,能獲得轉(zhuǎn)動角速度小時章動角的變化趨勢可定性地得出立地面的概率差異。當(dāng)然,更嚴(yán)謹(jǐn)?shù)挠懻搼?yīng)該是直接在建模時就考慮阻力的情況下進行的運動分析,理論分析會復(fù)雜很多。