求解无限域动力刚度矩阵的双渐近算法

采用连分式算法可以有效地求解无限域动力刚度表示的比例边界有限元方程, 它具有收敛范围广、收敛速度快等优点. 该文在高频渐近连分式算法的基础上考虑了低频渐近, 发展了一种针对矢量波动方程的双渐近算法. 随着展开阶数的增加, 双渐近算法可以在全频域范围内快速逼近准确解. 引入了系数矩阵?X(i)来增强连分式算法的数值稳定性. 通过在高频极限、低频极限时满足动力刚度表示的比例边界有限元方程, 建立了递推关系以求得动力刚度矩阵. 通过二维半无限楔形体、三维均质弹性半空间数值算例表明, 双渐近算法比单渐近算法更稳定、优越.     

基于动力刚度法的体外预应力梁自振频率分析

用预应力梁-杆组合结构模拟体外预应力梁,以体外预应力简支梁为例,建立了预应力梁杆单元动力刚度矩阵。采用动力刚度法,进一步推导出预应力梁-杆组合结构的整体动力刚度矩阵,利用Williams-Wittrick算法求解频率。本文以理论推导为基础,引入了动力刚度法。最后通过算例,讨论了各种因素对梁横向的振动特性的影响,并与试验值比较。计算结果表明:动力刚度法能够精确有效的求解体外预应力混凝土梁的横向振动问题。     

拉压刚度不同桁架的动力参变量变分原理和保辛算法

对于一些展开结构,为达到其设计性能,必须采用特殊的索、膜结构,这些索、膜部件表现出不同的拉压性质。具有拉、压不同性质的材料或结构的力学分析,体现出较强的非线性特征,需要针对这类问题发展有效的求解算法。本文建立了由拉压刚度不同杆单元组成的桁架结构的动力学参变量变分原理,将拉压刚度不同桁架问题的非线性动力分析转换为线性互补问题求解。结合时间有限元方法构造了求解此问题的保辛数值积分方法,此方法不需要迭代和刚度矩阵更新,避免了迭代求解方法的收敛问题,计算过程稳定、高效。     

阻尼矩阵与刚度矩阵的一种直接修正方法

提出了一种利用复模态测量数据同时修正有限元阻尼与刚度矩阵的有效方法。借助于矩阵的奇异值分解得到了满足动力方程的最小修正矩阵。该方法有一个简洁的表达式,修正过程简单而且容易实现,数值试验表明修改后的结构参数能准确地同试验值吻合。     

基于温度体模型的动网格生成方法及在流固耦合振动中的应用

提出温度体模型动网格生成方法,并将其应用于流固耦合算法,生成物体振动过程中的动态网格。温度体模型动网格方法将运动边界的位移映射为求解域的温度边界条件,以流体能量方程或固体导热方程作为控制方程求解得到计算域内的温度分布,将求解得到的温度分布作为网格节点的动态位移。基于温度体模型的动网格方法物理意义清晰,算法实现简单,能够快速而有效地生成高品质的动态网格,特别在边界位移大的情况下与其他网格生成方法相比有较大的优势。最后采用流固耦合有限元算法求解了定浆式轴流泵强迫振动过程中连锁特性和柱体由于旋涡脱落诱发自激振动两个工程问题。其中流场采用基于特征线方程的分离算法进行求解,固体场采用Newmark方法进行求解,在计算过程中采用温度体模型生成动态网格。结果表明该发展的算法性能优异,能有效解决流固耦合中的振动问题。     

地震作用下水-轴对称柱体相互作用的子结构分析方法

针对水-轴对称柱体动力相互作用问题,提出了一种地震作用下水-结构相互作用的时域子结构分析方法.基于三维不可压缩水体的波动方程和边界条件,利用分离变量法将其转换为环向解析、竖向和径向数值的二维模型;基于比例边界有限元推导了截断边界处无限域水体的动力刚度方程,并将水体内域有限元方程和人工边界处的动水压力进行耦合,从而得到结...     

模拟P波向无穷远域传播的一致边界

将人工边界设置在半无穷层单元和内部有限元区域的交界面上,建立了半无穷层单元的刚度矩阵后,得到了边界节点的动力平衡方程。任意给定激励圆频率,将边界节点系统的动力平衡方程转化为特征值方程。求解特征值方程得出边界节点系统的特征值和特征模态,利用模态叠加原理得到体现左半无穷层单元和右半无穷层单元对内部有限元区域作用的边界矩阵,这就是该文的一致边界。将其与内部有限元区域的刚度矩阵进行组装来模拟无穷远域介质对波的传播作用。最后用数值算例来说明一致边界的精确性和可行性。     

基于连分式与有限元法的坝库耦合瞬态分析方法

为高效模拟地震激励下坝库耦合瞬态响应,建立了无限水库的连分式与有限元法的耦合公式。结合坝体有限元公式,利用坝库耦合项,发展了坝库耦合瞬态分析迭代算法。利用该算法分析了水平向地震激励下重力坝的瞬态响应。比较了基于连分式法、动态刚度矩阵法、动态质量矩阵法模拟坝库耦合问题的计算效率。数值算例表明该耦合算法模拟坝库耦合瞬态响应的正确性及高效性。该方法继承了比例边界有限元法的精度高、离散单元少等特点,又避免了其卷积积分,提升其计算效率,为坝库耦合瞬态响应提供了一种高效分析方法。     

三维冲击-动力接触问题数值分析的一种高效算法

本文要二维冲击－动力接触问题^[3]和三维弹性接触有限元混合法的基础上,给出了三维冲击－动力接触问题的有限元混合算法,该方法将整个系统的等效刚度议程凝缩到接触边界,建立关于接触面的等效柔度方程,在采用Newmark直接积分法求解动力接触问题时必须作适当修正,本文给出了修正方法,由于用直接积分法进行每一时间步动力接触爱代计算时,只震反复求解接触面有效柔度议程,因此具有很高的计算效率。     

基于连分式与有限元法的坝库耦合瞬态分析方法EI北大核心CSCD

为高效模拟地震激励下坝库耦合瞬态响应,建立了无限水库的连分式与有限元法的耦合公式。结合坝体有限元公式,利用坝库耦合项,发展了坝库耦合瞬态分析迭代算法。利用该算法分析了水平向地震激励下重力坝的瞬态响应。比较了基于连分式法、动态刚度矩阵法、动态质量矩阵法模拟坝库耦合问题的计算效率。数值算例表明该耦合算法模拟坝库耦合瞬态响应的正确性及高效性。该方法继承了比例边界有限元法的精度高、离散单元少等特点,又避免了其卷积积分,提升其计算效率,为坝库耦合瞬态响应提供了一种高效分析方法。     

考虑索间作用的斜拉桥耦合振动研究

为了研究几何非线性条件下斜拉桥索梁耦合振动与索间作用问题,以两条斜拉索与简支梁组合体系为简化模型,利用D’Alembert原理建立考虑初始垂度的索梁体系非线性偏微分方程,设定索的前两阶复合振动模态与梁的基本模态,运用Galerkin方法将其离散为二阶常微分方程,并使用四阶—五阶Runge-Kutta方法对索与梁的振动响应进行了数值分析。结果表明:在双索单梁组合结构中,特定频率条件下一阶模态与主梁强烈耦合,二阶模态与主梁小程度耦合;与单梁单索结构相比,多索导致主梁频率增大,索间作用使得索振幅增大、拍频降低,面内一阶模态对索梁变化更敏感;当索梁频率不变时,索间作用对耦合振动产生的索大幅振动有明显抑制作用,且索梁结构对主梁初位移变化更敏感。     

不同冲击速度条件下矩形巷道围岩变形-开裂过程数值模拟

利用自主开发的拉格朗日元与变形体离散元耦合的连续-非连续方法,开展了不同冲击速度条件下矩形巷道围岩变形-开裂过程的数值模拟研究。当岩石单元的应力满足剪切开裂判据时,考虑了应力脆性跌落效应,在此过程中,围压保持不变。为了检验该方法的正确性,开展了单轴压缩条件下岩样变形-开裂过程的数值模拟研究。针对矩形巷道围岩的计算表明:①冲击速度低时载荷-位移曲线呈单峰特点;冲击速度高时载荷-位移曲线呈多峰特点;②冲击速度低时围岩的开裂呈渐进性特点;冲击速度高时围岩的开裂呈间歇性特点,这与在冲击过程中围岩中过去形成的一种较为稳定的结构(拱)被打破,并在其外围形成了新的稳定结构(拱)有关;③在模型相同垂直方向位移时,冲击速度低时围岩两帮的开裂深度较大,这说明冲击速度低时应力的传递较为均匀,冲击速度低时的结果(例如,V形坑内岩石碎块涌入巷道)与围岩的片帮类似;冲击速度高时的结果(例如,岩石碎块的弹射现象)与岩爆类似。     

基于多重分形分析的核爆与雷电电磁脉冲识别

根据核爆和雷电电磁脉冲信号非平稳、非线性特点,对NMEP和LEMP信号进行了多重分形分析,计算了二者的多重分形广义维数 和广义维数的变化率,通过分析发现,广义维数的变化率更能反映核爆和雷电的不同产生机理,并且和广义维数相比,以广义维数的变化率作为特征,对NEMP和LEMP的识别率更高,更稳定,可以对核爆和雷电电磁脉冲进行有效的识别。     

基于卡箍优化布局的飞机液压管路减振分析

飞机液压管路振动过大是一个亟待解决的问题,考虑了主动消振方式较为复杂,难以在飞机液压系统实现,提出了基于系统特征阻抗,通过优化卡箍布局的被动消振方式。首先建立系统管道及各元件数学模型,利用传递矩阵法（Transfer Matrix Method, TMM）得到其频域特征阻抗变化规律。以激振源固有频率点的特征阻抗加权和为目标函数,利用粒子群优化算法（Particle Swarm Optimization, PSO）在一定范围内调整卡箍位置,使在激振源频率点的特征阻抗加权和降到最低。通过优化卡箍位置,系统特征阻抗加权和较优化前衰减28.13%,验证了其有效性,为飞机液压管路的优化布局提供了一定的理论依据。     

齿轮非线性动力系统的振动功率流分析

针对一对含侧隙、具有时变啮合刚度的非线性齿轮系统动力学模型,通过对该系统的振源功率进行分析,推导了系统振动功率流的时域仿真算法,探讨了决定系统功率流水平的主要因素。并综合运用点映射和胞映射方法,分析比较了系统对不同参数连续变化的分岔图和功率流图谱,研究了参数变化对功率流的影响,以期为将功率流理论应用于非线性齿轮啮合系统提供参考。     

结构减震的稳定时滞LSSVM-LQR智能控制算法

针对时滞LSSVM-LQR智能控制算法存在的稳定性问题,提出相关的稳定性控制算法,以确保时滞LSSVM-LQR智能控制算法的鲁棒性。该算法的主要思路为：在时滞LSSVM-LQR控制算法中,加入控制力限制条件。当满足控制力限制条件时,控制程序继续运行;当不满足控制力限制条件时,控制程序自动跳出,便执行稳定性控制算法（或称为稳定/鲁棒的时滞LSSVM-LQR智能控制算法）。稳定性控制算法主要是通过调整反馈来控制作动器运行,从而确保控制系统的稳定性。数值结果表明,稳定性控制算法能够有效地保证时滞LSSVM-LQR智能控制算法的稳定性/鲁棒性;与时滞LSSVM-LQR智能控制算法相辅相成。     

上三角肘杆式压力机主驱动机构的运动学分析

针对压力机主驱动机构的运动特性,提出了一种上三角肘杆式压力机主驱动机构。根据该机构的工作原理,建立了上三角肘杆式主驱动机构的工作原理示意图和仿真模型,根据各连杆的几何关系,运用机械系统运动学分析软件ADAMS依次对连杆1~7进行了运动学仿真计算,得到了肘杆机构的所有杆件在不同杆长条件下对应的位移、速度和加速度的变化曲线。获得了上三角肘杆机构各个杆长的运动特征,分析了每个连杆对主驱动机构工作状态的影响,揭示了上三角肘杆式压力机主驱动机构在工作过程中的运动规律,为实际工况下确定机构参数提供理论依据,为后续更深入和详细的进行肘杆机构的研究奠定了一定基础。     

非线性转子-轴承系统的动力学降维分析与试验研究

应用有限元法建立转子-轴承试验系统的非线性动力学模型,采用固定界面模态综合降维法将原高维系统转换为低维少自由度系统,采用Newmark-β法对降维后模型进行求解,在和全自由度模型对比满意的前提下,得出转子系统的三维谱图、分岔图、三维振幅图、轴心轨迹图以及Poincare截面图,并和试验结果进行对比,结果表明本文所建立的动力学模型,较为真实地反映了试验系统的非线性特性,为计算复杂转子-轴承系统的深层次动态设计提供理论依据。     

侧铣加工动力学建模及仿真

本文以两自由度的铣削加工为研究对象,建立了闭环的动态力～振动切削动力学仿真模型;应用模型仿真切削力与振动并与实测数据进行比较;证实该模型仿真切削力与振动的准确性。该结果对预测铣削力与振动,选择合适的切削参数具有实际意义。     

舰船低频水下辐射噪声的声固耦合数值计算方法

针对舰船低频域水下辐射噪声计算问题,指出采用严格遵循声固耦合动力学方程的耦合声学有限元与远场自动匹配层(FEM/AML)方法以及耦合声学间接边界元(IBEM)方法是计算精度较高的策略。以某小水线面双体船(SWATH)为研究对象,使用声功率作为评价指标,探讨了声场区域特征尺度选取对计算精度的影响,比较了上述两种方法与常规基于流固耦合的两种方法在计算特性方面的差异。研究表明,声固耦合模式较流固耦合模式声学响应计算结果偏小,对于SWATH船的合成总声功率级两者偏差达到1dB~3dB,前者计算结果更为精确,基于声固耦合模式的耦合声学IBEM方法是舰船水下辐射噪声预报的首选算法。