我們關(guān)注CAE中的結(jié)構(gòu)有限元，所以主要選擇了商用結(jié)構(gòu)有限元軟件中文檔相對(duì)較完備的Abaqus來(lái)研究?jī)?nèi)部實(shí)現(xiàn)方式，同時(shí)對(duì)某些問(wèn)題也會(huì)涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問(wèn)題在數(shù)學(xué)上其實(shí)并不嚴(yán)謹(jǐn)，同時(shí)由于水平有限可能有許多的理論錯(cuò)誤，歡迎交流討論，也期待有更多的合作機(jī)會(huì)。

1.2 關(guān)鍵計(jì)算公式 漢航NTS.LAB傳遞路徑分析模塊可解決結(jié)構(gòu)聲TPA和空氣聲TPA兩類傳遞路徑分析問(wèn)題，軟件模塊已封裝核心計(jì)算邏輯，計(jì)算方法包經(jīng)典的直接法、懸置剛度法、單路徑分析法、逆矩陣法、工況TPA和高階分析方法—組件TPA。下面以逆矩陣法介紹結(jié)構(gòu)聲TPA和空氣聲TPA兩類計(jì)算模型以及最新研發(fā)的組件TPA分析方法。

10 月 24 日（周五）下午，Ansys 總部院士朱永誼博士首次線下開講，帶來(lái)四大“黑科技”： 1 混合多點(diǎn)約束 “一個(gè)接觸對(duì)”自動(dòng)識(shí)別固體-殼任意組合，依局部幾何秒選最優(yōu)約束，無(wú)需手動(dòng)修正偏移或?qū)R法向，前處理更省力，結(jié)果更精準(zhǔn)。

本文的重點(diǎn)不是通過(guò)開發(fā)修正這些劣勢(shì)，而是在評(píng)估VOF和LS在數(shù)值計(jì)算過(guò)程中的表現(xiàn)時(shí)需要合理地重視并盡量規(guī)避這些弊端。國(guó)外商軟與VirtualFlow搭載的都是基于有限體積法求解多流體Navier-Stokes方程的壓力基求解器。當(dāng)然，使用的離散化方案和時(shí)間積分方法是不同的。

在Ansys Motion中模擬洗衣機(jī)筒的復(fù)雜運(yùn)動(dòng)狀態(tài)（力、位移、加速度等），在Ansys Rocky中，使用SPH法模擬了平衡環(huán)內(nèi)液體的復(fù)雜流動(dòng)形式。通過(guò)兩個(gè)模塊的耦合計(jì)算，使得先前難以從測(cè)試和傳統(tǒng)耦合仿真（Mechanical和Fluent）進(jìn)行研究的平衡環(huán)問(wèn)題，得到了新的探索路徑。

，單元通過(guò)共旋坐標(biāo)法分離剛體運(yùn)動(dòng)與彈性變形，結(jié)合 von Karman 非線性板理論，可精確模擬載荷 - 位移曲線中的 “階躍” 現(xiàn)象。即使在粗網(wǎng)格（4×4×2）下，單元計(jì)算結(jié)果與解析解的誤差仍小于 5%，顯著優(yōu)于傳統(tǒng) C3D8R/Solid45 單元。 將擬協(xié)調(diào)單元CSS8與 ANSYS 的 Solsh190、ABAQUS 的 SC8R進(jìn)行對(duì)比，從精度、效率、穩(wěn)定性三方面評(píng)估優(yōu)勢(shì)。

</p><p>最終，通過(guò)求解由能量原理和加權(quán)殘差法導(dǎo)出的代數(shù)方程組，獲得了有限元法的數(shù)值解。這個(gè)解是對(duì)原始連續(xù)體問(wèn)題的近似，其精度取決于網(wǎng)格剖分的細(xì)密程度和所采用的插值函數(shù)的類型。</p><p><br></p><p>2 Ansys workbench有限元分析軟件</p><p>在ANSYS 7.0版本問(wèn)世之前，ANSYS公司致力于研發(fā)其核心產(chǎn)品ANSYS。

將疊加后的PSD譜直接復(fù)制到Ansys Workbench中，再進(jìn)行輸入Improved fit后即可進(jìn)行正常隨機(jī)振動(dòng)仿真。 示例： 1.模態(tài)疊加法隨機(jī)振動(dòng)分析，計(jì)算結(jié)構(gòu)模態(tài)。

3.3.2.2 高斯消元法求解問(wèn)題 第二個(gè)問(wèn)題是結(jié)構(gòu)有限元一般采用直接法來(lái)求解矩陣，高斯消元法最常見的直接法，高斯消元法將矩陣分解成單位下三角矩陣L和上三角矩陣U，高斯消元法不是直接求K的逆矩陣來(lái)求解。

通常，選擇性減縮積分，Galerkin混合公式和穩(wěn)定方法可以緩解體積鎖死，其中選擇性減縮積分的一種也稱為B-bar方法，在ABAQUS的C3D8和ANSYS的SOLID185中都采取B-bar修正作為八節(jié)點(diǎn)線性插值實(shí)體單元的默認(rèn)公式。