紧急求助：温度影响下的柔度最小拓扑优化

ajiang

对于目标为柔度最小化，在受温度影响和不受温度影响情况下相同吗？
1、我现在在矩形的中心位置X方向施加作用力为1，这个在不受温度影响下我可以求最小柔度
2、在温度影响下，目前我采用整体一致均匀变化的温度场，每个节点都会受到热应力，是不是我把这些节点的热应力计算出来，同时在矩形的中心位置X方向这个节点上还需要再加上原来的作用力1带入U (freedofs,:) = K (freedofs,freedofs) \ F (freedofs,:);（K是刚度矩阵） 求出位移， 再带入目标函数c = c + x(ely,elx,m)^penal*Ue'*KE*Ue;求出的就是温度影响下的柔度最小吧？
可是matlab画出的图柔度值却是越来越大？？？什么原因呀？扣扣:358939206有志同道合拓扑优化研究者可以多认识交流，呵呵