固体力学计算方法的发展研究.pdfVIP

固体力学计算方法的发展 孙秀山 岑章志 刘应华 (清华大学工程力学系，北京100084) 摘要本文简要回顾了固体力学计算方法的发展过程。从早期通 过解析方法求解简单问题开始，固体力学的计算方法经历了一个从精确解 法到近似解法、从解析方法到数值方法的发展过程，这一过程可以依据其 历史阶段分为三种类型：传统解析方法、近似求解方法(古典数值方法)和 现代数值方法。文中分析了不同发展阶段中典型固体力学计算方法的形 成及其特点，探讨了这些方法对固体力学发展的作用以及影响，总结了这 些方法之间的关系。 关键词 固体力学，计算方法，发展过程，继承关系 1 引言 固体力学是在经典牛顿力学框架下最先发展起来的学科之一，主要 研究可变形体在各种外界因素作用下，其内部各个质点所产生的位移、运 动、应力、应变以及破坏等的规律，是力学中形成较早、理论性较强、应用 较广的一个分支[1]。固体力学的发展首先是建立在弹性理论基础之上 的，随后在工业发展的推动下，固体力学中有关塑性理论、强度理论以及 稳定理论等得到了进一步的发展[2门]。在传统的固体力学理论中，一般把 研究对象看作是由无限个假象的元素组合在一起的连续体，因此研究对 象(连续体)中的力学量(如位移、应变、应力等)就可以假设为空间或时间 的连续函数。这样，对于一个确定的固体力学问题，借助于数学方法最终 可以将其转化相应的偏微分方程(或方程组)在给定条件下的边值问题或 初值问题，如经典弹性理论中L-N方程或昏M方程的狄利克雷 ·】52· 程组)的求解一直贯穿着固体力学的整个发展阶段，成为固体力学的重要 研究内容之一。 从早期通过解析方法求解简单问题开始，固体力学的计算方法依据其 历史发展过程大致经历了如下三个阶段：传统的解析方法、近似求解方法 (古典数值方法)和现代数值方法，其中每个阶段里都出现了多种分析方法 和计算方法。在这些方法的发展中，尤以计算机技术的出现和应用为转折 点，标志着固体力学计算方法的一个飞跃，促使了固体力学无论在理论研究 方面还是在实际工程应用中都有了显著的进步E“5I。另一方面，随着计算方 法的不断改进和完善，以有限单元法为代表的现代数值方法不仅把固体力 学的应用广泛推向了诸如航空、土木、机械等工程领域，而且也把诸如固体 力学、一般力学、流体力学等诸多力学分支的问题的求解统一在一个框架体 系内[6]。计算方法的成熟及其工程应用的成功，促成了计算固体力学的诞 生，这一学科成为计算力学学科的一个重要组成部分。在今后的几十年中， 计算固体力学很可能将突破经典力学的框架，继而渗入到诸如量子力学、分 子动力学、材料科学、生物力学等领域中，形成新的交叉学科ET]。 在以下的几节里，将分别讨论不同发展阶段中典型固体力学计算方 法的形成及其特点，并分析这些方法对固体力学发展的作用以及影响，最 后一节将对本文做一小结。 2传统解析方法 传统的解析方法是固体力学发展初期所采用的一种计算求解方法， 这一方法所解决的主要是一些简单的弹性力学问题、稳定问题以及后来 出现极少的塑性力学问题。除了少数简单问题可以由基本方程和边界条 件直接求出解函数之外，多数能求得解析解的问题往往借助于一些专门 的数学手段和工具，具有代表性的包括应力函数法、试凑法(反逆法和半 逆法)以及复变函数法等。 2．1应力函数法 力分布问题的论文中首次用到了应力函数嘲。虽然艾里的研究并不完 ·】53· 全，但他开创了求解弹性问题的应力函数方法。随后麦克斯韦(Maxwell， 方法L1 0|。应力函数法的基本思想是，在以应力为未知量的B—M方程中 引入某些函数，这样就把待求问题转化成求解以这些函数所表示的调和 方程，这些引入的函数就是应力函数，它们自动满足平衡方程，其偏导数 的组合即为应力分量。 应力函数法的出现进一步推动了弹性理论的发展，同时也使研究对 Plates)中利用应力 F．)在1947年出版的专著《各向异性板》(Anisotropic 函数法求解了多种正交各向异性平面问题，包括矩形截面梁