TWI506258B - 齒輪的誤差量測方法 - Google Patents

齒輪的誤差量測方法 Download PDF

Info

Publication number
TWI506258B
TWI506258B TW102131227A TW102131227A TWI506258B TW I506258 B TWI506258 B TW I506258B TW 102131227 A TW102131227 A TW 102131227A TW 102131227 A TW102131227 A TW 102131227A TW I506258 B TWI506258 B TW I506258B
Authority
TW
Taiwan
Prior art keywords
error
function
pattern
gear
frequency
Prior art date
Application number
TW102131227A
Other languages
English (en)
Other versions
TW201508257A (zh
Inventor
Yi Cheng Chen
Ming Hung Chang
Original Assignee
Univ Nat Central
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Univ Nat Central filed Critical Univ Nat Central
Priority to TW102131227A priority Critical patent/TWI506258B/zh
Priority to US14/167,221 priority patent/US9588014B2/en
Publication of TW201508257A publication Critical patent/TW201508257A/zh
Application granted granted Critical
Publication of TWI506258B publication Critical patent/TWI506258B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Description

齒輪的誤差量測方法
本發明係與一種齒輪的誤差量測方法有關,特別是與一種能夠與單齒腹測試機配合使用的齒輪的誤差量測方法有關。
目前業界廣泛使用之齒輪量測儀器,是以微小直徑的球型側頭直接碰觸齒輪的齒面來測量齒輪的精度。但此種齒輪量測方法只能量測單一個齒輪的精度,且量測位置僅限於齒面上的特定位置,而且不適合量測兩個齒輪(或稱齒輪對)相嚙合時的傳動誤差。
因此,對於高精度且低噪音的齒輪,常使用單齒腹測試機來量測其精度。單齒腹測試機包括一主動齒輪及一被動齒輪組成一彼此嚙合的齒輪對,以進行齒輪傳動。齒輪傳動之主要功能,為高效率且平順地以其齒數比的關係將主動齒輪之動力傳遞給被動齒輪,而達到想要的輸出轉速與扭力。實際地齒輪對嚙合的過程會受到組裝的誤差與加工上的誤差之影響,造成一些不連續,而這種不連續定義為傳動誤差。透過分析齒形誤差、鄰接節距誤差與累積節距誤差,來定義齒輪對的嚙合狀況。此方法能夠快速的量測綜合的傳動誤差,所以適合作為產業界的品管方法。
為了解析單齒腹測試機上的訊號,通常使用快速傅利葉轉換(Fast Fourier Transform;FFT)配合嚙合頻率,將訊號分為高頻與低頻的部分,其中高頻的部分為齒形誤差,其主要與齒輪外形有關,而低頻的部分為累 積節距誤差,其通常與節緣偏擺有關,並將訊號作頻譜分析得到頻率與幅值關係,並從其中找到與齒輪精度及傳動噪音相關的特徵。
由於FFT牽涉到濾波的概念,所以會產生相位的偏差。在定義齒輪對的嚙合狀況時會有一定的困難,容易造成齒輪精度的誤判。因此,若能克服定義齒輪對的嚙合狀況的困難,則能更精確的判斷齒輪的精度。
本發明之一目的係在於提供一齒輪的誤差量測方法,能更精確的定義齒輪對的單一節距誤差、鄰接節距誤差及累積節距誤差。
為達上述之一或部份或全部目的或是其他目的,本發明提供一種齒輪的誤差量測方法包括:提供一齒輪對(gear pair)以及一單齒腹測試裝置,單齒腹測試裝置包括一運算單元,單齒腹測試裝置將齒輪對進行單齒腹測試進而得到一傳動誤差數據,並利用傳動誤差數據產生一測試訊號圖形;運算單元利用經驗模態分析(EMD)分解測試訊號圖形進而產生複數個固有模態函數圖形(Intrinsic Mode Function),該等固有模態函數圖形包括一第一函數圖形及一第二函數圖形,其中第一函數圖形的頻率包括一第一範圍,且該第一函數圖形的頻率最接近齒輪對的一嚙合頻率,第二函數圖形的頻率包括一第二範圍,且該第二函數圖形的頻率最接近齒輪對的一轉軸頻率,並且第一範圍大於第二範圍;運算單元從該等固有模態函數圖形中選擇複數個第三函數圖形,其中該等第三函數圖形的頻率包括一第三範圍,該第三範圍大於第二範圍,且小於或等於第一範圍;運算單元疊加該等第三函數圖形進而產生一疊加函數圖形,疊加函數圖形中具有複數個波峰;運算單元利用疊加函數圖形中每兩相鄰波峰的差值計算出複數個單一 節距誤差,再利用該等單一節距誤差計算出複數個鄰接節距誤差及一累積節距誤差,其中鄰接節距誤差為每兩相鄰單一節距誤差的差值,累積節距誤差為每一單一節距誤差的疊加。
在一實施例中,齒輪的誤差量測方法更包括:運算單元量測第一函數圖形的振幅而得到一齒形誤差;運算單元量測第二函數圖形的振幅而得到一偏擺誤差。
在一實施例中,齒輪的誤差量測方法更包括:運算單元根據齒形誤差、偏擺誤差、單一節距誤差、鄰接節距誤差及累積節距誤差來定義齒輪對中之一齒輪的一精度。
在一實施例中,運算單元利用經驗模態分析(EMD)分解測試訊號圖形進而產生複數個固有模態函數圖形之步驟包括:運算單元判斷該等固有模態函數圖形是否有模態混雜(mode mixing);若有模態混雜,運算單元再利用集成經驗模態分析(EEMD)分解測試訊號圖形。
在一實施例中,運算單元量測第一函數圖形的振幅而得到一齒形誤差的步驟,係為量測第一函數圖形的一上包落線及一下包落線平均的差值而得到齒型誤差。
在一實施例中,運算單元量測第二函數圖形的振幅而得到一偏擺誤差的步驟,係為量測第二函數圖形的一上包落線及一下包落線平均的差值而得到偏擺誤差。
100‧‧‧齒輪對
200‧‧‧單齒腹測試裝置
210‧‧‧運算單元
圖1A係為本發明之一實施例的齒輪對誤差量測示意圖。
圖1B係為本發明之一實施例的齒輪的誤差量測方法步驟 圖。
圖2係為本發明之一實施例的測試訊號圖形。
圖3A及3B係為本發明之一實施例的固有模態函數圖形。
圖4係為本發明之一實施例的齒形誤差函數圖形。
圖5係為本發明之一實施例的偏擺誤差函數圖形。
圖6係為本發明之一實施例的疊加函數圖形。
有關本發明之前述及其他技術內容、特點與功效,在以下配合參考圖式之一較佳實施例的詳細說明中,將可清楚的呈現。以下實施例中所提到的方向用語,例如:上、下、左、右、前或後等,僅是用於參照隨附圖式的方向。因此,該等方向用語僅是用於說明並非是用於限制本發明。
請參閱圖1A及1B,圖1A係為本發明之一實施例的齒輪對誤差量測示意圖,而圖1B係為本發明之一實施例的齒輪的誤差量測方法步驟圖。齒輪的誤差量測方法步驟如下:
步驟S100:提供一齒輪對(gear pair)100以及一單齒腹測試裝置200,單齒腹測試裝置200包括一運算單元210。單齒腹測試裝置200將齒輪對100進行單齒腹測試進而得到一傳動誤差數據,並利用該傳動誤差數據產生一如圖2之測試訊號圖形。
步驟S110:運算單元210利用經驗模態分析(EMD)分解測試訊號圖形進而產生如圖3A及3B的複數個固有模態函數圖形(Intrinsic Mode Function)。
在步驟S110中,該等固有模態函數圖形包括一第一函數圖形及一第二函數圖形,其中第一函數圖形的頻率包括一第一範圍,且第一函數圖形的頻率最接近齒輪對100的一嚙合頻率,第二函數圖形的頻率包括一第二範圍,且第二函數圖形的頻率最接近齒輪對100的一轉軸頻率,並且,第一範圍係大於第二範圍。此外,運算單元210利用經驗模態分析(EMD)分解測試訊號圖形來得到的固有模態函數圖形(IMF)可能會有模態混雜(mode mixing)的情況發生。若遇到這種情形則執行步驟S111及S112。
步驟S111:判斷該等固有模態函數圖形是否有模態混雜(mode mixing)。若是,則執行步驟S112,若否,則執行步驟S120。
步驟S112:若有模態混雜,運算單元210再利用集成經驗模態分析(EEMD)分解測試訊號圖形。
步驟S120:運算單元210從該等固有模態函數圖形中選擇複數個第三函數圖形。該等第三函數圖形的頻率包括一第三範圍,該第三範圍係大於第二範圍,且小於或等於第一範圍。依此原則,第一函數圖形同時為該等第三函數圖形之其一。
在一實施例中,測試訊號圖形的取樣頻率為2000Hz,運算單元210利用經驗模態分析(EMD)分解測試訊號圖形進而產生十六個固有模態函數圖形IMF C1~IMF C16。其中,IMF C1包含頻率大約是500~1000Hz,IMF C2為250~500Hz,IMF C3為125~250Hz,IMF C4為62.5~125Hz等以此類推。在本實施例中,IMF C1係例如為第一函數圖形,IMF C4係例如為第二函數圖形,而IMF C1~IMF C3則係例如為該等第三函數圖形。
步驟S130:運算單元210疊加該等第三函數圖形進而產生一 如圖4的疊加函數圖形,疊加函數圖形中具有複數個波峰。
步驟S131:運算單元量測第一函數圖形的振幅而得到一齒形誤差。
步驟S132:運算單元量測第二函數圖形的振幅而得到一偏擺誤差。
步驟S140:運算單元210利用疊加函數圖形中每兩相鄰波峰的差值計算出複數個單一節距誤差,再利用該等單一節距誤差計算出複數個鄰接節距誤差及一累積節距誤差,其中鄰接節距誤差為每兩相鄰單一節距誤差的差值,累積節距誤差為每一該單一節距誤差的疊加。
步驟S150:運算單元根據該齒形誤差、該偏擺誤差、該單一節距誤差、該鄰接節距誤差及該累積節距誤差來定義該齒輪對中之一齒輪的一精度。
如圖2,其係為本發明之步驟S10)0)中產生的測試訊號圖形。測試訊號圖形之橫軸為主動齒輪轉動的角度;縱軸為齒輪傳動誤差。測試訊號圖形包含週期性且複雜的訊號,將圖2之測試訊號圖形利用經驗模態分析(EMD)來分解(步驟S110),則會產生如圖3之複數個固有模態函數圖形(Intrinsic Mode Function)。
圖3A及3B中共包括有十六個固有模態函數圖形。齒輪提供一嚙合頻率及一轉軸頻率,該等固有模態函數圖形提供一第一波動頻率及一第二波動頻率。從該固有模態函數圖形中挑選與第一波動頻率最接近的即為第一函數圖形IMF C1,而與第二波動頻率最接近的即為第二函數圖形IMF C4。量測第一函數圖形IMF C1的上包落線及下包落線的平均值會得到 如圖4的齒形誤差,量測第二函數圖形IMF C4的上包落線及下包落線的平均值會得到如圖5的偏擺誤差。
在步驟S120中由多個固有模態函數圖形(IMF)中挑選出與上述嚙合頻率或轉軸頻率最接近的或甚至是相同波動頻率的固有模態函數圖形(IMF)來進行分析。如圖4,與嚙合頻率相符合的固有模態函數圖形IMF C1,該固有模態函數圖形IMF C1的上包落線及下包落線的平均值係被定義為齒形誤差。如圖5,與轉軸頻率相符合的固有模態函數圖形IMF C4,該固有模態函數圖形IMF C4的上包落線及下包落線的平均值係被定義為偏擺誤差。若將固有模態函數圖形IMF C1與固有模態函數圖形IMF C4之間的該等固有模態函數圖形疊加會得到一如圖6的疊加函數圖形。其中,在疊加固有模態函數圖形IMF C1與固有模態函數圖形IMF C4之間的該等固有模態函數圖形時,包括IMF C1但不包括IMF C4,也就是說,疊加固有模態函數圖形IMF C1到IMF C3。如圖6的疊加函數圖形中具有複數個波峰,定義單一節距誤差為每兩相鄰的波峰的差值,例如但不限於E1、E2。鄰接節距誤差為每兩相鄰的單一節距誤差的差值,例如但不限於| E1-E2 |。累積節距誤差則係疊加每一單一節距誤差,例如但不限於疊加| E1-E2 |、| E2-E3 |、| E3-E4 |等。
最後在步驟S150中,根據齒形誤差值、偏擺誤差、單一節距誤差、鄰接節距誤差及累積節距誤差來定義齒輪對中之齒輪的精度。
惟以上所述者,僅為本發明之較佳實施例而已,當不能以此限定本發明實施之範圍,即大凡依本發明申請專利範圍及發明說明內容所作之簡單的等效變化與修飾,皆仍屬本發明專利涵蓋之範圍內。另外本發 明的任一實施例或申請專利範圍不須達成本發明所揭露之全部目的或優點或特點。此外,摘要部分和標題僅是用來輔助專利文件搜尋之用,並非用來限制本發明之權利範圍。

Claims (4)

  1. 一種齒輪的誤差量測方法,包括:提供一齒輪對(gear pair)以及一單齒腹測試裝置,該單齒腹測試裝置、包括一運算單元,該單齒腹測試裝置將該齒輪對進行單齒腹測試進而得到一傳動誤差數據,並利用該傳動誤差數據產生一測試訊號圖形,該測試訊號圖形為一傳動誤差相對於轉動角度之關係圖;該運算單元利用經驗模態分析(EMD)分解該測試訊號圖形進而產生複數個固有模態函數圖形(Intrinsic Mode Function),每一該固有模態函數圖形為一傳動誤差相對於轉動角度之關係圖,該等固有模態函數圖形包括一第一函數圖形及一第二函數圖形,其中該第一函數圖形的頻率包括於一第一範圍內,且該第一函數圖形的頻率最接近該齒輪對的一嚙合頻率,該第二函數圖形的頻率包括於一第二範圍內,且該第二函數圖形的頻率最接近該齒輪對的一轉軸頻率,並且該第一範圍內的數值係大於該第二範圍內的數值;該運算單元從該等固有模態函數圖形中選擇複數個第三函數圖形,其中每一該第三函數圖形的頻率包括於一第三範圍內,該第三範圍內的數值大於或等於該第二範圍內的數值,且小於或等於該第一範圍內的數值,並且該第一函數圖形同時為該等第三函數圖形之其一;該運算單元疊加該等第三函數圖形進而產生一疊加函數圖形,該疊加函數圖形中具有複數個波峰;該運算單元利用該疊加函數圖形中每兩相鄰波峰的差值計算出複數個單一節距誤差,再利用該等單一節距誤差計算出複數個鄰接節距誤差 及一累積節距誤差,其中該鄰接節距誤差為每兩相鄰單一節距誤差的差值,該累積節距誤差為每一該單一節距誤差的疊加;該運算單元量測該第一函數圖形的振幅而得到一齒形誤差,並且量測該第二函數圖形的振幅而得到一偏擺誤差;以及該運算單元根據該齒形誤差、該偏擺誤差、該單一節距誤差、該鄰接節距誤差及該累積節距誤差來定義該齒輪對中之一齒輪的一精度。
  2. 如申請專利範圍第1項所述之齒輪的誤差量測方法,其中該運算單元利用經驗模態分析(EMD)分解該測試訊號圖形進而產生複數個固有模態函數圖形之步驟包括:該運算單元判斷該等固有模態函數圖形是否有模態混雜(mode mixing);以及若有模態混雜,該運算單元再利用集成經驗模態分析(EEMD)分解該測試訊號圖形。
  3. 如申請專利範圍第1項所述之齒輪的誤差量測方法,其中該運算單元量測該第一函數圖形的振幅而得到一齒形誤差的步驟,係為量測該第一函數圖形的一上包落線及一下包落線平均的差值而得到該齒型誤差。
  4. 如申請專利範圍第1項所述之齒輪的誤差量測方法,其中該運算單元量測該第二函數圖形的振幅而得到一偏擺誤差的步驟,係為量測該第二函數圖形的一上包落線及一下包落線平均的差值而得到該偏擺誤差。
TW102131227A 2013-08-30 2013-08-30 齒輪的誤差量測方法 TWI506258B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
TW102131227A TWI506258B (zh) 2013-08-30 2013-08-30 齒輪的誤差量測方法
US14/167,221 US9588014B2 (en) 2013-08-30 2014-01-29 Error measuring method of gear

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW102131227A TWI506258B (zh) 2013-08-30 2013-08-30 齒輪的誤差量測方法

Publications (2)

Publication Number Publication Date
TW201508257A TW201508257A (zh) 2015-03-01
TWI506258B true TWI506258B (zh) 2015-11-01

Family

ID=52584380

Family Applications (1)

Application Number Title Priority Date Filing Date
TW102131227A TWI506258B (zh) 2013-08-30 2013-08-30 齒輪的誤差量測方法

Country Status (2)

Country Link
US (1) US9588014B2 (zh)
TW (1) TWI506258B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI518310B (zh) * 2012-09-21 2016-01-21 國立中央大學 齒輪精度的檢測方法
US10228303B2 (en) 2014-03-31 2019-03-12 Automation Controls & Engineering, LLC Flexible automation cell for performing secondary operations in concert with a machining center and roll check operations
US9551628B2 (en) 2014-03-31 2017-01-24 Automation Controls & Engineering, LLC Flexible automation cell for performing secondary operations in concert with a machining center and roll check operations
US9599535B2 (en) * 2015-05-19 2017-03-21 GM Global Technology Operations LLC Method and apparatus for excitation of gear rattle
CN105258942A (zh) * 2015-09-30 2016-01-20 南阳理工学院 一种采煤机摇臂齿轮故障诊断系统及方法
CN105241663B (zh) * 2015-11-06 2017-08-25 济南大学 弯扭复合动态加载装置
CN105699072B (zh) * 2016-01-11 2018-05-01 石家庄铁道大学 一种基于级联经验模态分解齿轮故障诊断方法
CN105928701B (zh) * 2016-04-29 2017-07-18 石家庄铁道大学 基于相关分析的emd过程中有效imf的判定方法
CN106053034B (zh) * 2016-05-11 2018-03-13 西安交通大学 一种基于调频信息重构的旋转机械全息诊断方法
CN106596116A (zh) * 2016-11-29 2017-04-26 西安理工大学 一种风力发电机组振动故障诊断方法
JP6897785B2 (ja) * 2017-09-29 2021-07-07 新東工業株式会社 ギヤ位置決め装置、応力測定システム、ギヤ位置決め方法及び応力測定方法
JP7325409B2 (ja) * 2017-10-23 2023-08-14 グリーソン メトロロジー システムズ コーポレイション 非接触センサを利用した複合機械試験機プラットフォーム上での歯付き物品の測定
CN110398363B (zh) * 2018-04-19 2021-03-16 中国科学院沈阳自动化研究所 一种基于图形化变换的rv减速器性能状态监测方法
CN109000926A (zh) * 2018-10-25 2018-12-14 中国计量大学 基于emd和近似熵的滚动轴承声发射信号特征提取方法
CN109975012A (zh) * 2019-04-24 2019-07-05 西安交通大学 一种传动误差微分信号与eemd算法结合的齿轮裂纹诊断方法
CN110646199B (zh) * 2019-09-06 2022-01-07 天津工业大学 一种基于加权导数动态时间规整的齿轮微弱故障诊断方法
CN110569478B (zh) * 2019-09-12 2021-02-02 北京航空航天大学 用于编码器信号分析的改进变分模态分解方法
CN114216677B (zh) * 2021-12-15 2024-06-28 重庆青山工业有限责任公司 齿面波纹阶次检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101701869A (zh) * 2009-12-04 2010-05-05 北京工业大学 一种适于控制锥齿轮啮合质量的便携式传动噪声测试仪

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5578229A (en) * 1978-12-11 1980-06-12 Ono Sokki Co Ltd Gearing test unit
US6862558B2 (en) * 2001-02-14 2005-03-01 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Empirical mode decomposition for analyzing acoustical signals
TWI460416B (zh) * 2011-03-28 2014-11-11 Univ Nat Taiwan 機械系統狀態之判斷方法及判斷裝置
US8732113B2 (en) * 2011-09-23 2014-05-20 National Central University Method for quantifying and modeling degree of nonlinearity, combined nonlinearity, and nonstationarity

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101701869A (zh) * 2009-12-04 2010-05-05 北京工业大学 一种适于控制锥齿轮啮合质量的便携式传动噪声测试仪

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陳仁傑,希爾伯特黃變換(HHT)於變轉速之齒輪故障診斷之應用,國立中央大學機械工程學系碩士論文,2011年6月(國家圖書館電子論文上傳日2011年12月07日) 劉鈞豪,齒輪對之單齒腹檢測與運動誤差分析,國立中正大學機械工程研究所碩士論文,2009年7月22日(國家圖書館紙本論文上架日2010年03月11日) 楊智宇,希爾伯特-黃變換(Hilbert-Huang Transform)結合主成份分析與類神經網路在齒輪故障程度之診斷,國立中央大學機械工程學系碩士論文,2011年7月(國家圖書館電子論文上傳日2011年12月07日) *

Also Published As

Publication number Publication date
US9588014B2 (en) 2017-03-07
US20150066390A1 (en) 2015-03-05
TW201508257A (zh) 2015-03-01

Similar Documents

Publication Publication Date Title
TWI506258B (zh) 齒輪的誤差量測方法
TWI518310B (zh) 齒輪精度的檢測方法
Simon Optimization of face-hobbed hypoid gears
CN103090960B (zh) 用于获取振动数据并对振动数据去噪声的系统和方法
CN107436236A (zh) 车辆变速箱的噪声检测方法及检测系统
CN107704695B (zh) 一种滚动轴承外圈缺陷全尺寸定量诊断方法
US11927501B2 (en) Method and device for monitoring a gear system
JP6488645B2 (ja) 走行抵抗算出方法、走行抵抗計測方法、及び走行抵抗算出装置
CN105300691B (zh) 基于最佳安装距的锥齿轮传动误差测量方法
JP2013545098A5 (zh)
Lin et al. Pitch deviation measurement and analysis of curve-face gear pair
CN112257192B (zh) 一种基于实测齿廓偏差的一齿切向综合偏差计算方法
CN106574881A (zh) 旋转机的异常检测装置、旋转机的异常检测方法以及旋转机
JP2006284443A (ja) 燃料消費演算方法
CN108931204A (zh) 三维曲面上的凹凸形状测定方法
JP4407553B2 (ja) 振動計測方法、及び、振動計測装置
CN104063605A (zh) 一种用于齿轮渐开线齿形误差检验方法
Liang et al. Modelling method, simulation and experimental verification of hypoid gear involved tooth surface deviation under manufacturing process
Kuzio et al. Identification of fault frequency variation in the envelope spectrum in the vibration-based local damage detection in possible changing load/speed conditions
CN105138800A (zh) 基于分段拟合数据处理算法的荧光光谱数据噪声滤波方法
Cho et al. Study of the transmission error prediction of a five-speed manual transmission system
Zhao et al. Multiple transient extraction algorithm and its application in bearing fault diagnosis
JP2014119042A (ja) 歯車対の設計装置
JP2009250954A (ja) 両歯面噛合式歯車検査装置
CN109978993A (zh) 渐开线斜齿面齿轮三维建模方法