一種壓痕數據分析方法
【專利摘要】本發明涉及一種壓痕數據分析方法，該方法的核心是基于Levenberg-Marquardt算法的數值優化，優化的目標函數為通過有限元模擬得到的壓痕數據和實驗壓痕數據之間的差異，優化參數即為被測材料待求的力學參數。由于可在有限元模型中定義多種加-卸載方案和材料本構模型，因此本數據分析方法具有很強的通用性。
【專利說明】一種壓痕數據分析方法
【技術領域】
[0001]本發明涉及壓痕數據分析方法。
【背景技術】
[0002]作為最簡單易行的一種方法，壓痕方法可能是人類最早應用的力學性能測試方法。根據產生壓痕的難易或殘留壓痕面積的大小，人們定義了各種硬度指標以評估材料抵抗永久形變的能力，其中最常用的是壓入硬度，即用一定的載荷將形狀確定的壓頭壓入被測材料，根據壓入區域殘余變形的大小評估該材料的軟硬程度。根據壓頭、載荷以及載荷保持時間的不同，壓入硬度有布氏硬度、洛氏硬度、維氏硬度和顯微硬度等多種類型。
[0003]硬度測試簡單易行，數據可重復性好。然而，硬度只是衡量材料抵抗塑性變形能力的一種綜合指標，它和材料的彈性模量、屈服強度等本構參數之間缺乏明確的換算關系，這導致了硬度測試難以在材料本構關系的表征中獲得直接的應用。隨著精密測量和控制技術的發展，上世紀七八十年代誕生了一種脫胎于顯微硬度測試但精度更高并且可連續記錄壓入力和壓入深度的測試方法，即儀器化壓入(instrumented indentat1n)方法或深度敏感壓痕(depth-sensing indentat1n)方法。由于空間分辨率高(壓頭尺寸可小至微/納米級)，儀器化壓入方法種特別適合于小尺度材料、薄膜及涂層材料系統的力學性能測試。經過近二、三十年的快速發展，該方法已廣泛應用于力學量測量和材料變形機理研究的各個方面，大大拓展了壓痕方法的應用范圍。
[0004]由于壓入過程涉及復雜的三維變形，壓痕測試無法像單軸拉伸或單軸壓縮那樣，可從實驗曲線中直接獲得彈性模量、屈服強度等本構參數。對壓痕測試，由于本構參數和測試數據之間缺乏簡單的對應關系，通常需要應用合適的模型和方法對實驗數據進行分析，才能從測試數據中獲得所需的本構參數，這個過程稱之為壓痕數據的反分析(reverseanalysis)。
[0005]目前商業化的儀器化壓入設備通常采用昂貴的電容式測距傳感器和電磁式的作動裝置。為實現高頻的加載，通常還采用復雜的測量和控制技術。這導致儀器化壓入設備的成本居高不下，嚴重制約了這種結構簡單、用途多樣、空間分辨率高的設備在工業檢測中的應用。目前的壓入設備仍然主要應用實驗室研究。另外，待測量的本構參數不同，采用的壓痕測試方案和壓痕數據反分析方法也往往不同。由于壓痕問題的復雜性，目前還沒有一種較通用的壓痕數據反分析方法，人們往往需要根據待測材料的特點選擇或開發合適的數據分析方法，這無疑也限制了儀器化壓入方法在工業測量中的應用。
【發明內容】
[0006]本發明的目的在于提出一種壓痕數據分析方法，其能解決傳統的分析方法不通用的問題。
[0007]為了達到上述目的，本發明所采用的技術方案如下:
[0008]一種壓痕數據分析方法，其包括以下步驟:
[0009]步驟1、建立壓痕測試過程的有限元模型，其中，被測材料的力學參數定義為變量;
[0010]步驟2、將給定的初始的力學參數輸入至所述有限元模型，并對所述有限元模型進行計算，得到模擬壓痕數據；
[0011]步驟3、調用所述模擬壓痕數據以及讀取實驗壓痕數據，利用公式一計算模擬壓痕數據與實驗壓痕數據之間的差異；
[0012]
【權利要求】
1.一種壓痕數據分析方法，其特征在于，包括以下步驟: 步驟1、建立壓痕測試過程的有限元模型，其中，被測材料的力學參數定義為變量； 步驟2、將給定的初始的力學參數輸入至所述有限元模型，并對所述有限元模型進行計算，得到模擬壓痕數據； 步驟3、調用所述模擬壓痕數據以及讀取實驗壓痕數據，利用公式一計算模擬壓痕數據與實驗壓痕數據之間的差異；
其中，F(Pk)為目標函數的返回值；P°為力學參數的初始值，Pk為修正k次后的力學參數；N為實驗壓痕數據包含的數據個數；rxp(ti)為加載時刻等于\時的實驗壓痕測試結果；fcal(Pk, 為加載時刻等于\時的模擬壓痕測試結果； 步驟4、利用Levenberg-Marquardt算法對目標函數進行優化，當判斷到F (Pk)小于預設閾值時，輸出Pk作為最終的優化結果。
2.如權利要求1所述的壓痕數據分析方法，其特征在于，所述步驟4具體包括如下子步驟: (a)將P°作為Levenberg-Marquardt算法的初始參數； (b)利用有限差分方法計算敏感度矩陣或Jacobian矩陣A= ?(ΡΑ )/5Ρ* ；
(c)求解方程(ΑΤΑ+λI)gk = -AtF (Pk)，得到修正量gk，其中I為單位矩陣，λ為非負的標量參數； (d)計算Pk+1= Pk+gk，并判斷目標函數是否小于預設閾值若是，則輸出Pk作為最終的優化結果，若否，則重復步驟a至步驟d。
3.如權利要求1所述的壓痕數據分析方法，其特征在于，所述有限元模型為軸對稱有限元模型。
【文檔編號】G06F17/50GK104077444SQ201410302066
【公開日】2014年10月1日 申請日期:2014年6月26日 優先權日:2014年6月26日
【發明者】張純禹 申請人:中山大學