三维弹塑性边界元1/r~3积分奇异性的处理

本文将文[2]的基本思想与方法,推广应用于处理三维弹塑性边界元法塑性单元中的1/r~3积分奇异性问题,推导了有关的具体公式,并从理论上证明了该方法处理1/r~3奇异性的有效性。针对奇异性单元。本文提出旋转与坐标转换的分割技术,可大大减少工作量。并使计算机编程更简洁。数值算例研究表明,文中采用的方法是十分有效的。     

一种采用解析积分的直线单元的边界元法

本文提出一种采用解析的直线单元计算受有离心力的二维问题的边界元法,推导了相应直线单元的积分解析式。算例表明,本文方法能获得较高的数值精度,尤其在计算靠近边界的内点应力时,效果更好。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

基于反馈定理的IB-LB方法(英文)

格子波尔兹曼方法(Lattice Boltzmann method,LBM)和浸入边界方法(Immersed boundary method,IBM)皆为近年来发展的可替代Navier-Stokes(N-S)方程求解复杂流体力学问题的数值模拟方法。本文采用浸入边界-格子波尔兹曼方法(IB-LBM),将IBM作为一种边界处理格式应用到LBM中,模拟了静止单圆柱及并排双圆柱绕流,与文献对比结果吻合良好,证明该方法具有模拟复杂边界物体绕流的能力。与采用虎克定理、直接力法计算体积力项的传统方法不同的是本方法运用反馈定理,即采用流场和流场内物体之间的速度、位移反馈来计算体积力项,使得该方法易于实施并具有优良的并行性。     

劲条止裂性能的分析方法

对于含裂壁板，采用劲条进行局部加强，可以阻止或延缓裂纹的扩展。本文对此作了理论分析，提出了劲条跨过裂纹对称加强时，裂纹尖端、劲条与壁板胶结界点以及劲条端点处应力强度因子的计算方法。在此基础上，通过数值计算，考察了劲条的刚度参数以及问题的一些几何特征参数对劲条止裂性能的影响规律。结果表明，本文方法对于劲条止裂性能的分析是有效的，它不仅能够准确地处理其中的应力奇异性问题，而且运算快速省时，便于进行优化设计。     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

一种快速的有约束气动CFD优化方法

提出了一种快速应用型罚函数法来求解有约束气动优化问题。多设计变量的优化算例表明：新方法与传统罚函数的序列极小化(SUMT)技术相比，避免了SUMT技术选择罚因子序列的问题，能在约束边界和不可行域进行计算，鲁棒性好，且求得约束最优点的效率高于SUMT技术，能减少30％～40％的计算量。本文分别将之与高分辨率Euler方程和N—S方程数值解相结合，进行了两个气动问题的优化。在二维无粘跨声速翼型的优化设计中，得到了升阻比明显提高的新翼型；对某型无人机高亚声速三维进气道设计，也获得了性能优秀的设计方案。     

直升机旋翼计算流体力学的研究进展

依据直升机旋翼计算流体力学(CFD)的发展，分别介绍了小扰动位势方程、全位势方程、Euler方程和Navier—Stokes方程在旋翼流场计算中的研究现状和发展趋势，着重从方程离散、数值算法、网格生成、计算时间和计算精度等方面分析了不同旋翼CFD方法的特点，指出了旋翼的尾迹在旋翼流场计算中的重要性，并针对位势方程和Euler／N—S方程分别讨论了求解的边界条件。最后，对旋翼CFD的发展提出了几点展望。     

新电位移边界条件下的压电介质断裂问题

应用边界元法,并结合一种特殊的复变函数基本解,研究了压电介质的断裂问题。由于该基本解完全遵循本构方程没有对电位移边界作假设,所以结论同非穿透性裂纹假设的结果有较大差异。通过计算表明,外加力载荷不但会引起裂纹尖端应力奇异性,也会引起电场奇异性,同样,外加电压也会引起裂纹尖端应力奇异性。     

有进动的高速旋转体的应力分析

本文讨论了采用有限元法对具有进动的高速旋转的轴对称体进行应力分析的方法，首先，推导了轴对称体中哥氏惯性力的分布公式，在柱坐标系下，哥氏力中含有沿径向的零阶对称不载荷和沿轴向的一阶对称及一载荷，首次，采用了具有计算精度，对边界适应性强的八节点四边等参数环形单元，并将受哥氏惯性力的轴对称体看作非轴对称问题，对有限元分析中位移模式的建立和哥氏惯性力的处理作了讨论。最后，简要地说明了所研制的计算机程序，并     

Winkler弹性地基上自由矩形板问题

在经典薄板理论中,弹性地基上自由矩形板弯曲问题的精确解,长期以来被认为是一个相当困难的课题。根据对矩形板精确解的研究,构造了一个包括三角函数和多项式组成的近似挠度函数。它满足四个自由边上的全部边界条件和自由角点条件,利用伽辽金方法得到Winkler弹性地基上自由边矩形板弯曲、自由振动和稳定问题的解。还给出了数值算例。     

弹性体表面覆盖层对冲击应力波的响应

本文研究了弹性体表面覆盖层对弹性体内部由于瞬态加载或者材料突然断裂等冲击因素所引起的应力波的动力响应问题。提出了一种对该类问题进行初步分析的简化模型,在此基础上建立了应力波场的分析方法,并借助于积分变换技术求得了问题的解析解。最后通过算例分析讨论了问题的两种特征参数对覆盖层位移响应及界面应力响应的影响规律。结果表明,本文所提出的模型及应大波场的分析方法对于初步研究覆盖层对冲击应力波响应问题的动态特征是有效的。     

高速碰撞中Lagrange有限元方法及其应用

根据连续介质力学守恒方程和有限元离散方法,针对三维高速碰撞现象的特性,应用瞬时最小闰能原理,完整地推导了高速碰撞有限元基本列式和离散方程,并给出了高速碰撞问题的中心差分求解算法,讨论了非线性Ｌａｇｒａｎｇｅ动力有限元程序数值计算中应重视的问题,本文给出的典型算例说明本方法对高速碰撞问题分析的有效性。     

大切深数控电解铣削加工阴极设计及试验研究

难加工材料整体叶轮广泛应用于航空领域,采用传统切削加工存在刀具磨损快、加工效率低等问题。本文针对某型号复杂整体叶轮,提出大切深五轴数控电解铣削预加工方法。通过设计锥形螺旋刃阴极,分析不同旋转角下单、双螺旋刃出口流场分布,得到旋转角720°的单螺旋刃阴极出口压力和流速分布均匀。同时开展大切深数控电解铣削加工试验,结果表明:在选取的工艺参数范围内,加工平衡间隙和进给速度随着加工电压升高而增大;较低的电解液温度有利于实现小间隙加工,可显著提高加工精度;主轴转速达到1 500 r/min后对加工速度影响较小。得到大切深数控电解铣削整体叶轮加工叶片,一次最大切深可达65 mm,余量误差控制在0.5 mm范围之内,提高了整体叶轮加工效率。     

双参数弹性地基上四边自由矩形薄板精确解

将弹性地基视为Vlazov双参数模型,利用二维有限域积分变换的方法推导出了Vlazov双参数弹性地基上四边自由矩形薄板在任意荷载作用下挠度和内力的精确解。在求解过程中,不需要预先人为选取挠度函数,而是直接从弹性薄板的基本方程出发。通过集中荷载计算实例验证了该方法的正确性。该方法计算简便,适用于不同边界的薄板问题求解。     

非线性薄梁结构对冲击波载荷的大挠度响应分析

采用大挠度弹性梁理论和中心差分法分析计算了非线性薄梁结构对冲击波载荷作用的大挠度响应，模型包括了纵向变形和横向变形之间的弹性耦合。算例给出了拉－弯耦合、拉／弯刚度比和冲击波冲量对挠度响应的影响。     

直升机旋翼/机身耦合系统气动/机械稳定性分析新方法(英文)

研究了前飞状态下直升机旋翼 /机身耦合系统的气动 /机械稳定性问题。根据柔性多体系统动力学理论 ,通过构造一种 2 4自由度的刚柔混合单元得到旋翼 /机身耦合系统的周期时变运动方程 ,建模中考虑了桨叶预锥、后掠、中等弹性变形以及直升机机身和传动轴的弹性影响 ,体现出铰接式桨叶绕挥舞、摆振和变距铰的整体刚性运动与桨叶中等弹性变形之间的动力学耦合作用 ,推导中对桨叶挥舞、摆振和变距转角幅值未加任何限制。根据 Floquet理论对稳态周期解的稳定性进行研究 ,采用 Newmark直接数值积分方法得到转移矩阵。对某新型直升机的气动 /机械稳定性进行了分析 ,结果发现对于给定的前飞状态是稳定的 ,但是随着传动轴弯曲和扭转刚度的降低出现不稳定现象。     

三维接触问题的接触状态有限元解法及应用

在文献［１］的基础上论述了有摩擦、有间隙的三维弹性接触问题接触状态有限元解法，并通过理论考题进行了验证。最后，在应用部分对汽轮机汽缸法兰连接汽密性进行了计算。     

具有多点反馈的可控约束层阻尼

讨论了具有可控约束阻尼结构的多点杂交阻尼控制问题,建立了弹性,粘弹性和多个压电片组成的多层复合梁的偏微分方程线,通过模态转换和应用粘弹性材料的振子模型对模型进行减缩。以实际中可测量量作为反馈量进行次优控制,数值模拟说明,提出了控制方式抑制振动效果好,控制电压低,易于实现。