四边简支压电热弹性体层合开口壳的精确解

根据压电热弹性体的本构关系和热传导关系的相似性,通过增维方式将热传导关系写进压电热弹性体的本构关系中,直接得到了一个新的增维的压电热弹性体本构关系。依据此增维的本构关系,结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体层合开口壳机、电、热耦合问题的齐次状态方程。这类方法大大简化了人们在分析压电热弹性体耦合问题时通常要求解的非齐次微分方程和关于温度的平衡方程和导热方程的二阶微分方程的繁琐方法,同时,齐次状态方程的提出大大减少了数值计算的工作量。     

热弹性复合材料层合板的响应分析

利用混合变分原理,考虑了热平衡方程与导热方程中变量的对偶关系,通过增加非齐次正则方程的维数,成功地将热弹性体机一热耦合问题的非齐次Hamilton正则方程转化为能独立求解的齐次Hamilton正则方程．所以大大简化了问题的求解。并通过数值实例研究了在温度栽荷和动力栽荷作用下及考虑阻尼的四边简支复合材料层合板的响应问题。     

伸展柔性附件对简谐激励的响应分析

从Hamilton原理出发推导了柔性附件伸展过程的时变系数动力学方程。对附件在简谐激励下的响应进行了数值仿真,讨论了伸展柔性附件所具有的一些独特的现象：齐次解的非瞬态特性、特解的非稳态特性、瞬时共振的延续和加强。由于这些特性的存在,大大增加了由激励引发剧烈振动而导致结构破坏的可能性。     

潜入式喷管喉衬界面间隙优化设计

为了解决固体火箭发动机潜入式喷管喉衬在热载荷与内压联合作用下界面间隙优化问题,数值模拟方法求解了喷管热结构响应.采用流体仿真软件,计算了潜入式喷管流动过程的稳态流场,获取了高温燃气流动参数与固体壁面对流换热系数.并采用三维有限元结构计算平台,编写了非均布壁面压力载荷与非均布对流换热系数子程序,求解了喉衬前后搭接界面间隙...     

球锥组合体在热外压作用下的稳定性

球锥组合体作为飞行器头部,热外压稳定性是结构设计的主要问题。本文对薄壳结构进行稳定性分析,将临界外压的有限元计算值与试验修正的理论值进行对比,获得有限元分析的修正系数为0.8;然后采用数值和试验相结合的方法、并考虑材料力学性能的温度效应,对复合材料球锥组合体在热外压载荷下的稳定性进行研究,高温下材料性能退化对临界载荷的影响显著,加热不均匀对临界载荷也有一定影响,线胀系数对临界载荷的影响很小。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

含固支边的强厚度层合板解析解的改进方法

基于含固支边层合板的非齐次状态方程,应用增维方法,建立了齐次状态方程,给出了静力问题的解析解。这种方法对于程序的实现和数值运算稳定性的提高十分有利,在求解的过程中避免了矩阵求逆,且大大提高了计算效率。数值算例显示了本文的方法是正确的。     

热声载荷作用下薄壁结构非线性响应分析和试验验证

针对热声载荷作用下薄壁结构大挠度非线性响应问题,开展了固支约束金属薄壁板结构热声激励试验及数值模拟分析。通过计算结果与试验结果对比,表明两者结果存在一致性,进而验证了薄壁板在热声载荷作用下动态响应计算方法和数值模型的有效性。在此基础上,针对加筋板结构完成了多种热声载荷组合作用下的动力学响应计算,获得了时域位移响应。重点对该结构在后屈曲状态下的3种典型振动形式进行分析,总结出热载荷与声载荷之间的相对强度决定了板的不同跳变形式。采用统计分析方法建立了位移响应的概率谱密度函数（PDF）并绘图,清楚地显示了后屈曲板的PDF表现出双峰现象。使用功率谱密度（PSD）函数分析了响应频率和峰值随着温度升高的变化,并确定了板的软化和硬化区域。总结了屈曲前/后结构特定区域拉应力和压应力的变化规律,并阐述了造成这种变化的原因。本文工作可对热声载荷作用下薄壁结构响应分析和动强度设计提供参考依据。     

基于PHengLEI的非稳态电热除冰过程仿真

为了填补国产自主可控CFD软件风雷平台的防除冰功能开发的空白,本文建立了电热除冰计算模型、计算方法和非稳态导热模型,并在国家数值风洞风雷平台（PHengLEI）基础上,集成了非稳态电热除冰计算和非稳态导热计算功能。通过与主流商业CFD软件仿真结果和实验数据的对比,验证了导热和非稳态电热除冰程序的准确性。针对某飞行工况进行了电热除冰计算,通过对表面溢流水、表面温度、结冰量的计算结果分析,发现合理布局加热片、设计加热热流密度和电热除冰控制率,可实现电热除冰系统的安全运行和能源的高效利用。     

加筋板在热声载荷作用下随机疲劳寿命估算

薄壁结构在热声载荷作用下的应力响应对其疲劳寿命估算有着很大的影响.给出了热声载荷作用下加筋板结构的大挠度运动方程,运用有限元方法进行数值模拟,分别计算了四边简支薄壁板与加筋板在相同有限带宽高斯白噪声载荷和变化温度下的应力响应,得到了应力时间历程、应力概率密度分布和应力功率谱密度,利用改进的雨流计数法,结合线性累计损伤理论,采用Smith-Watson-Topper（SWT）平均应力模型估算了两种结构在热声载荷作用下的疲劳寿命,并进行了比较分析.     

八节点Hamiltonian等参元列式

为了使平面八节点等参元的优越性在弹性力学Hamilton正则方程的半解析法得到应用。结合弹性材料修正后的Hellinger—Reissner（H—R）变分原理和二次插值函数表达平面外应力和位移函数,建立了Hamilton正则方程的八节点等参元列式。首先简要地介绍了弹性材料修正后的H—R变分原理．然后用二次插值函数表示平面外应力和位移变量,并详细地推导了Hamilton正则方程的八节点等参元列式。数值实例结果证明了本文等参元列式的正确性。     

Hamilton正则方程半解析法的收敛和对称分析

近年来,Hamilton正则方程半解析法在工程问题上的应用越来越广泛,但至今未见有关这种方法收敛性和对称性问题研究的文献。基于Hamilton正则方程的半解析法理论,通过变分原理详细推导了Hamiltonian元素的固支和简支边界公式及对称边界公式。多个实例的数值研究表明：随着网格加密,Hamilton正则方程半解析法的收敛速度快于一般传统位移有限元法,对称解法的效率明显优于整体解法。     

热空腔-喷气流的组合辐射

讨论了热空腔(喷气发动机的加力筒)和喷气流的组合辐射问题, 推演了组合辐射的辐射传输微分方程, 导出了有边界辐射条件下的解。讨论了由多个温度不同的壁面组成且腔内充满高温燃气的热空腔辐射传输方程解的具体形式, 导出了它向后半球的光谱辐射强度的表达式;同时还推导了热空腔-喷气流组合辐射的光谱辐射强度的表达式。最后给出算例, 算例结果表明, 热空腔-喷气流组合辐射的空间分布具有与实测分布较一致的梨形包络, 其光谱分布也与实测的光谱分布的一致性较好。      

总压形式的能量方程与热阻问题

本文根据热力学第二定律的熵方程和稳定流的能量方程导出了用总压表示的一维定常流的能量方程。这是不同于稳定流能量方程和伯努利方程的另一种形式的能量方程。     

A hybrid lattice Boltzmann flux solver for integrated hypersonic fluid-thermal-structural analysis

In this paper, a hybrid Lattice Boltzmann Flux Solver (LBFS) with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems. In the solver, the macroscopic Navier–Stokes equations and structural heat transfer equation are discretized by the finite volume method, and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation. To compute the numerical fluxes, two equilibrium distribution functions are introduced. One is the D1Q4 discrete velocity model for calculating the inviscid flux across the cell interface of Navier–Stokes equations, and the other is the D2Q4 model for evaluating the flux of structural energy equation. In this work, a new dual thermal resistance model is proposed to calculate the thermal properties at the fluid–solid interface. The accuracy and stability of the present hybrid solver are validated by simulating several numerical examples, including the fluid-thermal-structural problem of cylindrical leading edge. Numerical results show that the present solver can accurately predict the thermal properties of hypersonic fluid-thermal-structural problems and has the great potential for solving fluid-thermal-structural problems of long-endurance high-speed vehicles.     

热平衡等离子体的二维复合流动数值计算

采用中性气体流动欧拉方程的近似因子分解法，耦合等离子体的热平衡模型方程，对用微波产生的等离子体在喷管内二维轴对称复合流动进行了数值计算。计算结果表明：微波等离子体电子数密度不超出临界值，1kW以下功率微波产生的是弱电离等离子体；气流的马赫数分布不受微波加热功率、喷管入口压强的影响，喷管内等离子体中电子数密度随微波功率的加大而增加，随入口压强的增加而减小。     

尖端具有非线性吸附接触的界面裂纹的动态特性研究

 研究一类尖端具有非线性吸附接触型裂纹模型的弹性波散射问题。利用积分变换和积分方程方法推导了确定这类问题基本变量的奇异积分方程组,采用围道积分技术和切比雪夫多项式展开技术得到了待定系数的非线性代数方程组。最后给出了裂纹尖端接触区的大小和接触应力的数值结果,揭示了这种接触型裂纹模型的动力学特征及物理上的合理性。