反舰导弹三维流场数值模拟

三维流场的数值模拟是国外军事技术发达国家在飞行器设计过程中所采用的重要方法.本文建立了导弹的真实三维几何外形,然后结合实验数据,采用有限体积法,对国产某典型反舰导弹在高亚音速、不同攻角时的三维流场进行了数值模拟.计算出了导弹在不同工况下的阻力系数、升力系数和力矩系数.把计算所得结果与实验值进行了比较,用理论公式对计算结果进行了对比分析,并简单地分析了偏差的原因和解决的可行措施.文中所用的残差光顺技术和网格自适应技术能有效地提高计算的收敛速度和精度.     

CISA 3D Navier—Stokes程序在压气机中的鉴别和应用

用压气机叶栅试验数据、多级压气机的第4级静子计算结果和跨声速转子试验结果对CISA 3D N-S程序的计算精度进行了鉴别研究。结果表明,该程序对设计点的落后角的计算精度较高,而对损失的计算还存在一定的分散度。采用 CISA 3D N-S程序和S_2-S_1程序对一台多级压气机试验件的第3级静子进行三维造型修改设计,取得了预期的效果。     

有冷却的涡轮叶片三维和准三维温度场计算

通过对二维、准三维和三维程序的开发和运用,对叶片的准三维、三维和二维温度场进行了计算,通过对各种计算结果的比较评价,为气冷叶片的传热设计提供了更有价值的数据。     

普惠E~3高压涡轮级气动性能验算

利用自主开发的涡轮气动设计体系对普惠公司(P%W)E3高压涡轮叶型进行反设计.并在考虑气膜冷却的前提下对反设计得到的涡轮级进行全三维气动计算验证.计算结果与试验数据对比表明.三维计算结果与试验数据吻合较好,通过叶型反设计和计算验证掌握了该涡轮的设计特点.并证明了该涡轮气动设计体系具有较高的可信度和精度,可用于大涵道比涡扇发动机涡轮部件的设计.     

发动机进气道短舱前缘结冰三维模拟研究

为计算发动机进气道短舱前缘的结冰冰形,在Messinger结冰热力学模型的基础上发展了一套纯三维的表面溢流水流动结冰模型。考虑空气摩擦力为表面溢流水流动的主要驱动力,用空气对进气道表面的剪切力来确定溢流水的流动方向及流量分配。为求解溢流水结冰模型,发展了一套重复查找表面控制体状态的方法,能够快速完成整个三维表面的计算。用该方法对某三维发动机进气道进行计算,得到了三维结冰冰形,并将计算结果与FENSAP-ICE计算结果进行了对比,结果显示两者的冰形轮廓基本一致,仅在冰角处存在差异,表明本文三维发动机结冰计算模型与计算方法是有效的,其计算精度与FENSAP-ICE结果相当。     

有侧壁干扰的机翼半模型风洞实验的N-S方程数值模拟

应用三维可压、雷诺平均Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。计算采用中心有限体积多步Runge-Kutta时间步长格式,将Baldwin-Lomax两层代数紊流模型扩展应用于三维拐角区流动。用本文方法计算了3个算例 ,并与国外的有关实验数据进行比较,计算和实验吻合良好。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

基于凸包算法的陶芯弯曲度与扭曲度误差计算

陶芯弯扭变形直接关系到空心叶片的壁厚尺寸分布,为克服当前陶芯弯扭变形的计算中测量数据与理论模型三维配准、陶芯截面轮廓线提取或拟合的过程算法复杂、收敛速度慢、效率低等问题,提出了一种通过测量数据点直接计算陶芯弯曲度和扭曲度误差的算法,该算法不需要三维配准和提取陶芯外轮廓线,通过距离权值法计算陶芯弯曲度,凸包算法计算陶芯扭曲度,能大幅提高计算效率。仿真与实验结果表明:该算法弯曲变形计算精度为99.55%,与二维配准算法相差±0.01mm;扭曲变形计算精度为99.98%,与二维配准相差±0.006°。      

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。     

真实火箭燃气流场的测量与数值仿真

对用激光－莫尔偏折法测出的真实火箭燃气射流流场进行了数值仿真,在计算中采用了二阶精度的ＭＵＳＣＬ格式,将三维问题转化为二维问题,结果表明,理论计算与实验测量结果符合较好。     

复合材料层合板弯曲振动高阶理论及有限元模型

本文提出了一种新的复合材料层合板高阶弯曲位移模式,此位移模式既考虑横向剪切的影响,又考虑了法向正应力的分布,进而建立的相当于三维弹性理论解的二维有限元模型,具有精度高、计算量并不大的优点。     

一种含分层复合材料层合板高阶理论及有限元模型

本文提出了一种新的复合材料县合板高阶弯曲位移模式。此模式考虑了横向剪切的影响,得到一种新的厚板有限元模型。对于含有分层的单元,通过改变其单元内横向剪切应力的分布,达到改变单元刚度的目的。因此,此模型适用于任何厚度的或含有分层的复合材料层合板。     

板宽对胶焊接头应力分布影响的数值分析

采用三维弹塑性有限元方法,研究了酚醛树脂和丙烯酸树脂两种胶粘剂的胶焊接头中,板宽对应力分布的影响。结果表明,酚醛树脂胶粘剂胶焊接头中,搭接区边缘应力值高于焊点边缘应力值,导致接头从搭接区边缘起裂,板宽增大时,焊点边缘处和搭接区边缘的峰应力都减小,对于提高胶焊接头的强度有利,但不会改变接头的断裂形式;丙烯酸树脂胶粘剂的胶焊接头,焊点边缘应力远高于搭接区边缘应力,因此接头从焊点边缘起裂,板宽增大时,焊点边缘应力基本不变。因此增加板宽对提高接头强度无明显作用,接头断裂形式也不发生改变。     

环形过渡元素及其在轮盘振动分析中的应用

 本文发展了联结环形等参数元素和环形厚板元素的一类环形过渡元素,它对于轴对称体的结构动力分析具有重要意义。环形元素采用了半解析的方法,把三维问题降为二维问题来处理,大大节省了计算机存贮量及计算时间。文中对某实际涡轮盘的振动特性进行了计算,其结果是令人满意的。     

复杂受载承力机匣的强度及疲劳寿命

本文首先论述了航空发动机承力机匣应力分布计算与试验,指出这类机匣的有限元计算,一般可分为:模型估算、真实机匣应力分布计算和局邮结构的详细计算。其次介绍了圆筒机匣热应力理论与试验研究的基本方法与主要成果。最后从机匣低循环疲劳寿命预测的初步探索的角度,说明了缺口构件在常温与高温条件下的低循环疲劳寿命预测工作的一些进展。     

冷挤压孔板残余应力场分析

 <正> 1.冷挤压孔板残余应力场 如图1所示,设理想弹塑性材料厚壁圆筒内半径为R_1,外半径为R_2,受内压P作用,卸载后,形成塑性半径ρ,此时,圆筒残余应力场为     

循环加载下的热弹塑蠕变有限元分析

本文考虑在循环加载条件下, 采用混合强化模型, 导出相应的热弹塑性蠕变增量本构方程及有限元格式。研制出平面应力程序, 用牛顿迭代法求弹塑性权, 并提出精度高的屈服吸面校正方法, 程序经考核, 可供工程设计使用。      

喷丸强化残余应力场三维数值仿真分析

建立了喷丸强化三维有限元模型,实现了喷丸强化处理残余应力场的数值仿真。根据显式动力分析各种能量变化过程,确定了显式分析求解时间选择方法;研究了喷丸速度、弹丸直径、入射角度及摩擦力等因素对喷丸残余应力场影响的一般规律,并与现有相关文献试验结论进行对比;分析了单次和多次冲击下材料内部的残余应力场分布及塑性应变分布特点;建立了不同数目的多丸粒喷丸强化三维有限元模型,研究不同覆盖率对残余应力场分布的影响规律。数值分析结果表明,喷丸速度、弹丸大小、入射角度对残余应力场分布有显著影响。冲击次数对残余应力幅值影响小,但等效塑性应变增加明显。喷丸覆盖率的增加会显著改变残余应力场,并使残余应力场分布更加均匀,但残余压应力幅值与喷丸覆盖率没有必然相关性。     

用边界元-数学规划法解弹塑性接触问题

 本文基于塑性增量理论导出一种不显含应变的矩阵增量物理关系及平面应变问题σ_z的表达式,提出解决无摩擦弹塑性接触问题的边界元-数学规划法。据此编写了用边界元法解二维问题的通用程序,并用之解算了两圆柱体的弹性和弹塑性接触问题,将数值解与现有结果进行了比较。     

Postbuckling elastoplastic state analysis of three-dimensional bodies taking into account finite strains

A technique for analyzing the stress strain state of elastoplastic bodies taking into account large displacements, rotations, finite strains, and buckling is presented. The step-by-step loading method with formulation of the resolvent variation equation in the current configuration is used. The governing equations for elastoplastic bodies that connect the Jaumann derivative of the Cauchy-Euler stress tensor with the velocity strain are applied. The spatial discretization is based on the finite element method (FEM) in the framework of the polylinear three-dimensional isoparameteric approximation. The numerical solutions of the problems are given.