丁伯倫 張力娟 周金倩 朱貝貝
(揚州工業(yè)職業(yè)技術(shù)學(xué)院,江蘇 揚州 225127)
以2020年全國大學(xué)生數(shù)學(xué)建模競賽D題第一問為例,根據(jù)某一工件1在水平狀態(tài)下的測量數(shù)據(jù)(數(shù)字信號),來標注出被測工件表面輪廓線的各項參數(shù)值。

圖1 某工件1在水平狀態(tài)下測量的輪廓線
假設(shè)被測工件的輪廓線是由直線和圓弧構(gòu)成的光滑曲線,首先根據(jù)工件1的水平測量數(shù)據(jù)(數(shù)字信號),運用Matlab軟件可得工件x軸與z軸數(shù)字信號的平面圖,如圖2所示:

圖2 工件1在水平狀態(tài)下的數(shù)字信號坐標圖
可設(shè)定z軸初始點為-1.770,那么近似水平狀態(tài)輪廓線的表達式為z0=-1.770,再以第一段傾斜線為例,運用Excel進行數(shù)據(jù)擬合,可得輪廓線D1的表達式為z1=-2.7468x+134.95。采用類似的方法可得每段傾斜線的表達式,如表1所示:

表1 每段傾斜線的表達式
結(jié)合“人字形”兩邊輪廓線的表達式,求出相應(yīng)交點坐標,即可得到“人字形”輪廓線的長度。如表1,已知“人字形”兩邊的輪廓線D7、D8的表達式,可求出輪廓線D7、D8的交點坐標為(76.810,-0.782),運用Excel進行數(shù)據(jù)擬合,可求得輪廓線D7、D8的長度,D7=5.235,D8=5.147,結(jié)合之前水平線z0=-1.770,得到人字線性的高度為0.998。

圖3 工件1在水平狀態(tài)下的數(shù)字信號坐標圖的局部圖
運用Excel軟件擬合出工件1數(shù)字信號的局部圖,利用水平線段的表達公式:z0=-1.770與輪廓線D1的表達公式z1=-2.7468x+134.95求出交點A1的橫坐標為49.772,同理利用水平線段的表達公式z0=-1.770與輪廓線D2的表達公式z2=2.779x-148.05,求出A2的橫坐標為52.676,再用A2的橫坐標減去A1的橫坐標,即可求出槽口寬度x1=2.904;利用水平線段的表達公式z0=-1.770與輪廓線D3的表達式z3=-3.7092x+212.1,求出交點A3的橫坐標為57.653,再利用A3的橫坐標減去A2的橫坐標,即可求出水平線段長度為x2=4.977。以此類推,可求出其他槽口寬度和水平線段長度如表2所示:

表2 槽口寬度和水平線段長度
根據(jù)上述已知條件,求出了各輪廓線的表達公式,由已知的水平線段x0與斜線D1的表達式,得出輪廓線D1斜率kD1=-2.7468,則輪廓線D1與水平線段x0的夾角∠1=arctan(-2.7468),根據(jù)角度的關(guān)系可得∠1=π=arctan(-2.7468)=110.005°;已知水平線段x2與斜線D2的表達式,可得輪廓線D2斜率kD1=2.7779,則輪廓線D2與水平線段x2的夾角∠2為arctan(2.7779),根據(jù)角度的關(guān)系可得∠2=π-arctan(2.7779)=109.798°,同理可得其他斜線與水平線之間的夾角,如表3所示:

表3 斜線與水平線之間的夾角


表4 圓的半徑和圓心距
本文針對接觸式輪廓儀的自動標注問題,利用Matlab軟件與數(shù)據(jù)處理技術(shù),建立特定的圓方程、直線方程模型,合理標注出工件輪廓線的各項參數(shù)值。首先對被測工件的輪廓線假設(shè)成光滑平面曲線以獲取主體輪廓,然后利用Matlab軟件對測量數(shù)據(jù)(數(shù)字信號)進行輪廓提取,最后通過傾斜線,直線和圓弧等信息得出元件各項參數(shù)值。在同一狀態(tài)下的多組測量數(shù)據(jù)的情況下,擬合后得到的結(jié)果會相對精確,具有一定的參考性。