付美鑫

（長(zhǎng)春汽車工業(yè)高等專科學(xué)校，吉林 長(zhǎng)春130000）

全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽創(chuàng)辦于1992 年，隨著賽事的發(fā)展，參加競(jìng)賽的學(xué)校越來(lái)越多，尤其是?？平M。本文將對(duì)2020年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽D 題進(jìn)行分析和求解。

2020 年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽D 題是接觸式輪廓儀的自動(dòng)標(biāo)注的問(wèn)題，主要應(yīng)用最小二乘法的思想，建立了接觸式輪廓儀自動(dòng)標(biāo)注的模型，包括斜線模型、圓弧模型、旋轉(zhuǎn)平移模型等，通過(guò)模型對(duì)工件需要標(biāo)注的數(shù)據(jù)進(jìn)行求解，并對(duì)工件進(jìn)行標(biāo)注。

1 問(wèn)題簡(jiǎn)述

在工業(yè)生產(chǎn)中，工件表面測(cè)量方法分類很多，常見(jiàn)的分為接觸式和非接觸式[1]，其中接觸式輪廓儀，廣泛應(yīng)用于光學(xué)元件、精密五金、機(jī)械加工等工件，在測(cè)量?jī)x器中有著舉足輕重的地位[2]。接觸式輪廓儀在測(cè)量時(shí)，探針隨著被測(cè)工件表面勻速滑行，上下移動(dòng)，在X 和Z 方向分別采樣，并通過(guò)傳感器轉(zhuǎn)化為電信號(hào)，經(jīng)過(guò)計(jì)算機(jī)處理得到被測(cè)工件的表面輪廓圖[3]。現(xiàn)要對(duì)表面輪廓圖進(jìn)行標(biāo)注，主要標(biāo)注的各項(xiàng)參數(shù)值包括：槽口寬度、圓弧半徑、圓心之間的距離、圓弧長(zhǎng)度、水平線段的長(zhǎng)度、斜線段的長(zhǎng)度、斜線與水平線的夾角和人形線的高度。有時(shí)在對(duì)工件做多次測(cè)量時(shí)，由于放置的水平位置和角度的不同導(dǎo)致測(cè)量結(jié)果有明顯的差異，同時(shí)給輪廓圖的標(biāo)注帶來(lái)很大的阻礙，建立相應(yīng)的數(shù)學(xué)模型，完成以下問(wèn)題：

1.1 標(biāo)注工件1 的各項(xiàng)參數(shù)值。

1.2 針對(duì)有角度偏移和水平偏移的工件1 數(shù)據(jù)進(jìn)行校正，對(duì)校正數(shù)據(jù)進(jìn)行重新標(biāo)注，并比較矯正后數(shù)據(jù)與原始數(shù)據(jù)標(biāo)注之間的偏差。

圖1 提取數(shù)據(jù)之后的輪廓圖

1.3 針對(duì)工件2 給出的10 次測(cè)量數(shù)據(jù)，給出每次測(cè)量時(shí)的傾斜角度、并標(biāo)注輪廓線的各項(xiàng)參數(shù)值，畫(huà)出工件2 的完整輪廓圖。

1.4 根據(jù)進(jìn)一步給出的圓和角的測(cè)量數(shù)據(jù)，對(duì)1.3 中的結(jié)論進(jìn)行修正，并完成輪廓線的修正。

2 問(wèn)題分析

由于測(cè)量?jī)x測(cè)量時(shí)在工件上勻速滑行，在X 和Z 方向進(jìn)行采樣，由于采樣間距極短，且是等時(shí)間間隔獲取，導(dǎo)致相鄰很大部分?jǐn)?shù)據(jù)差異不明顯，本文所有模型建立求解所使用的數(shù)據(jù)都是每隔15 個(gè)數(shù)據(jù)進(jìn)行提取一次，獲取新數(shù)據(jù)。通過(guò)EXCEL 畫(huà)圖，可以清晰看到提取的數(shù)據(jù)保持工件原貌，（如圖1 所示）。

3 模型建立與求解

針對(duì)問(wèn)題一：以工件x=-1.77 為水平線，建立坐標(biāo)系，要求解槽口寬度、斜線段長(zhǎng)度、弧長(zhǎng)、圓心之間的距離等，首先要建立斜線模型和圓弧所在圓的模型，然后利用交點(diǎn)，確定相應(yīng)的槽口寬度、斜線段長(zhǎng)度、人形線高度、弧長(zhǎng)、圓心之間的距離等標(biāo)注參數(shù)。其中，建立斜線模型：

zi= kxi+ b( k , b為 參數(shù))；

圖2 輪廓線標(biāo)注圖

在求解圓弧弧長(zhǎng)與半徑時(shí)，基于點(diǎn)到圓心距離等于半徑理論，在工件凹槽處取最低點(diǎn)，以該點(diǎn)為中心，向前向后分別取20個(gè)數(shù)據(jù)，共計(jì)41 個(gè)點(diǎn)，利用點(diǎn)到圓心距離等于半徑建立圓弧模型：

( xi- x0)2+ ( zi- z0)2= r2,( x0, z0, r 為參數(shù))

對(duì)該線性方程組利用最小二乘法求解，即可得到圓弧的圓心及半徑，從而確定圓弧的方程。

利用Mathematica 軟件對(duì)上述斜線模型和圓弧模型進(jìn)行求解，即可確定標(biāo)注參數(shù)。對(duì)工件標(biāo)注如圖2。

針對(duì)問(wèn)題二，根據(jù)旋轉(zhuǎn)角度，與測(cè)量工件放置的位置及角度有變化，要對(duì)給出數(shù)據(jù)進(jìn)行校正，首先要確定旋轉(zhuǎn)角度，建立旋轉(zhuǎn)平移校正模型，先旋轉(zhuǎn)校正，然后進(jìn)行平移校正，利用前面建立的數(shù)學(xué)模型進(jìn)行求解，對(duì)工件進(jìn)行參數(shù)標(biāo)注。

旋轉(zhuǎn)校正模型：

平移校正模型：

針對(duì)問(wèn)題二，做誤差分析如下：

對(duì)校對(duì)之后的工件輪廓圖標(biāo)注參數(shù)和工件1 的原始數(shù)據(jù)的工件輪廓圖標(biāo)注參數(shù)分類進(jìn)行先對(duì)偏差分析（校對(duì)之后數(shù)據(jù)記為*），如下：

通過(guò)上述相對(duì)誤差分析發(fā)現(xiàn)，水平距離的相對(duì)偏差較大，分析原因，主要是由于d9的實(shí)際距離較小，測(cè)量產(chǎn)生偏差屬于合理范圍之內(nèi)，以后測(cè)量需要慎重。圓心間的距離偏差最小，在求解圓弧長(zhǎng)度、圓心距離、半徑時(shí)，主要是基于求解圓弧模型，其處理方式主要采取優(yōu)化處理方法，產(chǎn)生的偏差非常理想，也說(shuō)明了模型的合理性。

針對(duì)問(wèn)題三：前兩問(wèn)基本解決工件自動(dòng)標(biāo)注問(wèn)題，并建立相應(yīng)的數(shù)學(xué)模型，要對(duì)工件2 進(jìn)行標(biāo)注，首先對(duì)所給工件2 測(cè)量的十組數(shù)據(jù)，利用問(wèn)題二的數(shù)學(xué)模型進(jìn)行旋轉(zhuǎn)平移變換之后，對(duì)數(shù)據(jù)進(jìn)行校正整合，擬合成一組數(shù)據(jù)，然后對(duì)獲得的新數(shù)據(jù)，利用問(wèn)題一的模型進(jìn)行參數(shù)標(biāo)注。在對(duì)工件2 進(jìn)行標(biāo)注時(shí)，10 次測(cè)量數(shù)據(jù)，每組數(shù)據(jù)都存在旋轉(zhuǎn)和水平偏差，所以對(duì)10 組數(shù)據(jù)利用問(wèn)題二的旋轉(zhuǎn)校正模型和水平校正模型，對(duì)數(shù)據(jù)進(jìn)行校正。首先確定每組數(shù)據(jù)旋轉(zhuǎn)角度。通過(guò)問(wèn)題二中旋轉(zhuǎn)角計(jì)算模型，得到10 組數(shù)據(jù)旋轉(zhuǎn)角如下表所示：

表1 工件2 的10 組數(shù)據(jù)旋轉(zhuǎn)角

通過(guò)對(duì)工件2 的10 組數(shù)據(jù)進(jìn)行分析，對(duì)旋轉(zhuǎn)平移校正得到的數(shù)據(jù)，進(jìn)行補(bǔ)充，主要補(bǔ)充缺失數(shù)據(jù)信息，得到工件2 的大致輪廓圖，并對(duì)工件做簡(jiǎn)單標(biāo)注。針對(duì)問(wèn)題四：工件2 在問(wèn)題三中已經(jīng)基本標(biāo)注完成，但是圓、角的標(biāo)注還很粗糙，所以針對(duì)問(wèn)題四，將對(duì)圓、角數(shù)據(jù)進(jìn)行進(jìn)一步的處理，然后對(duì)問(wèn)題三的輪廓線進(jìn)行修正。

4 結(jié)論

本文針對(duì)2020 年數(shù)學(xué)建模競(jìng)賽D 題進(jìn)行分析求解，從模型求解結(jié)果和誤差分析可以看出，該模型能夠很好的對(duì)接觸式輪廓儀測(cè)量的工件進(jìn)行自動(dòng)標(biāo)注，并標(biāo)注參數(shù)誤差較小，該模型是有效的。由于時(shí)間問(wèn)題以及技術(shù)硬件要求不足等問(wèn)題，對(duì)于問(wèn)題三和問(wèn)題四，只做了數(shù)據(jù)的相關(guān)處理，在后續(xù)會(huì)進(jìn)一步對(duì)其進(jìn)行標(biāo)注。

該模型不僅可以解決針對(duì)本文的問(wèn)題，還可以應(yīng)用于非球面接觸式輪廓儀的自動(dòng)標(biāo)注、物體表面樣貌的分析等其它測(cè)量的自動(dòng)標(biāo)注問(wèn)題。