吳 飛，李德華，莫玉林

（1.貴州大學(xué) 機(jī)械工程學(xué)院，貴州 貴陽(yáng)，550003；2.中鋁長(zhǎng)城建設(shè)有限公司 越南項(xiàng)目部，河南 鄭州，450041；3.西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院，陜西 西安，710048）

1765年，瑞士數(shù)學(xué)家歐拉對(duì)張緊纏繞于圓柱體上線繩的張力與纏繞角，推導(dǎo)了數(shù)學(xué)力學(xué)關(guān)系式，即歐拉公式，從而奠定了撓性體摩擦理論的基礎(chǔ)。在機(jī)械工業(yè)中，帶傳動(dòng)、制動(dòng)器、對(duì)物體的牽引、提升、搬運(yùn)等都會(huì)涉及到柔性體的摩擦理論。經(jīng)典歐拉公式對(duì)小角度纏繞于圓弧面上的線繩張緊力的計(jì)算很方便，若用經(jīng)典歐拉公式計(jì)算接觸面為非圓弧面情況下的張緊力，誤差就很會(huì)大，對(duì)我們的工程應(yīng)用會(huì)產(chǎn)生一定的影響。

本文介紹一種纏繞于一般曲面體上的廣義歐拉公式，運(yùn)用ANSYA WORKBENCH對(duì)其進(jìn)行有限元分析計(jì)算，對(duì)分析的結(jié)果和公式進(jìn)行比較，最后得出結(jié)論。

1 纏繞于一般曲面體上的廣義歐拉公式

設(shè)一根線繩纏繞于一般曲面體上，接觸線為一般曲線。假設(shè)線繩是理想的撓性體，在工作中不能伸長(zhǎng)，也無(wú)彎曲阻力存在，忽略線繩轉(zhuǎn)動(dòng)時(shí)離心力的影響（v<10m/s）。如圖1所示，取一段接觸帶AB進(jìn)行分析，并在AB弧上取一微小段線繩，其所對(duì)應(yīng)的微小圓心角為，在該微弧段處的曲率半徑為R(s)。

圖1 線繩纏繞于一般曲面體

圖2 微段線繩的力平衡關(guān)系

如圖2所示，對(duì)微弧段dl作受力分析，其兩端受力分別為F和F+dF，線繩對(duì)轉(zhuǎn)筒的正壓力為dFN，其摩擦力為μdFN，μ為線繩與轉(zhuǎn)筒間的摩擦系數(shù)。線繩在即將打滑的臨界狀態(tài)下，各個(gè)力在水平和豎直方向的平衡關(guān)系為：

此即為纏繞于一般曲面體上的歐拉公式。

2 有限元計(jì)算

本文通過(guò)對(duì)一具體實(shí)例作有限元計(jì)算，分析卷筒纏繞時(shí)的線繩張緊力的情況，與用上述介紹的廣義歐拉公式分析的結(jié)果進(jìn)行比較。以下采用solidworks建立三維模型，然后導(dǎo)入ANSYS WORKBENCH進(jìn)行有限元分析計(jì)算。

2.1 建模

在solidworks中建立一纏繞于一般曲面體上的卷筒纏繞三維模型，該卷筒分別由6個(gè)圓弧面和平面交替組成。首先創(chuàng)建一個(gè)內(nèi)接圓半徑為200mm的六邊形，再將其6個(gè)角用半徑為100mm的圓弧進(jìn)行圓角，然后通過(guò)實(shí)體拉伸生成模型。該卷筒內(nèi)徑為280mm，卷筒厚250mm，繞繩厚為30mm，寬為100mm，卷筒與繞繩間的摩擦系數(shù)為0.3，建立纏繞包角分別為 60°、120°、180°、240°、300°、360°的三維模型。如圖 3 所示纏繞包角為360°時(shí)的裝配體模型。在solidworks中將所建立的裝配體另存為Parasolid格式，再分別導(dǎo)入到ANSYS WORKBENCH中進(jìn)行有限元分析。

圖3 卷筒纏繞有限元模型

圖4 接觸曲線形狀參數(shù)

2.2 公式分析

根據(jù)模型參數(shù)，這些模型的接觸面為圓弧面和平面的組合面。如圖4所示的接觸曲線，這些接觸曲線由圓弧和直線段交替結(jié)合組成，且各段圓弧和直線的參數(shù)相同。

2.3 模型的有限元分析計(jì)算

繞繩與卷筒是柔—?jiǎng)傂悦婷娼佑|問(wèn)題，分別將各個(gè)導(dǎo)入的模型在ANSYS WORKBENCH中設(shè)置好接觸參數(shù)、摩擦系數(shù)，繞繩一端施加1000N拉力，另一端設(shè)置為固定約束，將轉(zhuǎn)筒固定，劃分好網(wǎng)格，進(jìn)行結(jié)構(gòu)非線性分析。求得纏繞包角為120°時(shí)繞繩截面上的等效應(yīng)力分布如圖5所示，用探針probe分別探測(cè)出纏繞線繩在各個(gè)不同纏繞包角時(shí)的緊邊等效應(yīng)力和松邊等效應(yīng)力，結(jié)果如圖5、圖7所示，并將所得的結(jié)果列在表1中。

圖5 卷筒纏繞等效應(yīng)力圖

圖6 纏繞包角為60°、120°、180°時(shí)線繩兩端的等效應(yīng)力值

圖7 纏繞包角為240°、300°、360°時(shí)線繩兩端的等效應(yīng)力值

2.4 計(jì)算結(jié)果分析

由上述知，張緊力的計(jì)算公式為：F2≈1000/e0.3α。用數(shù)學(xué)分析軟件Matlab繪制出此公式的曲線圖，如圖8所示。x軸方向表示接觸時(shí)圓弧段參與的包角，y軸方向是繞繩的等效應(yīng)力。將表1中所得的結(jié)果點(diǎn)描在圖8等效應(yīng)力曲線圖中，如圖所示圓圈點(diǎn)即為有限元分析的結(jié)果。

圖8 公式與有限元分析計(jì)算結(jié)果比較

從圖8分析可看出，有限元分析計(jì)算的結(jié)果和公式計(jì)算的結(jié)果在一定范圍內(nèi)有一定程度的相符性。由此可知上述介紹的廣義歐拉公式在一定范圍內(nèi)可以指導(dǎo)實(shí)際的工程應(yīng)用。

3 結(jié)論

本文運(yùn)用有限元法對(duì)模型進(jìn)行計(jì)算分析，與公式計(jì)算所得的結(jié)果對(duì)比，借用數(shù)學(xué)分析軟件Matlab將它們繪制在一張圖上，直觀地進(jìn)行比較。

由結(jié)果比較可知，運(yùn)用ANSYS WORKBENCH有限元計(jì)算的結(jié)果和公式計(jì)算所得的結(jié)果在接觸圓弧參與的包角較小時(shí)，兩者的結(jié)果比較相符；圓弧參與的包角較大時(shí)，兩者的結(jié)果之間有一定的誤差。故本文介紹的纏繞于一般曲面體上的廣義歐拉公式在一定范圍內(nèi)可以指導(dǎo)實(shí)際的工程應(yīng)用，可以減少我們?cè)谟龅酱祟悊?wèn)題時(shí)計(jì)算張緊力時(shí)的誤差，比運(yùn)用經(jīng)典歐拉公式計(jì)算的結(jié)果更接近實(shí)際。所以上述介紹的纏繞于一般曲面體上的廣義歐拉公式具有更普遍的實(shí)用意義。

[1]楊可楨.機(jī)械設(shè)計(jì)基礎(chǔ)[M].北京：高等教育出版社，1989.

[2]楊景慧，陸玉，唐蓉城.機(jī)械設(shè)計(jì)（機(jī)械類）[M].北京：機(jī)械工業(yè)出版社，1995.

[3]曹助家，顧占山.歐拉公式在機(jī)械設(shè)計(jì)中的應(yīng)用[J].機(jī)械制造，1994.