趙　振　劉才山　魯建東

1. 北京航空航天大學(xué)航空科學(xué)與工程學(xué)院， 北京 100191； 2. 北京大學(xué)工學(xué)院， 北京 100871； 3. 北京印刷學(xué)院， 北京 102600；? E-mail: bhzhaozhen@buaa.edu.cn

?

空間物體點(diǎn)接觸純滾動(dòng)的幾何意義

趙振1，?劉才山2魯建東3

1. 北京航空航天大學(xué)航空科學(xué)與工程學(xué)院， 北京 100191； 2. 北京大學(xué)工學(xué)院， 北京 100871； 3. 北京印刷學(xué)院， 北京 102600；? E-mail: bhzhaozhen@buaa.edu.cn

空間物體間點(diǎn)接觸純滾動(dòng)的相互作用一般包含非完整約束， 而約束所限制的虛位移通常采用速度水平的 Appell-Chetaev 條件給出， 因此點(diǎn)接觸純滾動(dòng)約束對(duì)應(yīng)的幾何意義并不直觀。作者從多體系統(tǒng)中兩物體沿其輪廓面做點(diǎn)接觸純滾動(dòng)的問(wèn)題出發(fā)， 探討此類非完整約束對(duì)應(yīng)的幾何意義。首先， 提出兩物體保持點(diǎn)接觸的充分必要條件， 并以球-面系統(tǒng)為例推導(dǎo)接觸時(shí)的約束方程。然后， 由空間物體點(diǎn)接觸純滾動(dòng)的幾何和速度約束， 推導(dǎo)此時(shí)滿足的兩種幾何限制條件。結(jié)果表明， 采用兩種幾何條件獲得的虛位移與速度約束的Appell-Chetaev 條件相同。因此， 可以認(rèn)為保持點(diǎn)接觸純滾動(dòng)的空間兩物體在位形空間受到兩種幾何條件的約束限制。

點(diǎn)接觸； 約束方程； 非完整約束； Appell-Chetaev條件

北京大學(xué)學(xué)報(bào)（自然科學(xué)版）第52卷第4期2016年7月

Acta Scientiarum Naturalium Universitatis Pekinensis， Vol. 52， No. 4 （July 2016）

的作用， 并且此約束的虛位移滿足下式［1］:

則式（2）稱為Appell-Chetaev條件。經(jīng)典分析力學(xué)指出， 如果式（1）為完整約束， 則式（2）的幾何意義為約束積分曲面時(shí)間凝固后的任一無(wú)窮小位移。但是， 如果式（1）為非完整的， 則式（2）的幾何意義就很不直觀。本文從空間物體點(diǎn)接觸純滾動(dòng)的問(wèn)題入手， 討論此時(shí)非完整約束虛位移（式（2））的幾何意義。

1　球在水平面上的純滾動(dòng)

半徑為R的球在粗糙水平面上做純滾動(dòng)， 在水平面建立慣性直角坐標(biāo)系Oxyz， 球中心在慣性直角坐標(biāo)系下的分量表示為歐拉角為粗糙水平面的位置0z=， 球在粗糙水平面上做純滾動(dòng)的約束方程可寫(xiě)為其中rτv為接觸點(diǎn)的速度在水平面上的分量。式（3）的第1個(gè)等式為完整約束， 而第2個(gè)等式一般為兩個(gè)非完整約束。對(duì)于非完整約束， 其虛位移可由Appell-Chetaev條件（式（2））給出。通過(guò)拉格朗日乘子把約束嵌入到動(dòng)力學(xué)方程中。

受到完整約束作用的系統(tǒng)只能沿著完整約束規(guī)定的幾何限制運(yùn)動(dòng)， 而受非完整約束作用的系統(tǒng)，其位形空間似乎并不受這些約束的任何限制。以球-面系統(tǒng)為例， 只要外力允許， 平面上的任何點(diǎn)都是可達(dá)的。但是根據(jù)Appell-Chetaev條件， 系統(tǒng)在每一位置的位形空間似乎不是完全自由的， 因?yàn)閺V義坐標(biāo)虛位移并不是完全獨(dú)立的。

本文從點(diǎn)接觸問(wèn)題入手， 探討一般情況下空間物體點(diǎn)接觸和純滾動(dòng)時(shí)約束（式（3））的由來(lái)。我們認(rèn)為， 非完整約束的最初定義是一些力學(xué)家對(duì)某類問(wèn)題（如雪橇問(wèn)題以及點(diǎn)接觸的純滾動(dòng)問(wèn)題）直觀的總結(jié)， 他們直接把約束形式定義在速度水平上， 從而掩蓋了形成非完整約束所需的幾何位形條件。本文認(rèn)為， 對(duì)于點(diǎn)接觸的非完整約束由位形空間上兩個(gè)完整的約束組成， 這兩個(gè)完整約束在速度水平表達(dá)相同， 同時(shí)， 這兩種完整約束限制的虛位移正好滿足Appell-Chetaev條件。

2　保持點(diǎn)接觸的幾何條件

兩個(gè)物體形成并保持點(diǎn)接觸運(yùn)動(dòng)需要滿足以下兩個(gè)條件。

C1: 兩接觸點(diǎn)具有相同的空間位置（或它們的距離一直為零）。

C2: 兩物體的外表面在接觸點(diǎn)相切。

只要兩物體保持點(diǎn)接觸， 除點(diǎn)與線或者點(diǎn)與面接觸時(shí)C2 不存在外， 上述兩個(gè)條件都需成立。因此， 對(duì)于具有一般外形的兩物體的點(diǎn)接觸問(wèn)題， 可從上述兩個(gè)條件出發(fā)， 推導(dǎo)此時(shí)系統(tǒng)的約束方程［2-3］。下面以球在固定水平面上運(yùn)動(dòng)為例， 從上述兩個(gè)條件獲得約束（式（3））的第一個(gè)方程q3=R。

與此類似， 利用兩個(gè)表面參數(shù)

如果x和y表示球和平面上的兩個(gè)接觸點(diǎn)， 根據(jù)點(diǎn)接觸條件C1， 它們具有相同的空間位置:

將式（4）和（5）代入式（6）中， 得到

根據(jù)點(diǎn)接觸條件 C2， 球與水平面在接觸點(diǎn)相切:

由式（4）可得到接觸點(diǎn)的切矢量:

將式（9）代入式（8）， 并考慮球-面系統(tǒng)的特殊情況:

或者即

為了消去式（7）和（11）中的參數(shù)， 把最終的約束方程只表示為系統(tǒng)的廣義坐標(biāo)和時(shí)間的函數(shù)形式，則將式（11）代入式（7）的第3個(gè)方程， 得到

式（13）就是球與水平面接觸的約束方程， 式中的正負(fù)號(hào)以及式（12）中k和n的取值由初始時(shí)刻的接觸點(diǎn)決定。

參數(shù)方程（式（6）和（8））由點(diǎn)接觸保持的兩個(gè)條件C1 和 C2 給出。對(duì)于具有一般形狀物體的接觸問(wèn)題， 得出的參數(shù)方程很難像球-面系統(tǒng)一樣能夠消去接觸面的參數(shù)， 得到如式（13）所示的只含系統(tǒng)廣義坐標(biāo)和時(shí)間的約束方程。這時(shí)， 我們可以采用數(shù)值的微分方法， 解出式（6）和（8）決定的接觸點(diǎn)和接觸約束方程。數(shù)值方法需要初始接觸點(diǎn)的參數(shù)位置， 能夠唯一確定如方程（12）和（13）中的待定符號(hào)［2］。

3　保持點(diǎn)接觸純滾動(dòng)的幾何條件

為兩物體表面輪廓曲紋坐標(biāo)表示的接觸點(diǎn)對(duì)， 根據(jù)接觸點(diǎn)局部曲紋坐標(biāo)到整體慣性坐標(biāo)的關(guān)系， 兩個(gè)接觸點(diǎn)在慣性坐標(biāo)性下的坐標(biāo)可以表示為和。兩物體保持點(diǎn)接觸意味著它們滿足C1的幾何條件， 即兩點(diǎn)在空間重合:

式（14）兩邊對(duì)時(shí)間求導(dǎo)數(shù)（接觸點(diǎn)也隨系統(tǒng)位形或時(shí)間的變化而變化）， 得到

其中vr定義為接觸點(diǎn)的相對(duì)速度:

兩物體除保持點(diǎn)接觸外， 還要滿足接觸點(diǎn)切向相對(duì)速度為零的純滾動(dòng)條件（對(duì)應(yīng)球-面系統(tǒng)的方程（3）中的第2個(gè)方程）:

這里定義n為接觸點(diǎn)處公法線方向的單位矢量。

因?yàn)?n 被設(shè)為接觸點(diǎn)處公法線方向的單位矢量， 根據(jù)點(diǎn)接觸條件C2， 即兩物體在接觸點(diǎn)相切，可得因此， 式（15）兩邊同時(shí)乘以公法線單位矢量n， 可得

如果兩物體保持點(diǎn)接觸并純滾動(dòng)， 則式（17）和（18）意味著接觸點(diǎn)的相對(duì)速度為零:

將式（19）代入（15）中， 得到

只要兩物體保持點(diǎn)接觸純滾動(dòng)， 則式（20）成立， 意味著

其中1s和2s為兩個(gè)接觸點(diǎn)在兩物體表面形成跡線的弧長(zhǎng)， 于是可用和表示兩個(gè)接觸點(diǎn)。因此， 從式（21）可以得到點(diǎn)接觸純滾動(dòng)的幾何意義。

1） 兩接觸物體繞著接觸點(diǎn)轉(zhuǎn)動(dòng):

2） 接觸點(diǎn)在物體表面形成的跡線弧長(zhǎng)相等并相切:

式（23）的第2個(gè)等式雖然是空間矢量方程， 但定義在公切平面上， 因此只限制了一個(gè)自由度。

點(diǎn)接觸純滾動(dòng)的這兩個(gè)幾何限制條件可以總結(jié)為: 點(diǎn)接觸純滾動(dòng)時(shí)， 兩物體滾過(guò)的跡線長(zhǎng)度相同（或?yàn)榱悖?并且方向相切（如果存在）。

下面研究空間兩物體點(diǎn)接觸時(shí)， 具有式（22）和（23）幾何約束特征系統(tǒng)的虛位移。由于系統(tǒng)保持點(diǎn)接觸， 即 x - y=0， 所以系統(tǒng)廣義坐標(biāo)虛位移滿足此方程， 這些廣義坐標(biāo)虛位移也會(huì)引起接觸點(diǎn)的虛位移， 用接觸點(diǎn)參數(shù)表示為因此得到

由于兩物體保持點(diǎn)接觸純滾動(dòng)， 所以接觸點(diǎn)的虛位移滿足式（22）或（23）的幾何限制條件:

把式（25）代入（24）， 具有點(diǎn)接觸純滾動(dòng)的兩物體的系統(tǒng)在位形空間上應(yīng)滿足的虛位移限制方程為

這與利用 Appell-Chetaev 條件（式（2））從速度水平約束（式（19）， 同時(shí)參考式（16）獲得的虛位移相同。因此可以認(rèn)為， 做純滾動(dòng)的兩個(gè)物體在位形空間上滿足式（22）或（23）的幾何限制條件。

4　結(jié)論

對(duì)于完整約束， Appell-Chetaev 條件具有明確的幾何意義。但是對(duì)于非完整約束， 此條件的幾何意義并不直觀?？臻g兩物體保持點(diǎn)接觸純滾動(dòng)時(shí)，其速度水平的約束方程是兩物體接觸點(diǎn)的相對(duì)速度為零， 其中接觸切向相對(duì)速度為零的約束方程一般為非完整的。本文通過(guò)推導(dǎo)發(fā)現(xiàn)， 空間兩物體點(diǎn)接觸純滾動(dòng)要滿足兩種幾何條件的限制: 兩接觸物體滾過(guò)跡線的長(zhǎng)度相同（或?yàn)榱悖?并且方向相切（如果存在）。兩種幾何條件限定的虛位移與此時(shí)速度約束方程的 Appell-Chetaev 條件所限制的虛位移相同。本文研究對(duì)于加強(qiáng)分析力學(xué)中虛位移原理的基本地位具有一定的意義。

［1］ 梅鳳翔， 尚玫. 理論力學(xué)Ⅱ: 專題教程. 北京: 高等教育出版社， 2012

［2］ Zhen Zhao， Liu Caishan. Contact constraints and dynamical equations in Lagrangian systems ［J/OL］. Multibody System dynamics， 2016 ［2016-02-13］. http://dx.doi.org/10.1007/s11044-016-9503-1

［3］ Pfeiffer F. Unilateral problems of dynamics. Archive of Applied Mechanics， 1999， 69（8）: 503-527

［4］ 陳濱. 分析力學(xué). 2版. 北京: 北京大學(xué)出版社，2012

On Nonholonomic Constraints about the Pure Rolling of Point Contact

ZHAO Zhen1，?， LIU Caishan2， LU Jiandong3
1. School of Aeronautic Science and Engineering， Beihang University， Beijing 100191； 2. School of Engineering，Peking University， Beijing 100871； 3. Beijing Institution of Communication， Beijing 102600；? E-mail: bhzhaozhen@buaa.edu.cn

Nonholomonic constraints are involved for 3D point-contact problems. The virtual displacements restrictedby the constraints are usually given by Appell-Chetaev's rule. It has not been very clear of the geometric meaning in configuration space for Appell-Chetaev's rule of nonholonomic constraints. The authors investigate point contact with pure rolling by two rigid bodies in a multibody system to discover its geometric sense. First， the sufficient and necessary conditions of point contact are given. A ball-plane system is presented to demonstrate the validation of the conditions by deducing the system's obvious contact constraint originating from them. Two geometric restrictions for pure rolling are obtained by the nonholonomic constraints of pure rolling as well as the contact constraint in velocity level. It proves that the virtual displacements of the two restrictions are same as those of the constraints of point contact with pure rolling obtained by Appell-Chetaev's rule. So， it is thought that the constraints of pure rolling are constructed by the two geometric restrictions.

point contact； constraint equation； nonholonomic constraint； Appell-Chetaev's rule

O316

10.13209/j.0479-8023.2016.084

國(guó)家自然科學(xué)基金（11572017）、北京市教育委員會(huì)科技計(jì)劃面上項(xiàng)目（KM201310015001）和北京市優(yōu)秀人才培養(yǎng)項(xiàng)目（2012D005004000002）資助

2015-11-23；

2016-02-13； 網(wǎng)絡(luò)出版日期: 2016-07-12