ansys中如何計算積分的案例
有限元計算過程中積分點應(yīng)力如何外插至節(jié)點處？【公式推導(dǎo)篇】
注：由于技術(shù)鄰排版風(fēng)格有限，故部分內(nèi)容顯示不全，感興趣的小伙伴可點擊原文進行閱覽：
有限元計算過程中積分點應(yīng)力如何外插至節(jié)點處？【公式推導(dǎo)篇】
https://mp.weixin.qq.com/s/47byQ3b3e5UpbUp7Krs2mQ
本次分享的是：有限元計算過程中，單元積分點應(yīng)力如何外推至節(jié)點？
有關(guān)積分點與節(jié)點的概念可點擊跳轉(zhuǎn)閱讀歷史推文：有限元基本概念-【節(jié)點和積分點】，現(xiàn)科普一下Q4單元、Q8單元、Q9單元的形函數(shù)和高斯積分方案。
Q4單元
Q8/9單元
應(yīng)力外插
核心理念：坐標系的轉(zhuǎn)換。
假設(shè)是母單元的自然坐標系，是由高斯積分點控制的坐標系（術(shù)語可能不專業(yè)），假設(shè)高斯積分方案為。坐標系轉(zhuǎn)換關(guān)系：
單元內(nèi)任一點的應(yīng)力，由4個高斯積分點應(yīng)力進行插值時，可表示為
其中，是基于高斯積分點的形函數(shù)，第一個積分點的坐標在母單元坐標系下為（-1，-1），根據(jù)上述的坐標系轉(zhuǎn)換的方式，在高斯積分點的坐標系下，第一個單元節(jié)點在高斯積分點坐標系下坐標為，將此坐標值代入第一個形函數(shù)，得，相同的道理，可推導(dǎo)至四個節(jié)點在4個形函數(shù)下的外插矩陣：
對于Q8、Q9單元，依然可采用高斯積分方案（減縮積分）。
展開 abaqus C3D8 單元 計算中采用了多少個積分點？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結(jié)果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發(fā)現(xiàn)事實遠非如此，采用全積分計算得到的結(jié)果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應(yīng)變，采用8個積分，對于球應(yīng)變，采用中心點積分。這樣計算得到的結(jié)果才能與abaqus 完全對標，亦可從abaqus 幫助文檔得到答案。
展開 在 COMSOL 中如何求解積分上下限
那么，你知道 COMSOL 也可以計算積分嗎？
求解有限元問題需要對函數(shù)進行積分，COMSOL 不僅可以計算積分，還可以求解未知積分限的問題！
下面讓我來介紹方法。
對函數(shù)進行積分
考慮一個求解二次函數(shù)積分的問題：
積分可以獲得陰影區(qū)域的面積。
我們可以在 COMSOL Multiphysics 中使用積分函數(shù)來計算這個積分，這個函數(shù)的語法為：integrate(u^2, u, 0, 2, 1e-3)。其中，第一個參數(shù)是表達式，第二個參數(shù)是要積分的變量，第三和第四個參數(shù)是積分的極限，可選擇的第五個參數(shù)是必須在 0 到 1 之間的的積分相對容差。如果省略第五個參數(shù)，就將使用默認值 1e-3。我們可以在模型設(shè)置中的任意位置調(diào)用這個函數(shù)。
在這里，我們將在全局方程接口中使用它：
用于積分的全局方程計算了指定極限之間的積分。
到目前為止，這里沒有什么太大的驚喜。我們可以在 COMSOL Multiphysics 中求解這個問題，或者手算得到結(jié)果。但假設(shè)我們把問題稍微想復(fù)雜一點，如果我們知道積分的計算結(jié)果，但不知道積分的上限怎么辦？
我們來看看如何求解下面這個
的上限問題：
我們可以通過改變?nèi)址匠虂砬蠼膺@個問題，這樣就可以得到積分上限：
u_b 的全局方程求解了積分的上限，計算結(jié)果為 6。
上面的全局方程有一些變化。其中的變量更改為 u_b，必須等于零的表達式變?yōu)椋?-integrate(u^2, u, 0, u_b)。因此，軟件將找到一個 u_b 值，使積分等于指定的值。
請注意，u_b 的初始值不能為零。
展開 ansys之——如何將分析中前一次計算結(jié)果?
Q:挖分析中前一次計算結(jié)果導(dǎo)入下一部分析中
A:如果用dyna計算，有兩個可能：
1）如果網(wǎng)格需要重劃分，將ANSYS/lsdyna的計算結(jié)果插值到新網(wǎng)格中后輸出到數(shù)據(jù)文件，再組裝到lsdyna的.k文件中。
2）如果不需要網(wǎng)格重劃分，在用lsdyna計算之前，可用*set_part和*interface_springback_dyna3d將應(yīng)力應(yīng)變數(shù)據(jù)直接輸出到
dynain文件中，再編輯新的.k數(shù)據(jù)文件
A:我不明白為什么不能在你的新模型的第一載荷步進行重力加載計算，在第二載荷步進行挖掘計算。
即使按你所說的那樣，分成兩個模型，在lsdyna中也可以實現(xiàn)。即先進行重力載荷步計算，然后把計算結(jié)果輸出到另外一個計算模型中進行挖掘計算。這要求你在進行重力載荷步計算時，生成.K文件后，在此.K文件中加入（假設(shè)土體材料號為1），
*set_part
1
1
*interface_springback_dyna3d
1
計算結(jié)束后，會生成一個dynain文件，該文件中記錄了計算終點時的應(yīng)力分量和等效塑性應(yīng)變數(shù)據(jù)。至于土體的變形后的幾何模型很容易生成，有幾種辦法，最簡單的辦法是利用upgeom命令實現(xiàn)，如 UPGEOM,1,LAST,LAST,'test','rst',' ' 。
展開 
ANSYS mechanical如何在Workbench環(huán)境中使用高性能計算
ANSYS mechanical屬于隱式結(jié)構(gòu)有限元分析求解器，一般完成一個有限元分析過程需要前處理、求解和后處理三個步驟。前處理一般在圖形工作站上完成，有限元求解可在工作站、集群及SMP 服務(wù)器上進行。
對于中小型問題（例如1000 萬節(jié)點以內(nèi)的ANSYS mechanical問題），一般認為在圖形工作站上就可以進行求解；對于中大型問題（例如1000 萬節(jié)點以上的ANSYS mechanical問題），建議還是在計算性能更高的集群或SMP 服務(wù)器上進行。對于中小型問題，可以在圖形工作站上運行有限元后處理程序，讀取計算結(jié)果進行結(jié)果的分析。
因此對于ANSYS mechanical在Workbench環(huán)境中使用高性能計算的方法共有兩種：一種是直接通過workbench界面進行設(shè)置并行計算求解，在本地的工作站進行求解計算；另外一種是在workbench界面中將文件保存為ANSYS mechanical經(jīng)典界面的求解文件格式，提交給高性能計算平臺進行計算。
1、ANSYS mechanical在Workbench界面設(shè)置方法
此種方法適合中小型問題在本地的工作站進行求解計算，設(shè)置方法簡單方便。在Workbench界面環(huán)境下，打開Model模塊，在菜單中依次選擇Tools＞Solve Process Settings＞Advanced，進行CPU設(shè)置選擇對應(yīng)的CPU核數(shù)（建議關(guān)閉超線程，設(shè)置的核數(shù)不能超過工作站的CPU物理核數(shù)），默認使用分布式求解選項。
2、保存為經(jīng)典界面的求解文件格式方法
此種方法適合中大型問題在高性能計算平臺進行計算，需要在Workbench界面中存儲為指定的格式，設(shè)置步驟稍微繁瑣些。
展開 有限元中的高斯點與積分詳解(下)_《數(shù)值計算與程序設(shè)計》系列課程之九 ￥599
在下半節(jié)課中，詳細地分析了扭曲單元與有限元精度之間的關(guān)系。我們常聽到單元網(wǎng)格質(zhì)量要劃好，不然精度會不行，甚至?xí)蠼馐。@是為什么呢？我們透過表面來看有限元方法的本質(zhì)，用簡潔易懂的數(shù)學(xué)推導(dǎo)來展現(xiàn)誤差與單元形狀之間的關(guān)系。
本課從實際問題出發(fā)，帶著問題去講解有限元中的高斯點與數(shù)值積分。一開始拋出了以下3個關(guān)鍵問題：
1.對于一個任意函數(shù)怎么去得到它的積分？
2.數(shù)值積分的本質(zhì)是什么？為什么簡單地取幾個點就可得到積分值？此種方法的立足點在哪？
3.很多資料上都說“有限元求解精度嚴重依賴于網(wǎng)格質(zhì)量，過度扭曲的單元會導(dǎo)致結(jié)果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個問題，本課分別講了一下三個內(nèi)容：
1. 數(shù)值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
希望有興趣的同學(xué)多多支持下，你們的支持是我更新的動力
展開 有限元中的高斯點與積分詳解(上)_《數(shù)值計算與程序設(shè)計》系列課程之五 ￥599
本課從實際問題出發(fā)，帶著問題去講解有限元中的高斯點與數(shù)值積分。一開始拋出了以下3個關(guān)鍵問題：
1.對于一個任意函數(shù)怎么去得到它的積分？
2.數(shù)值積分的本質(zhì)是什么？為什么簡單地取幾個點就可得到積分值？此種方法的立足點在哪？
3.很多資料上都說“有限元求解精度嚴重依賴于網(wǎng)格質(zhì)量，過度扭曲的單元會導(dǎo)致結(jié)果不收斂或者精度極度惡化”，這只是為什么呢？扭曲單元到底影響的是有限元方法中的哪一步？
圍繞這3個問題，本課分別講了一下三個內(nèi)容：
1. 數(shù)值積分基本方法。
2. 有限元單元積分。
3. 誤差分析。
本次課程分為上下兩課，第一課講了第一和第二個內(nèi)容。關(guān)鍵詞是：數(shù)值積分的本質(zhì)，有限元高斯積分（附件中包含1個小時的詳細課程視頻以及PPT）。
在第二課中，再繼續(xù)展開第三部分內(nèi)容，誤差分析，解決問題“扭曲單元到底影響的是有限元方法中的哪一步”。
希望有興趣的同學(xué)多多支持下，你們的支持是我更新的動力
展開 【JY】ANSYS Workbench在減隔震應(yīng)用分析中的單元積分技術(shù)筆記
【JY】各類有限元軟件計算功能賞析與探討
我們需要更清楚減隔震元件的破壞模式，對減隔震元件進行破壞分析，除了對減隔震元件在正常工況下的性能進行評估，有限元技術(shù)還可以用于研究元件在極端條件下的破壞行為。這有助于了解元件的破壞機理，并為設(shè)計提供更全面的數(shù)據(jù)支撐。
并且在多物理場耦合分析也需要運用在實際應(yīng)用中，因為減隔震元件可能會面臨復(fù)雜的物理環(huán)境，如溫度變化、流體流動等。有限元技術(shù)可以考慮這些多物理場耦合效應(yīng)，從而更準確地預(yù)測元件在實際工況下的性能。
黏滯阻尼器的固流耦合分析：
對于ABAQUS的單元介紹已經(jīng)做了詳盡，個人感覺固體力學(xué)上ABAQUS還是上手比較方便，而多場耦合、快速建模預(yù)估Workbench會方便一些，因人而異：
【JY】有限單元分析的常見問題及單元選擇
ANSYS Workbench就像一個科技界的“瑞士軍刀”，集合了各種強大的單元技術(shù)，為減隔震元件提供全面且準確的分析支持。近期對于ANSYS Workbench進行了學(xué)習(xí)，本文將對ANSYS Workbench 各類單元技術(shù)做一個筆記總結(jié)，便于為減隔震元件分析提供理論基礎(chǔ)。（畢竟Workbench大部分時候會自動匹配相應(yīng)所需技術(shù)）
B-bar方法完全積分
Workbench中的B-bar方法是一種常用于處理低階單元完全積分的技術(shù)，也被稱為選擇性減積分策略。它是針對有限元分析（FEA）中的一種改進方法，旨在提高計算效率和準確性。
在傳統(tǒng)的有限元分析中，低階單元（如線性單元）在處理不可壓縮材料或近似不可壓縮材料時，常常遇到體積鎖定問題。體積鎖定是指在近似不可壓縮材料的有限元模擬中，由于體積應(yīng)變被過度限制，導(dǎo)致計算結(jié)果偏離實際情況的現(xiàn)象。為了解決這個問題，B-bar方法被引入到ANSYS Workbench中。
展開 Ansys Zemax | 公差的標準怎么計算的，如何確認計算細節(jié)?
使用 “SAVE” 公差操作數(shù)紀錄靈敏度計算過程
在敏感度 (靈敏度) 分析時，OpticStudio會把每一個公差操作數(shù)的最大最小值代入，并且計算在這些公差的極限狀態(tài)下，標準的變化如何，然后回報在文字報告中。而當(dāng)使用者對于這個結(jié)果有疑問時，會需要知道OpticStudio實際上是如何調(diào)整系統(tǒng)，并得到這個結(jié)果的，此時便是用 “SAVE” 這個操作數(shù)的時候。
舉例來說，假設(shè)我們有如下的公差設(shè)定以及靈敏度分析結(jié)果：
可以看到當(dāng)公差TRAD = -0.2時，標準 = 0.04967675
而TRAD = 0.2時，標準 = 0.04875308
假設(shè)我們想知道TRAD=0.2時的標準是如何計算的，我們可以在TRAD的下面加上一行SAVE指令，如下圖：
請注意在文件#（File#）字段代表檔案編號，如果有多個SAVE指令，則需要把編號分開。此外編號等于0的話，這個指令將不運作，不會存檔。
現(xiàn)在再執(zhí)行一次公差分析，文字報表的結(jié)果應(yīng)該相同，但是使用者可以發(fā)現(xiàn)在文件夾中多出兩個檔案，如下圖：
其中TSAV_MAX_0001代表TRAD=-0.2時的系統(tǒng)狀態(tài)，而TSAV_MIN_0001代表TRAD=0.2時的系統(tǒng)狀態(tài)。
讓我們打開TSAV_MAX_0001這個檔案，并開啟評價函數(shù)編輯器，可以看到如下圖：
可以看出系統(tǒng)計算TRAD=0.2的標準時，是計算RMS 光斑半徑，參考點為質(zhì)心，使用GQ算法，取樣是4個環(huán)以及8個臂。這反應(yīng)到我們之前的標準以及采樣設(shè)定。
此外評價函數(shù)值為0.0487530834843748，與公差報告吻合。
利用蒙特卡羅存檔了解公差擾動如何被執(zhí)行
前面我們介紹如何把靈敏度計算時用到的系統(tǒng)設(shè)定儲存下來。
展開 什么是光學(xué)計算？如何在 COMSOL 中分析光學(xué)計算器件
光學(xué)計算是替代當(dāng)前電子計算機的另一種可能形式。在這篇文章中，我們將探討光學(xué)計算的概念，并解釋了光學(xué)矩陣乘法網(wǎng)絡(luò)是如何工作的。我們還討論了如何使用 COMSOL Multiphysics? 軟件及其附加產(chǎn)品——波動光學(xué)模塊對光學(xué)計算設(shè)備進行建模。結(jié)合這些產(chǎn)品的使用，展示了在模擬大型光學(xué)系統(tǒng)時應(yīng)用波束包絡(luò)法的優(yōu)勢。
光學(xué)計算簡介
摩爾定律
在過去的幾十年里，計算機的能力一直呈指數(shù)級增長。這種增長遵循摩爾定律，即集成電路中的晶體管數(shù)量每兩年翻一番，而計算機的成本將降低。這使得我們今天享有的大部分現(xiàn)代技術(shù)成為可能。例如，主流計算機芯片完全基于晶體管等電子元件，每塊芯片的晶體管數(shù)量幾乎每兩年就會翻一番。為了跟上這種增長，并在可控的功率效率下提高計算機芯片的性能，芯片上的電子元件(包括晶體管)的小型化既關(guān)鍵又不可避免。盡管工程師們在這方面做了出色的工作，將晶體管從厘米尺度縮小到納米尺度，但重要的是要認識到，最終基本的限制將阻礙這類設(shè)備的發(fā)展。例如，當(dāng)一個電子元件的尺寸接近原子水平時，量子效應(yīng)將導(dǎo)致其功能不穩(wěn)定。科學(xué)和工程界長期以來一直在考慮電子計算機的替代形式。最近引起廣泛關(guān)注的一種替代是光學(xué)計算——指用光(光子)而不是電流(電子)進行計算。
雖然光學(xué)計算是一項新興技術(shù)，但光學(xué)在信息技術(shù)中的應(yīng)用已經(jīng)有相當(dāng)長的一段時間了，特別是利用光進行信息傳輸。損耗極低的光纖可以以光速長距離傳輸信息。光纖網(wǎng)絡(luò)設(shè)備常用于數(shù)據(jù)中心甚至普通家庭。然而，在商業(yè)化方面，利用光進行計算仍處于起步階段。
光學(xué)中的數(shù)學(xué)計算
眾所周知，某些光學(xué)過程對應(yīng)于數(shù)學(xué)計算。例如，考慮光的衍射。當(dāng)光通過衍射介質(zhì)時，本質(zhì)上是在進行傅里葉變換積分。然而，光學(xué)系統(tǒng)是否可以像我們今天擁有的計算機一樣進行通用數(shù)學(xué)計算，可能還不是很清楚。目前，光學(xué)計算有許多不同的形式。
展開 在ADINA中如何實現(xiàn)在計算過程中模擬單元移動
在ADINA中如何實現(xiàn)在計算過程中模擬單元移動，單元是平動的
感應(yīng)電爐生產(chǎn)球鐵中爐料如何選擇？如何調(diào)整計算成分？
（4）硅鐵加入量計算 計算式如下：
硅鐵加入量＝〔爐料總量×（球墨鑄鐵要求含硅量－孕育增硅量－球化增硅量－分析含硅量）〕/硅鐵含硅量
（5）錳和其他合金元素加入量 計算式如下：
合金鐵加入量＝爐料總量×（球墨鑄鐵要求合金含量－分析含量）/合金鐵中合金元素含量
上述計算中均不計燒損。這是因為感應(yīng)爐中熔化球墨鑄鐵，各元素?zé)龘p不大，同時鐵也有一定燒損。如遇到特殊情況，則應(yīng)注意燒損。
成分調(diào)整應(yīng)用舉例：
使用1t感應(yīng)電爐，生產(chǎn)珠光體球墨鑄鐵QT800—2，鑄件重1000kg。
控制化學(xué)成分為：wC＝3.5%、wSi ＝2.5%、wMn＝0.5%、wP＜0.07%、wS＜0.02%、wMg＝0.03%～0.06%、wRE＝0.02%～0.04%。
使用爐料成分如附表所示。
增碳后爐內(nèi)“廢鋼液”含碳量為2.5%。
生鐵加入量＝300×（3.5％－2.5％）/（4％－3.5％）＝600（kg）
回爐球墨鑄鐵使用量＝1000―300―600＝100（kg）
估算含硅量＝（300×0.3％＋600×1.8％＋100×2.5％）×100％＝1.42％
如果球化劑增硅量為0.6%、孕育劑增硅量為0.4%，則原鐵液含硅量應(yīng)為1.5%，高于估算含硅量。對此可以使用Z18鑄造生鐵調(diào)整碳當(dāng)量，使用回爐球墨鑄鐵補足爐料總量。
出爐前分析含硅量為1.38%，含錳量為0.15%。
展開 如何從Ansys APDL中提取剛度矩陣與質(zhì)量矩陣？ ￥69
1.引論
經(jīng)常使用Ansys、Abaqus等一系列有限元分析軟件進行計算、學(xué)習(xí)的學(xué)生或工程師們都會知道在有限元分析建模與計算中剛度矩陣與質(zhì)量矩陣的重要性。但是由于軟件的黑盒性質(zhì)，大家往往在實際使用十分成熟的商業(yè)化軟件的過程中慢慢忽視了有限元及其衍生出的商業(yè)軟件背后的原理與方法。
這時，不管是在學(xué)習(xí)中還是在工程應(yīng)用中往往都會遇到一個同樣的問題，那么就是如何將Ansys APDL運行中的產(chǎn)生的各種數(shù)據(jù)（例如：剛度矩陣、質(zhì)量矩陣）導(dǎo)出成為我們熟悉的形式或文件格式，從而為我們所用，所分析。
因此我決定寫下此篇文章來幫助很多實際工作或?qū)W習(xí)中需要用到此類技能的同學(xué)、同事們，讓大家更了解Ansys APDL背后的工作原理與數(shù)據(jù)導(dǎo)出方式。
當(dāng)然，在社區(qū)中早就有大佬回答過了這個問題，并給大家制作了相應(yīng)的提取矩陣軟件，其軟件具備了簡單、便捷的操作方式，讓很多想要提取剛度矩陣與質(zhì)量矩陣的同僚們受益，那么我為什么還要寫一篇這樣的文章重新提起這樣一個話題呢？這就又回到了我開頭所說的“原理與方法”，我在此更希望面對想要進一步學(xué)習(xí)了解軟件背后機理的群體，并在此基礎(chǔ)上保留教學(xué)的簡潔性，提供導(dǎo)出矩陣與轉(zhuǎn)換、列式、求解的源代碼，使其既兼顧基本原理，又可以讓大家直接上手使用，非常的便捷，也避免了很多因為優(yōu)化不完全導(dǎo)致的運行bug。
2.有限元軟件導(dǎo)出剛度矩陣與質(zhì)量矩陣的方法
在使用APDL進行求解時，每次在求解完成后都會在工作路徑下生成一個.full文件，而這個文件十分關(guān)鍵，其正是剛度矩陣與質(zhì)量矩陣的所在之處。
展開 ADAMS仿真過程中如何提高計算效率，縮短計算時間，相應(yīng)其他軟件也可以類似操作。（原創(chuàng)）
大家再用軟件做仿真計算時，總是感覺很費時間，有時候一算幾十個小時還沒有正確結(jié)果。個人總結(jié)了一下ADAMS中設(shè)置仿真計算節(jié)約時間的一些小規(guī)律，請采納，其他的軟件類似也是如此。
1、ADAMS中 end time和steps設(shè)置
endtime是仿真時長,step是仿真步數(shù)
顧名思義，仿真時長就是運動終止時間，如果是周期運動，一般計算幾個周期就行了，周期重復(fù)得到的結(jié)果就是一樣的，得到的曲線在一個周期時候系統(tǒng)已經(jīng)平衡，所以你的仿真時長不管改多大，曲線都會是同樣的。比如圓周運動和往復(fù)運動，計算兩三個周期的時間就夠了；
再說仿真步數(shù)，步數(shù)越多，仿真越詳細，計算量越大，但是精度也就越高，因為迭代的次數(shù)多，在你整個兒系統(tǒng)的驅(qū)動和約束已經(jīng)確定的情況下，對你仿真的結(jié)果不會產(chǎn)生太明顯的影響，所以這里適當(dāng)即可，幾百到小幾千已經(jīng)很好了，別大幾千上萬，那就是浪費了。
2、默認算法設(shè)置
系統(tǒng)中默認的算法采用的采用的GSTIFF算法，雖然不太懂什么意思，但是改成HHT算法計算效率能提高30%以上的，結(jié)果并沒有什么影響的，本人已經(jīng)通過算例驗算過。具體操作改正如下:
ADAMS view--settings--solver--dynamics--integrator--HHT
3、計算機多核設(shè)置
一般默認計算機只設(shè)置了單核計算，效率很低，大家都不會去修改，如果計算機是雙核，四核八核呢，是不是快很多。操作如下：
-ADAMS view--settings--solver--executable--左下角more--把1直接改成2、4、8
現(xiàn)在就這么多，后期發(fā)現(xiàn)還有再給補上吧。
展開 如何用ANSYS_WB做一桿斯諾克，采用顯示動力學(xué)模塊計算臺球碰撞問題,私信郵箱獲取計算文件。
采用ANSYS_WB的顯示動力學(xué)模塊模擬臺球碰撞問題，對于臺球碰撞屬于短時間接觸，計算所需要的時間步長足夠小才能捕捉到短時間的接觸過程，并且我們希望每個時間步計算應(yīng)該足夠快，不然硬件吃不消的。
理論上ANSYS_WB 中
瞬態(tài)結(jié)構(gòu)模塊
和
顯示動力學(xué)模塊
都可以模擬這樣一個臺球碰撞過程，但是
瞬態(tài)結(jié)構(gòu)模塊是采用隱式積分算法
，隱式積分可以使得時間步長很大，但每個時間步需要多次迭代才能達到收斂，時間步過多，計算時間將非常大，
顯示動力學(xué)模塊采用顯示積分
，時間步可以非常小足以捕捉瞬間碰撞行為，且不需要在每個時間步上進行剛度矩陣總裝，每個時間步計算非常快。因此這里采用顯示動力學(xué)模塊進行模擬。
有感興趣的朋友們
私信郵箱獲取計算文件
哦，創(chuàng)作不易，歡迎大家點贊轉(zhuǎn)發(fā)支持筆者。
計算結(jié)果
教程：Step by Step
建模：
采用ANSYS自帶的建模軟件進行建模，不做介紹。
計算模塊建立：
拖動Explicit Dynamics模塊到WB工作區(qū)域(左邊是我已經(jīng)計算完的模塊，拖到一個獨立的區(qū)域了)。
材料定義：
雙擊Engineering Data,建立新材料，選擇各向同性材料，輸入密度，模量，泊松比。
模型導(dǎo)入：采用ANSYS自帶的建模軟件進行建模，并導(dǎo)入顯示動力學(xué)計算模塊中。
剛性體定義：將臺球和臺球桌面定義為剛性體
網(wǎng)格劃分：
相互作用定義：小球間接觸采用摩擦接觸。
展開 