三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

多孔陶瓷等效导热系数的近似计算

提出了一个近似计算多孔陶瓷等效导效系数的公式。该公式综合考虑了孔隙大小，孔隙度以及气相辐射与导热的影响，可用于在任意温度下具有任意孔隙度的多孔陶瓷的等效导热系数的计算，用文献数据检验表明，该公式在高温区有较高的计算精度。     

涂层材料瞬态传热实验与导热系数辨识

针对涂层-基体一体化的双层结构,为测试评估其中涂层材料的导热性能,提出基于瞬态平面热源法（transient plane source, TPS）的涂层材料导热系数反演辨识方法。根据Hot-Disk实验测试原理,建立基体-涂层-探头整体的二维非稳态传热模型;结合测量过程中的瞬时温升数据信息,采用粒子群优化算法反演辨识获得涂层材料的导热系数;并通过实验和数值模拟论证了上述方法的可靠性。结果表明：该测量方法能够有效获得涂层导热系数,测试反演的数值偏差小于4.0%。最后,实际测量和反演辨识获得了一种涂层材料常温至773 K的导热系数,随温度提高呈现增大趋势,数值范围为0.18～0.29 W/(m·K)。     

基于高温瞬态热响应的石英窗口导热系数反演

提出一种利用高温加热下瞬态热响应信息来同时反演石英窗口在不同温度下导热系数的方法。结合实验测试过程,首先建立石英窗口-石墨板-石英窗口三明治结构高温非灰体辐射-导热耦合传热模型。再通过实际测试温度响应结果和材料光谱辐射特性参数,构建石英窗口导热系数的遗传算法反演辨识模型。采用1 100 K加热条件下实验数据,反演辨识获得了常温至1 100 K石英窗口导热系数介于1.35~2.58 $\mathrm{W}/(\mathrm{m}?\mathrm{K})$,并拟合出高温导热系数关联式。结果表明:该材料辐射对导热系数影响随温度升高而增大,在1 100 K下辐射影响占比16.12%。最终通过变工况的高温传热过程测试,验证了方法的可靠性及导热系数数据的准确性。      

计算滑动轴承动特性系数的有限元法

提出了一种快速求解8个轴承动特性的系数有限元法。该方法不仅比以往的有限元法计算速度快, 而且与其他方法具有同样的计算精度。      

一种考虑结构瞬态互辐射的有限元法

研究结构元件之间瞬态辐射过程及特点 ,对于元件温度未知时的瞬态辐射提出了一种有限元计算模型。通过广义边界条件简化了元件互辐射的数学描述 ,最后得到了广义互辐射元。所建模型具有简单、方便、实用和精度高等特点 ,使复杂的互辐射计算化为统一模式 ,并通过试验与实例计算证明其有效性。     

TC11钛合金压气机叶片锻后冷却过程瞬态温度场的数值模拟与试验研究

运用二维有限元法和试验方法,模拟、分析了TC11钛合金压气机叶片锻后冷却过程的传热现象,得到该叶片各个翼形截面上瞬态温度场的时变特性和分布规律.模拟结果为解决该类叶片锻后冷却过程翘曲变形提供了可靠依据.     

用电热类比法解导热系数随温度而变的导热问题

本文介绍一种简单而有效的方法, 即利用我们所熟知的λ变换法在电模型上解非线性问题(简称λ场法)。该法见特辛勃和施纳特的著作[10,4]。他们应用该法, 得出了导热系数和温度成线性关系的稳定导热问题的数值解, 并正确地指出了该法仅适用于第Ⅰ类边界条件及表面无换热的情况。在本文中, 我们将该变换法推广到解热传导问题的电模拟领域中来; 并利用了电模型固有的特点提出了对工程上最常见的第Ⅲ类边界条件的变换方法,成功地解决了求解导热系数和温度成线性关系的第重类边界条件问题; 文中还进一步提出了非线性的不稳定场的解法, 阐明了应用λ变换法在纯电阻网模型上求解导热系数和温度成线性关系的不稳定场的试验原理及试验方法。应用该法时可减少解题所需的工作量。根据所提出的变换方法, 在模型上计算了一个简单的第Ⅲ类边界条件下的稳定导热问题, 以阐明方法的应用及解题步骤。对于复杂的导热问题(如轮盘、叶片) 在电阻网模型上求解时, 方法当然仍是有效的。考虑到便于进行方法的比较, 故本文首先介绍了在电阻网模型上解导热系数随温度而变的稳定温度场问题的逐次近似法; 然后提出了把温度场转化为λ场的求解方法; 最后简单阐明了用λ场法在纯电阻网模型上求解不稳定场的方法原理。     

基于能量有限元法和虚拟模态综合法的高频冲击响应分析方法

为了分析结构受到高频冲击载荷激励后的瞬态响应,提出了一种基于能量有限元法（EFEM）和虚拟模态综合（VMSS）法的高频冲击响应分析方法。通过能量有限元法进行高频稳态分析,获得结构频响函数（FRF）的频段平均值,然后结合虚拟模态综合法得到虚拟模态振型系数,最后通过Duhamel积分获得结构在高频冲击载荷作用下的瞬态响应。对一简支梁模型进行算例分析,将本文方法的结果与传统有限元法（FEM）和统计能量法（SEA）的分析结果进行对比,验证了所提方法的有效性,也表明该方法具有模型简单、分析速度快等优点。     

基于有限元法的平纹编织C/SiC复合材料等效导热系数预测方法

提出了一种采用有限元法计算平纹编织C/SiC复合材料等效导热系数的方法.首先研究了材料的细观结构,根据材料显微照片建立了带基体碳纤维束复合材料的单胞模型,采用基体中孔隙分布随机生成的单胞模型,计算了孔隙率对基体等效导热系数的影响,通过施加3组边界条件计算出带基体碳纤维束和复合材料的等效导热系数.最后应用提出的方法计算分析了碳纤维体积分数和孔隙率对复合材料等效导热系数的影响规律.结果显示:复合材料等效导热系数随碳纤维体积分数增大而线性下降,碳纤维体积分数从54%增加到78%的过程中,复合材料y轴方向的等效导热系数下降了12.8%,x与z轴方向的等效导热系数同时下降了8.6%;复合材料等效导热系数随孔隙率增大呈加速下降趋势,孔隙率从0增加到30%的过程中,材料的x与z轴方向的等效导热系数下降了22.91%,y轴方向的等效导热系数下降了34.66%.     

一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

时域推进与边界元混合瞬态热传导问题

对时间变量在小区间内应用多项式插值,对空间变量采用边界元求解,提出一种新的分析瞬态热传导问题的计算格式。与现行的利用含时间变量的基本解所得到的边界元计算格式相比较,具有基本解简单易得、系数矩阵容易计算以及无重复积分等优点;与向前差分消去时间因子的边界元计算格式相比较,则在算法稳定性、计算精度以及时间步长的选取等方面具有明显的优势,并且很容易推广到求解非线性问题。算例表明,采用时域推进与边界元混合方法分析瞬态热传导问题有效可行。     

挤压油膜阻尼器的热平衡分析

利用变分有限元法对挤压油膜阻尼器（简称ＳＦＤ）的热平衡问题进行了分析计算。文中根据雷诺边界条件，通过联立求解ＳＦＤ的雷诺方程、能量方程得到了ＳＦＤ的压力分布、滑油流量等参数。为ＳＦＤ的设计与计算提供了更为可靠的理论依据。     

粘塑性材料结构的有限元分析方法

本文针对Ｗａｌｋｅｒ粘塑性本构方程,通过发展一自适应变步长积分算法,用户子程序的方法把该本构方程结合进了国际著名的非线性有限元软件ＡＢＡＱＵＳ,并进行了算例验证。它包括了四种单元类型,可考虑材料变形的率相关,加载的循环、非等温、变温等情况。适用于发动机热端部件非线性结构分析。     

卡箍刚度的有限元计算与实验测定

准确确定卡箍刚度对管路系统的振动特性分析具有重要意义。本文针对航空发动机管路系统中不同卡箍的结构特点,建立了计算卡箍刚度的有限元模型。对单管接地卡箍及双管悬地卡箍的刚度系数进行了有限元计算,还对一典型双管悬地卡箍刚度进行了实验测定。通过卡箍刚度计算结果与实验结果的对比表明,本文建立的卡箍刚度计算方法可用于管路系统振动特性分析工程实际之中。      

三维药柱的热粘弹性有限元分析

采用三维模型,用有限元法分析固体火箭发动机的药柱在硫化降温、内压和轴向过载作用下的热粘弹性准静态响应,为实际结构设计提供参考。     

微型固体脉冲推力器结构强度有限元分析

对用于飞行姿态控制和弹道修正的微型固体脉冲推力器的结构强度进行了动静态有限元分析, 考虑了脉冲推力器的结构设计特点和加工工艺, 模拟其台架静压试验和实际工作的载荷历程。数值分析结果表明, 当静态台架试验内压采用动压峰值时, 计算所得的结构应力、变形及其分布与模拟推力器实际工作过程的动态峰值应力、变形及分布具有较好的比拟性; 燃烧室壳壁沿平行其轴线方向的破裂以及燃烧室与点火器壳体螺纹连接处退刀槽根部的开裂是推力器结构强度的主要破坏形式。计算预示与台架试验出现的相符     

用有限元法计算边界移动的喷管温度场

以固体火箭发动机复合结构喷管为对象，用有限元法进行了二维轴对称瞬态热分析。计算中考虑了因烧蚀引起的边界移动．且实现了网格的自适应，给出了考虑与不考虑烧蚀情况下的温度场计算结果，并将计算结果与试验结果做了比较，二者吻合较好，说明该预估方法既具有理论意义又有实用价值。