温度梯度自应力计算

实际工程哦：
无意中在复核计算书时发现《公预规》JTGD60-2004给出的公式，计算的顶底面的应力比理论计算的有较大出入，顶面比理论计算的小10%，底面却比理论的大10%

@大海的怒爆，@所有土木人，说说什么原因造成的。本桥是一个人行桥，这个差异结果是不会影响结构安全的

你不是说输入有误？最后杀结果，请再粘一次图片

这个贴很早就看到了，不过工作中一直没有接触设计计算的机会，所以也就没有做啥研究，膜拜！最近部门接到个项目，钢混组合梁的，所以为了计算钢混组合梁的温度梯度应力，细细研究了下规范的方法，感觉搞工程的人弄出来的规范公式就是有点误人子弟，很简单的原理非要弄的这么复杂，所以我就在loupoo2老师的基础上，稍作修改。
目前混凝土截面的已经验证通过了，误差基本忽略，还差钢混组合梁的温度梯度应力。midas用的也不太熟，所以钢混组合梁的温度梯度应力总是怪怪的，所以用了ansys来分析对比，基本验证了我改写的正确性，最近整理的就贴上来。
不过又遇到问题了，ansys按照规范的温度梯度，给实体单元加温度怎么搞啊，线性变化的温度按单元分层赋值的话感觉又很复杂。@loupoo2 老师有没有好的建议？我还是喜欢用ansys，midas这软件我不大信，仅约束***这么多，工程人做的软件又恶心了我一回。