一种基于共轭方程法求解黏性反问题的简化方法

 对于给定压力分布的黏性气动反问题,考虑到壁面微小扰动造成的压力变化主要由势流作用引起,因此可以简化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解,则其相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解,网格数减少,同时收敛速度加快,简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明,对于附着的边界层流动简化方法计算得到的敏感性导数具有较高的精度,能够有效完成反问题设计且减少总的计算耗时。     

复合材料厚板双轴非线性刚度特性分析

基于非线性本构关系与厚板理论,研究了一种采用增量步迭代法对复合材料厚层合板渐进失效过程进行双轴刚度特性分析的方法。以三维Hashin失效准则为单层板失效判据,当层合板中有一层或多层被检测到失效时,对失效层进行刚度折减,并在增量过程中不断更新刚度矩阵。采用有限元软件MSC.Patran自带的编程语言PCL,编写了计算程序。应用该程序分别对不同材料和不同铺层的层合板算例进行非线性刚度计算,给出了不同的铺层和材料组合下层合板在不同载荷比下的应力-应变曲线,展示了层合板渐进失效过程中应力-应变的变化关系,为复合材料厚板失效分析过程中的刚度变化提供理论与计算依据。     

一种基于微观组织演化的率形式本构模型

为建立塑性变形中微观组织演化与宏观响应的关系,依据定义的率扩散势导出了考虑温度影响的率形式流动应力方程,修正了Chaboche非线性各向同性硬化准则以考虑位错密度演化对硬化的贡献,以位错密度为内变量建立了位错密度演化模型,提出了一种基于微观组织演化机理的率形式本构模型。以927℃钛合金Ti-6Al-4V热单拉实验数据为基础,通过遗传算法确定了模型材料常数。将计算结果与实验数据进行比较,最大应力误差在12%以内,最大晶粒尺寸误差在4%以内,证明了模型的有效性。     

一种新PSO混合算法在直升机配平中的应用

直升机配平计算是动力学分析的基础,其实质是求解高维复杂的非线性方程组。针对经典算法与智能算法的特点与不足,提出了一种求解非线性方程组的新粒子群方法。在粒子群(Particle swarm optimization, PSO)算法的基础上,根据模拟退火(Simulated annealing, SA)思想,引入了嵌入式LM (Levenberg-marquardt)优化 算子。该方法充分发挥了3种算法的优势,克服了LM算法初值敏感性,PSO算法易陷入局部极值等问题。通过UH-60A直升机实例配平计算,验证了本文算例模型的准确性。在此基础 上,针对某一前飞状态下的配平算例,在收敛可靠性和计算效率上通过与其他算法进行对比,表明该算法具有可靠的收敛性和较高的计算效率,进一步验证了该算法在配平问题上的可信度与实用性,为直升机飞行动力学问题的处理提供了一种新的有效方法。     

风力机翼型阵风响应计算与分析

基于网格速度方法,对阵风条件下风力机翼型运动的气动特性进行数值计算。中心格式有限体积法、双时间推进法求解预处理后的非定常N-S方程。计算结果表明:在弦向阵风作用下,翼型的运动能有效降低响应波动峰值,即弦向阵风对运动影响较小;而在法向阵风作用下,气动响应波动很大。计算翼型在不同参数阵风作用下的响应过程,阵风的频率变化对响应影响不大。弦向阵风的幅值增大,叶片运动对阵风响应峰值减缓效果增加;法向阵风的幅值增大1倍,响应峰值也增加1倍。     

压力环境下固体推进剂含损伤的统一本构研究

为了建立固体推进剂压力环境下的统一的损伤本构方程,利用改进的万能材料试验机对一种固体推进剂在常温下、不同的压力环境进行了准静态拉伸实验。利用ABAQUS软件的二次开发用户材料子程序UMAT接口,根据实验得到的结果及其应力应变响应特征,确定了本构方程中的系数。计算结果说明本构方程能够很好地描述该推进剂在压力环境下的力学行为。     

二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

数值模拟二维跨音速叶栅可压粘性流动

用有限体积法数值模拟了任意回转流面跨音速压气机叶栅中的粘性流动，四步龙格－库塔法用于N－S方程的时间推进。给出了一种抑制叶栅计算中“数值失速”的方法，并采用隐式残差光顺技术加速收敛。用所发展的方法计算了三种叶栅在宽广来流Ma数条件下的流动，计算结果与实验结果吻合良好，流场激波具有较高的分辨率。     

FGH95粉末冶金材料变形及破坏特点的试验研究与数值模拟及分析

进行了650℃下不同应变率的拉伸试验和应变率为10^-3／s的应变控制循环试验，用以研究FGH95材料的变形特征；进行了不同保载形式的疲劳试验，用以研究FGH95材料的破坏特征。采用Chaboche本构对材料的变形特征进行了数值模拟，同时也对其寿命进行丁评估，得到了较为理想的结果，为粉末高温合金构什的应力一应变分析及寿命预测打下了基础。同国外相近牌号Rene’95相对照，得出了粉末材料FGH95一些特有的性能特点，对其工程应用具有一定的参考价值。     

大型稀疏复线性方程组双共轭梯度法

有限元复线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改.在一般双共轭梯度法的基础上,本文利用广义变分原理对内积进行了重新定义,使双共轭梯度法求解复线性方程组更为有效.数值算例表明这种双共轭梯度法结合全稀疏存贮方案的求解算法在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.