輸出時間間隔的案例
LS-DYNA常用控制卡片詳解
*DATABASE_BINARY_OPTION（二進制文件的輸出設置）
DT/CYCLLCDT/NR BEAM NPLTC PSETID ISTATS TSTART LAVG
【DT/CYCL】輸出時間間隔/時間步長輸出間隔。
【LCDT/NR】dump文件之間的時間間隔/重啟文件數量（默認為1）
【BEAM】把彈簧阻尼單元的坐標和合力寫進D3PLOT和D3PART。
【NPLTC】DT=ENDTIME/NPLTC，覆蓋前面指定的DT。
【PSETID】SET_PARTID,只適用于D3PART。
【ISTATS】設置統計量級別，只適用于D3MENA，默認為0，不收集統計量。
【TSTART】設置收集統計量時間從哪里開始，只適用于D3MENA，默認為0。
【LAVG】設置統計量寫出時間間隔，只適用于D3MENA，默認為100。
展開 HyperMesh&LS-DYNA 控制卡片詳解（非原創）
【TSTAPT】僅用于D3BEAM，設定模擬開始的時間。默認為0。
【IAVG】設定寫出平均數據的間隔，僅用于D3BEAM。默認為100。
【IOOPT】僅用于D3PLOT的選擇。
EQ.1：在這時刻每個plot產生，載荷曲線的值也被加進來到當前的時刻，來決定下一個plot的時間。這個為默認的。
EQ.2：在這時刻每個plot產生，下一個plot的時間T被算出來，T=當前的時間+載荷曲線值在T時刻。
EQ.3：載荷曲線里每個縱坐標都產生一個plot。曲線準確值被忽略。
10.*DATABASE_BINARY_D3THDT（單元子集的時間歷程數據輸出控制）
【DT】輸出的時間間隔。
【LCDT】指定輸出時間間隔的曲線。
11.*DATABASE_BINARY_INTFOR（接觸面二進制數據輸出控制）
【DT】輸出的時間間隔。
【LCDT】指定輸出時間間隔的曲線。
12.*DATABASE_EXTENT_BINARY（輸出數據控制）
指定要輸入到D3PLOT、D3PART、D3THDT文件中的二進制數據。
【NEIPH】寫入二進制數據的實體單元額外積分點時間變量的數目。
【NEIPS】寫入二進制數據的殼單元和厚殼單元每個積分點處額外積分點時間變量的數目。
【MAXINT】寫入二進制數據的殼單元積分點數。如果不是默認值3，則得不到中面的結果。
【STRFLAG】設為1會輸出實體單元、殼單元、厚殼單元的應變張量，用于后處理繪圖。對于殼單元和厚殼單元，會輸出最外和最內兩個積分點處的張量，對于實體單元，只輸出一個應變張量。
【SIGFLG】殼單元數據是否包括應力張量。
EQ.1：包括。（默認）
EQ.2：不包括。
展開 Ansys/Ls-Dyna進行沖擊器活塞撞擊鉆頭的動力分析
接觸時間0-100000s范圍內
接觸剛度=0.1.剛性盾罰因子=0，接觸面初始穿透檢查
接觸不考慮單元厚度，接觸剛度選擇
不考慮接觸過程殼厚改變，接觸表面節點號自動定位
接觸面穿透深度增量=4，不穿透檢查
7.定義求解參數
終止時間、.RST文件的輸出時間間隔、.HIS文件的輸出時間間隔、輸出節點接觸力、
輸出合成界面接觸力、輸出材料能量計算結果、保存數據庫文件等.
8.進入求解器，輸出K文件.
例如:本段APDL文件
/SOLU
!進入求解器
EDWRITE,LSDYNA,'hzI','K',''
!輸出關鍵字K文件"hzI-a"
FINISH
/EXIT,NOSAV
二.
進入Ansys
Product
Lauhcher
點Run入計算界面，輸入APDL文件，輸出'hzI','K'，等等文件。
展開 輸出方法：abaqus最值隨時間變化
最近做項目，需要輸出abaqus中最大應變值隨時間變化的曲線，但是翻了很多帖子都沒找到相關的，所以后來自己寫了個Python代碼，用以提取最大值。需要用到這個功能的同學把代碼里的注釋刪掉直接運行就行。
viewportName = session.currentViewportName
j=getInput('how much frame','0') #要提取多少幀的最值，默認為0
#這一步也可以用在step上，我的需求始終是step-1，所以沒寫那部分代碼，和幀數的代碼同理
a=int(j)+1
i=1
while i<a:
cobject = session.viewports[viewportName].odbDisplay.setFrame(step=0, frame=i)
#有興趣的朋友可以print一下cobject1，看看odbDisplay里面包含了哪些函數
i+=1
cobject = session.viewports[viewportName].odbDisplay.contourOptions
#同上，想學習的朋友可以print一下cobject,
print("step=1", "frame=", i-1)
print(cobject.autoMaxValue)#輸出最大值
print(cobject.autoMinValue)#輸出最小值
值得注意的是，代碼里step和frame的值都比實際值要小1，比如我要提取本來是step-1的數值，但是代碼里要寫成step=0。當然我的代碼并不完整，這樣輸出的max/min只顯示在abaqus底部信息欄，還需要在Excel中處理數據。
展開 
RADIOSS 整車碰撞模型轉換方法
邊界條件的加載
整車模型邊界條件的加載包含整車基本模型、壁障的加載及相關設置、整車自接觸以及與壁障/地面之間的接觸 設置、整車模型需要重點研究的截面力、加速度、應力、應變以及零部件能量變化等信息要求的設置和輸出。不同軟 件均有其獨自的設置方式，且相關的控制參數均具有其軟件包含的獨自的特性。
模型計算的控制卡片
模型轉換時，HyperCrash 會自動生成一個求解計算的 engine 控制文件，根據不同的要求對求解計算控制參數 進行詳細設定，主要包括計算時間、計算步長控制、輸出信息的類別設定等。LS-DYNA 與 RADIOSS 的不同之處在 于，前者可以單獨給予各項時間歷程數據信息不同的輸出方式和輸出時間間隔，而后者對所有的時間歷程數據信息只 能采用相同的輸出方式和輸出時間間隔。RADIOSS 可以通過不同的 engine 文件對模型進行各個時間段的獨立求解。
在此以常用工況 50FFB、64ODB、50MDB 為例，結合相應的試驗數據，經過校核有限元模型中零部件的幾何 形狀、材料和厚度，調整臺車與整車的相對位置，利用 RADIOSS 有限元軟件進行計算，得到對標分析結果。
壁障臺車信息
RADIOSS 格式的偏置碰固定壁障為 solid 壁障，總質量為 26.1kg，節點總數為 49180，單元總數為 48311，其 中殼單元總數為 9310，實體單元總數為 39000，并包括一個剛體單元，如圖 2 所示。
RADIOSS 格式的側碰可移動變形壁障總質量為 950kg，其中分為前面碰撞塊 24kg，臺車 926kg。模型的基本 信息中，節點總數為 8629，殼單元總數為 348，實體單元總數為 6300，如圖 3 所示。
展開 Dyna中模擬材料失穩的GISSMO失效模型 ￥20
材料失穩時應變，損傷，應力三軸度急劇增加
損傷閥值DCRIT設定為1.0時計算結果如下：
材料不發生失穩時單元幾乎同時失效，有一定隨機性
材料不發生失穩時無法構建應力下降段
結果對比
當材料出現失穩時，塑性應變，損傷，應力三軸度在很短的時間內急劇增長，需在失穩時間內加密輸出頻率才能捕捉到峰值（損傷閥值0.5的例子在最小時間步長的輸出時間間隔內依然未捕捉到損傷為1的時刻點）。
Dyna中GISSMO和JC失效模型比較
JC失效模型無法模擬材料失穩過程，損傷計算只能是線性累積，失效應變與應力三軸度只能是單調的關系； GISSMO失效模型可以模擬材料失穩過程，損傷計算為非線性指數累積，失效應變與應力三軸度可以不是單調的關系。是一種更符合實際的失效模型。
展開 LS-DYNA中的操作及設置（一）（精度控制，單位制，截面，雙精度）
假如用戶非常保守，大大降低了時間步長（一個或兩個數量級），那么除非使用雙精度計算，否則會影響到計算精度。
輸出plot文件時注意
使用雙精度運行計算任務時，程序輸出的文件也將是雙精度的。由于輸出文件的大小是原來的兩倍，所以在使用FEMZIP 或 plotcprs等進行后處理時會出現一些問題。因此使用雙精度程序計算時，應該以32位格式輸出相關文件。
設置方法：
1.修改K文件: *DATABASE_FORMAT, IBINARY=1
2.修改環境變量: export LSTC_BINARY=32ieee
3.使用命令行: e. g. ls971 i=input 32ieee=yes
展開 圓形銅柱Taylor沖擊測試仿真的EFG算法實現
圖2計算模型
3.2使用的相關關鍵字
$EFG計算控制，采用默認值，核函數選項為立方樣條函數
0 0 12
$輸出控制
0 0 0 0 0 2 1000
$計算結束時間
0.00090000
$能量信息輸出
*DATABASE_GLSTAT
0.500000
$材料信息輸出
*DAIABASE_NODOUT
0.500000
$結果數據輸出時間間隔
*DATABASE_BINARY_D3PLOT
5
0
$BOX方式定義初始速度
*INITIAL_VELOCITY
0.000 0.000 1
0.000 0.00 -200
$定義BOX1，包括除底部節點外的所有節點
*DEFINE_BOX
$PART定義
*PART
1 1 1
*SECTION_SOLID_EFG
1 1 41
1.3 1.3 1.3
$銅柱采用彈塑性材料本構
*MAT_PLASTIC_KINEMATION
1 8930.000 1.170000e11 0.350000e6 0.004000e6 0.001000e6 1.000000
0 0 0 0
*ELEMENT_SOLID
*END
3.3后處理顯示
本質上，EFG和FEM兩種算法共享同樣的數據庫，因此采用LSPP讀入2.k結果文件后，將FEM網格顯示為EFG質點時，可以通過Page1的Appear按鈕選擇Pick Part，之后更改顯示方式為Sphere，點擊AIIVis，Done就可以顯示EFG粒子，此外，可以使用Setting按鈕調整質點的顯示效果，如圖3所示。
展開 顯式動態分析采用的命令
EDRD
：剛體和變形體之間的相互轉換
EDREAD
：把LS-DYNA的ASCII輸出文件讀入到POST26的變量中
EDRI
：為變形體轉換成剛體時產生的剛體定義慣性特性
EDRST
：定義輸出RST文件的時間間隔
EDSHELL
：定義殼單元的計算控制
EDSOLV
：把“顯式動態分析”作為下一個狀態主題
EDSP
：定義接觸實體的小穿透檢查
EDSTART
：定義分析狀態（新分析或是重啟動分析）
EDTERM
：定義中斷標準
EDTP
：按照時間步長大小繪制單元
EDVEL
：給節點或節點組元施加初始速度
EDWELD
：定義無質量焊點或一般焊點
EDWRITE
：將顯式動態輸入寫成LS-DYNA輸入文件
PARTSEL
：選擇部件集合
RIMPORT
：把一個顯式分析得到的初始應力輸入到ANSYS
REXPORT
：把一個隱式分析得到的位移輸出到ANSYS/LS-DYNA
UPGEOM
：相加以前分析得到的位移，更新幾何模型為變形構型
展開 【推薦】手把手教你用ABAQUS玩轉魔術 ￥1.98
將場變量輸出改為等時間間隔輸出100次，如下圖。
圖為場輸出變量設置
五、相互作用
在戒指中心的圓心位置創建RP（參考點），并對戒指創建剛體約束，剛體約束參考點為創建的參考點。
接觸屬性選擇切向接觸，摩擦系數penalty取0.05。
創建通用接觸general contact，接觸屬性選擇剛才創建的接觸屬性，特別注意的是，在surface property中進行設置如圖設置。
圖為通用接觸屬性界面
為了使項鏈重量更重，模擬效果更好，在整個項鏈上創建非結構質量。
圖為非結構質量設置
六、載荷與邊界條件
在戒指上端施加固定約束條件。
整個模型施加豎直向下的重力載荷。
在戒指剛體約束參考點上施加預定義角速度場。
邊界條件示意圖
角速度預定義場參數設置，此處特別注意轉軸的坐標位置通過戒指剛體約束參考點。
七、網格
使用C3D8R對戒指進行網格劃分，使用T3D2單元對項鏈進行網格劃分。
八、提交運算
新建job，提交分析。
九、問題思考
1、項鏈材料參數中的dampling對計算結果是否有影響？有何影響？
2、非結構質量參數的作用是什么？對計算結果的影響？
3、通用接觸設置是的surface property應如何設置？
通過修改模型參數，進行對比分析很容易理解上述問題，歡迎下載本例源cae文件試一試吧
本例cae文件的網格可根據計算機能力進行修改，可提高運算效率，更方便對比分析
如有問題歡迎留言討論吧
展開 模擬流體中的粒子運動時，選擇合適的公式以提升計算效率
粒子追蹤模型中的數值剛度處理
有兩種主要方法可以求解流體中的粒子運動的數值剛度模型，即輸出時間間隔比 τ_p 大幾個數量級的模型。
第一種是我們所說的“強力”方法：只需告訴求解器采取更小的時步即可。如果不想產生大量的輸出（可能會創建大量文件大小），那么可以不考慮輸出時間，而是在求解器序列的瞬態求解器設置中指定一個較小的步長或最大步長。
另外一種從 COMSOL Multiphysics? 5.6 版本開始提供的方法，可能是從方程4 中刪除慣性項。首先，我們把方程4 改寫成一對耦合的一階方程
然后，僅假設曳力始終與其他作用力處于動態平衡，而不是在 τ_p 最初一段時間完全解析粒子運動，
(6)
換句話說，我們僅假設粒子立即達到其自由沉降速度。如果達到自由沉降速度所需的時間比總仿真時間小許多數量級，那么這是一個合理的近似值。方程6 可用于求解 v，
或者，一般而言，
(7)
其中，Fother 是除了曳力以外的其他所有作用力的總和。
然后，我們要做的就是把 v 的表達式對時間進行積分以獲得粒子位置 q。
我們可以在粒子釋放和傳播部分，選流體流動顆粒跟蹤 接口求解的方程組。
展開 
基于LS-DYNA的平頭彈沖擊間隙雙層靶數值模擬
彈靶之間采用如下圖所示的三維面對面侵蝕接觸算法
除此之外，還要設置模型的終止計算時間以及結果輸出的時間間隔。
模型設置完成之后即可導出K文件，并利用ANSYS中的LS-DYNA求解器進行求解。
3、數值結果分析
以上為雙層靶板受平彈沖擊后米塞斯應力隨時間的變化展示。可以看出，頂層靶板首先被彈體沖擊發生擊穿現象，應力響應隨時間向周圍擴展，并且與平彈底部接觸的靶板單元發生整體脫落進而作用到底層靶板，底層靶板首先因受到脫落碎屑的沖擊而產生應力，但是這種沖擊載荷較小未使靶板失效，隨后平彈的繼續沖擊作用于底層靶板，最終導致靶板發生擊穿。本模型能夠較為準確地模擬出平彈沖擊靶板的真實效果，能夠為相關內容提供思路和參考。
展開 重啟動－－－for LS-DYNA
1,簡單重啟動；當ls-dyna運行過程中還沒有到達終止時間就被人為或其它原因中斷，需要重啟動繼續進行計算，在這種情況下，不需要對輸入文件作任何改動。
2，小型重啟動；當重啟動分析時，希望對關鍵字進行一些修改，如：重新設置求解終止時間，重新設置各種輸出文件的時間間隔，刪除接觸界面，把變形體轉換為剛體，改變阻尼選項等。SimWe仿真論壇.o c~ Y gL jp
3，完全重啟動；當要對關鍵字文件做出大量的修改時，如增加其它的part或接觸定義等，此時實際上是進行另一個全新的分析。在重啟動關鍵字文件中包含的關鍵字是在上次求解的基礎上對所關心的part進行變形和應力的更新。
展開 手把手教會你如何完成S7-200SMART與PID的應用
0
2
PID主要參數
采樣時間：CPU必須按照一定的時間間隔對反饋進行采樣，才能進行PID控制的計算。采樣時間就是對反饋進行采樣的間隔。短于采樣時間間隔的信號變化是不能測量到的。過短的采樣時間沒有必要，過長的采樣間隔顯然不能滿足擾動變化比較快、或者速度響應要求高的場合。
增益（Gain，放大系數，比例常數）
增益與偏差（給定與反饋的差值）的乘積作為控制器輸出中的比例部分。提高響應速度，減少誤差，但不能消除穩態誤差，當比例作用過大時，系統的穩定性下降。
積分時間：
偏差值恒定時，積分時間決定了控制器輸出的變化速率。積分時間越短，偏差得到的修正越快。過短的積分時間有可能造成不穩定。積分時間的長度相當于在階躍給定下，增益為“1”的時候，輸出的變化量與偏差值相等所需要的時間，也就是輸出變化到二倍于初始階躍偏差的時間。如果將積分時間設為最大值，則相當于沒有積分作用。
微分時間：
偏差值發生改變時，微分作用將增加一個尖峰到輸出中，隨著時間流逝減小。微分時間越長，輸出的變化越大。
展開 Workbench Ls-Dyna 參數化霍普金森桿(SHPB)沖擊仿真-2D ￥30
在分析設置中，求解時間設置為大于應力波在1.5走過1.5倍桿長所需的時間，建議設置為走過2倍桿長所需時間（本例設置為8e-4s）；根據電腦配置設置CPU個數及內存信息。全局沙漏控制類型為Flanagan-Belytschko Stiffness Form，數值為0.1。根據自己需求設置所需輸出的動畫個數（默認為20個，本例設置為40），并在時間歷史輸出中增加Elemate Data輸出，輸出時間間隔2e-7s。</p><p>E. 后處理中添加輸出應變片單元的總機械應變（或y方向機械應變），并采用圖表功能把入射桿應變片應變和透射桿應變片應變顯示在一張圖中。提交計算，本算例網格較粗（實際可足夠細化，2DS軸對稱計算很快），計算時間不到半分鐘。</p><p>3. 結果展示</p><p>A. Workbench后處理結果，數據提取自D3plot，只有40個點。入射桿、透射桿應變及能量曲線圖如下。
展開 