TATANICYELLOW的博客

原创 利用有限元法求解一维热传导问题

首先是进行网格划分，网格个数当然是越多越好，但是求解的时间也会相应大大增加，比如在下方程序中因为我把网格数设置为了500，相应的刚度矩阵就是501行501列，再加上我将时间步长设置较小，往往运行一次需要几分钟时间（不过我也有在考虑利用pytorch进行GPU加速来减小运行时间，希望可行吧）。实际上，运用有限元更简单的例子可能是稳态情况下的一维热传导问题，但是，笔者受外界不可抗力的影响，先上手的是非稳态的问题（不过稳态也没什么难度）。先写这些吧，接下来的一周要着手完成二维热传导的问题，祝我好运！