SolidWorks截面草圖曲線約束優(yōu)化技術(shù)(四)

2013-07-18 by:廣州Solidworks培訓(xùn)中心來源:仿真在線

2 截面草圖曲線約束優(yōu)化技術(shù)

截面草圖曲線集為STR={str1,str2,…,strQ},其中:strq表示第q段輪廓特征點集;Q表示截面特征單元的分段數(shù)。在確定了各截面特征之間的約束關(guān)系后,即可對各段截面草圖曲線進行幾何約束條件下的整體擬合。以曲線擬合的數(shù)學(xué)表達和約束表達為基礎(chǔ),可以方便地建立截面草圖上全部數(shù)據(jù)點的約束優(yōu)化數(shù)學(xué)模型,實現(xiàn)對任意多段曲線段的整體優(yōu)化求解。

2.1 數(shù)學(xué)模型的建立

1)用strw(w=1,2,…,W)表示第w,段截面輪廓數(shù)據(jù)所對應(yīng)的目標曲線;點pwj表示第w段數(shù)據(jù)中的第j個測量數(shù)據(jù)點(j=1,2,…,J),測量數(shù)據(jù)點pwj到目標曲線strw的距離為d(pwj,strw)。

2)這W段曲線的s維向量用X=[x1,x2,…,xs]表示,它是這W段曲線所有特征參數(shù)的集合。

3)這W段曲線之間滿足:

Ce(X)=0 e=l,2,…,E

式中:Ce(X)為約束集。

則約束優(yōu)化模型可表示為:

SolidWorks截面草圖曲線約束優(yōu)化技術(shù)(四) solidworks仿真分析圖片圖片1

式中:F(x)為測量數(shù)據(jù)點到目標曲線的距離平方和的最小值。

截面輪廓草圖在只有直線段和圓弧段的情況下,幾何約束類型用代數(shù)方程表示比較簡單,當存在樣條曲線段時,其幾何約束類型很難用簡單的代數(shù)方程表示,由于擬合出的B樣條曲線可以是任意的形狀,并且在隨后的約束添加與優(yōu)化求解過程中,B樣條曲線的控制多邊形頂點全部都作為優(yōu)化變量參與優(yōu)化計算,因此截面輪廓草圖含有B樣條曲線時,參與迭代計算的參數(shù)較多,迭代求解的收斂速度也會相應(yīng)地變慢。

2.2 整體優(yōu)化模型的求解

目前,已有優(yōu)化求解的方法多種多樣,各有其特點和適用范圍,而在實際應(yīng)用中,廣泛采用懲罰函數(shù)法將約束優(yōu)化問題變?yōu)榉蔷€性無約束優(yōu)化問題進行求解。本文將式(1)做如下處理,令:

SolidWorks截面草圖曲線約束優(yōu)化技術(shù)(四) solidworks simulation技術(shù)圖片2

式中:G(X,λ)為數(shù)據(jù)點逼近誤差和所有約束值的平方和的最小值;λe為懲罰因子,且λe>0。

式(2)采用Levenberg-Marquardt(L-M)方法迭代求解。迭代求解的初值通過按各段曲線特征進行最小二乘擬合的方法來給定,直線段和圓弧段的數(shù)據(jù)點在識別出曲線類型后直接進行曲線擬合,樣條曲線段的數(shù)據(jù)點則需要對給定數(shù)據(jù)點的最小二乘逼近得到。與Werghi的方法不同,對于式(2)用奇異值分解法(SVD)求解迭代過程中的方程A△X=-Q,其中A=δ2G(X、λ)/δ2X,Q=δG(X、λ)/δX,△X為X的增量。選取合適的初始懲罰因子λe,一般迭代20步左右即可得到滿足約束條件且擬合誤差在最大允許誤差范圍以內(nèi)的截面草圖曲線。具體的L-M迭代過程如圖3所示。增長因子v>1,取v=10,初始參數(shù)α0=0.12。迭代初值X(0)采用各分段曲線基于代數(shù)距離的最小二乘擬合得到,h為迭代次數(shù),△X為迭代步長,H表示與矩陣X相同維數(shù)的單位矩陣,K-T條件為迭代的終止條件。αh和X(h)分別表示迭代h次所對應(yīng)的參數(shù)值。迭代初值X(0)由各分段曲線基于代數(shù)距離的最小二乘擬合得到。

開放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才