MLPG混合配点法分析复合材料层合板的弯曲问题

提出了分析复合材料层合板弯曲问题的无网格局部彼得洛夫-伽辽金混合配点法,层合板理论采用经典理论。通过移动最小二乘近似构造形函数,采用配点法建立系统平衡方程,推导了系统刚度矩阵和载荷力向量。采用配点法和罚函数法分别施加自然边界条件和本质边界条件。计算实例表明,MLPG混合配点法具有较高的计算精度,并且相比其它的无网格法具有更高的计算效率。     

无网格法在旋翼空气动力学计算上的应用初探

本文介绍了无网格法与现有的旋翼空气动力学有限元计算方法的区别,无网格法的发展历史和基本原理,并通过在直升机悬停状态下的旋翼流场计算初步验证了该方法的可行性。     

基于有限元法的无刷交流励磁机空载特性分析及额定励磁电流计算

针对无刷交流励磁机铁心结构的复杂性和铁磁材料饱和特性的影响,采用非线性有限元法,利用有限元软件Maxwell,建立无刷交流励磁机的二维模型,对空载特性及额定励磁电流进行研究分析。以1台28 kW无刷交流励磁机为例,对电机气隙磁场仿真分析并求取空载特性。计算所得的基波电势与实测结果相符,为此类电机电磁场的设计提供理论依据。在此基础上通过场路耦合的方法,用外电路来模拟额定负载工况,计算出额定工况下的励磁电流。该方法具有实用性,可为一般的交流励磁机励磁电流计算提供参考。     

电磁场作用下包晶合金定向凝固过程传输模型的建立

建立了二元包晶合金在电磁场作用下定向凝固过程中多场耦合的数学模型,并根据包晶转变不同阶段的特点,提出了包晶转变中固相分数的计算办法。该模型可以计算电磁场作用下包晶合金定向凝固中的流动、传热、传质的过程,为电磁场作用下包晶合金定向凝固过程的数值计算及对实验结果的分析提供理论依据。     

临近空间高速飞行器抛罩过程的动网格方法对比

飞行器抛罩是一个复杂的动力学、流体力学耦合问题,涉及六自由度运动方程与N-S方程的耦合求解,其中动网格非定常计算是关键技术。针对二维临近空间高速飞行器,将耦合过程简化为匀速旋转抛罩,对比分析了光顺重构法、重叠网格法、域动分层法三种工程易行的动网格方法的仿真结果,并得到如下结论:三种方法在整流罩关闭状态下的定常流场结果一致,均能捕捉到非定常开罩过程的典型特征,得到正确的起动结果;在非定常过程中光顺重构法的分离区吞入速度慢于重叠网格法和域动分层法;重叠网格法的计算通用性最好,域动分层法的计算速度最快,光顺重构法的理论精度最高;由于光顺重构法的网格更新对于复杂模型容易失败,域动分层法只能处理运动轨迹已知的问题,三维动力学耦合计算建议采用重叠网格法;在进行整流罩的三维设计时,应考虑溢流效应以缩小整流罩前方的大分离区,降低飞行器的控制及热防护上的困难。     

背绕式有槽永磁同步电机的电磁场解析模型

以1台5 kW背绕式高速永磁同步电机为研究对象,建立其电磁场解析模型。将电磁场求解域划分为气隙子域、永磁体子域、槽口子域和槽子域,求解相应的拉普拉斯方程或泊松方程,解析模型计及电枢反应场、永磁场和定子开槽的影响。计算了该电机的气隙磁密、绕组磁链、绕组反电动势、齿槽转矩和电磁转矩,并将结果与二维有限元法计算结果和试验数据比较,比较结果说明了解析模型的准确性。最后以槽口开度为变量,研究其对气隙磁密分布和齿槽转矩的影响。     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

基于指数插值的浸没边界法在可压缩流模拟中的应用研究

由于浸没边界法在处理复杂外形和动态问题时具有极强的灵活性和便利性,因而受到广泛的关注。在当前浸没边界法的研究和应用中,主要工作集中在关于中/低雷诺数和不可压流动的模拟。本文引入 Choi 等提出的指数插值法对目标点速度进行插值,并根据多项式插值和 Crocce-Busemann 关系式等确定目标点压力和温度,并借鉴直角网格法中的虚拟网格法计算物面内的流动变量,建立了适用于可压缩流动的浸没边界法。通过亚、跨声速段翼型绕流和超声速圆柱绕流算例对方法进行验证,结果显示在 k 值取值合适时模拟结果与实验数据较为吻合。通过对相关算例网格收敛性分析,证明了本文发展的浸没边界法具有较高的精度。文中各算例的模拟结果表明,此方法能准确捕捉亚、跨、超声速绕流问题中的各类流动结构,适用于可压缩流动的模拟。     

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

气动加热下飞行器油箱温度场耦合分析

气动加热与结构热传递耦合问题在航空航天领域非常重要.分析了飞行器持续受热过程中的耦合性,指出了耦合性在分析持续气动加热问题中的必要性.针对存在的耦合性,提出了一种处理方法,即在“飞机结构三维温度场分析”程序中增加气动热流工程算法,同时采用对流传热来模拟油箱内燃油与油箱内壁面传热.对流换热系数以及出口油温通过与油箱壁面结构温度耦合计算来确定,即在每个时间点上,将当前结构内壁面温度咒代入对流换热系数α以及油温L的计算公式,求得该时间点的α和t,以此实现气动/温度场/流场三个物理场的实时耦合.计算结果表明,耦合计算结果与试验结果吻合得很好.     

考虑磁滞效应的电抗器二维电磁场计算

将JA模型与电磁场数值分析方法相结合,在考虑铁心的磁滞效应的情况下,推导出电磁场的有限元计算模型。分别采用基本磁化曲线和磁滞回线两种不同的方法对电抗器的电磁场进行仿真计算。研究结果表明磁滞现象对电抗器电磁场的计算影响较大,因此考虑磁滞效应的电磁场数值计算模型,为进一步研究电抗器的发热问题奠定了理论基础。     

基于波形重构法的永磁电机电磁振动噪声计算研究

提出了基于波形重构法的径向气隙磁密和径向电磁力波解析计算方法,简化了永磁电机电磁振动和噪声(EVAN)的计算过程,解决了永磁电机径向电磁力波计算较为复杂的难题。首先介绍了波形重构法计算永磁电机EVAN的理论依据。然后,以1台永磁电机为例,详细介绍了运用波形重构法计算永磁电机电磁振动噪声的具体步骤。最后,对电机EVAN进行多物理场有限元仿真,将波形重构法计算结果与有限元仿真结果对比,在误差允许范围内二者的计算结果具有一致性,由此验证了波形重构法计算电机EVAN的准确性。     

基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

一种基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

本文发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型。该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在二者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解。这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率。该模型先得到了Levine and Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响。     

温度对约束阻尼板振动性能影响分析

基于复模量本构关系,推导各向同性约束阻尼板的动力学控制方程,求出了四边简支约束阻尼板的固有频率以及损耗因子的解析解,与利用有限元方法数值计算结果进行了对比验证。根据粘弹性材料特性曲线,插入取样温度点,分析不同温度下约束阻尼结构的减振性能,与数值结果进行了对比验证。结果表明,存在最佳温度,使约束阻尼结构减振效果最好。     

齿轮耦合复杂转子系统弯扭耦合振动分析的轴单元法

针对存在非平行轴的多分支齿轮耦合复杂转子系统,基于子结构分析理论,将系统中的每一根轴视为单元,将齿轮啮合力视为单元的激振力,在局部坐标系中采用集总质量法建立该单元的振动方程;在系统总体坐标系中,将各单元的振动方程联立,形成齿轮耦合复杂转子系统弯扭耦合振动方程,通过数值求解,获得系统固有频率和动态响应。采用所提出的轴单元法,对某航空发动机启动系统进行了模态和强迫振动特性分析。      

基于IPDG方法的超声速混合层流动数值模拟

为了满足超声速混合层高精度模拟需求,实现了基于内罚方法的间断伽辽金(IPDG)方法数值模拟。通过将黏性通量作为辅助变量使得Navier-Stokes方程降阶,并利用间断伽辽金方法进行空间离散,最后采用Newton-Krylov隐式方法对空间半离散方程进行时间推进。相对于有限体积法数值精度提高到了3阶。将该方法应用于对流马赫数为0.2的二维平面超声速混合层的数值模拟中,通过与实验数据的对比验证了方法的可靠性。数值结果显示了混合层中层流到湍流的发展过程。在此基础上,将基于数值解误差分布的自适应网格技术与IPDG方法结合起来,对比发现自适应网格数量减少了9倍,计算时间减少了8倍,大幅提高了方法的计算效率。      

三维机翼亚临界定常流的局部线性化有限元法

本文使用了两种类型等参元,通过变分法,计算了三维机翼亚临界定常位势流中的压力分布,计算结果与试验结果符合较好。文中还使用了有限元法中的局部线性化理论,插值出元素中心点密度去处理迭代过程中元素系数阵的计算,其结果与通常的有限元计算相吻合,但计算时间却明显减少了。文中还对有限元法局部线性化理论和局部线性化有限元法做了阐述。     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。