《外文翻譯--關(guān)于評價特征去除所導致的工程分析錯誤的規(guī)范理論中文版》由會員分享,可在線閱讀,更多相關(guān)《外文翻譯--關(guān)于評價特征去除所導致的工程分析錯誤的規(guī)范理論中文版(10頁珍藏版)》請在裝配圖網(wǎng)上搜索。
1、關(guān)于評價特征去除所導致的工程分析錯誤的規(guī)范理論
SankaraHariGopalakrishnan,KrishnanSuresh
機械工程系,威斯康辛大學,麥迪遜分校,2006年9月30日
摘要:幾何分析是著名的計算機輔助設(shè)計/計算機輔助工藝簡化 “小或無關(guān)特征”在CAD模型中的程序,如有限元分析。然而,幾何分析不可避免地會產(chǎn)生分析錯誤,在目前的理論框架實在不容易量化。
本文中,我們對快速計算處理這些幾何分析錯誤提供了嚴謹?shù)睦碚?。尤其,我們集中力量解決地方的特點,被簡化的任意形狀和大小的區(qū)域。提出的理論采用伴隨矩陣制定邊值問題抵達嚴格界限幾何分析性分析錯誤。該理論通過數(shù)值例子說明。
2、關(guān)鍵詞:幾何分析;工程分析;誤差估計;計算機輔助設(shè)計/計算機輔助教學
1. 介紹
機械零件通常包含了許多幾何特征。不過,在工程分析中并不是所有的特征都是至關(guān)重要的。以前的分析中無關(guān)特征往往被忽略,從而提高自動化及運算速度。
舉例來說,考慮一個剎車轉(zhuǎn)子,如圖1(a)。轉(zhuǎn)子包含50多個不同的特征,但所有這些特征并不是都是相關(guān)的。就拿一個幾何化的剎車轉(zhuǎn)子的熱量分析來說,如圖1(b)。有限元分析的全功能的模型如圖1(a),需要超過150,000度的自由度,幾何模型圖1(b)項要求小于25,000個自由度,從而導致非常緩慢的運算速度。
圖1(a)剎車轉(zhuǎn)子
3、 圖1(b)其幾何分析版本
除了提高速度,通常還能增加自動化水平,這比較容易實現(xiàn)自動化的有限元網(wǎng)格幾何分析組成。內(nèi)存要求也跟著降低,而且條件數(shù)離散系統(tǒng)將得以改善;后者起著重要作用迭代線性系統(tǒng)。
但是,幾何分析還不是很普及。不穩(wěn)定性到底是“小而局部化”還是“大而擴展化”,這取決于各種因素。例如,對于一個熱問題,想刪除其中的一個特征,不穩(wěn)定性是一個局部問題:(1)凈熱通量邊界的特點是零。(2)特征簡化時沒有新的熱源產(chǎn)生; [4]對上述規(guī)則則例外。展示這些物理特征被稱為自我平衡。結(jié)果,同樣存在結(jié)構(gòu)上的問題。
從幾何分析角度看,如果特征遠離該區(qū)域,則這種自我平衡的特征可以忽略。但是,如果功能
4、接近該區(qū)域我們必須謹慎,。
從另一個角度看,非自我平衡的特征應(yīng)值得重視。這些特征的簡化理論上可以在系統(tǒng)任意位置被施用,但是會在系統(tǒng)分析上構(gòu)成重大的挑戰(zhàn)。
目前,尚無任何系統(tǒng)性的程序去估算幾何分析對上述兩個案例的潛在影響。這就必須依靠工程判斷和經(jīng)驗。
在這篇文章中,我們制定了理論估計幾何分析影響工程分析自動化的方式。任意形狀和大小的形體如何被簡化是本文重點要解決的地方。伴隨矩陣和單調(diào)分析這兩個數(shù)學概念被合并成一個統(tǒng)一的理論來解決雙方的自我平衡和非自我平衡的特點。數(shù)值例子涉及二階scalar偏微分方程,以證實他的理論。
本文還包含以下內(nèi)容。第二節(jié)中,我們就幾何分析總結(jié)以往的工作。在第三節(jié)中
5、,我們解決幾何分析引起的錯誤分析,并討論了擬議的方法。第四部分從數(shù)值試驗提供結(jié)果。第五部分討論如何加快設(shè)計開發(fā)進度。
2. 前期工作
幾何分析過程可分為三個階段:
識別:哪些特征應(yīng)該被簡化;
簡化:如何在一個自動化和幾何一致的方式中簡化特征;
分析:簡化的結(jié)果。
第一個階段的相關(guān)文獻已經(jīng)很多。例如,企業(yè)的規(guī)模和相對位置這個特點,經(jīng)常被用來作為度量鑒定。此外,也有人提議以有意義的力學判據(jù)確定這種特征。
自動化幾何分析過程,事實上,已成熟到一個商業(yè)化幾何分析的地步。但我們注意到,這些商業(yè)軟件包僅提供一個純粹的幾何解決。因為沒有保證隨后進行的分析錯誤,所以必須十分小心使用。另外,固有的
6、幾何問題依然存在,并且還在研究當中。
本文的重點是放在第三階段,即快速幾何分析。建立一個有系統(tǒng)的方法,通過幾何分析引起的誤差是可以計算出來的。再分析的目的是迅速估計改良系統(tǒng)的反應(yīng)。其中最著名的再分析理論是著名的謝爾曼-Morrison和woodbury公式。對于兩種有著相似的網(wǎng)狀結(jié)構(gòu)和剛度矩陣設(shè)計,再分析這種技術(shù)特別有效。然而,過程幾何分析在網(wǎng)狀結(jié)構(gòu)的剛度矩陣會導致一個戲劇性的變化,這與再分析技術(shù)不太相關(guān)。
3. 擬議的方法
3.1問題闡述
我們把注意力放在這個文件中的工程問題,標量二階偏微分方程式(pde):
許多工程技術(shù)問題,如熱,流體靜磁等問題,可能簡化為上述公式。
作為一
7、個說明性例子,考慮散熱問題的二維模塊Ω如圖2所示。
圖2二維熱座裝配
熱量q從一個線圈置于下方位置列為Ωcoil。半導體裝置位于Ωdevice。這兩個地方都屬于Ω,有相同的材料屬性,其余Ω將在后面討論。特別令人感興趣的是數(shù)量,加權(quán)溫度Tdevice內(nèi)Ωdevice(見圖2)。一個時段,認定為Ωslot縮進如圖2,會受到抑制,其對Tdevice將予以研究。邊界的時段稱為Γslot其余的界線將稱為Γ。邊界溫度Γ假定為零。兩種可能的邊界條件Γslot被認為是:(a)固定熱源,即(-kt)?n=q,(b)有一定溫度,即T=Tslot。兩種情況會導致兩種不同幾何分析引起的誤差的結(jié)果。
設(shè)T(x
8、,y)是未知的溫度場和K導熱。然后,散熱問題可以通過泊松方程式表示:
其中H(x,y)是一些加權(quán)內(nèi)核?,F(xiàn)在考慮的問題是幾何分析簡化的插槽是簡化之前分析,如圖3所示。
圖3defeatured二維熱傳導裝配模塊
現(xiàn)在有一個不同的邊值問題,不同領(lǐng)域t(x,y):
觀察到的插槽的邊界條件為t(x,y)已經(jīng)消失了,因為槽已經(jīng)不存在了(關(guān)鍵性變化)!
解決的問題是:
設(shè)定tdevice和t(x,y)的值,估計Tdevice。
這是一個較難的問題,是我們尚未解決的。在這篇文章中,我們將從上限和下限分析Tdevice。這些方向是明確被俘引理3、4和3、6。至于其余的這一節(jié),
9、我們將發(fā)展基本概念和理論,建立這兩個引理。值得注意的是,只要它不重疊,定位槽與相關(guān)的裝置或熱源沒有任何限制。上下界的Tdevice將取決于它們的相對位置。
3.2伴隨矩陣方法
我們需要的第一個概念是,伴隨矩陣公式表達法。應(yīng)用伴隨矩陣論點的微分積分方程,包括其應(yīng)用的控制理論,形狀優(yōu)化,拓撲優(yōu)化等。我們對這一概念歸納如下。
相關(guān)的問題都可以定義為一個伴隨矩陣的問題,控制伴隨矩陣t_(x,y),必須符合下列公式計算〔23〕:
伴隨場t_(x,y)基本上是一個預定量,即加權(quán)裝置溫度控制的應(yīng)用熱源??梢杂^察到,伴隨問題的解決是復雜的原始問題;控制方程是相同的;這些問題就是所謂的自身伴隨矩陣。
10、大部分工程技術(shù)問題的實際利益,是自身伴隨矩陣,就很容易計算伴隨矩陣。
另一方面,在幾何分析問題中,伴隨矩陣發(fā)揮著關(guān)鍵作用。表現(xiàn)為以下引理綜述:
引理3.1已知和未知裝置溫度的區(qū)別,即(Tdevice-tdevice)可以歸納為以下的邊界積分比幾何分析插槽:
在上述引理中有兩點值得注意:
1、積分只牽涉到邊界гslot;這是令人鼓舞的?;蛟S,處理剛剛過去的被簡化信息特點可以計算誤差。
2、右側(cè)牽涉到的未知區(qū)域T(x,y)的全功能的問題。特別是第一周期涉及的差異,在正常的梯度,即涉及[-k(T-t)] ?n;這是一個已知數(shù)量邊界條件[-kt]?n所指定的時段,未知狄里克萊條件作出規(guī)定
11、[-kt]?n可以評估。在另一方面,在第二個周期內(nèi)涉及的差異,在這兩個領(lǐng)域,即T管; 因為t可以評價,這是一個已知數(shù)量邊界條件T指定的時段。因此。
引理3.2、差額(tdevice-tdevice)不等式
然而,伴隨矩陣技術(shù)不能完全消除未知區(qū)域T(x,y)。為了消除T(x,y)我們把重點轉(zhuǎn)向單調(diào)分析。
3.3單調(diào)性分析
單調(diào)性分析是由數(shù)學家在19世紀和20世紀前建立的各種邊值問題。例如,一個單調(diào)定理:
"添加幾何約束到一個結(jié)構(gòu)性問題,是指在位移(某些)邊界不減少"。
觀察發(fā)現(xiàn),上述理論提供了一個定性的措施以解決邊值問題。
后來,工程師利用之前的“計算機時代”上限或下限同樣的定
12、理,解決了具有挑戰(zhàn)性的問題。當然,隨著計算機時代的到來,這些相當復雜的直接求解方法已經(jīng)不為人所用。但是,在當前的幾何分析,我們證明這些定理采取更為有力的作用,尤其應(yīng)當配合使用伴隨理論。
我們現(xiàn)在利用一些單調(diào)定理,以消除上述引理T(x,y)。遵守先前規(guī)定,右邊是區(qū)別已知和未知的領(lǐng)域,即T(x,y)-t(x,y)。因此,讓我們在界定一個領(lǐng)域E(x,y)在區(qū)域為:
e(x,y)=t(x,y)-t(x,y)。
據(jù)悉,T(x,y)和T(x,y)都是明確的界定,所以是e(x,y)。事實上,從公式(1)和(3),我們可以推斷,e(x,y)的正式滿足邊值問題:
解決上述問題就能解決所有問題。但是,
13、如果我們能計算區(qū)域e(x,y)與正常的坡度超過插槽,以有效的方式,然后(Tdevice-tdevice),就評價表示e(X,Y)的效率,我們現(xiàn)在考慮在上述方程兩種可能的情況如(a)及(b)。
例(a)邊界條件較第一插槽,審議本案時槽原本指定一個邊界條件。為了估算e(x,y),考慮以下問題:
因為只取決于縫隙,不討論域,以上問題計算較簡單。經(jīng)典邊界積分/邊界元方法可以引用。關(guān)鍵是計算機領(lǐng)域e1(x,y)和未知領(lǐng)域的e(x,y)透過引理3.3。這兩個領(lǐng)域e1(x,y)和e(x,y)滿足以下單調(diào)關(guān)系:
把它們綜合在一起,我們有以下結(jié)論引理。
引理3.4未知的裝置溫度Tdevice,當
14、插槽具有邊界條件,東至以下限額的計算,只要求:(1)原始及伴隨場T和隔熱與幾何分析域(2)解決e1的一項問題涉及插槽:
觀察到兩個方向的右側(cè),雙方都是獨立的未知區(qū)域T(x,y)。
例(b) 插槽Dirichlet邊界條件
我們假定插槽都維持在定溫Tslot??紤]任何領(lǐng)域,即包含域和插槽。界定一個區(qū)域e(x,y)在滿足:
現(xiàn)在建立一個結(jié)果與e-(x,y)及e(x,y)。
引理3.5
注意到,公式(7)的計算較為簡單。這是我們最終要的結(jié)果。
引理3.6 未知的裝置溫度Tdevice,當插槽有Dirichlet邊界條件,東至以下限額的計算,只要求:(1)原始及伴隨場T和
15、隔熱與幾何分析。(2) 圍繞插槽解決失敗了的邊界問題,:
再次觀察這兩個方向都是獨立的未知領(lǐng)域T(x,y)。
4. 數(shù)值例子說明
我們的理論發(fā)展,在上一節(jié)中,通過數(shù)值例子。設(shè)
k = 5W/m?C, Q = 10 W/m3 and H = 。
表1:結(jié)果表
表1給出了不同時段的邊界條件。第一裝置溫度欄的共同溫度為所有幾何分析模式(這不取決于插槽邊界條件及插槽幾何分析)。接下來兩欄的上下界說明引理3.4和3.6。最后一欄是實際的裝置溫度所得的全功能模式(前幾何分析),是列在這里比較前列的。在全部例子中,我們可以看到最后一欄則是介于第二和第三列。T Tdevice
16、T
對于絕緣插槽來說,Dirichlet邊界條件指出,觀察到的各種預測為零。不同之處在于這個事實:在第一個例子,一個零Neumann邊界條件的時段,導致一個自我平衡的特點,因此,其對裝置基本沒什么影響。另一方面,有Dirichlet邊界條件的插槽結(jié)果在一個非自我平衡的特點,其缺失可能導致器件溫度的大變化在。
不過,固定非零槽溫度預測范圍為20度到0度。這可以歸因于插槽溫度接近于裝置的溫度,因此,將其刪除少了影響。
的確,人們不難計算上限和下限的不同Dirichlet條件插槽。圖4說明了變化的實際裝置的溫度和計算式。
預測的上限和下限的實際溫度裝置表明理論是正確的。另外,跟預期結(jié)果一樣,
17、限制槽溫度大約等于裝置的溫度。
5. 快速分析設(shè)計的情景
我們認為對所提出的理論分析"什么-如果"的設(shè)計方案,現(xiàn)在有著廣泛的影響。研究顯示設(shè)計如圖5,現(xiàn)在由兩個具有單一熱量能源的器件。如預期結(jié)果兩設(shè)備將不會有相同的平均溫度。由于其相對靠近熱源,該裝置的左邊將處在一個較高的溫度,。
圖4估計式versus插槽溫度圖
圖5雙熱器座
圖6正確特征可能性位置
為了消除這種不平衡狀況,加上一個小孔,固定直徑;五個可能的位置見圖6。兩者的平均溫度在這兩個地區(qū)最低。
強制進行有限元分析每個配置。這是一個耗時的過程。另一種方法是把該孔作為一個特征,并研究其影響,作為后處理步驟。換言之,這是一個特殊的“幾何分析”例子,而擬議的方法同樣適用于這種情況。我們可以解決原始和伴隨矩陣的問題,原來的配置(無孔)和使用的理論發(fā)展在前兩節(jié)學習效果加孔在每個位置是我們的目標。目的是在平均溫度兩個裝置最大限度的差異。表2概括了利用這個理論和實際的價值。
從上表可以看到,位置W是最佳地點,因為它有最低均值預期目標的功能。