变厚度圆柱壳的轴对称弯曲问题

本文详细地研究了厚度h=h_0ξ~的圆柱壳的轴对称弯曲问题.文中通过引入一个位移函数H(ξ),将该问题的方程组化成一个关于H(ξ)的6阶常微分方程,用广义超几何函数给出问题的精确解.     

结构动力响应的时间有限元全域算法

利用状态方程将二阶结构动力学方程组变成一阶线性微分方程组,并采用加权残值伽辽金法将一阶线性微分方程组离散成线性方程组.该方法是一种全域算法,是真正意义上的时间有限元算法.数值算例表明:该方法有很好的计算精度,与解析解吻合较好.     

基于行消元回代法的多域边界元分析方法

提出了一种用多域边界元技术求解大型工程问题的新算法. 首先, 采用三步变量凝聚技术, 将由内部点、边界点和公共结点表述的每一子域的基本边界元代数方程表述成只有公共结点变量为未知量的代数方程, 然后, 根据公共结点的平衡方程和协调条件组集具有稀疏系数特征的总体系统方程组. 为了有效求解该系统方程组, 首次在边界元法中引进一种能有效求解大型非对称稀疏系数矩阵方程组的行消元回代法(REBSM), 该方法可在方程的每一行组形成时进行消元和回代, 当方程组组集完毕后即可得到方程的解, 不需要最后的回代过程. 因为一些项的重复计算在每一行的处理中合并掉, 因此REBSM要比传统的高斯消元法需要较少的内存, 而且计算速度具有数量级的提高, 可为边界元法求解大型工程问题提供有力的方程求解器.      

凸肩叶片的非线性振动特性与运动分岔

以凸肩叶片作为研究模型, 建立了考虑凸肩摩擦力, 几何大变形与阻尼的非线性振动方程.采用Galerkin法对振动方程离散化, 应用平均法对离散后模态方程组的非线性响应进行解析分析, 得到了非线性幅频特性曲线, 与数值解比较验证了解析解, 并讨论了系统周期解的稳定性. 用非线性振动理论详细研究了平均方程组的运动分岔现象, 揭示了平均方程组周期解的变化过程及其具有的非线性动力学性质. 解析结果表明, 凸肩之间的摩擦对系统第二阶非线性振动特性影响很大. 由于凸肩之间摩擦力方向的不断改变, 系统第二阶非线性幅频特性曲线不连续, 出现两个共振频域. 随着时间的推移, 系统振动的幅值会以$T/ 4$为周期在两个频域的幅频曲线上来回跳动, 这会使叶片的振动响应大幅降低.      

奇异Hamilton系统矩阵的精细积分法

常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点.     

处理一种非线性问题迭代法的计算机方法

本文依非线性问题的积分方程组,按迭代法得出一般迭代表达式.由此,任意次迭代解的所有待定系数可用计算机计算,使研究复杂的非路性问题成为可能.     

方腔内自然对流与混合对流的有限差分格子Boltzmann研究

基于有限差分法,建立了贴体坐标系下求解流体流动和传热的双分布格子Boltzmann模型.在密度分布函数和温度分布函数对应的离散速度方程中,时间项采用四阶Runge-Kutta法离散,空间离散采用二阶迎风和二阶中心差分的混合形式.采用此模型分别对瑞利数为10~3、10~4、10~5、10~6的方腔自然对流以及理查森数为0.1、1、10的方腔混合对流进行了数值模拟,获得了流体速度与温度分布的典型特征,得到的努塞尔数也与基准解高度吻合.计算结果表明了本文采用的数值方法和计算程序的有效性.     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

非线性振动系统的多项式向量方法

对于常微分方程描述的非线性振动系统,当采用摄动方法求近似解时,先是给出满足各阶近似解的二阶常微分方程组,继而依次对每一个常微分方程进行求解,以致多自由度非线性振动系统的求解过程相当繁琐.文章针对常微分方程表示的非线性振动系统,提出了一种求解非线性振动系统近似解的多项式向量方法,该方法将二阶常微分方程组表示成一阶状态方程组,将非线性部分写成常数矩阵和多项式向量之积的形式.然后,采用直接摄动方法,获得每个幂次近似解所满足的一组状态方程,此时状态方程的非线性部分成为常数矩阵和前一幂次近似解作为元素组成的多项式向量的乘积.进一步,借助Toeplitz矩阵将多项式向量之乘法表示成矩阵形式,以解决多项式相乘带来的幂次方系数的确定问题,再根据一阶非齐次方程组的求解方法,获得状态方程组的全部近似解析解.多项式向量方法将二阶常微分描述的非线性振动求解过程转换为一阶非齐次状态方程组的求解问题,计算过程主要是矩阵和向量之间乘法运算,提高了计算效率和程序化水平.     

最小二乘配点法解薄板几何非线性弯曲问题

用最小二乘配点法分析薄板弯曲问题,国内于1978年首先由徐次达、施德芳用幂级数作为挠度试函数进行了解算.1979年何广乾、张维岳成功地以最小二乘边界配点法解算了壳体的线性弯曲问题.本文系用最小二乘内部配点法解算薄板几何非线性弯曲问题,所用挠度试函数为双三角级数,以边界可动简支方形柔韧板为例进行计算,在解非线性方程组时采用Levenberg-Marguardt法(阻尼最小二乘法),克服方程组的奇异性和病态.计算结果与Levy的成果进行了比较.