基于有限元分析的特征值反问题求解的逆摄动方法

本文研究特征值反问题的求解方法,根据广义特征值反问题理论和有限元法的特点,以转子系统平面梁单元有限元模型结构分析的特征值反问题求解为例,给出一种新的逆摄动方法,给出了本逆摄动法较完整的理论基础,给出了其逆摄动参数的显式计算公式及相应的取值方法,本逆摄动法也可推广到其他单元类型的有限元模型特征值反问题的求解。     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题,推导了相应的离散方程. 与 其它基于网格的方法相比,有限点法采用移动最小二乘法构造形函数,只需要节点信息,不 需要划分网格,用配点法离散控制方程,可以直接施加边界条件,不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

以Kirchhoff应力张量-Green应变张量表示本构关系的粘弹性大变形平面问题的有限元法

本文基于Updated Lagrangian增量叠加方法,采用以现时Kirchhoff应力增量和现时Green应变增量表示的虚功方程和Kirchhoff应力-Green应变的积分本构关系,导出粘弹性大变形的变分方程.依此采用八节点二维等参数元编制相应的平面问题的计算程序,算例的计算结果与以往的工作完全符合.     

结构后屈曲分析中的混合Newton-Lanczos算法

本文提出一种适于结构非线性后屈曲分析的混合Newton-Lanczos算法。与当前流行的弧长法不同,本文提出的算法采用传统的载荷增量法进行逐步求解,可求出给定载荷下的结构变形且适于任意外加载荷。对于临界载荷附近的迭代应用了Lanczos法求解方程及相应变载技巧。文中给出的若干数值计算结果表明了该算法在结构非线性后屈曲分析中的适用性。     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

粘弹性大变形动力响应的有限元分析

本文基于Ｔｏｔａｌ Ｌａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

弹塑性问题的边界元—流动因子线性互补解

1.前言尽管求解弹塑性问题已有许多方法,但一般来说,如不采用迭代法则不可避免地会产生有效应力偏离屈服面的“漂移”现象.从近代发展起来的变分不等式和参变量变分原理出发所建立起来的有限元方程,较成功地把数学规划理论应用到非线性本构关系的应力分析中.本文方法与以往方法的最大区别在于,将边界元方程与弹塑性屈服准则联立,导出了按增量求解的线性互补方程.因此对任一荷载增量步,通过求解一次线性互补方程,将使边界元方程和屈服准则同时满足,从而既无须迭代又避免了有效应力的漂移.并且也适用于非法向的屈服流动、随机强化模型及Drucker-Prager 模型等.与文献[1-4]的     

求解二点边值问题打靶法的一种改进方法

本文讨论两点边值问题的数值解法,利用最小二乘法修正打靶法所需初始参数,将边值问题转化为相应的初值问题与左梯度控制方程合并,引入精细积分法进行求解。通过数值算例表明,该方法是一种求解两点边值问题的有效方法。     

高超声速三维化学非平衡流动简化NS方程数值解

本文研究用时间稳定法求解高超声速三维化学非平衡流动简化ＮＳ方程，把整个非平衡流体动力学问题分裂为流体力学问题和化学问题。流体力学与化学问题之间的耦合通过迭代方式加以处理。考虑高温空气分解与电离，七个化学组元，并认为化学反应是以有限速率进行的．差分格式采用无自由参数的混合反扩散格式，激波及壁面处理应用特征线理论．对来流攻角为１０°，马赫数为２６，飞行高度为８５ｋｍ的钝锥体，给出了流场参数和物面加热率的数值计算结果．     

弹塑性有限元分析的有效数值方案

本文对弹塑性有限元增量分析的三个基本步骤(即线性化、方程求解、状态决定)分别提出了预测中点刚度和自动选择步长、顺序-逆序修正法求解方程、渐近积分弹塑性增量本构关系等新的数值方案,并在计算机上通过实际计算表明这些方案在计算精度和效率方面改进了现行的数值方案。     

剪切波作用下埋藏刚性椭圆柱与周围介质部分脱胶时的动力分析

本文利用波函数展开法和奇异积分方程技术研究了ＳＨ型反平面剪切波作用下埋藏刚性椭圆柱与周围介质部分脱胶时的动力特性．将脱胶区看作表面不相接触的椭圆弧形界面裂纹，利用波函数（Ｍａｔｈｉｅｕ函数）展开法，并引人裂纹面的位错密度函数为未知量，将问题归结为奇异积分方程，通过数值求解积分方程获得了远场和近场物理参量，并讨论了共振特性和各参数对共振的影响．     

基于Updated Lagrangian法的三维粘弹性大变形问题的有限元分析

本文基于Updated Lagrangian增量迭加方法,采用以现时Kirchhoff应力增量和现时Green应变增量表示的虚功方程和Kirchhoff应力张量-Green应变张量的积分本构关系,导出粘弹性大变形的虚功方程。依此采用二十节点三维立方等参数元编制相应的计算程序。三个算例结果与以往一维、二维的计算结果完全符合。     

具损伤压电智能层合中厚板的非线性动力稳定性分析

基于各向异性损伤理论和压电理论，采用能量法和Lagrange方程，建立了具损伤压电智能层合中厚板在参数激励下的非线性运动控制方程，应用增量谐波平衡法，求解了具损伤压电智能层合中厚板的非线性动力稳定性问题。算例中，分析了损伤参数和压电效应对压电智能层合板非线性动力稳定性的影响，揭示了力电耦合效应的内在特征，并与有关文献的结果进行了比较。     

基于单源模糊数的模糊随机动态有限元方程的解法

本文提出一种模糊随机动态有限元方程的解法,指出利用单源模糊数和它的运算法则,可以把一个不含阻尼项的模糊随机动态有限元平衡方程转化为两类不同集合下的方程组,一种是模糊数方程,另一种是普通的动态有限元平衡方程。前者可用模糊数运算法则求解。通常这类方程的表达式非常简单,故很容易求解,后者可利用现有的求解随机动态随机有限元平衡方程的方法计算,这时求解该方程的计算量几乎等同于求解相应的普通随机动态有限元平衡方程的计算量。最后的算例表明,本文提出的方法与通常所用的γ截集法计算结果基本相同,而且所用的计算量远远小于用γ截集法所用的计算量。     

微分方程反问题的梯度正则化法

本文在微分方程反问题研究当前发展状况的基础上,提出了一种求解微分方程反问题的比较统一的方法——梯度正则化法。该方法将反问题的求解分为两个主要部分,一是展开以求补充条件对未知参数的梯度矩阵,二是正则化以解不适定的线性方程组。梯度正则化法从反问题的共性入手,未附加任何特殊的约束,所以可以适用各种类型的反问题,且在求解时不受空间维数的限制。     

一般叠层圆柱厚壳的非线性动力稳定性分析

基于Ｔｉｍｏｓｈｅｎｋｏ－Ｍｉｎｄｌｉｎ假设及Ｈａｍｉｌｔｏｎ原理，建立了一般纤维叠层圆柱厚壳在参数激励下的非线性振动方程；应用多模态近似和增量谐波平衡法求解了叠层圆柱厚壳的非线性动力稳定性问题。横向剪切变形、端部支承条件等因素的影响被讨论。     

一般叠层圆柱厚壳的非线性动力稳定性分析

基于Ｔｉｍｏｓｈｅｎｋｏ－Ｍｉｎｄｌｉｎ假设及Ｈａｍｉｌｔｏｎ原理，建立了一般纤维叠层圆柱厚壳在参数激励下的非线性振动方程；应用多模态近似和增量谐波平衡法求解了叠层圆柱厚壳的非线性动力稳定性问题。横向剪切变形、端部支承条件等因素的影响被讨论。     

低速冲击载荷下加筋板弹塑动力响应分析

给出了一种半解析方法求解加筋板在横向低速冲击载荷下的非线性弹塑动力响应分析。考虑大变形的影响，采用增量型本构关系，得到膜力、弯矩增量与应变增量的关系，忽略面内惯性力的影响，运用Hamilton变分原理，得到加筋板增量形式动响应控制方程，最后采用增量形式的迦辽金法来求解。     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

轴向冲击载荷下圆柱壳的弹塑性屈曲分析

研究时变轴向冲击载荷作用下的圆柱壳弹塑性动力屈曲问题。本构关系采用增量理论 ,借助增量数值计算方法对Karman Donnell运动方程进行求解。计算表明 :基于B R准则的屈曲判断方法和采用Southwell方法可以获得一致的临界屈曲载荷 ;分别讨论了应力波对屈曲的影响以及材料参数、几何参数、载荷峰值与持续时间和动力屈曲的关系。