基于动态本构模型的整流罩铰链系统自冲击分析

整流罩铰链系统结构分析是大型整流罩分离分析的重要环节,也是国内新一代运载火箭研制过程中的关键技术之一。基于30CrMnSiA和ZTC4-1的力学性能试验研究,得到材料的动态本构模型,并建立整流罩及动态本构模型下的铰链系统有限元模型。通过非线性显式动力学分析方法,对铰链系统自冲击过程进行仿真,得到冲击作用下铰链系统的应力和应变情况,并通过内力分析得到自冲击过程中铰链系统的主要内力形式。最后,通过比较铰链系统动态本构与准静态本构下的自冲击响应,验证并说明了动态本构模型对铰链系统自冲击仿真的影响。计算结果表明,自冲击作用会使铰链系统发生塑性变形,且动态本构模型有利于反映真实结果,为新一代运载火箭的研制提供一定的参考价值。     

石墨烯碳纳米管复合结构渗透特性的分子动力学研究

采用经典分子动力学方法研究了压力驱动作用下水在石墨烯碳纳米管复合结构中的渗透特性.研究结果表明,水分子渗透通过石墨烯碳纳米管复合结构的渗透率明显高于石墨烯碳纳米管组合结构.水在石墨烯碳纳米管复合结构中的渗透率随着压强的升高而增大,随着电场强度的增大而减小.考虑了温度和复合结构中双碳管轴心距对水渗透性的影响规律.系统温度越高,水的渗透率越高;随着双碳管轴心距的增加,水的渗透率逐渐降低.通过计算分析水流沿渗透方向的能障分布,解释了各参数变化对水在石墨烯碳管复合结构中渗透特性的影响机理.研究结果将为基于石墨烯碳管复合结构的新型纳米水泵设计提供一定的理论依据.     

基于超二次曲面的颗粒材料缓冲性能离散元分析

自然界或工业中普遍是由非球形颗粒组成的复杂体系,与球形颗粒相比,非球形颗粒间的高离散和咬合互锁可使冲击载荷引起的能量有效衰减实现缓冲作用.基于连续函数包络的超二次曲面单元能准确地描述非球形颗粒的几何形态,并可精确地计算单元间的接触碰撞作用.本文采用离散元方法对冲击载荷作用下非球形颗粒物质的缓冲性能进行数值分析,并与圆柱体冲击的理论结果和球体冲击的实验结果进行对比验证.在此基础之上,进一步研究了筒底作用力在不同颗粒层厚度和形状等因素影响下的变化规律.计算结果表明:不同颗粒形状都存在一个临界厚度H_c.当HH_c时,缓冲率随H的增加而增加;当HH_c时,缓冲率的变化不再显著并趋于稳定值.此外,减小颗粒表面尖锐度和增加或减小圆柱形和长方形颗粒的长宽比都会提高颗粒材料的缓冲效果.     

考虑弹簧阻尼作动器解析雅可比矩阵的多刚体动力学分析

弹簧-阻尼-作动器(spring-damper-actuator,SDA)是多体系统中常见的力元,在工程领域中有着广泛的应用.采用绝对坐标方法建立的多体系统动力学控制方程通常是复杂的非线性微分-代数方程组.为了保证数值解的精度和稳定性,通常需要采用隐式算法求解动力学方程,而雅可比矩阵的计算在隐式数值求解过程中至关重要.对于含有SDA的多体系统,SDA造成的附加雅可比矩阵是与广义坐标和广义速度相关的高度非线性函数.目前的很多研究工作专注于广义力向量的计算,然而对附加雅克比矩阵的计算则少有关注.针对含SDA的多刚体系统进行动力学分析,首先基于Newmark算法研究其在动力学方程求解中的雅可比矩阵的构成形式;然后推导SDA的广义力向量对应的附加雅可比矩阵,其中包括广义力向量对广义坐标和对广义速度的偏导数矩阵.最后通过两个数值算例研究附加雅可比矩阵对动力学分析收敛性的影响;数值分析表明:当SDA的刚度、阻尼和作动力数值较大时,SDA导致的附加雅可比矩阵对数值解的收敛性有重要影响;当考虑SDA对应的附加雅可比矩阵时,动力学分析可以以较少的迭代步实现收敛,从而减少分析时间.     

自适应一致性高阶无单元伽辽金法

近来提出的一致性高阶无单元伽辽金法通过导数修正技术大幅度减少了所需积分点数目,并能够精确地通过线性和二次分片试验,显著改善标准无单元伽辽金法的计算效率、精度和收敛性.本文在此基础之上,充分利用无单元法易于在局部区域添加节点的优势,发展了一致性高阶无单元伽辽金法的h型自适应分析方法.根据应变能密度梯度该方法自适应地确定需节点加密的区域,基于背景积分网格的局部多层细化要求生成新的计算节点,同时考虑了节点分布由密到疏渐进过渡的情形.采用相邻两次计算的应变能的相对误差作为自适应过程的停止准则,将所发展自适应无网格法应用于由几何外形、边界外载和体力等因素造成的应力集中问题的计算分析.数值结果表明,所发展方法能够自适应地对高应力梯度区域进行节点加密,自动给出合理的计算节点分布.与已有的标准无网格法的自适应分析相比,所发展方法在计算效率、精度和应力场光滑性等方面均展现出显著优势.与采用节点均匀分布的一致性高阶无单元伽辽金法相比,它大幅度地减少了计算节点数目,有效提高了一致性高阶无单元伽辽金法在分析应力集中等存在局部高梯度问题时的计算效率和求解精度.     

裂纹问题的一致性高阶无网格法

一致性高阶无网格法能高效精确地求解连续体问题,尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题（即非连续体问题）的数值分析。采用背景积分网格描述裂纹几何,基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断,提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展,采用相互作用积分方法计算应力强度因子,裂纹的扩展方向由最大周向应力准则确定。数值结果表明,本文发展方法能够精确地通过间断分片试验;相较于标准的高阶无网格法和低阶一致性无网格法,本文的一致性高阶无网格法显著改善了应力强度因子的计算精度,能够准确预测裂纹扩展路径。     

核主泵叶轮间隙的轴向振动特性研究

具有径向流的间隙结构广泛存在于轴承结构与旋转机械中,间隙中的流固耦合作用可能影响整体结构的运动稳定性。基于理论间隙模型和核主泵的实际结构,本文对径向间隙流引起的轴向振动进行了多方面的研究。当间隙的一个壁面产生轴向振动并处于旋转状态时,壁面受到由径向流引起的时变轴向力,因此间隙为叶轮提供附加的轴向刚度和阻尼。通过研究以水为介质的理论间隙模型,发现径向间隙流会引起负的等效轴向动力系数（刚度和阻尼）,并且流道形状是影响间隙轴向动力特性的重要因素。扩张流道和平行流道会产生负的轴向动力系数,特别是负阻尼会引起结构振动发散;而收缩流道间隙具有稳定的轴向动力特性。最后,对AP1000核主泵原型叶轮间隙模型进行分析,结果表明,间隙会引起轴向负刚度,并且在一定工况下出现负阻尼,此时系统轴向稳定性及结构安全运行将受到严重的不良影响。     

一种考虑瞬态效应的散热结构导热路径设计的拓扑优化模型

实际工程中,热载荷多数具有短时和周期性特点,瞬态效应显著。目前的散热结构导热路径设计多基于稳态热传导模型,未考虑瞬态效应。本文提出了一种以区域温度控制函数作为设计目标的瞬态热传导问题的拓扑优化模型,能够实现在整个时间历程上特定区域内最高温度最小。使用伴随变量法,推导了目标函数关于设计变量的敏度计算格式。算例表明,基于本文优化模型获得的散热路径设计与基于稳态热传导模型的结果有明显差别,具有更优的散热性能。因此,时变热荷载下的散热结构构型设计需要考虑瞬态响应的影响。     

铁路道床有砟-无砟过渡段动力特性的DEM-FEM耦合分析

针对铁路道床有砟-无砟过渡段的结构特点,采用离散元-有限元耦合模型分析散体道砟和无砟道床间过渡段的动力特性。散体道砟道床和无砟道床分别采用离散元方法 DEM和有限元方法 FEM模拟,而在过渡段将道砟颗粒嵌入无砟道床以增加道砟颗粒与无砟道床间的咬合力,并在离散元和有限元耦合区域实现了力学参数的传递。采用以上DEM-FEM耦合方法对有砟-无砟道床及其过渡段在列车荷载作用下的沉降过程进行了数值分析。计算结果表明,离散元方法中道砟颗粒间的力链呈现非对称梯形分布,其与有限元方法中的应力分布趋势一致;采用嵌入式道砟颗粒的方法可以增加有砟-无砟过渡段道砟间的咬合力,有效约束道砟颗粒的位移,减少有砟-无砟道床间的沉降差异。本文计算模型可以合理地分析有砟道床的力链分布以及无砟道床的应力分布,确定列车荷载下道床有砟-无砟过渡段的动力学行为。     

基于坐标映射的类周期电磁金属阵列布局设计的灵敏度分析方法

特定性能的类周期电磁金属阵列已成功应用于多种类型的电磁器件中,以实现其电磁辐射或散射性能的提升设计.类周期电磁金属阵列可描述为由多个相似形状的电磁金属单元按照特定的布局形式组合形成.布局形式由阵列的布局参数确定,通常包括每个单元的尺寸、转角以及位置参数.布局参数往往规模较大且相互独立,需要经过设计优化以满足阵列的性能要...