旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

胶接修补复合材料层合板失效分析的PDA-CZM方法

建立了一种预测胶接修补复合材料层合板的损伤演变与剩余强度的PDA-CZM方法。该方法应用三维渐进损伤分析（PDA）方法和粘聚区模型（CZM）分别模拟复合材料层合板和修补胶层的失效过程。对修补层合板的纤维断裂、基体开裂、层间分层以及胶层脱胶等损伤的萌生和扩展以及它们之间的耦合作用进行了研究。计算得到了修补结构的载荷--位移曲线,并预测了其极限强度。计算结果和试验数据吻合良好,验证了PDA-CZM方法的有效性。最后,讨论了修补参数对剩余强度的影响规律。     

航空发动机旋转帽罩结冰表面换热系数研究

旋转帽罩表面的对流换热系数是结冰、防冰研究的重要参数,旋转帽罩结冰模拟中对流换热量和蒸发量都需要已知对流换热系数。为了获取旋转帽罩表面的对流换热系数,采用数值模拟对旋转帽罩结冰表面的对流换热系数进行了研究。首先验证了网格与计算方法的合理性和可靠性,在此基础上,对影响表面换热的因素进行了分析,计算了不同转速、来流速度、来流温度对锥角40°,锥高176mm的旋转帽罩表面换热系数的影响。结果表明:锥尖区域,来流速度的影响与转速影响相比占主导;除锥尖外的其他区域,转速对换热系数的影响占据了主导地位;来流温度对整个表面的换热系数均有影响。采用不同工况下的数值模拟结果,在锥尖区域建立了来流雷诺数与努赛尔数之间的关联式;在除锥尖外的其他区域建立了旋转雷诺数与努赛尔数之间的关联式。     

发动机旋转整流帽罩的水滴撞击特性分析

基于计算流体力学软件FLUENT,采用拉格朗日法,提出了一种求取旋转部件表面水滴收集系数的方法,并采用用户自定义函（UDF）予以实现。将该方法运用于发动机旋转整流帽罩的水滴撞击特性研究,结果表明：随着帽罩锥体锥角及转速的增大,帽罩表面的局部水收集系数沿其表面外形线下降较快,且撞击区域变小。     

旋转帽罩结冰相似准则的研究

为了研究航空发动机旋转帽罩结冰相似理论,在静止部件结冰相似准则的基础上,忽略次要因素的影响,建立了一种适用于旋转帽罩的结冰试验相似准则。利用建立的旋转结冰相似准则,通过数值模拟方法,对旋转帽罩缩比前后的水滴收集特性和结冰冰形进行了误差分析,并对此准则的可靠性进行评估。计算结果表明:用建立的相似准则所确定的参数对于全尺寸旋转帽罩转速分别为1kr/min和5kr/min时进行数值模拟。在全尺寸旋转帽罩及其1/2缩比模型上的对应区域得到的水滴收集特性和结冰冰形,误差约为2.7%和5.3%,表明选取缩比参数的方法是可靠的,并且提出的相似准则是有效的,对结冰风洞试验的理论研究和参数选取有一定的参考价值。     

基于探针试验的预浸料黏性内聚力模型

以预浸料黏性为研究对象,通过探针试验,研究了预浸料铺放过程中预浸料与模具之间的脱粘失效行为及其随铺放工艺参数的变化规律。试验中发现铺放速率、压力和温度对预浸料黏性行为影响显著,且存在界面失效（低温）和内聚失效（高温）两种失效模式。针对这一现象,采用指数内聚力模型（CZM）表征了预浸料的黏性行为,利用该模型描述了不同铺放温度下粘结层应力快速上升、损伤扩展和最终失效的脱粘过程,定量地给出了黏性CZM参数与铺放工艺参数之间的关系。研究结果表明,预浸料粘结强度和特征位移均随铺放速率增加而近似线性减小,均随铺放压力增加而近似线性增大。随铺放温度的增加,黏性CZM参数先增大后减小,近似呈二次关系。为制定特定铺放条件下的铺放工艺规划提供参考。     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

旋转整流罩积冰过程中水膜流动和脱离的数值研究

为了更准确地模拟旋转整流罩积冰过程,在现有三维积冰与冰层表面薄水膜流动耦合模型基础上,基于功平衡分析的方法引入了旋转部件表面水膜脱离模型,并发展了相应的计算方法,给出了水膜脱离的判定依据:当气流曳力做的功和由于离心力使水膜增加的潜能之和大于黏附功时整流罩表面的水膜会发生脱离。对旋转整流罩积冰进行数值模拟,计算结果与实验结果吻合得较好,验证了该模型的合理性和计算方法的可行性。之后分析了转速和来流速度对整流罩表面水膜脱离和积冰的影响,结果表明:转速和来流速度越大,水膜发生脱离的比例越大。在研究范围内,转速为3000r/min和6000r/min时,因水膜脱离导致积冰总量分别减少13.4%和15.8%;来流速度为40、50m/s和60m/s时,因水膜脱离导致积冰总量分别减少为12.2%、13.4%、14.2%。      

离心力对于旋转叶片表面积冰影响的数值模拟研究

为了加深对旋转叶片表面积冰现象的认识,发展了一种适用于旋转表面积冰预测的数学模型,并采用该模型研究了离心力对于旋转叶片表面积冰的影响。通过对积冰控制体内的由于撞击水滴形成的薄水膜流动的质量、动量和能量守恒进行分析,发展了旋转贴体非正交曲线坐标系下的积冰模型并给出了相应的计算方法。所发展模型的动量方程中考虑了离心力对水膜流动的影响,并且在能量方程中考虑了净流出水带走的能量。采用文献中旋翼表面的积冰试验结果对所发展的模型和计算方法进行了验证,所计算得到的冰厚与试验结果吻合较好,计算得到的冰厚精度相比于LWICE软件有所提高。计算结果还表明,离心力增大会导致水膜速度在叶片展向上的分量变大,致使流出控制体的水质量增加,最终导致冰厚在驻点附近略微变小。     

航空弧齿锥齿轮风阻损失机理与影响因素分析

基于计算流体动力学（CFD）理论,采用剪切应力输运（SST） k-ω湍流模型及多参考系旋转模型（MRF）技术和Coupled解耦算法,仿真分析了单个弧齿锥齿轮周围的气流特性及齿轮风阻损失机理,研究了齿轮转速、旋转方向及挡风罩配置对风阻功率损失的影响规律。然后利用正交试验分析方法进行数值模拟仿真分析,研究挡风罩与齿轮齿面、大端及小端之间不同间隙值时的风阻力矩,对数值仿真数据进行方差分析与拟合得到齿轮风阻功率损失与间隙值的函数关系,结果表明:挡风罩间隙值均为1 mm时,风阻力矩最小,降低风阻功率损失的效果最好。利用最优化分析方法给出了挡风罩与齿轮间隙值的单目标最优参考值,为进一步研究挡风罩间隙值多约束下的多目标优化提供了理论基础及研究方法,也为挡风罩的工程应用提供了参考。      

含缺陷复合材料T型接头失效数值分析

 针对复合材料T型接头优化设计问题,利用粘聚区模型(CZM)建立了可预测复合材料接头拉伸破坏过程的分析方法.将填充物缺陷加入该模型,填充区裂纹的起始是随机的,重点分析填充物缺陷的位置、尺寸大小以及填充物刚度对接头承载能力的影响.计算结果表明:模拟结果与实验结果吻合较好;无论缺陷位置如何,缺陷尺寸越大,结构承载能力越弱;结构承载能力对上角缺陷最为敏感;填充物刚度在3~50 MPa时,填充物刚度较大的接头损伤容限特性更好;不同的缺陷位置,破坏模式差异较大,但起裂点或最终失效处都在缺陷附近的基体中.     

电热除冰的热力耦合特性及其对冰层的影响研究

采用计算和实验相结合的方法,初步研究了电热除冰过程中的热力耦合特性及其对冰层的影响。在电加热条件下,耦合外部气动力载荷的作用,采用有限元方法计算了不同热流密度下表面冰层和蒙皮间界面法向和切向应力的分布,比较了加热/不加热条件下界面法向和切向应力分布的差别,研究了冰层最大主应力随热流密度的变化规律。研究发现,电加热条件下,在表面冰层融化前,热力耦合特性将造成冰层内部应力的显著增加,从而造成冰层局部区域的破裂,加速冰层的破坏。同时,设计的原理性实验结果验证了热力耦合特性对冰层的破坏影响。研究结果对于电热除冰理论和除冰技术的发展有现实意义。     

三元乙丙橡胶薄膜黏接界面温度相关性力学性能

采用双悬臂夹层梁(DCSB)试件对三元乙丙橡胶(EPDM)薄膜黏接界面Ⅰ型断裂的温度相关性进行了研究;采用线弹性断裂力学(LEFM)方法获取了黏接界面断裂能,将其与单轴拉伸所得内聚强度作为双线性内聚力模型(CZM)参数,对不同温度条件下双悬臂夹层梁试件的Ⅰ型断裂行为进行了数值仿真.结果表明:仿真与实验曲线偏差较大,偏差出现的原因主要是未经修正的线弹性断裂力学方法所求断裂能存在较大的误差,需对其进行修正.调整模型参数使仿真曲线与实验曲线重合,获取了黏接界面的准确力学性能参数,采用此参数得到的仿真结果与实验结果具有较高的吻合度,证明了模型的可适用性.     

锥角对旋转整流罩积冰影响的模拟实验

采用缩比模型对4种不同锥角的旋转整流罩进行了冰风洞积冰模拟实验.推导出旋转表面积冰相似参数,并根据重要参数的匹配确定缩比整流罩模型的积冰实验参数.在实验中采用高速摄像系统记录冰生长过程及最终冰形.结果表明:4种锥角的整流罩表面积冰均由初期生成、分布连续的明冰和后期增长迅速的白色霜冰构成;锥角小于等于74°的整流罩表面霜冰为针状或粒状,积冰厚度较小;锥角大于80°的整流罩表面霜冰为羽毛状,积冰厚度较大,并伴随冰脱落现象,锥角较大的整流罩冰脱落位置向下游移动.      

双面贴补层合板压缩屈曲渐进损伤分析

贴补复合材料层合板在压缩载荷作用下的屈曲破坏强度及其损伤演化过程对于复合材料结构修理具有重要意义。本文基于应变和黏聚区模型(CZM)建立了贴补复合材料层合板的渐进损伤分析模型,引入复合材料与胶层的损伤判据和刚度退化方案,计算了结构屈曲强度。数值仿真结果和实验数据吻合较好,验证了模型的有效性。基于该模型,采用非线性有限元方法研究了压缩载荷下双面贴补复合材料层合板的屈曲损伤演化过程,并讨论了补片参数对结构屈曲强度的影响。研究结果表明:双面贴补复合材料层合板屈曲后,处于拉伸和压缩状态下的铺层中的损伤程度存在差异;增大补片直径与厚度可以在一定程度上提高双面贴补复合材料层合板的屈曲强度。     

基于CTOAc评估结构断裂特性的方法研究

介绍了一种基于临界裂纹尖端张开角(CTOAc)评估破损结构延性断裂特性的方法。引入了评估裂纹稳态扩展和止裂的CTOA准则,介绍了CTOAc的试验测定方法,重点推荐准静态的方法。通过对两类典型M(T)试件进行有限元分析,应用应力失效准则及节点力释放技术模拟了结构的裂纹扩展过程,得到了对应裂纹长度的CTOA值。结果表明,CTOA是模拟延性断裂的一个重要表征参数,CTOA准则在高韧性的航空铝合金结构中具有广泛的适用性和广阔的应用前景。     

四自由度飞机结冰脱落模拟研究

飞机表面的冰粒脱落可能会对安装在尾部的发动机和其他飞机部件造成破坏,直接影响到飞机的安全飞行。在对实际结冰脱落过程进行合理简化的基础上,建立了四自由度结冰脱落轨迹模型,应用3阶Runge-Kutta方法求解冰粒运动方程,并嵌入了气动阻尼对冰粒轨迹的影响;同时将计算得到的冰粒轨迹及转角与文献中的计算和试验结果进行了对比验证;最后模拟并分析了不同形状、尺寸、初始俯仰角及气动阻尼对于冰粒运动轨迹的影响,为飞机表面结冰脱落的进一步研究工作提供了参考。     

双轴应力下含裂纹加筋曲板剩余强度的弹塑性有限元计算

 采用考虑了大应变的弹塑性有限元方法,以裂纹尖端张开角为断裂准则,实现了双向受载条件下加筋曲板中裂纹起裂到失稳整个扩展过程的数值模拟计算;还对与加筋曲板有同样尺寸的加筋平板进行了分析和计算。结果表明,所提出的分析方法和所采用的断裂准则适用于对含裂纹加筋结构所进行的剩余强度分析。     

基于CDM-CZM的复合材料补片补强参数分析

 复合材料开口补强设计参数的确定对于结构设计具有重要的意义。针对复合材料层合板开口区补片补强结构,采用各向异性材料连续介质损伤力学模型(CDM)对复合材料层合板的损伤演化进行描述,采用粘聚区模型(CZM)对补片与母板间界面材料的分层损伤进行模拟,建立了复合材料开口区补片补强结构三维非线性渐进损伤模型,模型可预测补强结构强度和损伤演化过程。应用本文模型分析了补片铺层方式、补片厚度和补片半径3个主要设计参数对补强效果的影响,明确了补片与母板间界面材料分层损伤破坏是导致补强结构最终失效的主要原因。