超级计算机在结构力学领域的应用,计算结构力学

《计算结构力学》由会员分享，可在线阅读，更多相关《计算结构力学(80页珍藏版)》请在技术文库上搜索。

1、第1章 前言 1.1 工程结构分析 1.2 有限元法概述 1.3 弹性力学的基本方程 1.4 最小势能原理与弹性结构近似解 1.5 加权残数法构 1.6 结构数值分析的工程应用和未来发展 1计算结构力学 第1章 前言 有限元方法 第1章 前言2 1.1 工程结构分析 structural analysis 目的和意义： 从力学的观点，对 工程结构进行分析计算，从而充 分地了解结构的力学性态(定量 化)，保证结构能够完成预先的 功能，并使得结构形式尽可能的 合理，甚至最优。 2计算结构力学 第1章 前言 有限元方法 第1章 前言3 结构分析(一组控制方程)， 比如：梁、板的方程 解析解答 半解析。

2、 数值解 差分解，对微分 方程进行离散 有限元法或边界 元，对结构域进 行离散 结构分析示意 3计算结构力学 第1章 前言 解析法：问题的求解可以得到具体的数学 表达式，可以计算任意位置的问题的精确 解答。如悬臂梁求挠度等。 数值法：得不到具体表达式，只能得到某 些离散点处的近似值。该方法包括：有限 差分、有限元法、边界元法、变分法等.本 课程重点介绍有限元方法。 4计算结构力学 第1章 前言 有限元方法 第1章 前言5 结构分析的内容 应力分布位移分布 结构的屈曲和稳定性 结构动力学，频率与动力响应 热变形与热应力 应力集中 疲劳与断裂 结构优化分析 流体与结构的相互作用 5计算结构力学 第。

3、1章 前言 1.2 有限元方法概述 有限元方法的实质 Finite element method (FEM) 将复杂的连续体划分为有限多个简单的单 元体，化无限自由度问题为有限自由度问 题，将连续场函数的(偏)微分方程的求 解问题转化成有限个参数的代数方程组的 求解问题。 6计算结构力学 第1章 前言 左图所示,为分析齿轮上一个齿内的应 力分布,可分析图中所示的一个平面截 面内位移分布.作为近似解,可以先求出 图中各三角形顶点的位移.这里的 三角形就是单元,其顶点就是节点。 从物理角度理解从物理角度理解, 可把一个连续的齿形截面单元之间在节 点处以铰链相链接,由单元组合而成的结构近似代替原连续。

4、结 构,在一定的约束条件下,在给定的载荷作用下,就可以求出各 节点的位移,进而求出应力. 从数学角度理解从数学角度理解, 把这个求解区域剖分成许多三角形子域 ,子域内的位移可用相应各节点的待定位移合理插值来表示. 例如例如: 7计算结构力学 第1章 前言 有限元法的基本步骤 (1)结构的离散，将整个结构分为若干个小 的单元(单元的形状：三角形、四边形等) 。 (2)确定单元节点力和节点位移的关系(单 元位移和节点位移的关系，单元平衡方程) 。 (3)形成整体结构分析方程。 (4)引入边界条件，求解方程。 (5)进行后处理分析。 8计算结构力学 第1章 前言 有限元方法 第1章 前言9 有限元方。

5、法实际上是一种数值方 法。 因此也离不开计算机的发展。此 外作为一种数值方法，在数学表 示上，采用矩阵的方式。 9计算结构力学 第1章 前言 10 静力与动力 static and dynamic 线性与非线性 linear and nonlinear 确定性与概率与随机 determine problem and random 有限元方法进行求解 10计算结构力学 第1章 前言 有限元分类 基本未知量分类：位移法、力法、混合法 位移法应用最广，最早是位移协调单元 ，后来逐步提出非协调元、拟协调元、 广义协调元等。 混合法有杂交元，此外分区混合元在断 裂力学中得到了广泛应用。 11计算结构力学。

6、 第1章 前言 有限元方程的推导分类 ： 直接法(只适用于简单问题) 变分法(有限元有了坚实的数学基础) 加权残数(直接从方程出发，对于不存 在变分泛函的问题可也建立有限元方程 ，扩大了有限元的应用范围) 12计算结构力学 第1章 前言 有限元的发展历史 1956年MJ.Turner, R.W.Clough, H.C.Martin, L.J.Topp在纽约举行的航 空学会年会上介绍了一种新的计算方法 ，将矩阵位移法推广到求解平面应力问 题。他们把结构划分成一个个三角形和 矩形的“单元”，利用单元中近似位移 函数，求得单元节点力与节点位移关系 的单元刚度矩阵。 13计算结构力学 第1章 前言 1。

7、954-1955年，J.H.Argyris在航空工程杂 志上发表了一组能量原理和结构分析论 文 1960年，Clough在他的名为“The finite element in plane stress analysis”的论文中首次提出了有限元 (finite element)这一术语。 14计算结构力学 第1章 前言 在1963年前后，经过J.F.Besseling, R.J.Melosh, R.E.Jones, R.H.Gallaher, T.H.H.Pian(卞学磺)等许多人的工作，认识 到有限元法就是变分原理中Ritz近似法的一种 变形，发展了用各种不同变分原理导出的有限 元计算公式。。

8、 1965年O.C.Zienkiewicz和Y.K.Cheung(张佑 启)发现只要能写成变分形式的所有场问题， 都可以用与固体力学有限元法的相同步骤求解 。 1969年B.A.Szabo和G.C.Lee指出可以用加权余 量法特别是Galerkin法，导出标准的有限元过 程来求解非结构问题。 15计算结构力学 第1章 前言 我国的力学工作者为有限元方法的初期 发展做出了许多贡献，其中比较著名的 有：陈伯屏(结构矩阵方法)，钱令希 (余能原理)，钱伟长(广义变分原理 )，胡海昌(广义变分原理)，冯康( 有限单元法理论)。 16计算结构力学 第1章 前言 有限元法已经成功地应用在以下一些领域： 固。

9、体力学， 包括强度、稳定性、震动和瞬态问题的分析 ； 传热学；电磁场； 流体力学与声学。 土木结构、船舶结构、海洋结构、汽车、航 天与航空、医学等等 17计算结构力学 第1章 前言 目前发展了很多商用程序 18计算结构力学 第1章 前言 现代计算机的发展为FEM的发展提供了 物质基础，计算规模越来越大，最大有上 千万自由度计算,甚至还有在超级计算机 上进行的并行计算。 模型的输入、输出(即有限元的前后处理) 实现了图形化,典型的有限元前后处理软 件如:MSC/PATRAN. 19计算结构力学 第1章 前言 结构： 杆系，板壳，一般的固体 在弹性范围 内研究 分析方法 ： 材料力学 ， 板壳力学。

10、， 船 舶结构力学， 更一般的 弹性力学 1.3 弹性力学问题基本方程 Basic governing equation 1.3.1弹性力学的基本变量 20计算结构力学 第1章 前言 1.3.1 一点的应力状态 21计算结构力学 第1章 前言 有9个应力分量 或： 22计算结构力学 第1章 前言 1.3.2 平衡方程 equilibrium equation 23计算结构力学 第1章 前言 24计算结构力学 第1章 前言 上述结果极易推广到三微微元体的平衡： 25计算结构力学 第1章 前言 由对原点的弯矩为领， 导致： 称之为剪应力互等定理。 由此进一步得到， 微元体独立的应力变量为6 个： 。

11、26计算结构力学 第1章 前言 1.3.3 几何方程 Geometrical equation 27计算结构力学 第1章 前言 图(a) 28计算结构力学 第1章 前言 图(b) 29计算结构力学 第1章 前言 最后叠加, 并将其定义为应变 其中： 更一般的三微微元体： 30计算结构力学 第1章 前言 给出了位移 (displacement) 和应 变 (strain) 的关系 , 称之为几何方程， 注意它和结构是否 为弹性体并不相关 。 31计算结构力学 第1章 前言 32 弹性力学的基本变量-15个 应力 位移 应变 32计算结构力学 第1章 前言 1.3.4 应力应变关系 剪切弹性模量 。

12、Shear elastic modulus 广义虎克定律 Hookes Law 33计算结构力学 第1章 前言 34计算结构力学 第1章 前言 1.3.5 边界条件 位移边界 35计算结构力学 第1章 前言 力边界 力边界 36计算结构力学 第1章 前言 或针对一般的三维问题： 37计算结构力学 第1章 前言 整个弹性力学三维问题归结上述15个微分 方程的求解， 并由此给出15个变量。 适当的边界条件。 位移边界： 给定位移 力边界： 给定边界力 解析解答-(弹性力学课程) 数值解答-(其中之一就是有限元方法) 1.3.6 弹性力学边值问题的提法 38计算结构力学 第1章 前言 1.4最小势能。

13、原理与弹性结构近似解 The principle for minimum potential energy and ritz method 39计算结构力学 第1章 前言 意义： (1) 势能取极小值时， 对应的是真实解 (2) 将解方程的问题转化为求极小值问题。 (3) 通过尽可能使得势能极小，可以得到近 似解。此即所谓的RITZ方法。 40计算结构力学 第1章 前言 最小势能原理的证明 41计算结构力学 第1章 前言 注意到： 42计算结构力学 第1章 前言 上式右边第一个积分可以进一步写为 43计算结构力学 第1章 前言 44计算结构力学 第1章 前言 由于 必然有上述的平衡方程和力边界。

14、条件得以满足。 所以 (1) 满足位移边界的势能极小等价于平衡方程 和力边界条件 (2) 解微分方程的问题转化为极小问题 对于的极小 45计算结构力学 第1章 前言 经典弹性力学的近似解- RITZ方法 46计算结构力学 第1章 前言 举例1 简支梁的计算 一旦假定了位移的形式， 则极小化问题变为 对待定参数的极小化 只是对所选定位移函数形式，导致的极小化， 不是所有的满足位移边界条件的位移函数， 因 此是近似。 整个位移函数是在整个结构区域选取的。 47计算结构力学 第1章 前言 部分满足位移边界的 势能极小, 必然不 是满足全部位移边界位移势能极小, 于是 得到近似解. 求弹性力学的微分方。

15、程组可以转化为求 极值问题。 所有满足位移边界的 势能极小 真实解 不是极小 近似解 48计算结构力学 第1章 前言 特别说明 从解微分方程的角度讲: 完全的求解(找到满足边界条件的全部位移)很困难 但找到一组满足位移边界的位移还是可以的, 于是 得到近似解. 49计算结构力学 第1章 前言 讨论 50计算结构力学 第1章 前言 位移的界限 51计算结构力学 第1章 前言 1.5 加权残数法 52计算结构力学 第1章 前言 53计算结构力学 第1章 前言 应满足位移边界条件 54计算结构力学 第1章 前言 (1) 配点法 相当于简单地强迫在域内若干点微分方程 的残差值为零. 55计算结构力学 。

16、第1章 前言 (2) 子域法 在若干个子域内权为单位1, 之外为零. 56计算结构力学 第1章 前言 (3) 最小二乘 于是, 权函数为: 57计算结构力学 第1章 前言 (4) 力矩法 58计算结构力学 第1章 前言 (5) 伽辽金法 59计算结构力学 第1章 前言 关于加权残数法的简单讨论 以后只针对最小势能原理介绍有限元 60计算结构力学 第1章 前言 1.6 结构数值分析的工程应用和未来发展 海洋平台 1.6.1 各类工程问题的应用 61计算结构力学 第1章 前言 卫星结构 62计算结构力学 第1章 前言 CCD照相机 63计算结构力学 第1章 前言 Global Hull Struc。

17、ture analysis The global FEM model based PATRON PRO-PROCESS (SJTU) 64计算结构力学 第1章 前言 Global Hull Structure analysis Vibation analysis (SJTU) 65计算结构力学 第1章 前言 Sound prediction based nastran and sysnoise (sjtu) Dynamic response is obtained from NASTRAN and sound field is gotten from SYSNOISE 66计算结构力学 第1章。

18、 前言 Fatigue and Life Prediction Local mesh test result 67计算结构力学 第1章 前言 The results includes: Damage residual deformation; Fracture and tearing Collision force history Energy absorption 68计算结构力学 第1章 前言 6969 Ultimate Strength (lll.1) PROGRESSIVE HULL COLLAPSE Between Trans. frames Between Trans. bulkh。

19、eads Under combined hull girder actions 69计算结构力学 第1章 前言 1.6.2 有限元分析的优越性 (1) 对复杂几何构件的适应性 70计算结构力学 第1章 前言 (2)对各种复杂材料细观结构的分析 71计算结构力学 第1章 前言 单元脱层过程的模拟 脱层断裂准则 网格重划分技术 72计算结构力学 第1章 前言 复杂复合材料的损伤模拟 73计算结构力学 第1章 前言 解析单元嵌入有限元中 跨尺度计算或称多尺度计算 74计算结构力学 第1章 前言 (3)对于各种物理问题的可应用性 有限元法求解的是物理问题的控制方程 ， 对结构的线弹性，弹塑性问题，粘弹。

20、塑 性问题， 动力问题，屈曲问题， 热传导问题， 电磁场 75计算结构力学 第1章 前言 (4)建立于严格理论基础的可靠性 76计算结构力学 第1章 前言 (5) 适合计算机实现的高效性 单元、形函数等的规范化 求解方程矩阵形式的规范化 有限元法： 适合计算机程序的编制和执行 计算机软、硬件的飞速发展 77计算结构力学 第1章 前言 1.6.3 结构数值分析的未来发展 1. 为真实模拟新材料、新结构的行为，发展单元类 型、新材料本构。 2. 为真实模拟新材料、新结构的行为、模拟各类形 式的结构在复杂工矿和环境作用下的全寿命过程的 响应。 3. 各类工程分析的标准化问题(包括单元的选取网 格的大小后处理分析),流程再造 3. 有限元软件和CAD/CAM/CAE等软件系统共同集成 完整的虚拟产品发展系统 78计算结构力学 第1章 前言 本课程的目标 有限元方法的基本思想 线性静力学问题 线性动力响应、模态计算 非线性分析初步 借助于大型商用软件(如NASTRAN) 实 现典型结构的有限元分析 79计算结构力学 第1章 前言 80计算结构力学 第1章 前言 。