基于摄动原理的复杂土层地震反应分析的子结构法

把约束子结构模态综合法与直接模态摄动法相结合,建立复杂场地三维地震反应等效线性化分析计算方法.应用直接模态摄动原理,可简化各子结构模态分析过程,将特征值求解问题转化为线性代数方程组的求解,从而可有效提高计算效率.算例表明,该方法在提高大规模复杂场地地震反应分析计算效率方面优势明显.     

轴对称压电圆板的非线性动力响应

基于Von Karman板理论和压电材料力学,考虑横向剪切变形,建立了轴对称压电圆板的非线性运动方程,提出了相应的力学与电学边界条件.求解时,首先应用Galerkin方法,将非线性偏微分运动方程转化为仅含时间变量的非线性常微分方程.然后,应用Newmark-β方法将时间函数离散,整个问题应用Newtoni迭代法求解.算例中,求得了压电圆板线性振动基频,验证了方程和求解方法的可靠性,讨论了压电效应、几何非线性、结构尺寸、力学和电学荷载等因素对板非线性动力响应的影响.     

基于小波包技术的爆破地震效应计算模型及安全判据研究

采用小波包分解和重构的方法,将复杂的实测爆破地震波速度信号转化为多个简谐波的叠加,将 爆破速度荷载作用下结构的动态响应问题转化为一系列简谐荷载作用下的动态响应问题。通过提出速度因 子的新概念并反映在爆破地震效应计算模型中,考虑了结构动态响应中瞬态响应的影响。同时爆破地震效应 计算模型中又融入了归一化的能量比例,在考虑爆破荷载频率的影响时仅需考虑占有相当能量比例的优势 频率的综合作用。然后,在该计算模型基础上提出了一个新的爆破地震效应安全判据。该判据能反映出爆破 激励荷载作用下结构速度响应大小与结构特性、爆破荷载幅值、频率(包括多个优势频率)、持续时间及能量比 例等参数的关系。最后结合实际工程案例,通过使用基于小波包技术的爆破地震效应计算模型与时程分析 法,分别求解出速度响应幅值并将结果进行对比,验证了所建模型的可行性。     

一种求解局部非线性结构稳态响应的方法

为了降低求解局部非线性结构稳态响应的计算量,基于子结构和阻抗缩聚提出了一种用于求解局部非线性结构稳态响应的计算方法.将局部非线性结构分解为线性子结构和非线性子结构,利用谐波平衡构造各个子结构的阻抗方程,对线性子结构进行缩聚,将局部非线性动力学方程转化为求解一组非线性代数方程组问题,通过迭代求解非线性代数方程组,求解系统的稳态响应.     

非线性轨迹优化问题的保辛自适应求解方法

非线性轨迹优化问题一般是一个非线性最优控制问题。将非线性系统的最优控制问题导入到哈密顿体系的辛几何空间当中,基于对偶变量变分原理提出了求解非线性最优控制问题的一种保辛自适应方法。以时间区段两端协态作为独立变量,在时间区段内采用拉格朗日插值近似状态和协态变量,并利用对偶变量变分原理将非线性最优控制问题转化为非线性方程组的求解,保持了哈密顿系统的辛几何结构。并进一步,提出了基于多层次迭代的自适应算法,提高了非线性最优控制问题的求解效率。数值实验验证了该算法在求解非线性轨迹优化问题中的有效性。     

单层平面索网幕墙结构的几何非线性问题研究

单层平面索网支承式玻璃幕墙结构是近年来在国内外应用较为广泛的一种新型幕墙结构形式,由于单层平面索网不具有负高斯曲面形式,结构在平面外方向的刚度偏柔,表现出较明显的几何非线性特征.本文采用连续化方法建立了单层平面索网结构考虑几何非线性影响的静力平衡方程和振动方程,得到了结构刚度的解析表达式,并采用谐波平衡法求得非线性频率的简化解析表达式,以此为基础研究了单层平面索网结构的静力非线性和动力非线性问题.研究结果表明:结构的非线性和结构的初始位置密切相关;结构的非线性频率主要取决于索的初始应变、结构振动幅值与跨度的比值,几何非线性对于结构动力性能的影响要小于对结构静力性能的影响;本文得到的结构在地震荷载和平均风荷载作用下的非线性振动方程和非线性频率为结构在地震荷载和脉动风荷载作用下动力响应的求解奠定了基础.     

薄板几何非线性弯曲分析的深度能量法

发展了一种增量形式的深度能量法求解薄板几何非线性弯曲问题。根据最小势能原理和Von-Karman非线性理论,构建以薄板势能为驱动的增量式深度神经网络模型。首先用网格离散薄板求解域,通过Python读取网格数据计算Hammer积分点,并以此作为训练集代入网络模型预测板的弯曲位移,再将荷载分成一系列的荷载增量,每个增量步中计算薄板势能作为神经网络的损失函数,以最小化势能为目标,结合Adam优化算法更新网络模型参数,待势能取驻值后再继续下一个荷载步的计算。本文求解了不同形状、不同边界条件下薄板的几何非线性弯曲问题,并将计算结果与文献解或有限元Abaqus解进行对比,研究表明,本文方法在求解薄板的几何非线性弯曲问题上具备有效性和准确性,且增量式的神经网络模型能够减小计算内存,有效提高计算效率和模型的稳定性。     

线性随机结构的非平稳随机响应变异性分析

对于具有随机参数的结构受到非平稳随机激励的问题，给出了结构随机响应变异系数的虚拟激励摄动算法。它应用虚拟激励法先将随机荷载转化为确定性荷载，以使随机问题精确地转化为仅结构参数具有随机性的问题，从而将问题归结为应用随机摄动法求解单随机问题。求解过程简单高效，且有较高的精度。     

转换位移自由度法求解滑动接触问题

有限元分析中,一般都是以节点沿坐标轴方向的位移组成位移自由度空间,以单个节点上沿坐标轴方向的位移和等效节点荷载作为边界条件。而接触问题的位移边界条件,一般都采用多个节点位移线性组合的形式,其相应的荷载边界条件则为接触面上的法向接触力和切向摩擦力。本文从广义自由度的概念出发,通过对总刚矩阵进行处理,直接用接触点的位移自由度的组合表征相关节点沿法向和切向的位移自由度,其荷载边界条件亦相应转化为切向和法向荷载,从而可以方便地输入接触边界条件,完成接触分析。采用本文提到的方法,不必增加位移自由度,并可以保证总刚矩阵的对称性和各接触条件的独立性。试验证明,本方法是正确有效的。     

结构非线性随机地震反应的概率密度演化方法

运用随机过程的正交展开方法,将地震动加速度过程展开为由少量独立随机变量所调制的确定性函数的线性组合形式.结合概率密度演化理论,建立了结构非线性随机地震反应分析的密度演化方法.以滞回结构非线性反应分析为例,考察了一个具有10个自由度的剪切型框架结构的随机地震反应分析问题.分析表明,结构非线性地震反应具有明显的随机涨落现象,且概率密度曲线与常见的正态等分布相差甚远.     

旋转弹性体的有限元线法分析

本文将有限元线法应用于一般荷载作用下的旋转弹性体的分析。文中将任意荷载沿环向展开为Ｆｏｕｒｉｅｒ级数，利用正交性，将问题转变为一系列旋转子午面上二维问题的叠加。文中对任一环向谐波建立了旋转面上的曲线曲边单元，导出了相应的常微分方程体系，并对由柱坐标引起的ｒ＝０处的奇异性问题做了完备合理的处理。文中给出了具有代表性的数值算例，用以展示本法的出色的表现     

基于框-剪结构协同工作的抗振墙数量优化研究

针对当前抗震墙数量优化中多是定性分析的现状,根据框架-剪力墙结构整体受力特点,采用横向分类方法,将空间框-剪结构分解量化为榀结构形式。基于有限元思想及框剪结构空间协同分析的离散化原理建立了其力学分析模型及剪力墙优化模型,该力学分析模型能较好地反映结构整体工作性能的主要特征及考虑剪力墙剪切变形的影响,动力分析采用振型分解法计算地震作用,以遗传算法作为优化方法确定出了剪力墙的最优数量。计算结果表明该方法具有较高的精度,可为结构的设计提供参考。     

爆破地震波三要素对多层砌体结构弹塑性地震响应的影响

利用Matlab编写了一套Newmark时程分析法与刚度退化二线型恢复力模型相结合求解多层砌体结构体系弹塑性地震响应的程序.以一个具体的4层砌体结构为分析对象,输入荷载用实测及人工模拟爆破地震波作用的形式,并按某原则变换输入爆破地震波的特性参数,根据所求出的结构各弹塑性地震响应幅值结果,分别讨论了爆破地震波幅值、主频和...     

水下爆炸作用下基于声学的爆源子结构输入方法

无限域吸收边界和爆炸波输入是水下爆炸载荷数值模拟的两个关键问题.本文借鉴基于内部子结构的地震波动输入方法和爆源子结构多尺度分析方法,考虑水下爆炸载荷与地震载荷同属于波动问题,提出一种水下爆炸作用下爆源子结构的爆炸波输入方法,该方法首先通过对爆源区域进行水下爆炸波自由场的波场分解,并利用地震波的内部子结构输入方法,将该自由波场运动转换为等效爆炸波输入载荷,从而实现了水下爆炸问题中冲击波的输入;本文采用圆形爆源子结构区域,并采用AUTODYN软件中一维自由场模型计算水下爆炸作用下该区域的自由波场压力时程.进一步,基于连分式近似方法提出了一种模拟无限域水体辐射效应的高精度时域人工边界条件,该吸收边界可设置于离结构和爆源子结构较近处.本文提出的分析方法通过圆形爆源子结构进行水下爆炸载荷的转化,并采用高精度吸收边界大大减少计算区域,既保证了计算精度,又降低了单元数量,具有较高的计算效率和较强的实用性.通过数值算例分析,验证了本文模型与方法的准确性,模拟了水下爆炸作用下的冲击波和气泡脉动阶段测点的压力时程曲线,并研究了圆形结构对水下近场爆炸波散射效应的影响规律.     

小当量地下爆炸的质点速度模型

在第四纪砂砾层中进行了系列公斤至百公斤级的地下封闭爆炸,研究了小当量化学爆炸地震波传播规律。在近场,地震波持时很短,水平振幅强于垂向振幅。质点速度随当量呈现指数增加,水平向指数为1.09,垂直向指数为0.77,质点速度随距离呈现指数衰减,水平向和垂直向指数分别为2.07与1.57。Sadauskas模型、显函模型和双极模型都能定量描述小当量地下爆炸地震波质点速度变化,但他们之间的反演精度存在差异,双极模型的残差最小。换言之,采用双极模型反演的数据更接近实际。     

土-结构动力相互作用分析中基于人工边界子结构的地震波动输入方法

数值模拟是解决土-结构动力相互作用问题的重要手段,而合理地实现地震波动输入直接影响地震作用下土-结构动力相互作用问题数值模拟的精度。波动法是目前常用的地震动输入方法之一,该方法将输入地震动转化为人工边界上的等效荷载,相较于其他地震动输入方法,波动法模拟精度高,但实施上相对复杂。从有限元模型入手,推导了采用波动法确定等效输入地震荷载的另一种形式,以此为基础,提出了一种在人工边界上实现地震动输入的新方法。新方法通过对土-结构有限元模型中由包含人工边界节点的单元组成的子结构施加自由场位移时程并进行动力分析,直接获得可实现地震波动有效输入的等效荷载,然后将等效输入地震荷载施加在土-结构模型的人工边界节点上,从而完成土-结构动力相互作用问题分析中地震动输入和地震反应计算。与原有波动法相比,新方法避免了需分别计算人工边界上自由场应力和由引入人工边界条件带来的附加力,以及需要根据不同人工边界面确定荷载的作用方向等较为复杂的处理过程,具有等效地震荷载计算简便、地震动输入过程更易于实施的特点。采用竖直入射和斜入射地震波动作用下的弹性半空间和成层半空间地震反应算例验证了新方法的有效性。      

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.      

一种耗能剪力墙结构地震反应的简化计算方法

本文用哈密顿原理建立了一种带竖缝耗能剪力墙结构在地震作用下的动力平衡方程，用Ritz法求解该方程，采用基于割线刚度和阻尼比等概率 幅值平均的等效线性方法求解竖缝中的耗能装置的刚度和阻尼比，并通过对一耗能剪力墙结构模型振动台试验的计算分析表明了该方法计算简单，且有较好的计算精度，为便于工程应用，对由耗能装置提供给结构的附加阻尼进行强行解耦，用振型分解反应谱的方法进行耗能结构地震反应的简化计算，并通过对一算例的计算分析表明了该简化方法具有一定的精度，该简化方法使对耗能结构的抗震计算与普通结构比较一致，该方法主要适用于小震作用下的耗能结构的抗震计算。     

具有稳定数值解的三维谐振子

谐振子广泛应用于物理系统的描述和物理现象的数值模拟。由于二维或三维谐振子对于系统参数、初始条件和边界条件的高度敏感性,很多物理过程的动力学模拟都会出现数值解不稳定的现象。近年来发展的无网格法、物质点法和近场动力学法等数值模拟方法均绕开了对固体材料固有构形的量化描述。本文引入了定常耗散项和弹簧耗散项,考虑随机微扰效应,提出了一种三维耗散谐振子,构建了基于蛙跳法和边界松弛技术的数值积分算法。应用三维谐振子构建了耗散型弹簧摆、简化弦和简化梁三个模型,设定了13个定解问题进行动力学模拟。数值试验结果表明,三维谐振子是稳定的。基于简化弦模型,模拟了拉弦、放弦和重弦三个有界弦振动问题;其中,拉弦和放弦问题成功模拟了有界弦的三维振形;重弦问题模拟再现了悬链线在水平向的微幅振荡现象。基于简化梁模型,模拟了三维梁的拉伸、剪切和扭转行为,验证了三维谐振子对于非线性大变形问题动力学模拟的描述能力,及其对外部作用的高速响应能力。本文方法可以为弦振动问题和材料力学非线性大变形问题的动力学模拟提供一条可行的实现途径。     

SH波对直角平面区域内圆形衬砌的散射与地震动

研究了以任意方向入射的平面SH 波对直角平面区域内圆形衬砌的散射与地震动问题．利用复变函数法、多极坐标移动技术及镜像法,将直角平面区域内的波场延拓于半无限空间,根据边界条件将该问题的解答可归结为对一组无穷代数方程组的求解问题,通过算例考虑衬砌的动应力集中和地表位移．结果表明,衬砌的动应力集中和地表位移幅值取决于入射波频率、角度、衬砌的位置和厚度比．