瞬态随机温度场的Neumann展开Monte-Carlo解法

王涛¹，周国庆^1,2

1. 中国矿业大学深部岩土力学与地下工程国家重点实验室2. 中国矿业大学力学与建筑工程学院

摘 要：利用Neumann展开Monte-Carlo随机有限元法对瞬态随机温度场求解进行了分析,根据计算流程框图编写了求解瞬态随机温度场的Matlab随机有限元计算实施程序.对于每一次随机抽样,文中方法只需形成刚度矩阵,进行前代、回代以及矩阵乘和矩阵加减,便可求得各时刻温度场.因计算中无需矩阵求逆,可大大减少计算量,有限元节点数目越大,计算速度优势越明显.为说明解法的可行性及高效性,算例考虑了导热系数、比热容、换热系数、热流密度、环境温度、介质密度以及内热源等物理参数单独和同时具有随机性的瞬态随机温度场情况.结果表明:Neumann展开Monte-Carlo随机有限元法能有效解决瞬态随机温度场问题,且效率较高;编制的Matlab随机有限元计算程序可直接输出各时刻温度场的统计结果;各随机参数对温度场响应随机性的影响有差异.

关键词：Neumann展开;Monte-Carlo有限元法;随机有限元;瞬态温度场;