有限元離散元耦合的案例
快速了解離散元仿真軟件Altair EDEM（與多體/有限元/流體軟件實現耦合）
EDEM可以聯合主流的CAE工具軟件進行顆粒系統與流體、機械結構及電磁場的耦合模擬仿真。
1、與 MBD（多體動力學） 軟件耦合
EDEM與 MBD（多體動力學） 進行耦合，可以仿真設備的動態力學響應，不僅可獲取固體散料對機械設備的真實載荷大小及其對設備性能產生的影響，同時可通過分析固體散料的力學響應，為機械設備作業質量評估提供依據。
Altair HyperWorks? 2019.1實現了Altair MotionSolve和EDEM之間的實時雙向耦合。
在EDEM中創建散料的模型，設定顆粒的形狀和質量等屬性，創建顆粒間的接觸。在MotionSolve創建系統的多體動力學模型，與EDEM共享相關的幾何。耦合仿真同時計算，每個時間步交換數據：MotionSolve計算設備部件的位置和速度，共享數據給EDEM，EDEM計算散料顆粒之間的接觸力，以及與設備部件之間的相互作用，共享各部件上的合力和力矩給MotionSolve。兩者耦合計算整個系統的運動狀態。
借助EDEM與MotionSolve的雙向耦合，可以分析挖掘機的鏟斗在不同操作工況下的載荷，評估挖掘深度、鏟斗裝載率、結構件載荷分布、動力系統匹配等。
2、與FEA（有限元分析）軟件耦合
EDEM可以與有限元分析軟件耦合，從而實現對施加在機器零件的載荷進行仿真分析，并將結果直接導出到所選的結構分析工具中。
鏟斗應力分析
3、與 CFD（計算流體動力學） 進行耦合
EDEM可以與 CFD（計算流體動力學） 進行耦合，用于顆粒級的固-液相系統的建模。
展開 有限元與離散元。
目前世界上結構計算方法一般分為有限元（FEM finite element method）、離散元（DEM discrete element method）、還有邊界元(EEM)。
離散元方法是由分析離散單元的塊間接觸入手找出其接觸的本構關系建立接觸的物理力學模型并根據牛頓第二定律對非連續、離散的單元進行模擬仿真。而有限元方法是將介質復雜幾何區域離散為具有簡單幾何形狀的單元通過單元集成、外載和約束條件的處理得到方程組再求解該方程組就可以得到該介質行為的近似表達。
離散元方法的基本概念
離散元方法也被稱為散體單元法，最早是1971年由Cundall 提出的一種不連續數值方法模型離散元理論是由分析離散單元的塊間接觸入手找出其接觸的本構關系建立接觸的物理力學模型并根據牛頓第二定律建立力、加速度、速度 及其位移之間的關系對非連續、離散的單元進行模擬仿真。
離散元法是專門用來解決不連續介質問題的數值模擬方法。該方法把節理巖體視為由離散的巖塊和巖塊間的節理面所組成,允許巖塊平移、轉動和變形,而節理面可被壓縮、分離或滑動。因此,巖體被看作一種不連續的離散介質。其內部可存在大位移、旋轉和滑動乃至塊體的分離,從而可以較真實地模擬節理巖體中的非線性大變形特征。離散元法的一般求解過程為:將求解空間離散為離散元單元陣,并根據實際問題用合理的連接元件將相鄰兩單元連接起來;單元間相對位移是基本變量,由力與相對位移的關系可得到兩單元間法向和切向的作用力;對單元在各個方向上與其它單元間的作用力以及其它物理場對單元作用所引起的外力求合力和合力矩,根據牛頓運動第二定律可以求得單元的加速度;對其進行時間積分,進而得到單元的速度和位移。從而得到所有單元在任意時刻的速度、加速度、角速度、線位移和轉角等物理量。
展開 離散元-邊界元動力耦合模型
本文提出了一種二維變形體離散元與時域邊界元的耦合模型，這一模型可以將非連續體的模擬與無限域的模擬統一在一個模型中，可用于在地震波動輸入條件下，考慮輻射阻尼的巖體邊坡或地下結構等的動力穩定和變形分析，拓寬了離散元動力分析的領域。算例分析表明本耦合分析模型具有較高的精度
261281--.doc
離散元與有限元優缺點（附inp文件）
來源：ABAQUS大世界 （ABAQUSworld）
文末有離散元法模擬攪拌過程案例文件
目前世界上結構計算方法一般分為有限元（FEM finite element method）、離散元（DEM discrete element method）、還有邊界元(EEM)。
離散元方法是由分析離散單元的塊間接觸入手找出其接觸的本構關系建立接觸的物理力學模型并根據牛頓第二定律對非連續、離散的單元進行模擬仿真。而有限元方法是將介質復雜幾何區域離散為具有簡單幾何形狀的單元通過單元集成、外載和約束條件的處理得到方程組再求解該方程組就可以得到該介質行為的近似表達。
離散元方法的基本概念
離散元方法也被稱為散體單元法，最早是1971年由Cundall 提出的一種不連續數值方法模型離散元理論是由分析離散單元的塊間接觸入手找出其接觸的本構關系建立接觸的物理力學模型并根據牛頓第二定律建立力、加速度、速度 及其位移之間的關系對非連續、離散的單元進行模擬仿真。
離散元法是專門用來解決不連續介質問題的數值模擬方法。該方法把節理巖體視為由離散的巖塊和巖塊間的節理面所組成,允許巖塊平移、轉動和變形,而節理面可被壓縮、分離或滑動。
因此，巖體被看作一種不連續的離散介質。其內部可存在大位移、旋轉和滑動乃至塊體的分離,從而可以較真實地模擬節理巖體中的非線性大變形特征。
展開 
有限離散元FDEM
目前在學習FDEM，想用有限離散元方法做些計算，并在此基礎上做研究，希望學習過FDEM的大佬談談學習經驗，互相交流，共同學習~
歡迎有FDEM學習經驗的大佬、想要學習或了解FDEM的同學加QQ群：821022127
淺談有限元方法的核心思想：數值近似和離散化
有限元法是基于近代計算機的快速發展而發展起來的一種近似數值方法，用來解決力學，數學中的帶有特定邊界條件的偏微分方程問題。而這些偏微分方程是工程實踐中常見的固體力學和流體力學問題的基礎。
有限元和計算機發展共同構成了現代計算力學的基礎。有限元法的核心思想是“數值近似”和“離散化”，所以它在歷史上的發展也是圍繞著這兩個點進行的。
1. 數值近似
由于在有限元法被發明之前，所有的力學問題和工程問題中出現的偏微分方程只能依靠單純的解析解得到解答。這種方法對數學要求很高，而且非常依賴于一些理想化的假定。比如在土木工程中梁柱計算中出現的平截面假定，小應變假定，理想塑性假定。這些假定其實是和實際工程問題有很大偏差的，而且一旦工程問題稍微復雜一些我們就不能直接得到解析解，或者解析解的答案誤差過大。
而有限元法把復雜結構離散到有限個單元，再把這種理想化的假定和力學控制方程施加于結構內部的每一個單元，然后通過單元分析組裝得到結構總剛度方程，通過邊界條件和其他約束解得每個單元的反應，這樣就可以避免直接建立復雜結構的力學和數學模型了。其總過程可以描述為：
總結構離散化 — 單元力學分析 — 單元組裝 — 總結構分析 — 施加邊界條件 — 得到結構總反應 — 單元內部反應分析
在進行單元分析和單元內部反應分析的時候，形函數插值和高斯數值積分被用來近似表達單元內部任意一點的反應，這就是有限元數值近似的重要體現。一般來說，形函數階數越高，近似精度也就越高，但其要求的單元控制點數量和高斯積分點數量也更多。另外單元劃分的越精細，其近似結果也更加精確。但是以上兩種提高有限元精度的代價就是計算量幾何倍數增加。
為了提高數值近似精度同時盡量較少地提高計算量，有限元法經歷了很多發展和改良。
展開 接巖土工程領域有限元、離散元代做，擅長軟件abaqus，liggghts等
本碩985，巖土工程頂尖院校畢業不滿一年，工作較清閑，接代做，可直接私信或聯系qq1223558239
淺談有限元方法的核心思想：數值近似和離散化
有限元法是基于近代計算機的快速發展而發展起來的一種近似數值方法，用來解決力學，數學中的帶有特定邊界條件的偏微分方程問題。而這些偏微分方程是工程實踐中常見的固體力學和流體力學問題的基礎。
有限元和計算機發展共同構成了現代計算力學的基礎。有限元法的核心思想是“數值近似”和“離散化”，所以它在歷史上的發展也是圍繞著這兩個點進行的。
1. 數值近似
由于在有限元法被發明之前，所有的力學問題和工程問題中出現的偏微分方程只能依靠單純的解析解得到解答。這種方法對數學要求很高，而且非常依賴于一些理想化的假定。比如在土木工程中梁柱計算中出現的平截面假定，小應變假定，理想塑性假定。這些假定其實是和實際工程問題有很大偏差的，而且一旦工程問題稍微復雜一些我們就不能直接得到解析解，或者解析解的答案誤差過大。
而有限元法把復雜結構離散到有限個單元，再把這種理想化的假定和力學控制方程施加于結構內部的每一個單元，然后通過單元分析組裝得到結構總剛度方程，通過邊界條件和其他約束解得每個單元的反應，這樣就可以避免直接建立復雜結構的力學和數學模型了。其總過程可以描述為：
總結構離散化 — 單元力學分析 — 單元組裝 — 總結構分析 — 施加邊界條件 — 得到結構總反應 — 單元內部反應分析
在進行單元分析和單元內部反應分析的時候，形函數插值和高斯數值積分被用來近似表達單元內部任意一點的反應，這就是有限元數值近似的重要體現。一般來說，形函數階數越高，近似精度也就越高，但其要求的單元控制點數量和高斯積分點數量也更多。另外單元劃分的越精細，其近似結果也更加精確。但是以上兩種提高有限元精度的代價就是計算量幾何倍數增加。
為了提高數值近似精度同時盡量較少地提高計算量，有限元法經歷了很多發展和改良。
展開 ABAQUS中泰森多邊形Voronoi和有限離散元FDEM結合的晶體斷裂仿真
《ABAQUS中泰森多邊形Voronoi和有限離散元FDEM結合的晶體斷裂仿真》
作者：星辰北極星
這個專題是依托于POLARIS_Voronoi插件制作的一套仿真案例視頻，講述Voronoi多邊形結合FDEM在晶體仿真中的一些應用；FDEM是FEM和DEM的一個組合縮寫，也就是“有限離散元方法”，結合了有限元和離散元的特征，在ABAQUS中主要通過大量嵌入Cohesive單元來實現，這一方法目前廣泛應用于巖石、玻璃、陶瓷等脆性材料的破碎仿真。
【課程內容】
第1章：課程概述
第2章：POLARIS插件
2.1 POLARIS_Voronoi插件介紹
2.2 POLARIS_InsertCohElem插件介紹
第3章：ABAQUS-Standard隱式分析案例
3.1 基于Cohesive單元的彈塑性斷裂仿真基礎
3.2 平面二維晶體試件的彈塑性拉伸斷裂仿真（二維多邊形）
第4章：ABAQUS-Explicit顯式分析案例
4.1 晶體試件的切削仿真（三維多棱柱）
4.2 圓柱多晶體試件的壓縮破碎仿真（三維多面體）
【案例：晶體拉伸斷裂仿真】
本例采用ABAQUS/Standard隱式計算方法，模型為平面二維多邊形，Voronoi控制點的分布是非均勻的，兩邊密，中間稀疏，類似于金屬材料經過表面處理后的晶粒細化，這種模型需要人為指定晶體控制點位置才能實現；此外，模型中的實體單元采用彈塑性材料的，因此是一種基于Cohesive方法的彈塑性斷裂分析的案例（本案例已經添加到Cohesive專題中）。
【案例：晶體切削仿真】
本例采用ABAQUS/Explicit顯式動力學分析方法。
展開 在ABAQUS中使用Python腳本將有限元離散單元轉化成SPH單元方法介紹
以將CAE中的C3D8R單元轉換為PC3D為例：
使用ABAQUS建模離散為C3D8R單元，然后生成input文件。之后用Python腳本將進行處理轉換即可（腳本見附件）。
這里需要注意的是Python腳本轉換后的input文件只是將有限元離散單元轉換為sph單元的文件，還需要自己去修改后才能用。
SPH.zip
腳本運行方法：
abaqus python solidtosph.py -inp <inputFileName> -part
有限元熱力耦合分析？
Abaqus運行后顯示User subroutine utemp is missing Abaqus/Standard Analysis exite 應該怎么解決
有限元程序-熱力耦合彈性動力學 ￥19.89
摘要
熱力耦合的應用在科學技術中有重要的意義。熱應力和它所引起的強度、剛度問題，在航空、航天和核反應堆工程的設備和構件上的重要性是不言而喻的。所以我們要對其進行研究和求解。
本文采用線性有限元建模技術對熱環境下的梁結構建模，求解一個線性熱彈性問題。在熱彈性狀態下，溫度場與機械場不耦合，而機械場取決于溫度，因為熱彈性本構關系中存在熱應變。這種情況可以描述為弱熱力耦合。本報告將討論瞬態演化問題的完全熱力耦合。在給出溫度場的基礎上，給定彈性力學的邊界條件和初始條件后求解熱彈性運動微分方程，得到熱位移場。然后，再由溫度場和熱位移場，根據應力、應變和溫度關系的本構方程，求出熱應力 場。通過分析得出，由于左右橫向邊界ΔT=+50 的均勻溫升，隨著溫度的增加機械場中的形變量增大，進而使應力增加。
關鍵詞 耦合熱彈性;線性有限元建模;本構方程
1.1課題背景
隨著人類文明的進步和科學技術的迅速發展，傳統的單一功能材料已經不能滿足科學技術和工程實際的需求。20 世紀以來，許多高性能的新型材料開始
扮演著越來越重要的角色。它們具有輕質、高強、耐久、智能等多重優點而 且，一般而言，材料和結構通常都是在高溫和有限制的環境中使用，在這種
情況下必須考慮材料和結構的熱力學性能。顯然，對這類材料和結構的研究不能完全套用經典的連續介質力學理論，而需要發展相關的理論來合理描述材料的力學性能。
熱彈性力學的應用，在科學技術中有重要的意義。熱應力和它所引起的強度、剛度問題，在航空、航天和核反應堆工程的設備和構件上的重要性是不言而喻的。
展開 有限元熱力耦合調試
User subroutine utemp is missing Abaqus/Standard Analysis exited with an error - Please see the message file for possible error messages if the file exists. 數據提交檢查可以正常，但是提交以后會中斷顯示上面的東西，求教大神指點
有限元熱力耦合，顯示下面警告怎么解決？
Whenever a translation (rotation) dof at a node is constrained by a kinematic coupling definition the translation (rotation) dofs for that node cannot be included in any other constraint including mpcs, rigid bodies, etc. The option *temperature is used but the option *initial conditions,type=temperature is not. The initial temperature values are assumed to be zero. MPCS (EXTERNAL or INTERNAL, including those generated from rigid body definitions), KINEMATIC COUPLINGS, AND/OR EQUATIONS WILL ACTIVATE ADDITIONAL DEGREES OF FREEDOM Node set assembly_m_set-10 has no members and will be ignored. Nodes defined in this set may have been deleted because they were not connected to any elements. Node set assembly_m_set-11 has no members and will be ignored. Nodes defined in this set may have been deleted
展開 流固耦合（FESIM有限元分析）
將所有之流固耦合相關的參數、邊界與荷載等均融入流場與固體所共享之控制方程式組內，再采用數值計算法，求解耦合的(coupled)聯立方程式，故作較嚴格定義的文獻認為此法才是真正的流固耦合分析。此種方式雖最完整，但難度也最高，故使用者最少。
