燃气轮机叶轮非线性瞬态温度场的边界元分析

首先采用克希荷夫变换将包含变物性及混合和辐射边界条件的热扩散方程转换成线性方程，进一步由加权余量法和离散化将方程转变成边界元方程组。用所开发的二维线性—非线性瞬态温度场边界元程序分析计算了燃气轮机叶轮起动工况下的非稳态温度场，与有限元数值计算结果比较表明，边界元法应用于这类非线性工程问题的分析是可靠、初始精度高，数值稳定性好及计算输入数据简单等优点。     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

加权变系数的瞬态导热有限元法

研究了两点时间差分格式的瞬态导热温度场有限元法数值解的稳定性和振荡现象,推导出消除振荡的方法:将时间差分格式的加权系数σ与时间步长△t分别和温度振荡联系起来,得出变σ和变△t的函数关系式。      

喷管喉部沉积层温度场与沉积速率分析

以传热学的观点,喷管喉部发生沉积时的传热分析属于移动边界不稳定导热问题.本文建立了传热分析模型与传热微分方程.通过动座标变换,求得了微分方程的近似解析解,即温度场随时间变化的规律.由此进一步分析了金属氧化物沉积速率随时间变化的规律.     

硅基内衬喷管扩张段烧蚀和温度场的耦合计算方法

本文叙述固体火箭喷管扩张段硅基内衬烧蚀和温度场的耦合计算方法.烧蚀采用非定常的液体层烧蚀模型;温度场计算采用圆柱坐标系下的一维径向瞬态导热模型,考虑材料在高温下形成的炭化层和热解层,而由烧蚀产生的边界移动问题用坐标变换的方法处理.利用烧蚀和导热计算方程组中二者的结合部,通过充分耦合,获得了更准确的预示计算结果.     

固液相变问题的边界元法数值解

相变问题的数值计算方法有多种,本文用边界元法解应用工程中常见的相变导热问题。     

MTPS蜂窝夹芯结构传热性能及热应力分析

对于金属防热结构的蜂窝夹芯结构的稳态情况,基于材料的全灰体假设,同时考虑热传导和热辐射两种传热形式对温度场的耦合作用,利用热流量守恒建立了蜂窝芯层温度场的非线性积分方程。离散化后利用数值方法得到方程组的数值解。对于美国兰利研究中心的实验结果,与本文方法的对比计算结果基本吻合。进一步,利用计算结果讨论了给定面板温度边界情况下,下面板、柱体层的灰度、蜂窝结构长径比对夹芯温度场的影响。并根据温度场和近似的应力分析模型,用半解析结果讨论了稳态情况下蜂窝芯层上的热应力。     

转子轴系耦合振动模型一般性解析式推导与分析

基于转子受到的线性和非线性作用机理，引入线性尺度因子、非线性尺度因子和耦合强度比作为特征参量，分别构建了阻尼力和刚度力的一般表达式，进而推导出转子轴系一般形式的线性与非线性耦合的振动模型.基于仅考虑质量不平衡激励的单圆盘转子轴系，应用多尺度法求解出该方程的稳态解及瞬态解.通过对解析解的分析，揭示了线性与非线性耦合作用及响应的振动机理.通过数值计算分析了瞬态时间尺度因子对响应的影响，瞬态时间尺度因子越大，瞬态解衰减越快，越接近于稳态解.对稳态解短时傅里叶变换的幅频特性进行分析，可发现1倍频出现了非线性特征，3倍频存在双峰特征，进一步阐述了非线性尺度因子对转子轴系耦合振动的影响.      

边界元法在非稳态导热问题中的应用

本文用边界元法分析与计算了航空发动机涡轮工作叶片温度场的瞬态分布。     

飞机整体瞬态热状况的数值仿真研究

 在分析飞机内外热环境的耦合传热机制基础上，建立了描述飞机整体瞬态热状况的物理数学模型。引入壁面热流函数，将飞机蒙皮外的气动对流与辐射换热作用转化为舱内热分析的浮动热边界条件，实现了飞机内外耦合热作用的解耦计算，避免了难以实现的直接耦合求解。采用区域分解技术，将飞机整体分解为特性不同的多个求解区域。采用蒙特卡罗法求解舱内设备之间的辐射换热，采用热网络法求解设备内部及各设备之间的辐射-导热-对流耦合换热。利用该方法对某型飞机飞行过程中的瞬态热状况进行了数值模拟，获得了飞机舱内的瞬态温度场，分析了相关因素的影响，为飞机设计及相关研究提供了重要依据。     

燕尾形榫联接修形的数值研究

将边界元法用于榫联接修形的数值分析研究。对分析具有摩擦接触问题的边界元法的理论及实现方法作了讨论;将边界元分析应用于航空发动机轮盘榫槽与叶片样头联接的计算,并与激光散斑法的实测结果进行比较;对修形的榫联接进行了计算分析,初步揭示了修形榫联接接触边界上力学参数的变化规律,与试验结果相符合。     

二维可压粘流N—S方程的边界元解法

本文给出了一种求解二维可压Navier-Stokes(简称N-S)方程的边界元法。使用线化技术把控制方程变成变系数的线性偏微分方程,其基本解由组合代数的方法构造,从而给出了N-S方程的积分表达式。本文计算了几例低R_e数的跨声速翼型绕流,结果表明边界元法可以应用到可压粘流中。     

用边界元法计算变截面变轴力柱的屈曲

 1.变截面变轴力柱屈曲计算的边界元法 变截面变轴力柱的挠曲方程的一般形式为式中 d²/dx²[EI(x)dw/dx²]+p(x)d²x/dx²=0,E——弹性模量;W——挠度;I（x）——截面处的横截面惯性矩;P（x）——x截面处的轴向压力;x——轴向坐标。 （1）式是非线性常微分方程,本文采用微分方程分解法,将（1）式分解成两个二阶变系数常微分方程。     

燕尾形榫联接抗微动损伤修形研究

为了预防和减轻榫联接中微动磨损,提高其抗微动疲劳能力,主要措施有三:第一,在结构设计中,使联接的方式有利于减轻磨损或尽可能地消除接合面;第二,对接触表面作适当的处理,改善表面质量,如表面涂层、离子注入和表面冷加工等;第三,在接触面上加填充物。本文提出用修形方法改变榫联接接合面的形状,改善接触处的有关力学参数,达到减轻联接的微动磨损和提高抗微动疲劳能力的目的。      

高阶有限元格式研究的一个新途径

 在作者提出的有限差分－有限元混合方法的基础上，将差分法中的三阶无波动、无自由参数的ENN格式较成功地推广应用到有限元法中，为发展高精度有限元格式提供了一个新途径。     

三维铸件定向凝固过程模拟Ⅱ——非稳态定向凝固过程边界元解

定向凝固过程可分为两个阶段即为抽拉前的稳态过程和抽拉后的非稳态过程,此部分用边界元法对三维瞬态过程进行了数值模拟,给出了定向过程中的热源及温度分布特性。     

Theoretical aspects of selecting repeated unit cell model in micromechanical analysis using displacement-based finite element method

Repeated Unit Cell(RUC)is a useful tool in micromechanical analysis of composites using Displacement-based Finite Element(DFE)method,and merely applying Periodic Displacement Boundary Conditions(PDBCs)to RUC is almost a standard practice to conduct such analysis.Two basic questions arising from this practice are whether Periodic Traction Boundary Conditions(PTBCs,also known as traction continuity conditions)are guaranteed and whether the solution is independent of selection of RUCs.This paper presents the theoretical aspects to tackle these questions,which unify the strong form,weak form and DFE method of the micromechanical problem together.Specifically,the solution’s independence of selection of RUCs is dealt with on the strong form side,PTBCs are derived from the weak form as natural boundary conditions,and the validity of merely applying PDBCs in micromechanical Finite Element(FE)analysis is proved by referring to its intrinsic connection to the strong form and weak form.Key points in the theoretical aspects are demonstrated by illustrative examples,and the merits of setting micromechanical FE analysis under the background of a clear theoretical framework are highlighted in the efficient selection of RUCs for Uni Directional(UD)fiber-reinforced composites.     

用边界元法研究非保守力系下杆的稳定性

 本文采用边界元法,研究了完全弹性支持杆,在均布的和集中的跟随压力作用下的稳定性,给出了结构失稳的临界载荷和两种失稳形式（颤振失稳和发散失稳）的转变随边界条件变化的一般规律,以及加载方式（载荷分布情况）对临界载荷和失稳形式的影响,说明了边界元法对解决非保守稳定问题的有效性。     

随机有限元分析定向结晶涡轮叶片在随机压力场中的可靠性

本文采用随机有限元法和二阶矩Hasofer-Lind方法,分析随机压力场对于定向结晶叶片可靠性的影响。利用有解析表达式的、做为随机量采用Monte-Carlo方法计算可靠性指标,与本文编制的随机有限元程序计算的可靠性指标对比,以验证后者的正确性。      

Neumann 随机有限元法计算涡轮叶片的可靠度

推导了具有随机性时Ｎｅｕｍａｎｎ随机有限元法的基本方程,用级数来拟合结构响应的概率密度函数,并通过数值积分计算结构的可靠度。对涡轮叶片的可靠度分析表明,该方法可利用较小的样本容量获得较高的计算精度,并具有较高的计算效率。