浮點溢出的案例
LS—DYNA求解問題及解決方案
錯誤名稱
原因及修改措施
輸入格式不正確
找到錯誤位置，修改數據格式
浮點溢出，核心代碼崩潰
（1） 材料和單元特性定義錯誤，如密度、零厚度
（2） 單元定義錯誤
（3） 各種曲線定義有問題，如求解時間比定義曲線的時間長導致數據外推失敗
（4） 單位不一致
（5） 時間步太大，導致求解不穩定
程序一開始就崩潰
（1） 單位不一致
（2） 邊界條件、載荷曲線定義錯誤
（3） 某一材料未定義
節點速度無限大
（1） 應用LS-PREPOST檢查最后一幀動畫，找出崩潰的單元和PART
（2） PART的材料定義有問題，如某個參數設置過大或與該PART相關的接觸定義有問題，檢查接觸的設置
負體積
（1） 材料參數設置錯誤，選擇合適材料模式
（2） 沙漏模式的變形積累，改為全積分單元
（3） 太高的局部接觸力，需調整間隙、降低接觸剛度或降低時間步
展開 LS-DYNA求解過程中出現問題如何排查？
<p><span style="color: rgba(0, 0, 0, 0.88); background-color: rgb(246, 246, 246);">在 </span><strong style="color: rgba(0, 0, 0, 0.88); background-color: rgb(246, 246, 246);"><em>LS-DYNA求解 </em></strong><span style="color: rgba(0, 0, 0, 0.88); background-color: rgb(246, 246, 246);">過程中，</span>經常會出現奇奇怪怪的問題，例如，求解時間過長、求解時間在慢慢變長（與初始預計時間相比）、節點速度無窮大、負體積（畸變）、浮點溢出等等。要解決上述問題，就需要耐心檢查<span style="background-color: rgba(18, 18, 18, 0);">d3hsp和messag（mes000）文件，仔細編輯調試關鍵字文件，</span><span style="background-color: rgb(246, 246, 246); color: rgba(0, 0, 0, 0.88);">必要時還需要借助LS-PrePost等后處理工具加載d3plot動畫文件（檢查最后一步動畫出現的問題，網格畸變）、glstat文件（查看各類輸出算法正常）。
展開 Matlab中幾個數值積分函數的比較和優缺點（轉載Matlab技術論壇）
Infinite or Not-a-Number function value encountered'
意味著在積分計算時，區間內出現了浮點數溢出或者被零除。
例2 計算積分1/(x^3-2*x-p)，其中參數p=5，積分區間為[0,2]
%by dynamic
%all rights reserved by www.matlabsky.com
>>F = @(x,n)1./(x.^3-2*x-n);
>>Q = quad(@(x)F(x,5),0,2)%或者使用 quad(F,0,2,[],[],5)效果是一樣的，只是前者使用的函數嵌套
Q =
-0.4605
>>quad(F,0,2,[],[],5)
ans =
-0.4605
三、[q,fcnt] = quadl(fun,a,b,tol,trace,p1,p2...)
自適應Lobatto數值積分，適用于精度要求高，被積函數曲線比較平滑的數值積分
注意事項：
同quad
可能警告：
同quad
例3 計算積分1/(x^3-2*x-p)，其中參數p=5，積分區間為[0,2]
%by dynamic
%all rights reserved by www.matlabsky.com
>>F=@(x,p)1./(x.^3-2*x-p);
>>Q = quadl(F,0,2,[],[],5)%或者Q = quadl(@(x)F(x,5),0,2)
Q =
-0.4605
四、[q,errbnd] = quadgk(fun,a,b,param1,val1,param2,val2,...)
展開 