整体弹性结构法及在叶片流固耦合分析中的应用

针对三维叶片时域流固耦合振动响应计算普遍耗时的问题,采用一种假设整体结构的模态位移,求解固体时域响应和实现高效网格变形,发展了一种3维时域流固耦合分析的整体弹性结构方法,并应用于压气机叶片0°和180°相位角的气动稳定性分析中。结果表明,所发展方法的计算结果与传统双向时域算法和文献的结果较为吻合,而计算效率相比于传统算法显著提升;在所分析的两个相位角下,叶片振动的气动阻尼均随流量减小先增大后降低,相比于0°相位角,180°下叶片的气动稳定性大幅提高,表明该方法能有效应用于叶片的工程流固耦合研究。     

燃烧室结构高温变形非接触测量试验研究

为获得地面试验中主动冷却燃烧室在高温状态下的结构总变形和局部应变,基于数字图像相关法开展了非接触测量试验研究。分别采用加速鲁棒特征算法和改进的圆柱曲面数字图像相关法完成了燃烧室结构总变形和局部应变测量及分析。测量结果表明,燃烧室稳定工作条件下,结构轴向总变形约4.8 mm、局部平均应变约0.0049。测量结果与工程粗略估算结果比较吻合,说明所使用的非接触测量方法有效,能够支撑燃烧室热结构设计,测量数据能够用于三维结构强度数值计算的验证。     

C/SiC复合材料螺钉拉伸强度分布模型

C/SiC复合材料螺钉是在高超声速飞行器上应用越来越广泛的一类重要紧固件,但其拉伸性能存在较大散差,分布规律尚不明确,给材料选用和结构设计带来了很大困难。本文采用电子万能试验机对M8、M10、M12三种规格的平头C/SiC复合材料螺钉进行力学性能试验,并分析了拉伸强度的分布规律。在此基础上应用双参数Weibull模型对统计数据进行拟合,并对拟合结果进行了柯尔莫哥洛夫检验。结果表明:C/SiC复合材料螺钉的拉伸性能分布满足双参数Weibull模型,其特征强度β为212 MPa,形状参数α为9. 45,可以据此进行复合材料许用强度设计。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

考虑弯扭变形的叶片模型配准方法

叶片测量数据与其CAD模型的配准定位是叶片几何形状分析的关键步骤。针对涡轮叶片的弯扭变形对叶片模型配准定位的不利影响,提出一种考虑弯扭变形的叶片模型配准方法,以提高叶片模型配准的可靠性与鲁棒性。在迭代最近点算法的目标函数中引入预变换函数,建立考虑弯扭变形的叶片模型配准目标函数。对叶片模型做切片化处理,进行弯扭变形分析生成叶片弯扭变形曲线,使用叶片弯扭变形曲线指导叶片测量数据配准。使用仿真生成带弯扭变形的叶片数据与通过坐标测量机获取的叶片实测数据对配准方法进行了验证。结果表明,所讨论方法能提高模型配准的可靠性,避免弯扭变形可能导致的叶片模型配准失败。     

TC11合金高温变形行为及其机理

采用Gleeble-1500热模拟机研究了TC11合金在800～1 050℃、应变速率0.005～5/s条件下的高温变形行为.根据动力学分析,确定了不同温度区间的热激活能和热变形方程.结合变形微观组织观察确定了TC11合金的高温变形机制.结果显示:TC11合金在(α β)两相区和β相区的热变形激活能分别为285.38和141.98 kJ/mol,表明不同温度区间的热变形机理不同;在两相区变形主要发生片状组织的球化,在β相区变形时低应变速率下(0.005～0.05/s)主要发生β相的动态再结晶,高应变速率下(0.05～5/s)主要发生动态回复.研究结果为确定该合金的最佳变形工艺参数提供了理论依据.     

RFI工艺成型两种环氧树脂基复合材料性能比较

采用RFI工艺分别成型了648和5228A环氧树脂基复合材料层合板,其增强材料为碳纤维无屈曲织物,铺层方式为[(0,90)/( 45)]s;测试了两组层舍板的拉伸性能、弯曲性能和层间剪切性能并做了比较分析;对破坏形式和机理进行了探讨.结果表明:5228A相对于648环氧树脂膜有较宽的低黏度区域,较长的凝胶时间;5228A与648层合板相比,拉伸强度高106%,拉伸模量、泊松比接近;弯曲强度高58%,弯曲模量高16%;层间剪切强度高62%.     

预拉伸条件下铝合金焊接残余应力的数值模拟

利用ANSYS有限元分析软件,通过APDL语言编程,实现了不同预拉伸应力条件下铝合金平板对接焊过程中,温度场及应力场的模拟分析.计算结果表明,预拉伸焊接法可以有效地控制焊接残余应力,随着预拉伸应力的增大,其焊后残余应力值逐渐减小.当预拉伸应力σP从0增加到90%σ0.2时,纵向残余应力降低了85.6%.模拟分析结果与实验测试结果基本吻合.     

关于变形网格“几何守恒律”概念的讨论

在国内外文献中应用包含有变形的动网格进行流动数值模拟时,要满足所谓的几何守恒律条件;通过对这一概念最早提出推导过程和近年来应用发展情况分析以后发现,逻辑上存在不适当之处.从流体力学基本理论出发,对包含有变形的动网格技术的理论基础进行讨论,证实几何守恒律是流体力学控制方程的伴随方程或退化方程;通过简单模型分析有限体积方法离散过程,发现计算过程中不满足几何守恒律所引起的非物理现象本质上是目前有限体积离散处理中微元体积计算方法不符合物理定律所致;因此在变形动网格计算方法中不存在必须要满足的几何守恒律.最后根据理论提出物理背景清晰的微元体积计算方法,数值验证可行.     

基于有限元模拟的三维型材拉弯轨迹设计

针对三维型材拉弯,提出一种基于变形控制有限元模拟的轨迹设计方法。首先提取模具长度方向特征线并将其离散成诸多线段单元,这样毛料逐步包覆模具的过程就变成毛料包覆这些线段单元的过程。给定毛料依次包覆线段单元发生的变形,根据切线接触条件（拉弯成形过程中毛料离开模具的位置处两者相切）计算出每步毛料末端位移。将这些位移作为边界条件输入有限元模型中计算毛料应力、应变和回弹。根据成形极限图和规定最大回弹超差量来调整变形模式,找到优化的变形模式和拉弯轨迹。以中空矩形截面型材三维拉弯为例,给出了轨迹设计的详细流程。