计算固体力学

绪论 固体力学主要研究在各种外界因素作用下可变形体内部各点所产生的位移、皮力。应变以及破坏等的媒律：假设研究对象中的位移、应变、应力等为空间或时间的连续函数，借助于数学方法将其研究问题转化成相应的偏微分方程边值向题过初边值问题。用微分方程来描述工程技术问题是科学的一大成果，其求**直贯穿于固体力学的发展阶段。 固体力学遇到的数理模型是复杂多样的，其计算方法已经历了三个发展时期解析方然，古典数值方然和现代数值方法。 在固体力学发展初期，科学家针对基本方程和边界条件的定解问题提出了许多解析方法，如应力函数法，试凑法(反道和半逆法)及复变函数法等，这些方法所解决的主要是一些简单的弹性力学问题.稳定问题及后来出现得极少的整性力学问题。除了少数简单图体力学问题外，解析方法是不可行的。随着固体力学白身的发展及实际工程问题的出现，许多复杂的问题求解开始逐渐引入近似的求解方法。 与传统解析方法对数学的完美要求相比，近似解法更注重在工程问题中的实用性。古典数值方法在数学形式上就是利用近似解代替精确解，近似解不一定严格满足基本方程成边界条件，即放松了对解的限制。历史上*早采用的数值方法是有限差分法，从微分方程出发，将连续的定解区城用有限个离散点构成的网格来代替，用奉勒展开式将原力程和定解条件中的微商用不同时间成空间点差商来近似把原微分方程和边界条件的求解转变为求解一个线性代数方程组，从而得到原问题在离散点上的近似解，再利用插值方法便可以得到整个区城上的近似解。另一种近似方法是某于等效和分的数值方法。倒如，瑞士科学家里兹于1908年尊里兹法作为一种有效方法提出.基于变分法(积分方程)中*小势能原理或虚位移原理，选择一个满足位移边界条件的近似函数，对泛画求驻值，得到一组线性代数方程，从而获得问题的近似解。另外，苏联数学家伽辽叠于1915年发明了伽辽金法，采用微分方程对应的等效积分弱形式，选择满足位移边界条件(与力边界条件)的近似函数，并把近似函数中基所数或形函数为权函数，得到一组线性代做方程，可得问题的近似解。这种方法属于加权余量法。里兹法和伽辽金法均用有限白山座体系近似代替了无限自由度体系，莴者在某个特定的条件下是等效的。有限差分法.里旅法，伽辽金法等近似解法的出现表明固体力学从初期的单纯理论研究逐新转入到实际工程应用之中。但是，还存在不满意之处：有限差分法局限于规则的差分网格，如正方形矩形或正三角形网格等；里兹法和伽辽金法选择的近似函数必须调足整个求解区域。当研究对象是一个复杂的结构或具有复杂几何形状，近似函数不易得到满足；所以古典数值方法对于模拟复杂边界的二维或三维问题有一定难度。 现代数值方法则抛弃了这种在整个求解域上选取近似函数的思想。求解模型中采用了“离散化"的思想，近似函数仅雷在局部满足微分方程。基于变分法成加权佘量去，提出了流行的有限单元法，边界单元法，无网格法等。 20世纪60年代，计算机技术的出现和应用为标志着国体力学计算的一个飞跃，使图体力学的现代数值方法进人了前所未有的深度与广度。电子计算机应用的飞速发展，以及计算方法的不断改进和完善，促成了计算固体力学学科的诞生。 计算固体力学是采用离散化的数值方法，并以电子计算机为工具，求解固体力学中各类问题的学科。借助于计算机，有限元法与有限差分钟相辅相成，已成为现代工程计算中不可缺少的强有力工具。但是.有限差分法只看到了节点的作用，没有注意到连接节点的单元所起到的作用；有限元法吸取了有限差分法中离敏化处理的内核，又继承了空分计算中选择插值函数井对区城进行积分的合理方法，并且充分估计了单元对节点参数的贡献，使计算结果更为精确。因而。在工程计算中，以有限元法为核心的现代数值方法得到了广泛的成用。 计算固体力学应用到的工程问题及其求解的特点： (1)静力学问题。离散化后归结为求解线性代数方程组。常见于求解结构的应力和变形。 (2)特征值问题。离散化后归增为求解矩阵的特征值和特征向量问题，常见于求解结构成系统的频率和振皇，稳定极限载荷和屈曲形状。 (3)动态响应问题，离散化后得到一常微分方座组，可直接数值机分成利用先求得特征向量将它转换为一组互不耦合的常微分方程，再进行时同积分求解。常见于求解结构的振动和弹性波的传播。 (4)非线性问题。例如，黏弹(塑)性等物理非线性问题，大变形和后屈曲等几何非线性问题，一般采用增量算法将它们转化为一系到线性问题求解。 (5)含裂纹的非连续问题。可采用奇异单元模拟裂纹尖端的应力场。 (6)复合材料和结构的非均质问题。目前，对此类问题求解还不完着，正在发展之中。 (7)多场耦合问题。对此类问题求解也在发展之中。 计算固体力学研究和应用的领城不断扩大，随计算机技术的发展，解题能力成数量级地提高。例如，借助计算机，已能对整个“鸟巢”、整艘航空母舰，或整架飞机等工程问题进行详细分析，并得到满意结果。 计算固体力学的发展，既有其学科自身的要求，也有实际工程问题的推动。1997年9月，钱学森先生给予清华大学力学系赠言：“随着力学计算能力的提高，用力学理论解决设计何题成为主要途径，而试验手段成为次要的了。由此展望21世纪，力学加电子计算机将成为工程设计的主要手段。”计算固体力学的发展方向是：在数值方法方面，利用多种数值方法的优点，取长补短，提高大型系统的非线性分析.随机分析.*合分析等算法的精度和效率，改进其稳定性和收效性；在应用方面，充分利用计算机图像、数据库，人工智能等技术，并可与优化设计，可靠性设计等相结合，发展多功能，自动化的道用成专用工程软件系统，将突破经典力学的植架，继而渗人到诸如生物力学，量子力学等领城，形成新的交叉学科。 本书将讨论线性，材料非线性，几何非线性.动力学问题和热分析向题；涉及加权余量铁，有限元法，边界元法，无网格法和离散元法等，主要介绍这些方法的基本原理和概念。鉴于研究生曾学过“教值分析"，“线性代数"与“弹性理论"，为了避免重复，不再整述有限益分法，矩阵算法和变分法。 著名的有限元法，边界元法既可以从变分原理推出，也可以用加权余量法推出。已经证明，加权余量法用于存在泛函极值的微分方程与泛函级值是等价的，遗憾的是，至今还有一部分散分方程没有找到对应的泛函。换句话说.直接针对原始微分方程推导出来的加权余量法比变分法重有优势，这是因为它也适用于不能给出近(需对其求板小值)的事选间以，既然加权余量法包容了壹分原理中泛函极值.有限元法.边界元法等的*普遍原理，它更容易推广应用到不同微分方程的其他向题，所以本书以加权余量法开篇。有限元法是当代计算固体力学应用的核心，着量*多。边界元法是对有限元法的补充，二者取长补姬，其期合计算方法将来也许是个发展方向，故将边界元法作为单独章节进行核心技术的介绍。目前，无网格法是计算图体力学研究领城的前沿热点，出生较晚，有特成长，甚至有些概念的“名字”还在争议之中。无网格法有许多优点，甚至有专家预测无网格法将成为继有限元法之后新一代的数值方法。所以，安排在边界元法章节之后，单独介留无网格然，便于启迪这方面的研究。 不论有限元法.边界元法，还是无网格法，均是基于连续介质力学基础之上，在非连续介团力学领城的计算又如何呢？离散元法就是该领域败值计算方法的典型代表，主要用来模拟大量颗粒在给定条件下如何运功。在计算机的辅助下。离散元法甚至可以完成模拟“介于流体和固体之间的颗粒或者粉末”受力与运动分析。离散元法的应用已扩展到连续介质及连续介质向非连统介质转化的力学问题。例如。冲击、侵彻等动载荷作用下材料的破坏。为了知识点的全面性，在计算固体力学的简尾，把离散元法作为其扩充部分。 思考 1.有限益分法的优缺点各是什么？ 2.计算固体力学的生命力如何？ 3.回忆虚位移原理，虚功原理，*小势能原理，里兹法。 4.什么是自然变分原理和广义变分原理？ 弹性力学*小势能原理和*小余能原理都属于自然变分原理。在自然变分原理中试探函数事先应满是规定的条件。例如，*小势的原理中位移试探函数应事先满足几何方程和给定的位移边界条件；*小余能原理中应力试探函数应事先满足平衡方程和给定的外力边界条件。如果这些条件未事先满足，则需要利用一定的方法将它们引入泛函。这类变分原理称为的束变分原理，或广义变分原理。利用广义变分原理可以扩大选择试探函数的范圈，从而提高利用变分原理求解数学物理问题的能力。 第1章加权余量法 1.1微分方程的等效积分 应用科学和工程问题往往可以归结为：在一定边界条件，初始条件下，求解问题的控制微分方程(组)。微分方程(组)可以是常微分方程，偏微分方程，线性的或非线性的。例如，某一应用科学问题中的控制微分方程式及边界条件分别为 式中，为待求的函数；A、B为微分算子；为不的已知函数。 微分方程组(1-1-1)的等效积分形式 式中为为任意函数。也称为权函数。由于前响为任意函数上式与微分形式(1-1-1)是完全等价的。 假设微分方程组(1-1-2)中A的微分算子为x阶，对微分方程的等效积分形式(1-1-2)进行次分部积分，得到微分方程组(1-1-1)等效积分弱形式 分部积分后，微分算子EF为阶，微分算子C，D为m阶。将微分方程转化为弱形式，这个弱并不是满化对方程解的结果，而是弱化对解方程的要求，具体是弱化待求函数的连续性，当然这种购化是以提高权函数的连续性为代价的。权函数为选择的已知函数，能够满足分部职分方法对权函数连续性要求。这种弱化摸来了以下好处： (1)降低对未知函数的连续性的要求，从而可以在更广泛的范国内选择试深函数； (2)对连续介质问题，便于有限元构造单元和插值函数； (3)在物理上更符合实际问题对来如函数连续性的要求。 如果在微分方程的等效机分弱形式中，对场函数和任意权函数的连续性要求是相同的，则称为微分方程的对称等效机分弱形式；如果对场函数和任意权函数的连续性要求是不相同的，则称为微分方程的非对称等效积分场形式。 以二维稳态热传导问题的微分方程和边界条件等效的积分弱形式(1-1-4)进行说明。