本發(fā)明涉及計算機輔助工程分析，更具體地說，涉及用于數(shù)值模擬產(chǎn)品(例如汽車、飛機等)的結(jié)構(gòu)性能的方法和系統(tǒng)，該方法和系統(tǒng)采用具有時間步長控制方案的顯式有限元分析(FEA)。

背景技術(shù)：

有限元分析(FEA)是一種計算機實現(xiàn)方法，廣泛用于工業(yè)建模和解決工程問題，涉及復雜系統(tǒng)，如三維非線性結(jié)構(gòu)設(shè)計和分析。FEA的名稱來自于所考慮的對象的幾何形狀的指定方式。隨著現(xiàn)代數(shù)字計算機的出現(xiàn)，F(xiàn)EA已經(jīng)被實現(xiàn)為FEA軟件。基本上，F(xiàn)EA軟件提供幾何描述的模型以及在模型中每個點的相關(guān)的材料屬性。在這模型中，要分析的系統(tǒng)的幾何形狀由各種尺寸的實體、殼體和梁表示，稱之為單元。單元的頂點被稱為節(jié)點。模型包括有限數(shù)量的單元，其被分配材料標識以關(guān)聯(lián)其材料屬性。模型因此表示由要分析的對象沿其周圍的環(huán)境所占據(jù)的物理空間。FEA軟件然后引用一表格，其中每種材料的屬性(例如應力-應變本構(gòu)方程、楊氏模量、泊松比、熱導率)被制成表格。此外，對象的邊界條件(即載荷、物理約束等)被指定。對象及其環(huán)境的模型通過這種方式創(chuàng)建。

FEA有兩種求解技術(shù)：隱式有限元分析(“隱式方法”)和顯式有限元分析(“顯式方法”)。兩種方法都用來求解運動的瞬態(tài)動力學方程，從而獲得方程的均衡解。該方法從時間(t)經(jīng)過離散時間間隔或時間步長Δt推進到時間(t+Δt)。時間步長Δt是在每個求解周期中確定(即計算)的，用于下一求解周期。這種方法有時被稱為時間推進模擬，它包括了若干連續(xù)的時間步長或求解周期。

本發(fā)明涉及顯式方法，其僅在時間步長非常小時是穩(wěn)定的，特別是，時間步長必須小于彈性波從單元的一側(cè)傳播到另一側(cè)所用的時間。顯式方法中維持穩(wěn)定求解的最大時間步長稱之為臨界時間步長Δt_cr。彈性波的速度是材料質(zhì)量和結(jié)構(gòu)剛度的函數(shù)，由有限元和單元的大小或尺寸表示。對于具有基本相似的材料的FEA模型，最小的單元通常控制著臨界時間步長。

在FEA模型中即使一個較小的單元也能導致臨界時間步長對于FEA模型中的多數(shù)單元而言不必要的小。其結(jié)果是，需要將一個非常小的時間步長用于其余的數(shù)值模擬。模擬不僅將變得非常耗時，而且也不切實際。例如，它需要很長時間才能得到工程師用來及時做出滿足商業(yè)需求的設(shè)計決策的模擬結(jié)果。

當采用具有一較小尺寸的厚殼有限元或單層實體/結(jié)構(gòu)有限元時，這個問題會變得更糟。厚殼有限元的厚度將控制時間步長，這會導致時間推進數(shù)值模擬消耗太多的計算機資源從而不符合實際生產(chǎn)。

這將因此需要一個時間步長控制方案以確保高效的顯式FEA以獲得產(chǎn)品的模擬結(jié)構(gòu)性能。

技術(shù)實現(xiàn)要素：

本發(fā)明公開了一種產(chǎn)品的具有時間步長控制方案的高效顯式有限元分析。根據(jù)一示例性實施例，在安裝有顯式FEA應用模塊的計算機系統(tǒng)中定義并接收表示產(chǎn)品或結(jié)構(gòu)的有限元分析(FEA)模型。FEA模型包括至少第一層和第二層相應節(jié)點以形成一組節(jié)點對。每一節(jié)點對包括位于第一層的一節(jié)點和位于第二層的一相應的節(jié)點。在每一節(jié)點對中的節(jié)點位于與離FEA模型中其他節(jié)點的距離相比基本上距離彼此更近的位置。一組候選有限元被確定用于控制時間步長。每一候選有限元由一個或多個節(jié)點對限定。

通過在多個求解周期中采用FEA模型運行時間推進數(shù)值模擬以獲得產(chǎn)品的模擬結(jié)構(gòu)性能。在每一求解周期中，包括所有候選有限元的FEA模型的節(jié)點力根據(jù)顯式FEA方案得到?？s放尺寸由候選有限元的最小尺寸乘以加速縮放因子(例如大于1)計算出。然后，縮放尺寸與其他尺寸相比較以確定用于控制下一求解周期的時間步長的臨界尺寸。根據(jù)基于節(jié)點質(zhì)量、節(jié)點力和加速縮放因子的公式，為了維持穩(wěn)定求解而重新分配的相應的節(jié)點力的分量被計算并重新分配。

用于下一求解周期的新的時間步長采用臨界尺寸(例如可以是或不是縮放尺寸)計算出來。模擬時間增加計算出的新時間步長以在新的求解周期內(nèi)重復上述步驟直到達到時間推進數(shù)值模擬的結(jié)束條件。

本發(fā)明的其他目的、特征和優(yōu)點是顯而易見的，在仔細閱讀以下實施例的詳細描述，并結(jié)合附圖。

附圖說明

本發(fā)明的這些和其他特征、方面以及優(yōu)點會被更好的理解，對于以下敘述、所附權(quán)利要求以及如下附圖：

圖1A-1B共同示出了采用具有時間步長控制方案的顯式有限元分析(FEA)的產(chǎn)品的數(shù)值模擬結(jié)構(gòu)性能的方法的流程圖，根據(jù)本發(fā)明的一實施例；

圖2是圖示了根據(jù)本發(fā)明一實施例的第一層節(jié)點和第二層節(jié)點的示意圖；

圖3A是示出了根據(jù)本發(fā)明一實施例的厚殼有限元或單層實體有限元的示意圖；

圖3B是示出了根據(jù)本發(fā)明一實施例的梁/桁架有限元的示意圖；

圖4A是示出了根據(jù)本發(fā)明一實施例的相比距離其他節(jié)點的距離更近距離的節(jié)點對的示意圖；

圖4B是示出了根據(jù)本發(fā)明一實施例的采用加速縮放因子之后用于時間步長計算的確定臨界尺寸的比較的示意圖；

圖5是示出了根據(jù)本發(fā)明一實施例的在時間步長控制方案中節(jié)點力重新分配的示意圖；

圖6是示出了本發(fā)明的一個實施例在其中實現(xiàn)的計算機示例的主要部分的功能圖。

具體實施方式

圖1A-1B共同示出了采用具有時間步長控制方案的顯式有限元分析(FEA)的產(chǎn)品的數(shù)值模擬結(jié)構(gòu)性能的方法100的流程圖，根據(jù)本發(fā)明的一實施例。方法100優(yōu)選結(jié)合其他附圖理解并在軟件中實現(xiàn)。

方法100由步驟102開始，在計算機系統(tǒng)(例如圖6所示的計算系統(tǒng)600)中接收表示產(chǎn)品或結(jié)構(gòu)(例如汽車、飛機、汽車的面板等)的有限元分析(FEA)模型。FEA模型包括多個節(jié)點，節(jié)點包括至少第一層和第二層相應節(jié)點，以形成一組節(jié)點對。每一節(jié)點對包括位于第一層的一節(jié)點和位于第二層的一相應的節(jié)點。在每一節(jié)點對中的節(jié)點位于與離FEA模型中其他節(jié)點的距離相比基本上距離彼此更近的位置。在一例子中，較近距離比其他距離小5-10倍。圖2示出了可以在根據(jù)本發(fā)明一實施例的FEA模型中使用的節(jié)點的兩個示例層。節(jié)點的第一層211(如空心圓所示)和節(jié)點的第二層212(如實心圓所示)包括相應的節(jié)點以形成節(jié)點對。每一節(jié)點對包括兩個節(jié)點，例如節(jié)點N1 201(空心圓)和N2 202(實心圓)。

接著，在步驟104，一組候選有限元被確定用于控制時間步長(即用于計算下一求解周期的時間步長的最小尺寸)。候選有限元由一個或多個節(jié)點對限定。圖3A示出了三個厚殼有限元或單層實體有限元301-303，根據(jù)本發(fā)明一實施例。每一厚殼或?qū)嶓w有限元由四個節(jié)點對限定。圖3B示出了一對梁/桁架有限元311-312，每個由一個節(jié)點對限定，根據(jù)本發(fā)明一實施例。

有許多技術(shù)來確定這樣的組。例如，用戶可以在FEA的輸入指定一個或多個組。它也可以由通過檢查每個有限元的初始尺寸和/或拓撲結(jié)構(gòu)來確定的自動化過程完成。

接著，在步驟108，通過采用具有顯式FEA應用模塊的FEA模型運行時間推進數(shù)值模擬以獲得產(chǎn)品的模擬結(jié)構(gòu)性能。

時間推進模擬包括多個求解周期或時間步驟，通過在步驟112中設(shè)置模擬時間為零來啟動。

在步驟114，包括所有候選有限元的FEA模型的節(jié)點力通過顯式FEA方案得到。在步驟116，每一候選有限元的縮放尺寸由最小尺寸乘以加速縮放因子(例如大于1的因子)計算出。通常，最小尺寸是每個節(jié)點對之間的較近的距離。

在步驟118，縮放尺寸與其他尺寸(例如厚殼有限元或?qū)嶓w有限元的長度和寬度)相比較以確定用于控制下一求解周期的時間步長的臨界尺寸。圖4B為示出了這種比較的示意圖。等于厚度T(例如T 410)乘以加速縮放因子S的縮放尺寸411(S×T)與其他候選有限元的尺寸(例如寬度W 412和長度L 414)相比較。在該例子中，縮放尺寸仍然是臨界尺寸，因為在所有尺寸中它是最小尺寸。

此外，圖4A為示出了由節(jié)點N1 401和節(jié)點N2 402形成的節(jié)點對的示意圖，根據(jù)本發(fā)明的一實施例。W 412和L 414為FEA模型中距離其他節(jié)點的距離，而T 410是節(jié)點對(N1-N2)的節(jié)點401-402之間的較近的距離。

候選有限元的每個節(jié)點對(例如節(jié)點對N1-N2)的節(jié)點力的分量(ΔF)根據(jù)基于加速縮放因子、節(jié)點力和節(jié)點質(zhì)量的以下公式計算。節(jié)點力的分量在每一節(jié)點對中的節(jié)點之間重新分配以確保在下一求解周期可以實現(xiàn)穩(wěn)定求解。值得注意的是，下一求解周期由更大的時間步長增加，取決于縮放尺寸，當縮放尺寸為臨界尺寸時。

ΔF＝β₁F₁-β₂F₂

F₁^final＝F₁-ΔF

F₂^final＝F₂+ΔF

其中，S是加速縮放因子，m₁和m₂分別是每一節(jié)點對的節(jié)點質(zhì)量，F(xiàn)₁和F₂分別是每一節(jié)點對的節(jié)點力，ΔF是計算出的相應的節(jié)點力的分量，F(xiàn)₁^final和F₂^final分別是每一節(jié)點對的重新分配后的節(jié)點力。

在許多情況下，顯式FEA方案需要節(jié)點質(zhì)量的倒數(shù)用于計算。當節(jié)點質(zhì)量倒數(shù)都存儲在計算機系統(tǒng)中時，兩個系數(shù)可以由如下的倒數(shù)量計算：

合并中間步驟，相應的節(jié)點力的分量可以寫為如下:

節(jié)點力的重新分配方案示例如圖5所示。節(jié)點力F₁ 511和節(jié)點力F₂ 512分別通過用于節(jié)點N1 511和節(jié)點N2 512的顯式FEA獲得。在節(jié)點力重新分配之后，各節(jié)點力分別變?yōu)镕₁^final＝F₁-ΔF 521和F₂^final＝F₂+ΔF 522。

接著，在步驟122，新的臨界時間步長用臨界尺寸計算出。縮放尺寸應當是最小的，時間步長將由加速縮放因子縮放。例如，當5作為加速縮放因子，時間步長將比正常顯式FEA方案大5倍。為了補償更大和可能不穩(wěn)定的時間步長，候選有限元的節(jié)點對中的節(jié)點力相應地重新分配。

接著，在步驟124，確定時間推進數(shù)值模擬是否達到結(jié)束條件。如果沒有，在步驟126中，模擬時間增加計算出的時間步長。方法100重復步驟114-124直到結(jié)束條件達到。結(jié)束條件可以是由用戶指定的總的模擬時間。

根據(jù)一方面，本發(fā)明涉及一個或多個能夠執(zhí)行在此描述的功能的計算機系統(tǒng)。計算機系統(tǒng)600的例子在圖6中示出。計算機系統(tǒng)600包括一個或多個處理器，例如處理器604。處理器604連接到計算機系統(tǒng)內(nèi)部通信總線602。關(guān)于該示范性的計算機系統(tǒng)，有各種軟件實現(xiàn)的描述。在讀完這一描述后，相關(guān)技術(shù)領(lǐng)域的人員將會明白如何使用其它計算機系統(tǒng)和/或計算機架構(gòu)來實施本發(fā)明。

計算機系統(tǒng)600還包括主存儲器608，優(yōu)選隨機存取存儲器(RAM)，還可包括輔助存儲器610。輔助存儲器610包括例如一個或多個硬盤驅(qū)動器612和/或一個或多個可移除存儲驅(qū)動器614，它們代表軟盤驅(qū)動器、磁帶驅(qū)動器、光盤驅(qū)動器等。可移除存儲驅(qū)動器614用已知的方式從可移除存儲單元618中讀取和/或向可移除存儲單元618中寫入。可移除存儲單元618代表可以由可移除存儲驅(qū)動器614讀取和寫入的軟盤、磁帶、光盤等?？梢岳斫猓梢瞥鎯卧?18包括其上存儲有計算機軟件和/或數(shù)據(jù)的計算機可讀媒介。

在可選實施例中，輔助存儲器610可包括其它類似的機制，允許計算機程序或者其它指令被裝載到計算機系統(tǒng)600。這樣的機制包括例如可移除存儲單元622和接口620。這樣的例子可包括程序盒式存儲器和盒式存儲器接口(例如，視頻游戲設(shè)備中的那些)、可移除存儲芯片(例如可擦除可編程只讀存儲器(EPROM))、通用串行總線(USB)閃存、或者PROM)以及相關(guān)的插槽、以及其它可移除存儲單元622和允許軟件和數(shù)據(jù)從可移除存儲單元622傳遞到計算機系統(tǒng)600的接口620。通常，計算機系統(tǒng)600由操作系統(tǒng)(OS)軟件控制和管理，操作系統(tǒng)執(zhí)行例如進程調(diào)度、存儲器管理、網(wǎng)絡連接和I/O服務。

可能還設(shè)有連接到總線602的通信接口624。通信接口624允許軟件和數(shù)據(jù)在計算機系統(tǒng)600和外部設(shè)備之間傳遞。通信接口624的例子包括調(diào)制解調(diào)器、網(wǎng)絡接口(例如以太網(wǎng)卡)、通信端口、個人計算機存儲卡國際協(xié)會(PCMCIA)插槽和卡等等。計算機600基于一組特定的規(guī)則(也就是，協(xié)議)通過數(shù)據(jù)網(wǎng)絡與其它計算設(shè)備通信。通用協(xié)議的其中一種是在互聯(lián)網(wǎng)中通用的TCP/IP(傳輸控制協(xié)議/互聯(lián)網(wǎng)協(xié)議)。通常，通信接口624將數(shù)據(jù)文件組合處理成較小的數(shù)據(jù)包以通過數(shù)據(jù)網(wǎng)絡傳輸，或?qū)⒔邮盏降臄?shù)據(jù)包重新組合成原始的數(shù)據(jù)文件。此外，通信接口624處理每個數(shù)據(jù)包的地址部分以使其到達正確的目的地，或者中途截取發(fā)往計算機600的數(shù)據(jù)包。在這份文件中，術(shù)語“計算機程序媒介”和“計算機可用媒介”都用來指代媒介，例如可移除存儲驅(qū)動器614和/或設(shè)置在硬盤驅(qū)動器612中的硬盤。這些計算機程序產(chǎn)品是用于將軟件提供給計算機系統(tǒng)600的手段。本發(fā)明涉及這樣的計算機程序產(chǎn)品。

計算機系統(tǒng)600還包括輸入/輸出(I/O)接口630，它使得計算機系統(tǒng)600能夠接入顯示器、鍵盤、鼠標、打印機、掃描儀、繪圖儀、以及類似設(shè)備。

計算機程序(也被稱為計算機控制邏輯)作為應用模塊606存儲在主存儲器608和/或輔助存儲器610中。也可通過通信接口624接收計算機程序。這樣的計算機程序被執(zhí)行時，使得計算機系統(tǒng)600執(zhí)行如在此所討論的本發(fā)明的特征。特別地，當執(zhí)行該計算機程序時，使得處理器604執(zhí)行本發(fā)明的特征。因此，這樣的計算機程序代表計算機系統(tǒng)600的控制器。

在本發(fā)明采用軟件實現(xiàn)的實施例中，該軟件可存儲在計算機程序產(chǎn)品中，并可使用可移除存儲驅(qū)動器614、硬盤驅(qū)動器612、或者通信接口624加載到計算機系統(tǒng)600中。應用模塊606被處理器604執(zhí)行時，使得處理器604執(zhí)行如在此所述的本發(fā)明的功能。

主存儲器608可被加載有一個或多個應用模塊606，所述應用模塊606可被一個或多個處理器604執(zhí)行以實現(xiàn)期望的任務，所述處理器可具有或不具有通過I/O接口630輸入的用戶輸入。在運行中，當至少一個處理器604執(zhí)行一個應用模塊606時，結(jié)果被計算并存儲在輔助存儲器610(也就是，硬盤驅(qū)動器612)中。有限元分析的狀態(tài)通過I/O接口630報告給用戶。

雖然參照特定的實施例對本發(fā)明進行了描述，但是這些實施例僅僅是解釋性的，并不用于限制本發(fā)明。本技術(shù)領(lǐng)域的人員可得到暗示，對具體公開的示范性實施例做出各種修改和改變。雖然加速縮放因子在圖4B中已經(jīng)示出和描述為5，其他因子(大于1)也可以用于實現(xiàn)相同的目的。此外，雖然在圖2中節(jié)點的第一層和節(jié)點的第二層被示為平面層，其他形狀的層可以用于實現(xiàn)相同的目的?？傊?，本發(fā)明的范圍不限于在此公開的特定示范性實施例，對本技術(shù)領(lǐng)域人員來說暗含的所有修改都將被包括在本申請的精神和范圍以及后附權(quán)利要求的范圍內(nèi)。