導(dǎo)讀：動(dòng)力學(xué)仿真計(jì)算技術(shù)用來(lái)確定慣性、阻尼起重要作用時(shí)結(jié)構(gòu)動(dòng)力特性計(jì)算技術(shù)。結(jié)構(gòu)動(dòng)載計(jì)算分析滿(mǎn)足但動(dòng)力學(xué)分析要求不一定滿(mǎn)足，例如尾氣排氣管固有頻率與發(fā)動(dòng)機(jī)激勵(lì)頻率接近時(shí)可能會(huì)振散開(kāi)裂，因此動(dòng)力學(xué)分析計(jì)算是必要的。

動(dòng)力學(xué)計(jì)算一般基于模態(tài)疊加法后再進(jìn)行諧響和PSD計(jì)算研究，對(duì)于關(guān)注頻率范圍較寬高達(dá)上千Hz以上航天系統(tǒng)或?qū)τ陉P(guān)注頻率周邊簇點(diǎn)較多掃頻計(jì)算時(shí)，全零件級(jí)網(wǎng)格計(jì)算時(shí)長(zhǎng)和存儲(chǔ)量數(shù)據(jù)驚人。

本文中筆者鑒于動(dòng)力學(xué)計(jì)算工作中遇到相關(guān)問(wèn)題，結(jié)合ANSYS子結(jié)構(gòu)CMS方法幫助文檔，給出子結(jié)構(gòu)CMS法在動(dòng)力學(xué)計(jì)算中應(yīng)用簡(jiǎn)介，以能為工程企業(yè)同仁提供一種計(jì)算處理方法，若能同時(shí)為社會(huì)工程應(yīng)用做出技術(shù)貢獻(xiàn)則深感高興。

二、子結(jié)構(gòu)法簡(jiǎn)介

子結(jié)構(gòu)技術(shù)能過(guò)相對(duì)少的超單元主DOF去描述一組單元的等效質(zhì)量、剛度、阻尼矩陣等。使用子結(jié)構(gòu)技術(shù)能夠有效節(jié)省計(jì)算求解時(shí)長(zhǎng)和降低存儲(chǔ)文件規(guī)模，通常用于一個(gè)系統(tǒng)結(jié)構(gòu)包含許多重復(fù)組件時(shí)，整個(gè)結(jié)構(gòu)模型線(xiàn)性部分，子結(jié)構(gòu)防止重新計(jì)算單元矩陣，減少平衡迭代計(jì)算時(shí)間，不同設(shè)計(jì)組獨(dú)立工作“組裝”系統(tǒng)級(jí)模型，適合于大型動(dòng)態(tài)分析結(jié)構(gòu)組裝形成復(fù)雜系統(tǒng)水平模型計(jì)算，例如飛機(jī)、核電站、石油平臺(tái)等大型結(jié)構(gòu)/系統(tǒng)。

ANSYS Mechanical子結(jié)構(gòu)技術(shù)能夠利用CMS法進(jìn)行支持動(dòng)力學(xué)計(jì)算內(nèi)容，包括模態(tài)、模態(tài)疊加法諧響應(yīng)分析、隨機(jī)振動(dòng)PSD、響應(yīng)譜以及剛性動(dòng)力學(xué)分析計(jì)算模塊等。子結(jié)構(gòu)技術(shù)借助Condensed Part功能進(jìn)行實(shí)現(xiàn)，Condensed Part功能將幾何視為向量組成的超單元，其自由度遠(yuǎn)低于完整有限元網(wǎng)格模型。Condensed Part功能子結(jié)構(gòu)行為將幾何組合綜合考慮更復(fù)雜結(jié)構(gòu)響應(yīng)，也能在剛性動(dòng)力學(xué)求解器中建模柔性零件等。

CMS（Component Mode Synthesis）指系統(tǒng)矩陣被簡(jiǎn)化為主DOF集合和其他組件之間一系列接口，是一種子結(jié)構(gòu)技術(shù)表達(dá)方法。Mechanical提供了Condensed Part功能生成超單元的方法，使用Expansion Settings功能擴(kuò)展求解方案。

三、CMS方程推導(dǎo)

（1）動(dòng)力學(xué)基本方程

（2）超單元位移向量和矩陣組合方式

M：接觸節(jié)點(diǎn)的主DOF；

S：非主節(jié)點(diǎn)的全部DOF；

{u}：位移向量用主坐標(biāo)和廣義坐標(biāo)表示；

{yδ}：截?cái)鄰V義坐標(biāo)集；

[T]：變換矩陣。

三、CMS設(shè)計(jì)流程

1、生成子結(jié)構(gòu)

一組單元及其相關(guān)的接觸面節(jié)點(diǎn)被縮減到一個(gè)單一的超單元中，該超單元由被由縮減質(zhì)量、剛度、阻尼矩陣組成。Mechanical中實(shí)現(xiàn)生成子結(jié)構(gòu)的工具是Condensed Parts。

2、使用子結(jié)構(gòu)參與計(jì)算

使用超單元模型來(lái)代表結(jié)構(gòu)的一部分分析。生成Condensed Part后，Mechanical將在模態(tài)求解或諧響應(yīng)求解中自動(dòng)處理子結(jié)構(gòu)的計(jì)算“使用通道”。

3、擴(kuò)展子結(jié)構(gòu)求解結(jié)果

利用超單元的主DOF位移以及廣義坐標(biāo)和變換矩陣來(lái)計(jì)算超單元內(nèi)的位移和應(yīng)力。Condensed Part會(huì)在求解分支中創(chuàng)建擴(kuò)展設(shè)置對(duì)象Expansion Settings。子結(jié)構(gòu)位移的計(jì)算來(lái)自于主位移計(jì)算的延伸擴(kuò)展，Expansion Settings設(shè)置主要進(jìn)行子結(jié)構(gòu)位移的計(jì)算，默認(rèn)情況子結(jié)構(gòu)的計(jì)算結(jié)果不作為求解內(nèi)容的部分進(jìn)行展開(kāi)。

四、Condensed Parts基礎(chǔ)

Mechanical采用CMS法支持模態(tài)和模態(tài)疊加法諧響應(yīng)分析以及剛性動(dòng)力學(xué)， CMS自動(dòng)化程度較高，工作流程以捕獲在一個(gè)壓縮幾何對(duì)象中進(jìn)行裝載，即生成子結(jié)構(gòu)的過(guò)程。

圖2 Condensed Parts工具

Condensed Parts提供所有信息通道生成所需的數(shù)據(jù)的工具，包括矩陣縮減方法，接口交互方法、求解設(shè)置等。接口功能定義主自由度與模型其他部分的接口（這些接口通常以contact/jionts形式存在）以及包括Condensed Parts中的Load、Support、Mass、Remote points、Name selection等。

1、Condensed Parts菜單細(xì)節(jié)

（1）矩陣縮減方法：CMS。

（2）接口方法：Fixed。特征向量在全部fixed約束進(jìn)行計(jì)算。

（3）模態(tài)數(shù)量使用：用于描述超單元的特征向量的數(shù)量。

（4）接口：定義連接到其余模型的DOF，使用自動(dòng)檢測(cè)建立自動(dòng)接口。

圖3 Condensed Parts菜單

2、Condensed Parts生成流程

(1)將進(jìn)行子結(jié)構(gòu)創(chuàng)建零件建立設(shè)置關(guān)系，包括綁定接觸、載荷、約束、點(diǎn)質(zhì)量定義等。

(2)添加Condensed Parts項(xiàng)，將該零件指定為Condensed Parts。當(dāng)指定Condensed Parts后，該幾何所賦予載荷、約束等相關(guān)設(shè)置內(nèi)容臨時(shí)掛起。不可編輯，呈現(xiàn)？狀態(tài)。

右鍵Condensed Parts項(xiàng)快捷Detect Condensed Parts Interface進(jìn)行接口關(guān)系創(chuàng)建。子結(jié)構(gòu)“接口”在工作表視圖中可以瀏覽，呈現(xiàn)？的未定義項(xiàng)將再次完全定義。注意可以使用工作表進(jìn)行手動(dòng)接口內(nèi)容定義，需要Name selection進(jìn)行頂點(diǎn)、邊、面等命名選擇。

圖5 Condensed Parts Interface工作表

3、Condensed Parts應(yīng)用注意事項(xiàng)

（1）Condensed Parts計(jì)算過(guò)程中進(jìn)行非線(xiàn)性特性忽略。

（2）Condensed Parts非連接件共享拓?fù)洹eam連接、綁定接觸、分布質(zhì)量、固定關(guān)節(jié)、點(diǎn)質(zhì)量、彈簧連接等數(shù)據(jù)傳遞。

（3）Condensed Parts接口連接結(jié)構(gòu)模型有效性取決于幾何、遠(yuǎn)程點(diǎn)、連接、加載、約束、點(diǎn)質(zhì)量、命名選擇等的合理定義，接口生成過(guò)程偵測(cè)創(chuàng)建。

（4）Automatic contact detection接觸對(duì)創(chuàng)建時(shí)常導(dǎo)致冗余接觸表面，自動(dòng)接觸創(chuàng)建過(guò)程應(yīng)控制接觸容差并對(duì)接觸面、目標(biāo)面內(nèi)容確認(rèn)合理性。Condensed Parts封裝無(wú)效接觸面節(jié)點(diǎn)會(huì)增加求解時(shí)長(zhǎng)，一種好的方法是采用Fixed jiont進(jìn)行接觸對(duì)的取代。



圖6 自動(dòng)接觸創(chuàng)建封裝主節(jié)點(diǎn)數(shù)量

五、Condensed Part阻尼定義

對(duì)于動(dòng)力學(xué)分析來(lái)說(shuō)，阻尼特性是影響結(jié)構(gòu)幅值響應(yīng)表現(xiàn)的重要因素，Condensed Part阻尼在工程數(shù)據(jù)中進(jìn)行定義，不支持阻尼比定義，整機(jī)計(jì)算模型阻尼比可以設(shè)置。Condensed Part阻尼另外一種定義方式采用Bushing Joints定義進(jìn)行6自由度剛度、阻尼施加。

Condensed Part阻尼的求解需要使用阻尼模態(tài)求解器定義：

（1）Solver Type = Reduced Damped

（2）Store Complex Solution = No

圖7 阻尼模態(tài)求解器

圖8 Bushing Joints阻尼與剛度定義

六、子結(jié)構(gòu)求解擴(kuò)展

生成Condensed Part后，Mechanical在模態(tài)求解、諧響應(yīng)求解中等自動(dòng)處理子結(jié)構(gòu)計(jì)算“使用通道”。

求解后Condensed Part在求解分支創(chuàng)建擴(kuò)展設(shè)置對(duì)象Expansion Settings。子結(jié)構(gòu)位移計(jì)算來(lái)自于主位移計(jì)算的延伸擴(kuò)展，Expansion Settings設(shè)置主要進(jìn)行子結(jié)構(gòu)位移的計(jì)算，默認(rèn)情況非子結(jié)構(gòu)計(jì)算結(jié)果作為求解內(nèi)容部分進(jìn)行展開(kāi)。

圖9Expansion Settings

子結(jié)構(gòu)幾何完成計(jì)算求解后具備完整后處理結(jié)果，子通過(guò)對(duì)壓縮對(duì)象進(jìn)行通道擴(kuò)展，否則只有子結(jié)構(gòu)的主接口位移結(jié)果可用。Expansion Settings worksheet控制Condensed Part后處理結(jié)果，包括特定Condensed Part或多個(gè)Condensed Parts幾何；擴(kuò)展設(shè)置與輸出控制內(nèi)容類(lèi)型請(qǐng)求組合輸出子結(jié)構(gòu)位移和所有評(píng)價(jià)的求解計(jì)算結(jié)果。

圖10 控制輸出內(nèi)容

圖11 Condensed Parts 擴(kuò)展項(xiàng)選擇

圖12 運(yùn)行擴(kuò)展

七、CMS模態(tài)法和完全模態(tài)法頻率比對(duì)

某大型幾何結(jié)構(gòu)測(cè)試計(jì)算中，采用CMS法子結(jié)構(gòu)模態(tài)分析與一般全零件級(jí)模態(tài)分析計(jì)算比對(duì)結(jié)果如圖所示，計(jì)算穩(wěn)定性較高，且求解時(shí)間縮短較為明顯（依據(jù)計(jì)算機(jī)性能表現(xiàn)各異）。

圖13CMS子結(jié)構(gòu)與全零件模態(tài)計(jì)算比對(duì)

八、 案例應(yīng)用

某產(chǎn)品結(jié)構(gòu)處于路譜行駛顛簸環(huán)境，供應(yīng)商復(fù)現(xiàn)問(wèn)題不能還原產(chǎn)品失效原因，且具有明顯失效位置偏頗。

在缺少路譜數(shù)據(jù)轉(zhuǎn)化PSD條件下，筆者以顛簸環(huán)境簡(jiǎn)化為諧響應(yīng)分析為出發(fā)點(diǎn)確定最大激勵(lì)頻率對(duì)應(yīng)的結(jié)構(gòu)響應(yīng)（當(dāng)載荷激勵(lì)頻率與結(jié)構(gòu)固有頻率接近共振，激勵(lì)頻率共振下結(jié)構(gòu)應(yīng)力狀態(tài)可以是靜載數(shù)倍至幾十倍）；基于此仿真計(jì)算研究方法假設(shè)，從動(dòng)力學(xué)諧響應(yīng)分析入手，以焊縫振動(dòng)疲勞計(jì)算研究結(jié)束。

圖14 某結(jié)構(gòu)幾何模型

首先在仿真計(jì)算處理中將大型網(wǎng)格結(jié)構(gòu)進(jìn)行子結(jié)構(gòu)CMS法進(jìn)行處理，建立多個(gè)Condensed Part，降低計(jì)算時(shí)間和存儲(chǔ)空間。

其次完成結(jié)構(gòu)計(jì)算模型固有頻率、模態(tài)振型計(jì)算，基于模態(tài)疊加法諧響應(yīng)分析計(jì)算確定最大激勵(lì)頻率和應(yīng)力響應(yīng)。

圖15 模態(tài)和諧響應(yīng)分析計(jì)算

最后基于nCode DesignLife掃頻振動(dòng)疲勞完成焊縫振動(dòng)疲勞求解，預(yù)測(cè)焊縫位置狀態(tài)與開(kāi)裂狀態(tài)一致，振動(dòng)壽命和實(shí)際產(chǎn)品服役時(shí)長(zhǎng)表現(xiàn)接近，完成該故障計(jì)算研究工作。

圖16 焊縫振動(dòng)疲勞求解

審核編輯：劉清