三维低速NS方程的并行计算

本文在求解低速NS方程的串行程序基础上,针对分布式存贮网络异构的并行编程环境建立了隐式并行计算软件,并行计算的方法采用区域分解法.并行软件对高速列车进行了气动特性计算,并对并行效率和并行可靠性进行了研究,结果表明并行计算结果与串行计算吻合良好,并行改造后的软件可以应用于工程实际问题.     

双三角翼大迎角绕流的并行计算

本文在ＬＵ－ＳＧＳ方法求解可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的串行程序基础上,发展了适合ＭＰＰ系统类型的并行计算机的并行算法程序,并针对双三角翼大迎角绕流问题,就不同的并行计算方法进行了对比研究,探讨了优化并行策略,提高计算效率和计算加速比等问题。     

一种三维并行自适应计算方法

基于串行网格划分软件METIS与并行化消息传递编程接口（ MPICH2）对现有串行自适应程序进行简单的并行化改造,给出了一种三维可压缩无粘流数值模拟的并行自适应方法。首先利用单个进程调用METIS,串行划分网格;然后对所有进程并行计算以获得初始网格下的流场解;再次利用单个进程对整个流场运用自适应方法进行局部网格加密并调用METIS串行划分网格;最后全部进程在流场初始解的基础上继续并行计算,以获得自适应网格下的流场解。数值模拟算例验证了此方法的可靠性与高效性。     

飞行器亚跨超声速流气动力并行计算研究

隐格式并行求解是计算流体力学并行算法研究的难点,本文针对分布式存储并行处理系统,在异构编程PVM环境下,采用"局部隐"方法,使用发送/接收数据结构和自动剖分等技术,有效地实现了对飞行器亚跨超声速流气动力数值模拟串行程序的并行化改造,保证了并行计算结果与串行结果的一致性.该方法具有简单实用、并行效率高以及并行程序可移植性好和应用前景广泛等特点.     

有限元的OpenMp并行计算技术

传统有限元串行计算技术耗时长、效率低,已远远不能满足工程实践需要。通过分析串行计算的耗时分布、程序可并行性及数据相关性,在保持串行程序不变的基础上,采用基于多线程的OpenMP技术实现单元刚度矩阵形成总体刚度矩阵的程序并行化。针对有限元计算机翼实例进行了不同线程数目下并行加速比的研究,结果表明,OpenMP技术可有效提高有限元计算效率,且存在一个最佳的线程数,使得并行加速比达到最大,此外加速比还与问题的计算量等因素有关。     

真空羽流场的DSMC并行数值模拟

对DSMC方法进行小喷管真空羽流数值模拟的计算软件进行了并行化改造。计算表明在多处理器(节点)上实现并行计算能够大幅度地提高运算速度,缩短运行时间,增加计算规模,计算结果与串行计算一致。并对目前串行计算不能完成的真空羽流计算的算例成功地进行了DSMC并行计算,计算结果与理论分析一致。      

并行多重网格算法求解跨声速Euler方程

在工作站网络分布式存贮ＰＶＭ环境下组织了有效求解三维跨声速Ｅｕｌｅｒ方程多重网格算法的并行计算，讨论了如何将串行算法转变为并行算法，并有效提高并行效率的一些基本方法，最后分析了影响并行效率的一些瓶颈问题。     

喷管内流场并行计算方法研究

建立了基于区域分解的并行计算方法,将已有的后台阶喷管内流场三维ＴＶＤＭａｃ－Ｃｏｒｍａｃｋ显式格式的串行计算程序改造成为并行程序,获得了较高的并行加速比。并行计算结果同串行结果进行了对比,证明计算是成功的。     

复合材料螺栓连接结构参数化建模与并行计算

使用MSC.Marc对复合材料螺栓连接结构进行三维接触有限元分析,使用MSC.Patran的二次开发语言PCL开发了螺栓结构的界面化分析插件,实现了建模的参数化、自动化;提出了考虑接触的对螺钉连接结构进行求解域分解的新方法,大大提高了并行计算效率,减少分析时间。计算结果与之前研究吻合较好。参数化建模和并行计算的结合,为大规模连接结构的三维接触有限元分析乃至优化设计提供了可行的方法。     

战术导弹绕流流场并行计算

本文采用求解NS方程的方法对战术导弹的流场进行了数值模拟.应用区域分割技术和拼接网格的并行策略,发展了一种适合于分布式存贮并行环境的隐式有限体积并行算法,在PVM并行环境下,对战术导弹绕流流场实现了单块网格多机并行计算和多块网格多机并行计算,通过负载平衡等方法得到了较高的加速比.计算结果表明并行计算结果与串行计算结果完全一致.     

分布存储结构并行处理机上隐式有限差分算法并行尝试

通过对传统的Beam-Warming隐式格式的并行化尝试,实现了无粘和粘性空气动力学方程的并行求解。计算中采用了SPMD并行编程模式,并将静态非重叠网格区域分裂技术用于建立计算区域和并行处理系统处理器节点之间的映射。结果表明原等价串行算法的稳定性和收敛性被基本保留,获得的加速比和并行效率较为满意。     

一种计及静气动弹性变形影响的跨声速机翼气动优化设计方法研究

耦合流场控制方程和结构静平衡方程求解得到大展弦比跨声速弹性机翼气动性能,分析研究了结构静弹性变形对气动载荷的影响.在此基础上以机翼典型剖面外形和机翼型架外形为设计变量进行了基于控制理论的跨声速弹性机翼气动优化设计方法研究,实现了在真实飞行条件下考虑静气动弹性变形影响的大展弦比跨声速弹性机翼一体化优化设计.优化设计算例结果表明所发展的方法是成功的,优化得到的弹性机翼在满足升力系数约束条件下提高了升阻比.     

全复合材料主承力梁结构研制与分析验证

根据结构特点和受力特性,研制了一种全复合材料主承力梁结构。给出了主梁的铺层设计,详述了其成形工艺。基于MSC.Patran/Nastran建立了主梁有限元分析模型,采用线性求解器对多级载荷作用下的主梁进行了静强度分析,给出了结构应变与位移结果。同时设计了强度试验方案,对主梁结构进行了静强度试验验证,采集了应变与位移信息,确定了结构的响应。结果表明,全复合材料主承力梁结构满足强度设计要求,为复合材料应用于主承力结构提供了依据。     

复合材料大斜削结构压缩试验和损伤模拟研究

对复合材料斜削结构进行了压缩试验和相应的有限元模拟分析,编制了损伤子程序模拟结构内部的损伤.通过对试验数据和数值模拟结果的对比,发现数值模拟结果和试验数据吻合较好.试验和模拟结果显示,由于削层的存在,使得结构更复杂,尤其在过渡段内,各截面的性能互不相同.通过损伤模拟,预测了三种损伤形式的萌生位置和扩展方向,发现结构同时发生基体开裂损伤和纤维-基体剪切损伤,然后再发生纤维屈曲损伤.基体开裂损伤和纤维-基体剪切损伤的发生对结构的承载能力有严重影响,是灾难性的.而随后发生的纤维屈曲损伤加剧了这种影响,此时结构完全失去承载能力.     

复合材料机身壁板的强度分析与试验验证

复合材料在大型客机机身主承力结构应用是近年来的发展趋势，通过分析与大型壁板试验结合的方式研究了复合材料机身壁板的静力承载能力，以及疲劳与损伤容限特性。采用理论公式、半经验公式、有限元模态分析研究了蒙皮的屈曲载荷、壁板的承载能力。依靠创新的机身壁板多轴载荷试验系统，模拟机身壁板的实际受载情况，实现充压载荷、拉伸/压缩、剪切载荷的独立施加与组合施加。通过静力试验验证蒙皮屈曲的工程及有限元分析方法和壁板的剩余强度承载分析方法。引入预埋缺陷、BVID、VID冲击损伤，通过试验研究了损伤对应变分布的影响，并通过疲劳、损伤容限试验，验证了壁板的设计以及损伤的无扩展特性。     

民机中央翼加筋壁板承载能力的非线性有限元分析

应用非线性有限元分析技术,采用MSC.Patran为前后置处理器进行建模,以具有较强非线性功能的Marc软件为求解器,细致研究了某民机中央翼加筋板结构在受力状态下的非线性变形、后屈曲破坏过程、最终失稳破坏形式及载荷,与试验结果进行了比较,破坏过程及破坏形式与试验相吻合,计算的破坏载荷与试验的破坏载荷误差在5%以内,显示文中方法可以满足工程设计和分析要求。     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

NNW-FSI软件静气动弹性耦合加速策略设计与实现

在国家数值风洞（NNW）工程的资助下,依托NNW-FSI流固耦合模拟软件平台,从气动载荷作用下飞行器结构静变形大小与收敛过程无关的物理机制出发,基于变形增量叠加的方式,设计和实现了一种静气动弹性耦合加速策略,通过松弛因子对耦合迭代的收敛过程进行调整。结合超大展弦比无人机和CHN-T1模型两种不同外形,开展了不同松弛因子下的静气动弹性耦合数值模拟,对耦合加速策略的参数影响和加速效果进行了测试和评估。从计算误差控制角度对松弛因子加速耦合迭代收敛的作用机制进行了理论分析,弄清3种类型静气动弹性耦合模拟过程中松弛因子发挥的作用,并给出了松弛因子选取范围的建议。静气动弹性耦合模拟和理论分析结果表明,针对不同类型的静气动弹性耦合问题,选取合适的松弛因子,能够达到抑制振荡并加速收敛的效果。     

大型 CFD 软件自动化测试平台的初步设计与实现

针对大型通用 CFD 软件研制的需求，为了提升 CFD 软件的开发效率，初步设计与开发了一个 CFD 软件自动化测试平台。该平台基于 MySQL 数据库，通过前台界面和远程集群后台管理程序的网络信息交互，实现了自主研发的通用 CFD 软件平台(HyperFLOW)的自动测试及验证与确认。根据需要适当修改前置处理接口，该平台亦可推广应用于其他 CFD 软件的自动测试。本文重点介绍了平台的基本框架和总体设计思路、平台的数据库和前台管理 GUI 设计、网络交互和后台管理等，并利用典型算例对计算结果的验证与确认过程进行了简要介绍，最后对平台构建的后续工作进行了展望。     

Buckling and post-buckling behavior of titanium alloy stiffened panels under shear load

Titanium alloy has been increasingly applied in aviation industry due to its superior performance. However, the titanium alloy structures are less studied. This work investigates the structural behavior of Ti6Al4V titanium alloy stiffened panels under in-plane shear load by experiments and numerical analysis. After the shear tests, the buckling instability, the post-buckling process and the failure mechanism of the specimen were obtained. The Finite Element(FE) models were established with the subsequent validation verification. A parametric analysis was implemented to study the influence of stringer thickness and stringer height on the behavior of the stiffened panels. The results show that after the initial local buckling on the skin, the buckling mode jumps several times with the increase of load. The stringers twist when the load reaches a certain level, and finally the structure damages due to the plastic deformation and the global buckling. The shear clip has little effect on the buckling and failure loads. Compared to the relatively large effect on the buckling load, the influence of the stringer thickness and stringer height on the failure load is neglectable.According to the parametric analysis, the stringer thickness influences the final buckling mode and failure mode, while the stringer height affects the buckling mode transformation.