范祥彬，盧奕合，劉逍娜，王鐵成

(山西大學(xué) 物理電子工程學(xué)院，山西 太原 030006)

2022年CUPT第7題：硬幣落地時側(cè)面站立的情況十分罕見.為了使圓柱形骰子落下時立在側(cè)面、正面與反面概率相同，應(yīng)該具有怎樣的物理性質(zhì)和幾何特征？

目前已經(jīng)有研究者對三面骰子的問題展開過詳細(xì)的描述[1]，但是仍然存在一些細(xì)節(jié)沒有得到解答.例如，在投擲三面骰子時，如果給予它水平方向的速度，側(cè)面出現(xiàn)的概率會發(fā)生什么變化？ 三面骰子反復(fù)彈跳是否會有不同的結(jié)果？基于以上問題，本文進(jìn)行了相關(guān)研究和討論.

首先，把硬幣定義為圓柱體，當(dāng)它充分翻轉(zhuǎn)落地時，雖然側(cè)面出現(xiàn)的概率極小，但卻不能忽略其可能性.硬幣除了有1個平移自由度外，還有2個旋狀自由度，但是在實際情況下，由于硬幣的對稱性(側(cè)面高度趨于零)，平移自由度作用并不明顯，因此可以認(rèn)為理想硬幣落在正面或者反面的概率近似為1/2，但是一旦改變其側(cè)面的厚度，立在側(cè)面的概率就會改變，正面或者反面的概率將不再是1/2.

對于硬幣、正方體骰子和陀螺骰子而言，立在不同面的概率可以通過物體的對稱性來預(yù)測.但當(dāng)物體不再具有精確的對稱性時，立在不同面的概率就不可能只受幾何特征的影響，還需要考慮彈跳的影響.

1 三面骰子的定義

三面骰子的立體圖及截面圖如圖1所示，設(shè)圓柱半徑R與高度h的比為

(a)立體圖

(1)

則內(nèi)角為

β=2arctan(2η),

(2)

(3)

2 骰子的投擲方式

研究者們大都假設(shè)三面骰子豎直下落，而不具有水平方向的速度，下面給出簡單的理論解釋.

假設(shè)在t=t1和t=t2時刻，系統(tǒng)的位置由2組坐標(biāo)q(1)和q(2)確定，則系統(tǒng)在這2個位置之間的運(yùn)動積分

(4)

(5)

故

(6)

函數(shù)L為給定系統(tǒng)的拉格朗日函數(shù)或拉格朗日量，積分S為作用量.

當(dāng)僅考慮三面骰子落地時的狀態(tài)(落地前的運(yùn)動行為不考慮)，假設(shè)骰子在地面沿直線滾動.現(xiàn)以三面骰子向前滾動時的2種極端的運(yùn)動方式為例，即側(cè)面滾動和正反面輪流翻滾(圖2).利用式(1)求出這2種運(yùn)動方式的η值，其他運(yùn)動方式的η值介于這兩者之間.

(a)側(cè)面滾動 (b)正反面輪流翻滾

2.1 側(cè)面滾動

(7)

(8)

拉格朗日量

(9)

根據(jù)拉格朗日方程，則有

(10)

(11)

積分可得

x=C1t+C2.

(12)

2.2 正反面輪流翻滾

如圖2(b)所示，骰子勢能V不再是常量，并其轉(zhuǎn)動慣量為

(13)

故動能可表示為

(14)

拉格朗日量為

(15)

根據(jù)拉格朗日方程可得：

(16)

(17)

(18)

對比圖2(a)和圖2(b)，可以看到在圖2(b)情況下，骰子會受到額外的耗散力影響，因此圖2(a)消耗能量較少，圖2(b)消耗能量較多，其他運(yùn)動狀態(tài)則介于二者之間.由于圓形物體滾動時，其重心高度相對于地面接觸點不變,滾動時不需要克服重力做功，處于能量最低態(tài)，是最穩(wěn)定的運(yùn)動方式.因此骰子總是趨向于側(cè)面著地.

在現(xiàn)實情況中，首先骰子會受到地面的影響，使得骰子滾動時的軌跡為曲線，但曲線軌跡并不會影響結(jié)論；其次，滾動摩擦力比最大靜摩擦力小得多，側(cè)面滾動主要受到滾動摩擦，正反面翻滾主要受到最大靜摩擦，即前者相比后者更能保證運(yùn)動最大量，可以很大程度上減小克服摩擦力做功.因此，即使考慮摩擦影響，骰子也是趨向側(cè)面滾動狀態(tài).所以可得出以下結(jié)論：三面骰子如果具有水平方向的初速度，骰子總是趨向于側(cè)面著地向前滾動.性質(zhì)會影響骰子立地概率的公平性，從而導(dǎo)致實驗結(jié)果與實際情況不符.

那么，是否可以得出具有一定厚度的三面骰子都是趨向于側(cè)面著地的結(jié)論呢？對此，對不同尺寸的三面骰子的著地情況進(jìn)行統(tǒng)計，結(jié)果如表1所示.

表1 不同尺寸的三面骰子著地的情況

實驗采用的骰子直徑2R=2.500 cm，h=0.833 cm，如若不給初速度，正常投擲，大都會認(rèn)為骰子正反面著地的可能性大，然而表1卻是側(cè)面落地的次數(shù)多，表明骰子向前投擲會增加側(cè)面著地的概率，與理論猜想相同.

同樣利用硬幣進(jìn)行翻滾實驗時，大多數(shù)情況下也都是側(cè)面向前滾動的情況，但由于硬幣的側(cè)面尺度很小，立地條件不足，容易受到外界擾動的影響，所以不論怎么拋擲都不會影響公平性.

如果地面光滑，采用向前拋擲的方法，硬幣便不會受到外界擾動，將保持自身的慣性，從而有可能出現(xiàn)側(cè)面立地的情況.因此，在投擲圓柱形骰子時，盡量保持骰子在水平方向無初速度.

3 三面骰子實驗設(shè)計

圖3所示骰子是以ABS材質(zhì)為原材料的實心圓棒，利用游標(biāo)卡尺測量長度，通過手工鋸和銼刀切割磨皮，雖以上制作方法存在較大誤差(手工切割所帶來的正反兩面的不對稱)，但由于目前只是驗證模型，即只需證明側(cè)面出現(xiàn)的概率接近1/3，因此可以忽略誤差帶來的影響.骰子直徑固定2.500 cm，制作不同高度的骰子，實驗時1次性投擲5個相同高度的骰子，多個數(shù)據(jù)點可以避免單一骰子帶來的固有偏差[3].

圖3 手工骰子

第2章已證明三面骰子水平方向不應(yīng)具有初速度，為了避免這種情況，設(shè)計了如圖4所示的投擲裝置.

圖4 投擲裝置示意圖

將骰子從寬口丟下，讓骰子充分旋轉(zhuǎn)，在骰子接近窄口時，由于路口變窄，骰子劇烈震蕩，骰子的轉(zhuǎn)動動能部分損失，并轉(zhuǎn)化為豎直向下的動能，這樣就滿足了第2章中論述的限制條件.改變投擲裝置的高度，可改變骰子的彈跳情況.投擲裝置采用5 L礦泉水桶裁去桶底部分制作.

當(dāng)三面骰子豎直落地時，由于邊緣落地帶來的反彈，能量在豎直方向可能無法抵消，從而讓骰子具有水平方向的動能(在同一豎直線反復(fù)彈跳的結(jié)果不太可能出現(xiàn))，那么當(dāng)出現(xiàn)該情況時，需要舍棄該結(jié)果并重新投擲.此外，設(shè)計骰子分別在水泥和PVC材質(zhì)的接觸面上著地.

4 三面骰子模型

在完全不考慮反彈的情況下提出了3種三面骰子的模型，出發(fā)角度各不相同，得出的結(jié)果也有差異.為了驗證前文內(nèi)容，設(shè)置了3種實驗環(huán)境：在無初速度時水泥面和PVC面著地，以及有初速度時，骰子在水泥面著地的情況.

4.1 圓柱體表面積模型

通過借鑒標(biāo)準(zhǔn)正方體骰子出現(xiàn)等概率性的特性，猜想表面積可能是影響因素之一，于是得出：

(19)

表2 圓柱體表面積模型骰子著地情況

4.2 立體角模型

(a)立體圖

設(shè)球冠的表面積為A，利用微元法求解.在球冠上取細(xì)小的圓環(huán)，半徑為r，弧長為dl=Rdθ，于是面積為

dA=2πrdl=2πrRdθ,

(20)

(21)

因r=Rcosθ,則

(22)

(23)

由立體角公式，有

(24)

為使其平均，應(yīng)該使得立體角為

(25)

(26)

表3 立體角模型骰子下落情況

由表3中數(shù)據(jù)得：水泥接觸面時骰子的正面概率PA=0.443,反面概率PB=0.383,側(cè)面概率PS=0.174；PVC接觸面時骰子的正面概率PA=0.427,反面概率PB=0.363,側(cè)面概率PS=0.210.很明顯，側(cè)面概率與理想值0.333有較大偏差.故馮諾依曼假設(shè)僅在數(shù)學(xué)模型上合理，但與實驗結(jié)果不符.

4.3 重心模型

圖6為三面骰子落地時與接觸面接觸時的一種狀態(tài)，重心線的偏向決定了三面骰子倒地的結(jié)果[5]，θ為重心線與x軸的夾角，OH為重心線，決定骰子的偏向.

圖6 圓柱骰子截面圖

表4 重心模型下骰子下落情況

由表4中數(shù)據(jù)可得：接觸面為水泥且無初速的條件下，骰子正面概率PA=0.335,反面概率PB=0.322,側(cè)面概率PS=0.343,側(cè)面相對偏差Er=3.0%；接觸面為PVC且無初速的條件下骰子正面概率PA=0.320,反面概率PB=0.351,側(cè)面概率PS=0.328,側(cè)面相對偏差Er=1.5%.數(shù)據(jù)與理論吻合較好.

5 彈跳模型

上文僅考慮了三面骰子落地瞬間便損失全部動能的情況，現(xiàn)以簡單彈跳情況為例進(jìn)行分析.假設(shè)每次落地?fù)p失一部分動能，即三面骰子經(jīng)過多次反復(fù)彈跳才會停止，并且每次落地彈跳之后會經(jīng)過充分旋轉(zhuǎn)，再次落地時方向依舊隨機(jī).假設(shè)初始時，三面骰子具有能量E，第一次與地面碰撞后能量為γE，第n次碰撞后能量為γnE(n=1,2,3…)，其中γ<1為衰減系數(shù)[6].

定義2個能量值EAB和ES，它們分別表示底面和側(cè)面的能量勢阱，如圖7所示.假設(shè)底面的勢阱高，即ES>EAB(若ES

圖7 ES>EAB情況下對應(yīng)的能量勢阱圖

此時落在2個山谷的概率正比于兩側(cè)山峰的間距，則會出現(xiàn)如下情況：

1)假設(shè)經(jīng)過n次彈跳，能量為EAB<γnE,γn+1E

(27)

該假設(shè)還可以表述為：假設(shè)前m次碰撞，其能量為γmE>ES，而EAB<γm+1E

2)假設(shè)經(jīng)過n次碰撞后能量為γnE

(28)

(29)

(30)

(31)

現(xiàn)考慮式(28)，則

(32)

(33)

該結(jié)果與式(30)和式(31)相同.

當(dāng)ES=EAB時，對應(yīng)如圖8所示的能量勢阱.骰子經(jīng)過n次彈跳后，能量為γnE

圖8 ES=EAB情況下對應(yīng)的能量勢阱圖

(34)

此條件下得到的結(jié)果，與前面結(jié)果相同，但此種情況很特殊，等式成立的條件比不等式成立更加苛刻，因此不易實現(xiàn).

在本節(jié)中，得出結(jié)論有：三面骰子不允許多次彈跳，否則實驗結(jié)果會出現(xiàn)很大偏差.所以實驗過程中要求三面骰子落地后，動能必須立刻損失或彈跳1次后減為零，否則將不會出現(xiàn)等概率事件.

在構(gòu)造彈跳模型時，沒有給出E和r的數(shù)值，因此做了如下假設(shè)：

a.每次碰撞后的能量衰減參量不變，并且各個接觸點的衰減參量都相同；

b.每次彈跳落地前都經(jīng)過充分旋轉(zhuǎn)，落地后再次隨機(jī)；

c.E和γ不同取值不會導(dǎo)致其他影響；

d.從始至終只有骰子的豎直方向有動能，且反彈不會帶來方向的偏移.

上面4個假設(shè)的約束性較強(qiáng)，使得彈跳模型局限性較大.根據(jù)熱力學(xué)模型的建立過程可知，模型的限制越多，模型與實際情況就越可能畸變.三面骰子不像正方形骰子、異體形骰子，每個面都有精準(zhǔn)的對稱性，因此需做出一些條件假設(shè)求解三面骰子模型.

根據(jù)本節(jié)內(nèi)容，n=0的情況與4.3節(jié)相同，現(xiàn)給出n=1,2時的實驗數(shù)據(jù)，如表5～6所示.

表5 n=1,η=時骰子的下落情況

由表5可知：當(dāng)n=1，接觸面為水泥時，骰子的正面概率PA=0.330,反面概率PB=0.349,側(cè)面概率PS=0.321,側(cè)面相對偏差Er=3.6%.

猜想n=0，1情況下，實驗結(jié)果相差不大的原因只是由于碰撞1次帶來的效果不特別明顯，為了滿足骰子只碰撞1次，需要調(diào)節(jié)高度來實現(xiàn)，在實驗中，這種高度不足以產(chǎn)生特別明顯的誤差.圓柱形骰子是不完美對稱的，因此碰撞之后會使骰子的運(yùn)動方向改變.因此骰子具有水平方向的初速度，也會導(dǎo)致側(cè)面的概率發(fā)生變化.

由圖9可以看出，不完美對稱的物體反彈以后，會偏離原來的運(yùn)動方向，由能量守恒可知，骰子具有了水平方向的速度，根據(jù)第2章內(nèi)容，骰子側(cè)面著地的概率將會增加.

(a)高度對稱球的反彈

通過簡單的測試發(fā)現(xiàn)，骰子豎直落下，發(fā)生彈跳,側(cè)面著地的次數(shù)會大大增加(表7)，這樣的結(jié)果與第2節(jié)的內(nèi)容相符.所以給出的結(jié)論是：實驗過程中，保持n=0才最為合理.

表7 發(fā)生彈跳后的骰子下落情況

表6 n=2,η=0.391 2時骰子下落情況

綜上所述幾何和反彈模型共同說明:圓柱形骰子須無反彈、無水平初速度的投擲才能保證其等概率著地.

6 標(biāo)準(zhǔn)骰子實驗

圖10 3D打印骰子

表8 標(biāo)準(zhǔn)骰子的下落情況

接觸面為水泥面時，骰子正面概率PA=0.332,反面概率PB=0.336,側(cè)面概率PS=0.331，側(cè)面相對偏差Er=0.6%.

接觸面為PVC面時，骰子正面概率PA=0.327,反面概率PB=0.342,側(cè)面概率PS=0.330，側(cè)面相對偏差Er=0.9%.

利用Pycharm軟件進(jìn)行模擬如圖11所示，模擬投擲次數(shù)為1 000次，做了10次模擬的數(shù)據(jù)如表9所示.

表9 Pycharm軟件模擬投擲情況

圖11 Pycharm軟件模擬

從實驗數(shù)據(jù)上可以看出，實驗結(jié)果與理論分析較吻合，如果需要更接近333的數(shù)據(jù)，可以通過調(diào)節(jié)滾動摩擦力來實現(xiàn).

7 結(jié) 論