用二维N—S方程计算汽车流场的地面效应

采用有限体积方法求解了Navier-Stokes方程,计算了绕二维类汽车物体的绕流并研究了地面效应。在非正交网格系统上以强守恒形式构造出有限体积方程,采用压力校正算法来求解所得到的方程。对于湍流情况,采用了双方程模型的κ_ε模型。 为了验证采用本程序计算移动壁面效应的正确性,本文给出了Couette流动的计算结果。 本文也给出了汽车绕流的计算结果,它包括层流情况,从汽车头部开始由层流转为湍流的混合流动情况及湍流情况。对于每一种情况都比较了移动壁面的计算结果和固定壁面的计算结果。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

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

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

粉末高温合金粘塑性试验评定与本构模型参数估计

 开展了粉末高温合金 FGH95 550℃、600℃和 650℃等 3种温度下控制应变率单向拉伸试验和 550℃下循环加载试验研究,结果表明 :600℃以下,快、慢应变率时,5%的试验应变范围内应力—应变曲线都一直上升,不存在应力饱和现象,热恢复效应不显著;但 650℃下慢应变率时则存在较明显的应力饱和现象,反映出在此条件下必须考虑蠕变效应。温度越高应变率对 FGH95的拉伸力学性能影响越明显,但总的说来是一种应变率不甚敏感的循环硬化材料。最后,在试验的基础上建立了 FGH95的 Bonder-Partom统一弹-粘塑性本构模型,理论与试验吻合较好,表明该模型能够模拟 FGH95的应力-应变关系曲线、应变率响应特性以及循环硬化特性,从而为 FGH95粉末高温合金构件的高温应力分析打下了基础。     

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

直升机非线性运动方程及其数值分析

基于牛顿法建立了单旋翼带尾桨直升机的非线性运动方程，提出了一种求解直升机非线性运动方程的数值方法。建立方程时，用欧拉角描述直升机在空中的姿态，旋翼采用有挥舞铰且在铰上带有弹性约束的模型，桨盘处的入流采用线性模型，通过对一算例直升机的操纵响应进行数值仿真，将非线性与小扰动线性化的直升机运动方程进行了对比分析。结果表明：从操纵输入到随后的2s左右，非线性和小扰动线性化的直升机运动方程的计算结果非常接近，但在此之后，非线性运动方程的计算结果更能反映直升机的运动规律。     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

改性BMI／DPA和CTBN增韧环氧树脂的温度,力学性能及本构关系研究

研究了ＣＴＢＮ增韧环氧树脂和烯丙基双酚Ａ改性双马树脂（ＢＭＩ／ＤＰＡ）在不同温度和加载速率下的拉伸和压缩力学性能。对于ＣＴＢＮ 增韧环氧树脂,材料在动态拉伸情况下,随着应变率的增大,材料由韧变脆,模量升高,断裂伸长率下降;在压缩情况下,弹性模量和屈服应力均随着加载速率的增大而增大。对于改性双马树脂,在室温的动载压缩下,其模量Ｅ和屈服应力σｙ均基本上与ｌｎε呈直线关系增大,并给出了该树脂在不同加载速率下的压缩应力应变曲线。计算了两种材料的温度敏感度和加载速率敏感度以及温度和应变率对两种树脂材料屈服应力影响的本构关系。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。