复合材料广义热弹性问题格点型有限体积法研究  

Vertex-centerfinite element method for generalized thermoelastic analysis of composite materials

作　　者：刘琦 肖光明[1] 齐文亮 杜雁霞[1] 刘磊[1] LIU Qi;XIAO Guang-Ming;QI Wen-Liang;DU Yan-Xia;LIU Lei(State Key Laboratory of Aerodynamics,China Aerodynamics Research and Development Center,Mianyang 621000,China;College of Power and Energy Engineering,Harbin Engineering University,Harbin 150001,China;Xi’an Aeronautics Computing Technique Research Institute,AVIC,Xi’an 710068,China)

机构地区：[1]中国空气动力研究与发展中心空气动力学国家重点实验室,绵阳621000 [2]哈尔滨工程大学动力与能源工程学院,哈尔滨150001 [3]中国航空工业集团公司西安航空计算技术研究所,西安710068

出　　处：《四川大学学报（自然科学版）》2022年第2期111-122,共12页Journal of Sichuan University(Natural Science Edition)

基　　金：国家自然科学基金(11972359,11972361);国家数值风洞工程;陕西省自然科学基础研究计划(2022JQ-564)。

摘　　要：基于Lord-Shulman(L-S理论)、Green-Lindsay(G-L理论)及经典热弹性耦合理论(T-C理论),本文发展了一种适用于二维复合材料广义热弹性问题的格点型有限体积法(CV-FVM).运用交错网格技术,待解变量存储在单元节点,物性参数存储在单元中心,控制方程空间项采用双线性四边形单元进行离散,时间项采用欧拉隐式进行离散.针对均质无限大方板热冲击问题,本文对CV-FVM进行了数值验证,计算结果表明,本文发展的数值方法可以很好地捕捉热波波前和弹性波前的温度阶跃特性及热弹耦合特征.本文采用CV-FVM进一步研究了不同梯度系数p下钛合金/氮化硅复合板的热冲击问题,发现L-S理论复合板p=1时预测的功能梯度界面应力最小;T-C、G-L理论下,复合板p=10时预测的功能梯度界面应力最小,不同耦合理论选系数p对界面热弹性行为影响规律并不相同,并非材料线性分布(p=1)时界面应力幅值最小.本文发展的CV-FVM可作为复合材料热波问题和广义热弹性问题求解的一种备用工具.To analyze the generalized thermoelastic problem in two-dimensional composite material, a new two-dimensional vertex-center finite volume method(CV-FVM) has been developed based on Lord-Shulman(L-S)、Green-Lindsay(G-L) and traditional coupled theories. Using the staggered grid technique, the unknown variable is defined at the cell vertex, while the material property is defined at the cell center. The space terms of governing equations are discretized by bilinear quadrilateral element, and the time terms are discretized by Euler implicit formula. Thermal shock problem in infinite plate with homogeneous material is studied by CV-FVM, The results show that the present method can effectively capture the temperature jump and thermoelastic coupling characteristics at the front of the thermal wave and elastic wave. Then, the developed CV-FVM was used to study the thermal shock problem in composite with Ti-6 Al-4 V/ZrO;with different material constant p, the results show that the value p=1 minimizes the maximum(tensile) stress applied at the middle of the functionally graded layer under L-S theory, and the value p=10 minimizes the maximum(tensile) stress under G-L and T-C theories. The effects of p on interfacial thermoelastic response is different under different coupling theories, one cannot conclude that a linear variation of the properties minimizes the maximum stress. The developed method can be used as an alternative tool for solving thermal wave and generalized thermoelastic problems.

关 键 词：有限体积法 非傅里叶热传导 复合材料 热冲击 广义热弹性理论 

分 类 号：O341[理学—固体力学] O343.6[理学—力学]