边界元法中非对称满系数矩阵方程组的拟波阵解法及程序

本文提出求解边界元法中非对称满系数矩阵方程组的拟波阵法,即边形成方程边采用高斯-约当消去法进行消元,在计算机中只需存储整个系数矩阵的四分之一部分。文中给出了相应的FORTRAN子程序。经使用表明,该程序经济有效,简单实用,大量节省了计算机的容量,加倍提高了解题能力,为小机器解大题目提供了方便。此外,该法普遍适用于任何可以逐行形成方程的非对称满系数矩阵方程组的求解。     

常微分方程技术及其在固体力学计算中的应用

常微分方程(ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展推出了一批高质高效的通用软件,而工程中大量的ODE问题并非呈现为这些求解器(Solver)所接受的标准形式。然而,运用一些简单的ODE变换技巧可以将大量的不同类型的特殊问题转化为标准形式。本文列举了若干常用的变换技巧,并广泛地应用于各种固体力学问题的计算中,使大量的ODE问题在形式上得到统一,得以用标准的ODE Solver方便有效地求解。     

介绍一个常微分方程边值问题求解通用程序——COLSYS

1.引言 随着近年来常微分方程(简称:ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展,一批高质高效的通用求解软件程序相继问世。就笔者对一些软件的使用经验及比较来看,COLSYS程序是一个很值得推荐的ODE求解器。本文特对此程序作一介绍。     

空间几何构造分析的有限单元法

提出空间杆系几何构造分析的有限单元法，构造了两种单元(链杆单元和准梁单元)的几何约束矩阵，集成为整体矩阵并引入承条件后，通过对其阶数与秩的比较分析确定体系的几何可变性及静定性．本法原理简单，便于计算机实施，结果完备：对于几何不变体系，可指出多余约束的数目；对于几何可变体系，可给出体系的自由度数及相应的运动模态，并确定自由度的常变瞬变性质．     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

二维切口问题完备特征解的研究现状与进展

简要回顾和介绍了文献中常见的几种求解二维切口问题特征根的方法,评述这些方法 的特点及其求解范围．重点介绍了在求解二维切口问题完备特征根方面的一些进展,包括特征 根的求解方法、特征根的重根探讨以及相应的特征函数计算,这一问题基本上得到了解决．为 基于特征展开研究二维断裂问题打下了基础．      

P型有限元线法分析板弯曲问题

Ｐ型有限元线法分析板弯曲问题袁驷（北京清华大学，１０００８４）关键词有限元线法，板弯曲问题，中厚板理论互方法简介有限元线法（ＦＥＭＯＬ）是一种新型的以常微分方程（ＯＤＥ）求解器为支撑软件的半解析数值方法，其基本理论可参见文【１］【２］及有关文献，这里...     

极坐标有限条法解扁球壳问题

本文用有限条法解扁球壳问题,对于轴对称问题与位移和内力沿环向按cosnθ或sinnθ分布的非轴对称问题,分别推导了三类壳元(圆底壳元,圆孔外无限壳元,圆环底壳元)刚度矩阵的精确解。对于沿环向按cosnθ或sinnθ分布的环形面荷载和环形线荷载推导了单元等价节线荷载的表达式。编制了相应的计算程序,计算了典型例题。结果表明,此法具有计算工作量少而计算精度高的优点。     

变截面框架-剪力墙-薄壁筒结构弯扭耦连振动的常微分方程求解器解法

本文用沿高度方向分段连续化的方法,对沿高度方向为阶形变截面,在平面内为任意斜向布置的框架-剪力墙-薄壁筒协同工作体系,建立了弯扭耦连的振动方程,用常微分方程求解器(COLSYS)求解了自振频率及相应的振型。