法線和副法線螺旋彈簧勁度系數(shù)的研究

戴璐，臧濤成，葛麗娟

（蘇州科技學(xué)院數(shù)理學(xué)院，江蘇蘇州215009）

法線和副法線螺旋彈簧勁度系數(shù)的研究

戴璐，臧濤成，葛麗娟

（蘇州科技學(xué)院數(shù)理學(xué)院，江蘇蘇州215009）

借助于Cosserat曲線理論定量研究了法線和副法線螺旋彈簧的勁度系數(shù)。計(jì)算結(jié)果表明，在拉伸螺旋彈簧過程中勁度系數(shù)并非始終保持定值。并且副法線螺旋彈簧比法線螺旋彈簧有更好的線性彈性性，所以副法線螺旋彈簧更有利于應(yīng)用在機(jī)械系統(tǒng)。

法線和副法線螺旋彈簧；Cosserat曲線理論；勁度系數(shù)

在大學(xué)物理的教學(xué)研究中彈簧的勁度系數(shù)被認(rèn)為是常數(shù)[1-2]。然而筆者在研究微納米螺旋機(jī)械性質(zhì)的時(shí)候發(fā)現(xiàn)，螺旋結(jié)構(gòu)的胡克常量會(huì)隨著螺旋軸向拉伸而變化[3]。微納米螺旋具有超彈性，其在軸向高負(fù)載的作用下幾乎可以被拉伸成直線，當(dāng)撤去負(fù)載時(shí)，螺旋恢復(fù)原形，所以微納米螺旋可以被當(dāng)作彈簧[4-5]。筆者將借助考慮了拉伸和剪切形變的Cosserat曲線建立一套理論來(lái)定量研究整個(gè)拉伸和拉斷區(qū)域法線和副法線螺旋彈簧的勁度系數(shù)。

1 Cosserat曲線理論

Cosserat曲線這一概念首先在1909年由E·Cosserat和F·Cosserat提出，他們用3個(gè)方向矢量來(lái)定義一條曲線，并用4個(gè)矢量場(chǎng)來(lái)表示該曲線，即Cosserat曲線[6]。Whitman和DeSilva將該理論應(yīng)用于彈性細(xì)桿的研究，從而形成彈性細(xì)桿三維非線性平衡理論，即彈性細(xì)桿的Cosserat曲線理論[7-8]，其平衡方程為

其中τ，m分別為作用于整個(gè)細(xì)桿的力和力矩矢量，ε為輪轉(zhuǎn)算符，W和y分別為位置矢量和方向矢量。下標(biāo)α，β，γ取1，2，3。，拉伸量λ＝?s/?S，S和s分別為沿著固定結(jié)構(gòu)和形變結(jié)構(gòu)的弧長(zhǎng)。

圖1 法線螺旋彈簧結(jié)構(gòu)（a）及副法線螺旋彈簧受到軸向拉力F作用前后的結(jié)構(gòu)（b）、（c）

如圖1（a-c）所示，筆者考慮了兩種橢圓截面的螺旋形彈簧，分別稱為法線螺旋和副法線螺旋[9-10]。橢圓截面的兩個(gè)半軸長(zhǎng)分別為r1和r2（r2＞r1）。假設(shè)一N匝均勻螺旋HI，半徑為a0，螺距為b0，長(zhǎng)為l0。螺旋HI在軸向拉力F作用下被拉伸為螺旋HF，半徑為a，螺距為b，螺旋拉伸量為L(zhǎng)。在模型中HI為Cosserat曲線理論中固定參考結(jié)構(gòu)，而HF為形變結(jié)構(gòu)。Di（i=1，2，3）為螺旋HI的方向基矢，而di為螺旋HF的方向基矢，它們分別由兩組歐拉角φ0，θ0，ψ0和 φ，θ，ψ定義。對(duì)于螺旋結(jié)構(gòu)都為常數(shù)，不隨弧長(zhǎng)變化[11]。選擇HI的第三個(gè)方向矢量D3沿著螺旋線中心線的切向，拉力F沿著固定笛卡爾坐標(biāo)系的e3軸。D1、d1和D2、d2分別沿著彎曲剛度最大和最小的方向。

法線螺旋和副法線螺旋HI的方向變化矢量W（0）有以下形式

將螺旋HF的方向變化矢量代入平衡方程（2）得到法線螺旋的方程

和副法線螺旋的方程

其中E1=E2=KGtw，E3=Etw，A=EI1，B=EI2和C=4GI1I2/（I1+I2）[12]。K為Timoshenko切系數(shù)，可用泊松比ν來(lái)表示[12]。E和G=E/2（1+ν）分別是楊氏模量和切模量。I1=（πr23r1）/4和I2=（πr2r13）/4是截面的轉(zhuǎn)動(dòng)慣量。

因?yàn)樨?fù)載力決定了螺旋HF的形狀，所以方程（6）和（8）表示，對(duì)于法線螺旋來(lái)說(shuō)方向變化矢量W有如下形式

對(duì)于副法線來(lái)說(shuō)W有如下形式

以及副法線螺旋的變化矢量

根據(jù)方程（11）、（12），不僅可以得到拉伸量λ

還能得到法線和副法線螺旋HI和HF的半徑和螺旋角，用歐拉角表示為

其中Δ≡I2/I1，對(duì)于法線（或者副法線）螺旋來(lái)說(shuō)i=1（或者i=2），δi2為克羅內(nèi)克函數(shù)。并且力矩M沿著螺旋軸向，即負(fù)載力F方向的表達(dá)式為

根據(jù)方程（14）、（15）、（16）和胡克定律k=dF/dL，螺旋HF的勁度系數(shù)可表示為

通過測(cè)量螺旋HI的幾何參數(shù)a0，b0，r1，r2，N和負(fù)載或者力矩，就可以根據(jù)方程（13-18）以及螺旋線長(zhǎng)度l=λl0和體積守恒條件分別是螺旋線形變前后的長(zhǎng)度（如圖1（a）和1（b）））來(lái)計(jì)算螺旋HF的半徑a，螺距b和勁度系數(shù)k。

2 討論與結(jié)果

筆者將基于Cosserat曲線模型，利用Si3N4微米螺旋來(lái)定量研究在整個(gè)軸向拉伸過程中法線和副法線螺旋的勁度系數(shù)。Si3N4微米螺旋具有圓截面，其半徑為r0=2.58 μm[5]。在計(jì)算中假設(shè)螺旋線截面面積為定值πr02，然后令長(zhǎng)短半軸比為2來(lái)得到橢圓截面螺旋。Si3N4微米螺旋的其他幾何參數(shù)和材料參數(shù)為a0=80 μm，b0=200 μm，N=8，E=240 GPa[5]，ν=0.25[13]，其中幾何參量來(lái)自于SEM圖像。首先要確定螺旋能承受的最大負(fù)載，即拉斷螺旋的力Fbreak。根據(jù)楊氏模量和勁度系數(shù)的定義，在拉斷點(diǎn)有k=kS=EA/l0，其中kS為螺旋線彈性常數(shù)，A為拉伸后螺旋線截面面積（A=πr1r2/λ）。結(jié)合方程（14）-（18），得到Si3N4微米法線和副法線螺旋能承受的最大負(fù)載分別為Fbreak=120×10-3N和Fbreak=51×10-3N。

在已知材料的Fbreak后，根據(jù)方程（18）圖2（a）和2（b）分別給出了整個(gè)拉伸區(qū)域Si3N4法線和副法線微米螺旋的勁度系數(shù)。勁度系數(shù)—拉伸曲線可分為3個(gè)區(qū)域：區(qū)域I為螺旋線性彈性區(qū)域。在該區(qū)域內(nèi)，即負(fù)載初始階段，勁度系數(shù)幾乎不隨伸長(zhǎng)變化，法線螺旋平均值為0.39 N·m-1，副法線螺旋的平均值為0.28 N·m-1。區(qū)域II為螺旋非線性彈性區(qū)域。在該區(qū)域中繼續(xù)拉伸螺旋，螺旋被拉緊，勁度系數(shù)呈非線性增長(zhǎng)到高應(yīng)力區(qū)域：勁度系數(shù)隨拉伸迅速增加。區(qū)域III為螺旋線線性彈性區(qū)域，在該區(qū)域螺旋被更大負(fù)載力拉成一根緊繃直線直到拉斷，其勁度系數(shù)進(jìn)入一飽和值，即螺旋線彈性常數(shù)kS。

通過比較圖2（a）和2（b）可以得到法線和副法線螺旋的線性彈性區(qū)域分別為大約1.2 mm和2.1 mm處，這兩個(gè)螺旋拉伸量分別為螺旋原長(zhǎng)的600%和1 050%。顯然副法線螺旋的線性彈性區(qū)域更大，這表明副法線螺旋的線性彈性性更好，也就是機(jī)械穩(wěn)定性強(qiáng)。并且在線性彈性區(qū)域內(nèi)副法線螺旋的勁度系數(shù)比法線螺旋的勁度系數(shù)小，在相等的負(fù)載力作用下，副法線螺旋會(huì)有更大的形變量，即彈性更好。所以副法線螺旋更有利于應(yīng)用在機(jī)械系統(tǒng)。

3 結(jié)語(yǔ)

基于Cosserat曲線理論定量研究了大學(xué)物理中彈簧的勁度系數(shù)。以橢圓截面法線和副法線螺旋彈簧為研究對(duì)象，推導(dǎo)出以螺旋幾何參數(shù)、材料參數(shù)和負(fù)載力為自變量的勁度系數(shù)表達(dá)式，并利用Si3N4微米螺旋彈簧定量解釋了整個(gè)拉伸和拉斷區(qū)域的勁度系數(shù)。研究結(jié)果發(fā)現(xiàn)，在拉伸過程中螺旋彈簧的勁度系數(shù)會(huì)發(fā)生變化，而副法線螺旋彈簧比法線螺旋彈簧有更長(zhǎng)的線性彈性性區(qū)域，因此，更有利于應(yīng)用。

[1]馬文蔚，周雨青.物理學(xué)教程（上冊(cè)）[M].2版.北京：高等教育出版社，1999.

[2]程守洙，江之永.普通物理學(xué)（上冊(cè)）[M].6版.北京：高等教育出版社，2006.

[3]Dai L，Shen W Z.Cosserat curve model for superelasticity of helices[J].Nanotechnology，2009，20（46）：465707（1-6）.

[4]Chen X，Zhang S，Dikin D A，et al.Mechanics of a carbon nanocoil[J].Nano Lett，2003，3（9）：1299-1304.

[5]Cao C B，Du H L，Xu Y J，et al.Superelastic and spring properties of Si3N4microcoils[J].Adv Mater，2008，20（9）：1738-1743.

[6]Cosserat E，Cosserat F.Theorie des Corps Deformables[M].Paris：Hermann，1909.

[7]Whitman A B，DeSilva C N.Stability in a linear theory of elastic rods[J].Acta Mechanica，1972，15（3、4）：295-308.

[8]Whitman A B，DeSilva C N.An exact solution in a nonlinear theory of rods[J].J Elasticity，1974，4（4）：265-280.

[9]Fonseca A F，Malta C P，Galvao D S.Mechanical properties of amorphous nanosprings[J].Nanotechnology，2006，17（22）：5620-5626.

[10]Goriely A，Shipman P.Dynamics of helical strips[J].Phys Rev E，2000，61（4）：4508-4517.

[11]Zhou Z，Lai P Y，Joós B.Elasticity and stability of a helical filament[J].Phys Rev E，2005，71（5）：052801（1-4）.

[12]Love A E H.A Treatise on the Mathematical Theory of Elasticity[M].New York：Dover，1944.

[13]Fate W A.High-temperature elastic moduli of polycrystalline silicon nitride[J].J Appl Phys，1975，46（6）：2375-2377.

Stiffness coefficients of normal and binormal helical springs

DAI Lu，ZANG Taocheng，GE Lijuan
（School of Mathematics and Physics，SUST，Suzhou 215009，China）

With the aid of Cosserat curve theory，we investigated the stiffness coefficients of normal and binormal helical springs.The results show that the stiffness coefficients are not immune to the elongation in the entire stretching process.The binormal helical springs are quantitatively confirmed to excel over the normal ones in linear elasticity.Therefore，the binormal helical springs are more beneficial to the applications in engineering.

normal and binormal helical springs；Cosserat curve theory；stiffness coefficients

O469

A

1672－0687（2015）02－0037－04

責(zé)任編輯：李文杰

2014－10－04

國(guó)家自然科學(xué)基金資助項(xiàng)目（11347136）；江蘇省自然科學(xué)基金資助項(xiàng)目（BK20130265）

戴璐（1981-），女，江蘇蘇州人，講師，博士，研究方向：納米材料。