ABAQUS計算周期的案例
abaqus簡單立方體胞元周期性邊界條件施加計算腳本源代碼 ￥39.9
</p><p class="ql-align-justify">在abaqus中通過file→run script選擇該腳本運行，腳本會輸出X，Y，Z三個方向的計算結果ODB文件。</p><p><br></p>
Abaqus三維周期性邊界和一般周期性邊界施加
針對ABAQUS周期性邊界手動施加繁瑣，復雜的問題，開發了兩款腳本文件，用于施加周期性邊界和一般周期性邊界。其中，周期性邊界的單元類型沒有任何限制；一般周期性邊界的單元類型需為四節點，如C3D4、C3D4R等。這兩款代碼，實現的效率比較高，對于節點數量在10W的模型，其需要的時間在1分鐘內（一般筆記本電腦）；計算結果合理，其測試模型為100mm*100mm*100mm的立方體，材料彈性模型為2.1e5MPa，泊松比為0.3，施加x向為5mm的拉伸位移，用周期性或一般周期性代碼進行施加邊界，具體如下圖所示。
作者QQ：2812468512
展開 LAMMPS分子動力學基于周期擾動法的黏度計算
關鍵詞：黏度,周期擾動法，SPC/E水分子，分子動力學，lammps
目前分子動力學計算黏度主要有以下方法：（1）基于 Green - Kubo 關系的方法。從微觀角度出發，利用壓力張量自相關函數積分計算黏度。理論基礎強，能考慮復雜微觀因素，但計算量極大，對計算機性能和時間要求高，積分上限選擇需謹慎。（2）愛因斯坦關系法。通過分析粒子擴散行為間接求黏度，依據愛因斯坦關系，由粒子擴散系數計算。計算相對簡單，只需粒子運動軌跡信息。不過對于粒子間作用強、擴散不明顯的復雜體系誤差較大，粒子半徑選擇也會影響結果。（3）非平衡分子動力學方法。對系統施加外力場或速度梯度使其處于非平衡態，根據非平衡條件下的響應，如應力和速度梯度，依牛頓黏性定律計算。可模擬復雜流場中黏度，但施加外力和處理邊界條件要注意，理論較復雜。
本文將介紹非平衡分子動力學計算黏度的方法之一：周期擾動法。周期擾動法計算黏度的原理主要基于對流體系統施加周期性的外力擾動，然后通過分析系統的響應來計算黏度。具體如下：
(1)施加周期擾動：在流體系統中引入一個周期性變化的外力場，比如在模擬體系的某個方向上按照一定的頻率和振幅施加外力。這個外力會使流體中的粒子或分子產生相應的運動和位移。
(2)系統響應分析：流體系統在受到周期外力擾動后會產生響應。由于流體具有黏性，粒子或分子的運動不會立即跟上外力的變化，而是會有一定的滯后。這種滯后表現為流體內部產生剪切應力，并且剪切應力會隨著時間周期性地變化。
(3)計算黏度：通過測量流體系統在周期擾動下的應力響應和外力的關系，結合相關的理論公式來計算黏度。在周期擾動法中，通常會測量系統在多個周期內的響應，并對數據進行統計分析，以獲得準確的黏度值。
接下來我們以計算水的黏度為例子，用lammps中的周期擾動法計算其黏度。
展開 周期性瞬態導熱有限元計算網格剖分規則研究
流固耦合
周期性瞬態導熱有限元計算網格剖分規則研究.pdf
內燃機機體內冷卻水腔的三維精確建模.pdf
自定義模塊：根據光導的導光條件計算光柵周期
摘要
為了滿足光波導的導光條件，在VirtualLab Fusion中生成了一個計算耦合光柵周期范圍的模塊。為了輔助設計基于波導的顯示器件，給定某個視場(FOV)作為所需的輸入參數。在該模塊中，利用光波導的全內反射限制和傳播光限制來計算可能的光柵周期范圍。
編程任務：在k域中定義視場
任務：生成一個計算耦合光柵周期范圍的模塊，以滿足平面光波導的導光條件。
說明耦合過程的平面光導圖
定義入射光空間頻率矢量的x、y分量為
而笛卡爾角α和β則用于定義一組入射方向的特定視場(FOV)方向。角度和方向之間的關系描述如下
編程任務：定義導光條件
k域圖說明導光條件
導模必須滿足導光條件，包括全內反射條件和傳播模條件
光柵是一種優良的耦合元件，因為在考慮光柵矢量G的情況下，FOV在k域中發生位移，進而可將導光條件推廣到
編程任務：計算周期范圍
k域圖說明導光條件
在一維周期光柵的情況下，光柵矢量的一個分量變為零，并且FOV總是可以旋轉到光柵的內部坐標系中，使???? = 0不失一般性。
在導光條件之后，將一定FOV耦合到光導中的一維周期光柵的周期范圍可以通過以下公式計算：
VirtualLab模塊的輸入
VirtualLab模塊的輸出
文件信息
延伸閱讀
- 單入射方向光導耦合光柵的優化
- 期望視場上用于光導耦合的二元光柵優化
- 期望視場上用于光導耦合的斜光柵優化
展開 自定義模塊：根據光導的導光條件計算光柵周期
摘要 為了滿足光波導的導光條件，在VirtualLab Fusion中生成了一個計算耦合光柵周期范圍的模塊。為了輔助設計基于波導的顯示器件，給定某個視場(FOV)作為所需的輸入參數。在該模塊中，利用光波導的全內反射限制和傳播光限制來計算可能的光柵周期范圍。
編程任務：在k域中定義視場
任務：生成一個計算耦合光柵周期范圍的模塊，以滿足平面光波導的導光條件。
說明耦合過程的平面光導圖 定義入射光空間頻率矢量的x、y分量為
而笛卡爾角α和β則用于定義一組入射方向的特定視場(FOV)方向。角度和方向之間的關系描述如下
編程任務：定義導光條件
k域圖說明導光條件
導模必須滿足導光條件，包括全內反射條件和傳播模條件
光柵是一種優良的耦合元件，因為在考慮光柵矢量G的情況下，FOV在k域中發生位移，進而可將導光條件推廣到
編程任務：計算周期范圍
k域圖說明導光條件 在一維周期光柵的情況下，光柵矢量的一個分量變為零，并且FOV總是可以旋轉到光柵的內部坐標系中，使???? = 0不失一般性。在導光條件之后，將一定FOV耦合到光導中的一維周期光柵的周期范圍可以通過以下公式計算：
VirtualLab模塊的輸入
VirtualLab模塊的輸出
文件信息
展開 Abaqus混凝土周期性邊界代表體單元插件：Random Sphere RVE 3D (Mesh) - AbyssFish ￥698
插件介紹
Random Sphere RVE 3D (Mesh) - AbyssFish 插件可在Abaqus生成三維具備周期性邊界條件(Periodic Boundary Conditions, PBC)的隨機球體骨料及骨料-水泥界面過渡區(Interfacial Transition Zone, ITZ)模型。即采用周期性代表性體積單元法(Periodic Representative Volume Element，PRVE)，以代表體積單元(Representative Volume Element，RVE)或稱為表征單元體(Representative Elemental Volume, REV)微觀結構的計算來準確地模擬和預測混凝土材料的宏觀行為。插件采用體素網格方式，通過背景網格將砂漿、骨料、ITZ劃分為三個集（Set），并對單元映射三種空材料。
插件支持設置長方體部件的長度(Length)、寬度(Width)、高度(Height)，以及在網格劃分中單元的尺寸(Element size)。可設置生成球體的最小粒徑(D_min)及最大粒徑(D_max)，即球體尺寸的分布范圍，球體占整個長方體試件的比例(Ratio)，界面過渡區的厚度(ITZ)，以及超時終止參數(Time)。
模型可分為砂漿基體、界面層、球體骨料三相材料。
插件生成的模型均滿足周期性分布邊界條件。
可對每個集(Set) 批量插入嵌入0厚度cohesive粘結單元（注：需要自行添加，本插件不具備此功能）。
?
展開 comsol中施加周期性邊界條件計算任意橫截面介質的導波頻散曲線 ￥1
<p>計算結果</p><p><img src="https://img.jishulink.com/upload/202010/9150fe93517a46adb4170dd665759fc5.png" alt="mmexport146300ec987f46a175e2841788c25b20.png"></p><p>紅色為本方法結果，黑色為半解析有限元法結果。</p><p>本方法建模速度快，計算精確，能實現任意橫截面介質的頻散曲線計算。</p><p>個人wx29996883 注明來意</p><p><br></p><p><br></p>
展開 基于ANSYS裂紋擴展模擬和生命周期預測計算實例（原創，如轉載，請注明出處）
分析類型：斷裂力學
技術難點：斷裂 裂紋擴展 生命周期預測
完成人：技術鄰ANSYS專家
網址：http://www.yqgqt.org.cn/content/other/402981
模擬過程：
裂紋擴展模擬和生命周期預測
abaqus模擬周期性邊界條件（單向纖維復材單胞） ￥19.89
本實驗在邊界上只有12個節點，總共添加了14個約束方程，但考慮到實際問題有很對對稱的節點，這時用這種手動添加的方法會非常麻煩，這種情況下需要編寫自動識別對應節點并添加周期性邊界的腳本。</p><p><br></p>
Abaqus周期幅值曲線
在Abaqus的分析中，有時我們不能采用表格類型的幅值曲線定義載荷曲線，而需要采用周期幅值曲線。下面我們簡單介紹下Abaqus中的周期幅值曲線定義
Abaqus中是以傅里葉基數定義的周期幅值，如下公式：
對應于Abaqus中的定義參數如下：
當上面的參數都給出后，傅里葉基數就確定了，當然關于t的曲線也就確定了。當N不等于1的時候可以給多個A和B的值。例如：對于X=sint可以如下圖進行定義
基于ABAQUS的旋轉周期對稱結構振動仿真
ABAQUS是一款功能強大又方便操作的通用有限元仿真軟件。本文主要介紹ABAQUS在旋轉周期對稱結構仿真中的便捷性。在ABAQUS環境下，通常我們都對結構的強度和振動進行仿真時，都將整個結構模型進行網格劃分，然后進行整體分析。但對于一些結構如光盤、風扇、輪胎，甚至是汽輪機轉子等的旋轉周期對稱結構，我們則不必對整個模型進行建模，而是可以截取其中的一個扇區，將其作為計算模型，進行適當的設置便可進行整個模型的振動仿真。
以一個空心盤為例。如下圖所示：
若我們對這個模型進行強度與振動仿真，我們只需截取其中的一個扇區，如截取其中1/72（即5°）的扇區如下圖：
將其導出并劃分好網格，再導入ABAQUS中，設置旋轉周期對稱條件便能仿真整個盤的振動了。具體視頻操作見鏈接：http://www.yqgqt.org.cn/college/video/c10169
在這給出視頻中的相應結果：
一階一節徑振型
一階二節徑振型
………………………………
展開 abaqus幅值中周期函數（傅里葉級數）的應用
最近系統具體的學習了相關幅值曲線的知識，分享一個周期函數的知識。附件有具體的傅里葉級數的講解。
ABAQUS交流群：1063594113
隨便在這吐槽一下：竟然有人把這種簡單的知識收費，也是厲害了！
周期信號的傅里葉級數.pdf
ABAQUS周期性邊界條件插件-EasyPBC1.4
EasyPBC is an ABAQUS CAE plugin developed to estimate the homogenised effective elastic properties of user created periodic (RVE).
各位有條件的同學可以去youtube上看一下作者的視頻，我不知道怎么下載，希望對各位的科研有所幫助
https://www.youtube.com/watch?v=pyZXvl3C6fI&pbjreload=101&ab_channel=SadikOmairey
https://www.youtube.com/channel/UCI2n0gg_4hFJ5XZ6eGG1yZA
EasyPBC V.1.4_5.zip
展開 Abaqus質量-彈簧-杠桿系統振動周期
在本教程中，我們將利用Abaqus軟件計算質量和彈簧系統的振動周期。我找了一個簡單點的模型，以便我們將Abaqus軟件的計算結果與理論計算結果進行比較。
問題描述如下，一根長1米的鋼桿繞O點鉸接，并通過兩個剛度系數為K和3K的彈簧在兩端連接。在此例中，我們以每秒10弧度的初始速度旋轉整個鋼棒，此模擬的目的是繪制隨時間變化的角位移圖并獲得旋轉周期。
建模過程如下：
首先我們需創建一根梁長度1m，將材料密度及截面參數調成與已知條件一致，即質量m=5.549kg；接下來需要創建兩端的彈簧，使用special spring/dashpot建立彈簧，注意左端彈簧剛度15000N/m，右端彈簧剛度為5000N/m；最后需要創建鉸接，這里可以使用connector來模擬，注意鉸點位置為0.25m處。
模型搭建過程
建立Dynamic， Implicit分析步，分析時長0.5s，為了較準確的捕捉運動狀態，設置固定增量步長0.001及歷時輸出連接器轉動位移UR1；設置梁的初始初始轉速10rad/s （考考大家最后的轉動的周期與初速度有關嗎？)。
結果
進行求解后，提取connector轉動角度曲線如圖。
轉動角度結果
如上圖所示，Abaqus軟件計算的結果可知系統周期為0.0925秒，如果我們用理論方程式和公式解決求解這個問題，同樣可以得到0.092s秒的周期（此處不進行推導，有興趣的同學可以試試哦）。可見，Abaqus求解的答案與理論結果完全一致。
轉自公眾號：ABAQUS仿真世界，歡迎關注！
展開 