由于ALE梁單元只有一個積分點,我們在*DATABASE_EXTENT_BINARY中使用BEAMIP,以類似于NEIPH的方式告訴LS-DYNA有多少數據需要輸出到d3plot文件中。通常情況下,BEAMIP只是告訴LS-DYNA需要輸出多少個積分點的數據。通過這種方法,我們可以為元素中的每個 ALE 組獲取 6 個應力、1 個塑性應變和一些歷史變量,這些變量取決于 ALE 組的材料和狀態方程。
mises應力
for s_mises in s.values:
vtk_file.write("{}\n".format(s_mises.mises))
對于C3D4單元或C3D8R單元的積分點就只有一個,以上語句的添加不會有任何問題,但是對于C3D8單元,每個單元有8個積分點,寫起來就需要進一步操作了,我選擇的方法是,將單元中8個積分點的應力做一個平均,相當于每個單元有一個平均應力