板料冲压成形过程的一种数值模拟方法

在板料冲压成形过程中存在着大变形、大应变,且几何、物理及边界条件三重非线性相互耦合使得计算过程较为复杂.在虚功原理的基础上对板料冲压成形过程的弹塑性有限元分析方法进行了研究,并将有限变形弹塑性有限元法引入板料冲压成形过程,采用更新拉格朗日(ULM, Updated Lagrange Method) 方法描述,推导了有限元列式,建立了板料成形过程的有限元分析模型,开发了分析程序.通过对汽车油底壳的成形模拟说明了模型的正确性.      

基于多模型最优融合的双星定位系统一体化精密定轨方法

针对卫星动力学模型复杂且不准,考虑到卫星定轨中待估参数在时间和空间上的相关性,提出了一种基于多模型最优融合的双星定位系统一体化参数建模的近地卫星精密定轨新方法.利用节点自由分布B样条描述卫星运动,实现了对卫星粗略动力学模型的抑制作用;同时结合双星观测模型,使该方法转化为关于求解卫星轨道样条表示参数和定轨系统误差的多模型融合的非线性优化问题;通过引入模型结构确定最优融合权值的选取准则,在最小二乘准则下,采用非线性最优化方法搜寻样条的最优节点分布,得到了待估参数的最优估计,完成了近地卫星的精密定轨.理论分析和仿真计算表明,该方法确实有效,不仅提高了卫星的定轨精度,而且使状态估计的结构更加稳定.      

双级低滞后刷式密封级间不均衡性分析

针对一种结构参数的双级低滞后刷式密封,采用多孔介质模型进行流场计算,对刷丝安装角度以及第二级前板间隙高度和背板间隙高度与分压比的关系进行分析;基于梁弯曲理论和刷丝区域压力分布,采用有限元法计算得到刷丝的受力和变形情况,分析参数对转子所受摩擦扭矩的影响,讨论参数对双级低滞后刷式密封级间不均衡性的影响.分析表明,第一级刷丝安装角度的增大和第二级刷丝安装角的减小能在一定程度上使两级分压趋于均衡.第一级结构参数不变时,第二级的前板间隙高度变化对级间不均衡性影响不大,背板间隙高度在特定值时,两级可达到均衡状态.     

Reissner型矩形板分析的U变换-有限元法

将U变换法推广应用于Reissner矩形板的有限元分析中.对原结构进行等效变换,形成周期循环的板单元,使刚度矩阵成为循环矩阵,应用双重U变换解耦了有限元的矩阵方程,使有限元计算只须在一个板单元上进行,并且仍能方便分析整个板的一般分布载荷.所发展的U变换-有限元法不仅提高了计算效率,很快收敛于精确解,对于简支板还给出了精确分析的有限元解、准确的误差估计表达式和收敛速度,可以直接掌控计算精度,这是其它方法难以得到的.对简支和固支矩形板的数值算例及与其它方法的对比说明了U变换-有限元法的优点和重要的工程实用价值.      

反射器夹层面板精密成形原理

对高精度、大面积不可展双曲反射面的制造问题,提出了基于"离散钉模、真空负压、蜂窝夹层、应力释放、回弹补偿"的高精度反射面板成形原理;建立了四边自由面内多点约束的薄板大挠度弯曲解法及夹层板成形过程仿真模型,用于确定成形工艺参数;利用动力显式和静力隐式算法相结合的方法实现成形过程仿真;通过开缝和开槽方法释放面内压应力,确定了面板开缝和开槽方案;通过逆向修正离散模包络面,使其与成形模拟得到的夹层板型面误差最小,得到模具修正因子,补偿回弹误差.在此基础上开发了夹层板精密成形CAE系统.该方法制造的面板型面精度RMS达到25 μm以下,已成功地应用于多项大型紧缩场和毫米波天线的制造.      

四元数、Rodrigues参数在卫星姿态解算上的对比研究

详细阐述了罗德里格(Rodrigues)参数及改进型罗德里格(Modified Rodrigues)参数与四元数间的关系;推导了Rodrigues参数及改进型罗德里格参数微分方程和姿态阵,提出了后者避免奇异性的方法。通过仿真对以Rodrigues参数及改进型罗德里格参数描述的卫星姿态算法和以四元数描述的姿态算法在计算量和计算精度上进行了详细比较。仿真结果表明,Rodrigues参数和改进型罗德里格参数计算精度相当且计算量上优于四元数。     

复杂薄壁零件板多级充液成形及过程数值模拟

板充液成形已被证明是制造复杂零件一种有效的方法.采用板多级充液成形方法对难成形复杂薄壁结构的飞行器零件进行了研究,利用有限元数值模拟的方法,分析了预成形凸模形状、预成形深度,以及液室压力等关键工艺参数对成形的影响.对成形中出现的诸如破裂、起皱等失效形式进行了探讨,给出了终成形时实现成功成形的最高液室压力区域范围.提出了控制成形质量的措施,优化了工艺参数,并根据数值模拟的结果进行了试验验证.结果表明,模拟结果和试验结果达到了较好的一致.      

板料成形数值模拟接触搜索模型的研究

基于参数曲面描述法和有限元网格描述法,提出模具几何形状模型的曲面单元描述法,采用了逐级更新的接触搜索算法,推导出求交迭代初值优化选择方案及满足无穿透原理的接触判断准则,并建立了板料成形数值模拟接触搜索模型,此模型不仅改善了数值模拟的收敛性,而且提高了计算效率和稳定性.     

模糊变分原理在求解结构固有频率中的应用

在处理含模糊参数的结构固有频率问题时,现有的方法局限于先将模糊参数转化为区间数求解,再由区间解构造出相应的模糊解.由于现有方法计算量较大,所以提出了基于模糊变分原理的模糊里兹法和模糊有限元法.通过在瑞利商变分中直接引入模糊参数,得到了结构固有频率问题的模糊变分原理.作为模糊变分原理的应用,推导了模糊里兹法和模糊有限元法.这两种方法都可以直接得到问题的模糊解,与传统的区间解法相比,提高了计算效率.算例表明,这两种方法不但具有较高的计算精度,而且可以很大程度上减少计算量.     

双曲率开槽板精密成型有限元数值模拟

通过有限元数值模拟方法对双曲率开槽板的贴模性问题进行了系统研究.建立了开槽板在离散钉模上贴模性分析的有限元模型,研究了槽深、槽宽等开槽参数对面板型面精度的影响规律,其中相对槽深是影响开槽板贴模性的主要因素.通过分析计算确定了优化的开槽参数,为开槽方案的制定提供了可用的数据和重要依据.数值模拟和实验验证表明,工作板开槽可以解决紧缩场双曲率连续表面面板的胶接成型问题.      

板料冲压成形回弹的数值模拟

回弹是板料成形过程中不可避免的现象.回弹的数值模拟作为解决回弹问题的有效方法,越来越受到工业界和学术界的重视.完整的回弹数值模拟分为2步,即冲头加载成形过程模拟和冲头卸载回弹过程模拟.本文比较了静力隐式和动力显式2种算法在计算回弹问题时的优劣,采用静力隐式算法建立了计算冲头加载成形过程模拟和冲头卸载回弹过程模拟的有限元模型,并利用自行开发的软件给出了计算实例.      

闭截面型材的拉弯成形性

拉弯是型材的基本成形方法.截面畸变是型材拉弯成形中的常见缺陷,并对型材的回弹和变薄有交叉影响,型材拉弯成形性的分析评估远比板材成形复杂.提出了型材拉弯成形性的概念,建立了矩形截面型材拉弯成形性分析模型,基于有限元数值模拟研究了截面结构、相对高度和壁厚分布等几何参数对拉弯成形性的影响规律,从提高拉弯成形性出发提出了型材截面优化设计的原则.      

复杂结构振动控制设计与仿真方法研究

为了避免在建立结构数学模型过程中因简化和降阶引起的模型误差,解决复杂结构建模难的问题,提高控制器设计的准确性,提出了一种新的基于有限元法(FEM)的结构振动控制设计方法.该方法可借助于商用结构有限元分析软件来实现,简便易行,在PATRAN/NASTRAN环境中利用PATRAN命令语言PCL(Patran Command Language)建立了通用结构控制仿真平台.通过对压电自适应桁架结构的振动控制进行控制设计仿真和实验验证,证明了新方法的有效性;与Simulink使用简化模型的仿真结果相比,新方法的结果与实验结果更为接近,表明其准确性更高.该方法适用于对处在设计阶段的复杂结构进行控制设计与仿真评估.      

空间碎片超高速碰撞数值模拟的SPH方法

利用光滑质点动力学SPH(Smoothed Particle Hydrodynamics)方法对Whipple防护结构在空间碎片超高速碰撞下的物理过程进行了数值模拟.在数值模拟中,为了充分发挥SPH方法和有限元方法FEM(Finite Element Method)的优点,利用有限元单元和SPH节点混合建模,将有限元单元和SPH节点(SPH nodes)通过定义接触条件相结合,在大变形和飞溅区域采用SPH节点建模,而小变形区域则采用有限元单元建模,从而大大节省求解时间,提高计算效率.计算结果表明,弹丸在穿透前板后,形成二次碎片,碎片云经膨胀和拉长,对后板造成轻微的损伤,这和文献的相关试验数据是符合的.利用SPH方法对空间碎片的超高速碰撞过程进行数值模拟,不仅很好地预测了Whipple防护结构的破坏情况,而且对整个碰撞过程,包括碎片云的形成、膨胀和拉长过程都有形象的描述,符合超高速碰撞的试验测试结果.     

橡皮液压成形工艺的可靠性稳健优化设计方法

针对当前橡皮液压成形工艺过程中材料特性、工艺参数和其他随机因素都存在的固有波动问题,面向飞机精准制造,提出了一种考虑不确定性的可靠性稳健优化设计方法.该方法基于稳健优化设计和合理的可靠性约束,对比了可靠性稳健优化解和传统优化解相关度量指标,并对2B06-O态铝合金板材橡皮成形典型的曲弯边零件进行了工艺试验验证.结果表明该方法可有效降低质量评价指标对不确定性因素的敏感度,提高了橡皮成形工艺的稳健性和可靠性.     

2024-T3航空铝合金板材电磁V形弯曲应变分析

金属塑性成形后的应变特征影响零件的力学性能和使用寿命。电磁成形技术能提高金属的成形性,有助于克服铝合金室温成形性低的缺点,因而得到重视。本文以2024-T3航空铝合金板为研究对象,开展了电磁V弯成形和传统机械V弯成形的实验及数值模拟,研究了2种成形方法对弯曲成形试件V形弯曲圆角位置沿试件长度方向的应变特征的影响。数值模型表明,电磁V弯试件弯曲圆角处的拉应变峰值低于机械V弯试件;电磁V弯实验件的拉应变峰值比机械V弯试件低13.9%。同时,与机械V弯试件相比,电磁V弯试件有更大区域的金属材料参与了弯曲变形。