TW201822948A - 偵測工具機的線性滑軌預壓值變化的方法 - Google Patents

偵測工具機的線性滑軌預壓值變化的方法 Download PDF

Info

Publication number
TW201822948A
TW201822948A TW105142621A TW105142621A TW201822948A TW 201822948 A TW201822948 A TW 201822948A TW 105142621 A TW105142621 A TW 105142621A TW 105142621 A TW105142621 A TW 105142621A TW 201822948 A TW201822948 A TW 201822948A
Authority
TW
Taiwan
Prior art keywords
linear slide
natural frequency
machine tool
vibration
calculation module
Prior art date
Application number
TW105142621A
Other languages
English (en)
Other versions
TWI609737B (zh
Inventor
鄭志鈞
蔡秉均
程文男
郭育昕
鄭尹俊
鄭錫勳
余冠德
Original Assignee
上銀科技股份有限公司
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 上銀科技股份有限公司 filed Critical 上銀科技股份有限公司
Priority to TW105142621A priority Critical patent/TWI609737B/zh
Application granted granted Critical
Publication of TWI609737B publication Critical patent/TWI609737B/zh
Publication of TW201822948A publication Critical patent/TW201822948A/zh

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

一種偵測工具機的線性滑軌預壓值變化的方法,工具機包含線性滑軌及設置於線性滑軌上的工作載台,工作載台上設置有多個振動感測器,該方法包含以下步驟:每一振動感測器在第一時點產生一第一振動信號;根據至少一相關於工作載台的理論模態振型與該等第一振動信號計算出至少一相關於工作載台的第一自然頻率;每一振動感測器在第二時點產生一第二振動信號;根據至少一理論模態振型與該等第二振動信號計算出至少一相關於工作載台的第二自然頻率;及根據該至少一第一自然頻率與該至少一第二自然頻率判斷線性滑軌的預壓值變化。

Description

偵測工具機的線性滑軌預壓值變化的方法
本發明是有關於一種偵測線性滑軌預壓值變化的方法,特別是指一種偵測工具機的線性滑軌預壓值變化的方法。
工具機的線性滑軌在出廠時即對應有一預壓值,一般來說,線性滑軌的預壓值會隨著線性滑軌的使用損耗而漸漸變小;而當預壓值下降到不符使用需求時,則需汰換線性滑軌。
目前來說,若要量測線性滑軌的預壓值,需先將線性滑軌從工具機上拆卸下來,再利用量測工具測得預壓值的大小;因此,若要藉由實際量測線性滑軌的預壓值來觀察其變化,需一再地在工具機上裝卸線性滑軌,相當地不便與耗時。
因此,本發明之目的,即在提供一種不需從工具機上拆卸線性滑軌即可實施的偵測工具機的線性滑軌預壓值變化的方法。
於是,本發明偵測工具機的線性滑軌預壓值變化的方法由多個振動感測器及一計算模組實施。該工具機包含至少一線性滑軌及一設置於該至少一線性滑軌上的工作載台,該等振動感測器設置於該工作載台。該偵測工具機的線性滑軌預壓值變化的方法包含一步驟(a)、一步驟(b)、一步驟(c)、一步驟(d),及一步驟(e)。
該步驟(a)是每一振動感測器在一第一時點感測該工作載台的振動而產生一對應該第一時點的第一振動信號。
該步驟(b)是該計算模組根據至少一相關於該工作載台的理論模態振型與該等第一振動信號計算出至少一相關於該工作載台的第一自然頻率。
該步驟(c)是每一振動感測器在一第二時點感測該工作載台的振動而產生一對應該第二時點的第二振動信號。
該步驟(d)是該計算模組根據該至少一理論模態振型與該等第二振動信號計算出至少一相關於該工作載台的第二自然頻率。
該步驟(e)是該計算模組根據該至少一第一自然頻率與該至少一第二自然頻率判斷該至少一線性滑軌的預壓值變化。
本發明之功效在於:不需從工具機上拆卸線性滑軌即可有效地偵測出線性滑軌的預壓值的變化。
在本發明被詳細描述之前,應當注意在以下的說明內容中,類似的元件是以相同的編號來表示。
參閱圖1,本發明偵測工具機的線性滑軌預壓值變化的方法由多個振動感測器1及一電連接該等振動感測器1的計算模組(圖未示)實施。
該工具機2為一習知的工具機,並包含一馬達21與一進給單元22。該進給單元22包括一連接該馬達21的螺桿221、多個樞接該螺桿221的軸承222、二線性滑軌223,及一設置於該等線性滑軌223上的工作載台224,且該二線性滑軌223均具有相同的一初始預壓值。
該等振動感測器1設置於該工作載台224,並用於感測該工作載台224的振動;在此以三個振動感測器來說明,但並不限於此,振動感測器的個數只要是複數即可。該計算模組包含一具有一般計算功能的電腦。
參閱圖2,以下詳述本發明偵測工具機的線性滑軌預壓值變化的方法的一實施方式。
首先,在步驟31,該馬達21在一第一時點激振(excite)該進給單元22而使該工作載台224振動;每一振動感測器1感測該工作載台224的振動而產生一對應該第一時點的第一振動信號。雖然在此是利用該工具機2本身所具有的該馬達21來激振該進給單元22,但可不限於此;在另一實施方式中,也可利用一外部的激振器(vibration exciter)來激振該進給單元22而使該工作載台224振動。
接著,在步驟32,該計算模組利用一現有的操作模態分析法(operational modal analysis, OMA)根據該等第一振動信號計算出相關於該工作載台224的多個自然頻率(natural frequency)與多個分別對應該等自然頻率的操作模態振型(operational mode shape)。該操作模態分析法為本發明所屬技術領域的具有通常知識者所熟悉,然而,為了便於說明本發明的技術內容,以下簡述之。
令該等振動感測器1所感測到的該等第一振動信號所分別對應的頻域信號為,及,其中k 表示頻率。接著,根據該等頻域信號產生對應該工作載台224的多個位移傳遞率(displacement transmissibility),其中,且函數G 為功率頻譜密度函數(power spectrum density function)。接著,對於每一頻率k ,針對矩陣進行奇異值分解(singular value decomposition),即,其中奇異值矩陣,且
對於每一頻率k ,均可計算出一對應的奇異值,而所有的倒數,也就是所有的,所形成的曲線中的每一峰值所對應的頻率即為一自然頻率,而該自然頻率所對應的操作模態振型也就是該自然頻率對應的奇異值在奇異值分解過程中所對應的特徵向量(eigenvector)。例如,參閱圖3,其中的頻率、頻率,及頻率均為自然頻率;而所分別對應的特徵向量也就是頻率、頻率,及頻率所分別對應的操作模態振型。
接著,在步驟33,該計算模組根據預先藉由習知的有限元素法(finite element method, FEM)分析出來的最符合該工作載台224的振動方式的一理論模態振型,與每一操作模態振型的相似程度從該等自然頻率中選出一第一自然頻率,其中係選出與該理論模態振型有最大相似程度的操作模態振型所對應的自然頻率做為該第一自然頻率。
在此是利用現有的模態確認指標(modal assurance criterion, MAC)來計算該理論模態振型與該操作模態振型之間的相似程度。模態向量與模態向量的模態確認指標值的計算方式為,其中模態確認指標值愈高代表的正交性愈低,相似程度也就愈高。若向量表示該理論模態振型,向量,,…,表示步驟32所計算出的該等操作模態振型,則在此是計算 {},並以具有最大MAC值的操作模態振型所對應的自然頻率為該所選出的第一自然頻率。
接著,在步驟34,該馬達21在一晚於第一時點的第二時點激振該進給單元22而使該工作載台224振動;每一振動感測器1在該第二時點感測該工作載台224的振動而產生一對應該第二時點的第二振動信號。
接著,在步驟35,類似於針對該等第一振動信號的處理,該計算模組利用該操作模態分析法根據該等第二振動信號計算出相關於該工作載台224的多個自然頻率與多個分別對應該等自然頻率的操作模態振型。
接著,在步驟36,該計算模組根據步驟33所述的該理論模態振型,與步驟35所述的每一操作模態振型的相似程度從該等自然頻率中選出一第二自然頻率,其中係選出與該理論模態振型有最大相似程度的操作模態振型所對應的自然頻率做為該第二自然頻率。在此也是利用該模態確認指標來計算該理論模態振型與該操作模態振型之間的相似程度。
接著,在步驟37,該計算模組根據該第一自然頻率與該第二自然頻率判斷該等線性滑軌223的預壓值變化。當該第二自然頻率小於該第一自然頻率時,該計算模組判定該等線性滑軌223的預壓值變小,且判定該等線性滑軌223的預壓值的下降比率為由該第一自然頻率下降至該第二自然頻率的下降比率。
由上述可知,若該第一時點為該等線性滑軌223在裝設至該工具機2後開始運作的時點,且該第二時點為該等線性滑軌223運作一段時間後的時點,則該等線性滑軌223在該第一時點對應的預壓值為該等線性滑軌223在出廠時即具有的該初始預壓值,藉由上述自然頻率的變化估計出該等線性滑軌223的預壓值的下降比率後,能在不將該等線性滑軌223從該工具機2上拆卸下來的情況下,進一步根據該初始預壓值與該下降比率估計出該等線性滑軌223在該第二時點的預壓值。
雖然上述實施方式僅利用到單一個理論模態振型,但可不限於此。在另一實施方式中,可預先利用該有限元素法分析出多個相關於該工作載台224的理論模態振型,然後計算根據該等第一振動信號所分析出的每一操作模態振型與每一理論模態振型的模態確認指標值,並選出具有最大模態確認指標值的操作模態振型所對應的自然頻率做為該第一自然頻率;並計算根據該等第二振動信號所分析出的每一操作模態振型與每一理論模態振型的模態確認指標值,且選出具有最大模態確認指標值的操作模態振型所對應的自然頻率做為該第二自然頻率。
也就是說,若向量,,…,表示該等理論模態振型,向量,,…,表示根據該等第一振動信號所分析出的該等操作模態振型,則計算 {},並以具有最大MAC值的操作模態振型所對應的自然頻率做為該第一自然頻率;而求取該第二自然頻率的方式類似前述求取該第一自然頻率的方式,故不再贅述。
在另一實施方式中,也可利用根據有限元素法所分析出來的多個相關於該工作載台224的理論模態振型與根據該等第一振動信號所計算出來的多個操作模態振型來獲得對應該第一時點的多個第一自然頻率;並利用該等理論模態振型與根據該等第二振動信號所計算出來的多個操作模態振型來獲得對應該第二時點的多個第二自然頻率;然後,根據該等第一自然頻率與該等第二自然頻率來預測該等線性滑軌223的預壓值的下降比率。
舉例來說,可從根據該等第一振動信號所計算出的該等操作模態振型中選出分別具有最大模態確認指標值與次大模態確認指標值的兩個操作模態振型,並以該兩個操作模態振型所對應的自然頻率做為兩個第一自然頻率,在此以向量分別表示該兩個第一自然頻率所對應的操作模態振型,其中對應的自然頻率小於對應的自然頻率。
類似地,令向量分別表示根據該等第二振動信號所獲得的兩個第二自然頻率所對應的操作模態振型,其中對應的自然頻率小於對應的自然頻率。然後,將[,,,] T 輸入一預先訓練好的預測模型,例如類神經網路或線性迴歸模型,來預測該等線性滑軌223的預壓值的下降比率。
綜上所述,本發明偵測工具機的線性滑軌預壓值變化的方法,根據至少一理論模態振型與該等第一振動信號分析出對應該第一時點的相關於該工作載台的至少一第一自然頻率;並根據該至少一理論模態振型與該等第二振動信號分析出對應該第二時點的相關於該工作載台的至少一第二自然頻率,且根據該至少一第一自然頻率與該至少一第二自然頻率來估計/預測該等線性滑軌的預壓值的下降比率,能在不將該等線性滑軌從該工具機拆卸下來的情況下,判斷出該等線性滑軌的預壓值變化,故確實能達成本發明的目的。
惟以上所述者,僅為本發明的實施例而已,當不能以此限定本發明實施的範圍,凡是依本發明申請專利範圍及專利說明書內容所作的簡單的等效變化與修飾,皆仍屬本發明專利涵蓋的範圍內。
1‧‧‧振動感測器
2‧‧‧工具機
21‧‧‧馬達
22‧‧‧進給單元
221‧‧‧螺桿
222‧‧‧軸承
223‧‧‧線性滑軌
224‧‧‧工作載台
31~37‧‧‧步驟
本發明的其他的特徵及功效,將於參照圖式的實施方式中清楚地呈現,其中: 圖1是一俯視示意圖,說明一工具機與將多個振動感測器設置於該工具機的一工作載台的態樣; 圖2是一流程圖,說明本發明偵測工具機的線性滑軌預壓值變化的方法的一實施方式;及 圖3是一示意圖,說明利用奇異值的倒數的曲線求取相關於該工作載台的自然頻率。

Claims (7)

  1. 一種偵測工具機的線性滑軌預壓值變化的方法,由多個振動感測器及一計算模組實施,該工具機包含至少一線性滑軌及一設置於該至少一線性滑軌上的工作載台,該等振動感測器設置於該工作載台,該偵測工具機的線性滑軌預壓值變化的方法包含以下步驟: (a)每一振動感測器在一第一時點感測該工作載台的振動而產生一對應該第一時點的第一振動信號; (b)該計算模組根據至少一相關於該工作載台的理論模態振型與該等第一振動信號計算出至少一相關於該工作載台的第一自然頻率; (c)每一振動感測器在一第二時點感測該工作載台的振動而產生一對應該第二時點的第二振動信號; (d)該計算模組根據該至少一理論模態振型與該等第二振動信號計算出至少一相關於該工作載台的第二自然頻率;及 (e)該計算模組根據該至少一第一自然頻率與該至少一第二自然頻率判斷該至少一線性滑軌的預壓值變化。
  2. 如請求項1所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(b),該計算模組利用一操作模態分析法根據該等第一振動信號計算出相關於該工作載台的多個自然頻率與多個分別對應該等自然頻率的操作模態振型,並根據每一操作模態振型與每一理論模態振型的相似程度從該等自然頻率中選出該至少一第一自然頻率,且在該步驟(d),該計算模組利用該操作模態分析法根據該等第二振動信號計算出相關於該工作載台的多個自然頻率與多個分別對應該等自然頻率的操作模態振型,並根據每一操作模態振型與每一理論模態振型的相似程度從該等自然頻率中選出該至少一第二自然頻率。
  3. 如請求項2所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(b),該計算模組從該等自然頻率中選出一第一自然頻率,其中該第一自然頻率所對應的操作模態振型所對應的相似程度為該等相似程度中的一最大者,且在該步驟(d),該計算模組從該等自然頻率中選出一第二自然頻率,其中該第二自然頻率所對應的操作模態振型所對應的相似程度為該等相似程度中的一最大者。
  4. 如請求項3所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(b)與該步驟(d),該操作模態振型與該理論模態振型的相似程度對應該操作模態振型與該理論模態振型的一模態確認指標值。
  5. 如請求項3所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(e),當該第二自然頻率小於該第一自然頻率時,該計算模組判定該至少一線性滑軌的預壓值變小。
  6. 如請求項1所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(b),該計算模組計算出一第一自然頻率,在該步驟(d),該計算模組計算出一第二自然頻率,且在該步驟(e),當該第二自然頻率小於該第一自然頻率時,該計算模組判定該線性滑軌的預壓值變小。
  7. 如請求項5或6所述的偵測工具機的線性滑軌預壓值變化的方法,其中在該步驟(e),該至少一線性滑軌的預壓值的下降比率為由該第一自然頻率下降至該第二自然頻率的下降比率。
TW105142621A 2016-12-22 2016-12-22 偵測工具機的線性滑軌預壓值變化的方法 TWI609737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW105142621A TWI609737B (zh) 2016-12-22 2016-12-22 偵測工具機的線性滑軌預壓值變化的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW105142621A TWI609737B (zh) 2016-12-22 2016-12-22 偵測工具機的線性滑軌預壓值變化的方法

Publications (2)

Publication Number Publication Date
TWI609737B TWI609737B (zh) 2018-01-01
TW201822948A true TW201822948A (zh) 2018-07-01

Family

ID=61728491

Family Applications (1)

Application Number Title Priority Date Filing Date
TW105142621A TWI609737B (zh) 2016-12-22 2016-12-22 偵測工具機的線性滑軌預壓值變化的方法

Country Status (1)

Country Link
TW (1) TWI609737B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI689371B (zh) * 2019-05-24 2020-04-01 國立虎尾科技大學 硬軌工具機的即時檢測補償裝置及其方法
TWI821491B (zh) * 2018-12-27 2023-11-11 日商迪思科股份有限公司 具備自我診斷功能的加工裝置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI756780B (zh) * 2020-08-13 2022-03-01 國立屏東科技大學 車輛之車架及其設計方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7185546B2 (en) * 2004-06-14 2007-03-06 Ascenx Systems and methods for measuring belt tension
TWI470189B (zh) * 2011-04-27 2015-01-21 Yu Hsiang Pan 機械系統狀態之判斷方法及判斷裝置
CN102323033B (zh) * 2011-05-30 2013-08-21 北京工业大学 一种测试结合面法向动态特性的装置
TWI504478B (zh) * 2014-06-04 2015-10-21 Hiwin Tech Corp Monitoring System of Ball Screw Pre - pressure and Its Monitoring Method
CN103994812A (zh) * 2014-06-09 2014-08-20 北京航空航天大学 可施加预紧力的固有频率综合测试试验平台

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI821491B (zh) * 2018-12-27 2023-11-11 日商迪思科股份有限公司 具備自我診斷功能的加工裝置
TWI689371B (zh) * 2019-05-24 2020-04-01 國立虎尾科技大學 硬軌工具機的即時檢測補償裝置及其方法

Also Published As

Publication number Publication date
TWI609737B (zh) 2018-01-01

Similar Documents

Publication Publication Date Title
KR101935148B1 (ko) 선형 안내 홈에 가해진 선-하중의 편차를 추정하는 방법
JP6404893B2 (ja) 工具寿命推定装置
Zhou et al. Tool wear monitoring using acoustic emissions by dominant-feature identification
Shah et al. Tool wear prediction in face milling of stainless steel using singular generative adversarial network and LSTM deep learning models
Guo et al. An integrated wireless vibration sensing tool holder for milling tool condition monitoring with singularity analysis
JP6673216B2 (ja) 要因分析装置、要因分析方法とプログラム、及び、要因分析システム
CN107132064B (zh) 基于多传感器的旋转机械系统运行状态监测方法及系统
JP6673002B2 (ja) 情報処理装置、情報処理方法、情報処理プログラム、および情報処理システム
TW201822948A (zh) 偵測工具機的線性滑軌預壓值變化的方法
ZA201808645B (en) A method of operating a condition monitoring system of a vibrating machine and a condition monitoring system
JP6247627B2 (ja) 異常値検出装置およびその動作方法
WO2021090842A1 (ja) 工作機械および表示装置
JP2017502390A (ja) データに基づく関数モデルを定めるための方法及び装置
US11275005B2 (en) Fatigue limit stress specification system, fatigue limit stress specification device, and fatigue limit stress specification method
JP6511573B1 (ja) 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム
Huang et al. Automatic classification of magnetic tiles internal defects based on acoustic resonance analysis
Zhang et al. Modeling of tool wear for ball end milling cutter based on shape mapping
CN108318158B (zh) 侦测工具机的线性滑轨预压值变化的方法
US11880427B2 (en) Time-series data processing method, corresponding processing system, device and computer program product
WO2017002242A1 (ja) 電力対策支援システム、電力対策支援装置、電力対策支援方法及び電力対策支援プログラム
JP2005051269A5 (zh)
JP2021135513A (ja) 情報処理装置、情報処理方法、およびプログラム
US11921001B2 (en) System for inspecting deviation in dynamic characteristics of linear transmission device of feeding system
TWI789167B (zh) 檢測進給系統之動態特性偏差的方法及系統
JP2020067750A (ja) 学習方法、装置及びプログラム、並びに設備の異常診断方法