静力学的数值方法.docVIP

静力学的数值方法 Zheng Z.X. Wei Z. 理论力学作为工程力学的一个分支必然是要与现代计算技术相结合，在进行理论探究的同时还要能提供数值结果。在理论力学的学习中、对基本概念的深入理解及运算技巧的掌握是重要的，但作为其工程应用，在计算机上进行数值计算的基本技能也是必须掌握的。为此我们编写了本教材，希望它有助于同学们在掌握理论力学基础理论的同时也熟悉数值计算的基本方法。这里列出的计算程序均使用FORTRAN95语言。由于运动学与动力学部分问题的“个性化”比较突出，一时还未能建立起一个能涵盖在运动学或动力学中常见问题的通用程序,所以本文目前仅限于静力学部分,以后将逐步补充运动学与动力学部分并将更名为“理论力学的数值方法”。 第一章 矩阵与矢量的基本运算 矩阵是呈矩形排列的一组数。矩阵[A]的i行j列元素记为Aij。 n1行n2列的矩阵称为nl×n2型矩阵,记作[A]n1×n2 。行数与列数相等的矩阵(即此处n1=n2)称为方阵。如果方阵中，只有等于1的主对角元素，而其它元素均为0，则称之为单位阵，记为[I]。 只有一行的矩阵（n1=1）称为行阵；只有一列的矩阵（n2=1)则称为列阵。 矢量可用一个列阵表示。表示一个矢量所必须元素的数目称为该矢量的维数，即相应列阵的行数。矢量作为一个列阵除了如同一般