現在Abaqus、LS-DYNA、Ansys等結構商軟都說可以處理復雜的上萬零部件接觸的整車、整機等模型仿真，沒做過實際的這種仿真分析，很好奇，接觸分析算法往往涉及大變形、邊界不連續，只要輸入條件或者算法稍微變化一些，兩個零部件算出來的接觸結果就可能差異很大，更不用說上萬個零部件的接觸結果了，對這種大規模組裝模型的仿真結果不知如何來判斷它的可靠性，像普通的只校核一下材料的應力還是看一下動畫是否和試驗一致

現在Abaqus、LS-DYNA、Ansys等結構商軟都說可以處理復雜的上萬零部件接觸的整車、整機等模型仿真，沒做過實際的這種仿真分析，很好奇，接觸分析算法往往涉及大變形、邊界不連續，只要輸入條件或者算法稍微變化一些，兩個零部件算出來的接觸結果就可能差異很大，更不用說上萬個零部件的接觸結果了，對這種大規模組裝模型的仿真結果不知如何來判斷它的可靠性，像普通的只校核一下材料的應力還是看一下動畫是否和試驗一致

而Ansys、Abaqus暫時沒有虛擬質量方法，都是直接采用基于聲學有限元的流固耦合來求濕模態。 本章只介紹基于虛擬質量的濕模態計算。

圖6 受剪砌體試件的精細化模型 1砂漿與磚塊本構模型 在砌體精細化模型中，砂漿與磚塊的本構模型采用ABAQUS自帶的混凝土塑性損傷本構模型(CDP)進行模擬。

上一章講過聲學有限元只要加入聲學單元，求出類似的剛度陣和非平衡力，就很容易嵌入到基于增量迭代法的有限元結構流程中，但邊界元實際上融入這個流程還有相當的困難，按照最終的方程來說，我們可以把P(r)前的系數陣當成剛度陣，然后也可以采用迭代法來求非平衡力，正常來說也是一次平衡，但邊界元基本不這么做，我們理解困難點在于全局剛度陣的組裝，有限元中由于節點只與跟它相連的單元節點影響，可以先求出單元剛度陣得到該單元內節點之間關系

4 系統算例 以10T22.5M雙梁橋式起重機的箱形偏軌主梁為例，運用系統進行算例分析，系統主界面如圖5所示，先后點擊“Creo4.0鏈接”、“ANSYS19.0鏈接”,前處理界面如圖6所示，結構設計界面如圖7所示，ANSYS求解界面如圖8所示，輸入對應參數，即可完成設計，設計完成后查看求解結果和生成計算書，查看結果圖如圖9所示。

圖中第一條藍色虛線之前屬于螺栓的預緊過程，在后半段（累計收斂步 數 450~850），力收斂值一直小于力標準值，這一段計 算收斂速度特別快；在前半段（累計收斂步數 450 之 前）收斂速度較慢，因此對加載曲線適度調整使該段 對應的螺母旋轉角速度減緩，讓接觸面的接觸狀態的 改變不太劇烈，從而使計算更容易收斂。

回答當然是肯定的，主要有幾點原因: ①復雜裝配體動輒成百上千的螺栓連接，大量的螺栓連接直接導致費時的接觸對創建工作 有伙伴會說:現在很多軟件可以使用批處理，自動識別接觸對或者通用接觸大幅度縮減這部分工作量 但是， ②螺栓連接涉及接觸非線性問題，非線性的引入使得求解需要迭代，對于大型裝配體，其調試成本，計算時間不容小覷 有伙伴會說:公司電腦擱那放著，啥時候算完啥。

以下是我建立的一個三維模型，但是由于個人筆記本電腦算力不足，作為學習，我采用簡化的二維模型進行了計算，計算結果如下圖所示。 （1）首先是建立模型，拖拽geometry模塊進入操作界面即可建模，模型建立可以通過軟件自帶的Design model模塊，或者其他建模軟件，如solidworks等。主要原則是建立一個為大流場所包圍的固體模型，這里不詳細介紹。

我們不是什么編程天才，只能靠時間堆積慢慢烏龜式的前進，沒有每個人朝著同一個目標步調一致的前進，很多年的默默付出，哪來的真正能看得著摸得著的軟件實體，可以說團隊即軟件本身。