焊点热循环仿真分析中遇到的问题

最近在做焊点热循环仿真分析， 遇到以下问题，请大家多多帮助，谢谢了！1. 计算完成后如何选取应变能密度？ 我选择最后一个循环的应变能密度的差值，但是如果选择不同的载荷步得到的结果是不同的，对此不太理解。 例如， 如果选择倒数第二个载荷步的降温到-40度的deltaW作为Delta W1， 最后一个载荷步降温到-40度的DeltaW作为delta W2，计算得到的应变能差值与以125度计算的到的差值是不同的。 请一定给予帮助，多谢了2. 看了一些文献，在做焊点热循环时都没有考虑焊点的屈服， 这是为什么？  我在做的的时候，加屈服和不加屈服时的最大应力会相差近一倍

SAC 305， 采用稳态蠕变模型 Garofalo-arrhenius 模型如果增加屈服的采用双线性模型施加您也在做这方面的仿真吗

Anand 模型我也用过，个人觉得这两种模型模拟出来的结果相差不大你的能量密度的提取是怎么做的

Ananda模型包含两部分，一部分是流动方程，一部分时内部变量演化方程。 从你这个错误来看看应该还是材料的定义上出现了问题，你检查一下你的材料单位是不是统一。 Ananda模型中有几个参数的单位是MPa，和你建立的几何模型，其它材料模型是不是统一的？

你的热疲劳仿真做的怎么样了。 使用开氏温度计算的结果感到很奇怪，能不能帮我看一下哪里有问题？