(2) 模型文件：taiqiu.cae，為案例的Abaqus原生模型文件，可直接用Abaqus軟件打開，包含所有部件、材料、裝配、分析步等設置。 (3) 輸入文件：taiqiu.inp，Abaqus分析的輸入文件，可用于提交計算或二次修改模型參數。 (4) 結果文件：taiqiu.odb，案例計算完成后的結果文件，可直接用于后處理分析，查看應力、速度等關鍵物理量分布及曲線。

材料庫計算，右圖為vumat子程序計算結果）： 等效塑性應變對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）： 反力曲線對比如下： 塑性耗散曲線對比如下： 3.2 圓棒拉伸測試 對一圓棒骨料進行單軸拉伸，其邊界條件如下： von Mises應力對比結果如下（左圖為Abaqus材料庫計算，右圖為vumat子程序計算結果）

使用Everything軟件搜索Launcher.bat文件，注意需要是Abaqus下的這個文件，不是其他軟件的。打開文件路徑，使用記事本打開文件，進行修改Launcher.bat文件。

各部件的幾何形狀和尺寸取自Sawa等人(1991)，并稍加修改以簡化建模。這兩個輪轂和墊片的內壁半徑是25 mm。管道法蘭外壁半徑為82.5 mm，墊片外壁半徑為52.5 mm。墊片的厚度為2.5 mm。管法蘭上有8個螺栓孔，螺栓孔在半徑為65mm的節圓內等間距。在本分析中，螺栓孔半徑修改為與螺栓相同，為8mm。螺栓頭(承載面)假定為圓形，其半徑為12mm。

可見采用實驗模態法必須在實驗輪轂已經制造完成的前提下開展，且對結構進行優化修改更為困難，花費的資金與時間更多。 基于有限元分析的原理，有限元模態分析法可以通過計算機模擬，以質量矩陣、剛度矩陣和阻尼矩陣的形式，結合振動理論，對結構的力學性能進行有效的分析。

圓形給出了圓心坐標和圓上一點的坐標。我們可以對這行代碼稍加修改，完成批量隨機幾何建模。

為實現隨機顆粒模型的代碼編寫，我們分3步進行： ? 首先通過錄制獲得ABAQUS中建模的相關代碼。 ? 然后修改相應代碼，刪除無用代碼，實現代碼的參數化。 ? 最后在代碼中添加循環和干涉判斷，實現多個顆粒的隨機投遞。

當然楔形單元沒問題，也可以通過切分修改掃掠方向為軸向，此時的階梯軸必須從階梯位置斷開，這樣可以在軸的</span><span style="font-family:'Calibri';font-size:12.0pt;font-weight:bold;white-space:pre-wrap;">圓截面上劃分</span><span style="font-family:'Calibri';font-size

（29） 靜力問題依然自動輸出總質量 （30） 修改應變各分量和Abaqus一樣，但是max principal等值不一樣的bug （31） 修正主軸應變顯示和abaqus不一致的問題 （32） 支持聲學邊界元的輻射聲場計算 （33） 輸出Set集，支持后處理按set集查看結果 7.

s.VerticalDimension(vertex1=v[2], vertex2=v[4],textPoint=(0.115524396300316, 0.0262394621968269), value=0.0) s.ObliqueDimension(vertex1=v[5], vertex2=v[3],textPoint=(0.0519323498010635, 0.0), value=0.005) ?修改圓尺寸