陳彥 陳豐 邵紅才
(揚州市職業(yè)大學(xué)土木工程學(xué)院,江蘇揚州225000)
平衡微分方程和幾何方程是彈性力學(xué)課程中非常重要的內(nèi)容。教材[1]中平面問題極坐標(biāo)中的平衡微分方程是先取徑向和環(huán)向的微元體,然后按照靜力學(xué)的平衡條件導(dǎo)出;幾何方程是先取徑向和環(huán)向的微線段,然后按照線應(yīng)變和切應(yīng)變的定義導(dǎo)出。推導(dǎo)方法雖然和直角坐標(biāo)系一樣,但極坐標(biāo)中環(huán)向的坐標(biāo)線是曲線,導(dǎo)致推導(dǎo)過程的難度變大,容易出錯,最后方程的項增多,形式也不如直角坐標(biāo)中的方程那樣對稱,這種方法可以歸結(jié)為幾何定義法。事實上,對于極坐標(biāo)中的平衡微分方程和幾何方程,可以從直角坐標(biāo)中的方程直接導(dǎo)出,不需要作圖。根據(jù)直角坐標(biāo)和極坐標(biāo)中偏微分算子的轉(zhuǎn)換式、應(yīng)力和應(yīng)變分量的坐標(biāo)變換式以及基矢量變換矩陣,提出了推導(dǎo)極坐標(biāo)中平衡微分方程和幾何方程的解析法,將直角坐標(biāo)中的方程直接變換成極坐標(biāo)中的形式。
《彈性力學(xué)》[1]中平面直角坐標(biāo)系中的平衡微分方程為
《高等數(shù)學(xué)》[2]中直角坐標(biāo)系和極坐標(biāo)系中偏微分算子之間的轉(zhuǎn)換關(guān)系為
《彈性力學(xué)》[1]中極坐標(biāo)系中的應(yīng)力分量σρ,σφ,τρφ和直角坐標(biāo)系中的應(yīng)力分量σx,σy,τxy之間滿足坐標(biāo)變換式
體力矢量f在極坐標(biāo)系中的徑向分量fρ,環(huán)向分量fφ和直角坐標(biāo)系中的分量fx,fy如圖1所示,
圖1 體力分量坐標(biāo)變換
滿足M為直角坐標(biāo)系變換到極坐標(biāo)系的基矢量變換矩陣[3]。
將應(yīng)力分量坐標(biāo)變換式(3)中的σx和τxy分別代入偏微分算子轉(zhuǎn)換式(2)得到
將式(5)和式(6)以及式(4)中的fx代入平衡微分方程(1)的第一式,并化簡得
同樣的方法變換平衡微分方程(1)的第二式,化簡得
對比方程(7)和方程(8),設(shè)法消去sinφ和cosφ項。方程(7)×cosφ+方程(8)×sinφ得
方程(7)×(-sinφ)+方程(8)×cosφ得
方程(9)和方程(10)即為極坐標(biāo)中的平衡微分方程。
《彈性力學(xué)》[1]中平面直角坐標(biāo)系中的幾何方程
《彈性力學(xué)》[1]中雖然沒有給出兩坐標(biāo)系中應(yīng)變分量的坐標(biāo)變換式,但可以通過應(yīng)力分量的坐標(biāo)變換式和平面問題的物理方程導(dǎo)出,這里推導(dǎo)過程省略,
結(jié)論為
位移矢量u在極坐標(biāo)系中的徑向分量uρ,環(huán)向分量uφ和直角坐標(biāo)系中的分量ux,uy仍然滿足
變換幾何方程(11)中的第一式,用應(yīng)變分量坐標(biāo)變換式(12)中的εx代入左邊,將式(13)中的ux代入偏微分算子轉(zhuǎn)換式(2)的第一式,計算出?ux/?x然后代入右邊,得
同樣的方法計算出?uy/?y,?uy/?x和?ux/?y,這樣幾何方程(11)的第二式變換為
幾何方程(11)的第三式變換為
將方程(14)和方程(15)相加得
將方程(14)減去方程(15)得
對比方程(16)和方程(18),設(shè)法消去sin 2φ和cos 2φ項。方程(16)×sin 2φ+方程(18)×cos 2φ得
將方程(16)×cos 2φ+方程(18)×(-sin 2φ)得
方程(20)即為極坐標(biāo)中幾何方程第三式,由方程(17)和方程(19)得
方程(21)即為極坐標(biāo)中幾何方程前兩式。
在將直角坐標(biāo)系中的偏微分方程等價變換到極坐標(biāo)系中的方程時,出現(xiàn)了三角函數(shù)項,如果令φ=0,此時直角坐標(biāo)中的x軸和y軸分別與極坐標(biāo)中的ρ軸和φ軸方向重合,雖然能很快得出結(jié)論,但這只是一個特例。事實上,通過方程之間的消元法可以看到,不管φ的取值如何,對應(yīng)的三角函數(shù)項只是一個參數(shù),是可以消去的,這樣極坐標(biāo)中的平衡微分方程和幾何方程確實不含三角函數(shù)項,使得推導(dǎo)過程更加嚴謹,也更具有一般性。此方法與傳統(tǒng)的幾何定義法相比有異曲同工之處,在教學(xué)中可以對比使用。