有限元与数值方法-讲稿1

有限元与数值方法第一讲授课教师：刘书田Tel:84706149; Email:stliu@ 教室：综合教学楼 351

时间：2013年3月15日：8:00—10:50

教学内容计算固体力学的基本理论 固体力学(以弹性力学为主描述)的基本理论 能量、变分原理和变分法 特殊问题的数值计算方法 介绍各类方法的构造过程 有限差分法(Finite Different Method) 加权残数法(Weighted Residual Method) 有限元法(Finite Element Method) 无网格法(Meshless Method) 边界元方法(Boundary element Method)

计算固体力学的主要方法

有限元法的应用和前后处理

参考教材 R.D.Cook,有限元分析的概念与应用 (Concepts and Application of Finite Element Analysis)关正西等译，西安交大 出版社 王勖成等, 有限单元法基本原理和数值方法， 清华大学出版社 杨庆生, 现代计算固体力学，科学出版社 刘正兴等, 计算固体力学，上海交大出版社

第一章 前言一 计算固体力学的任务：1.力学的任务物体机械运动的规律 ——研究物体受到的力和物体发生的运动的关系 物体(流体，固体，气体) 力(热，电，磁等环境) 运动

流体力学：研究对象是流体 2.固体力学的任务

(水，空气等);

固体力学

研究固体(结构)在外部作用(外力，温度等变化)下的变形 和应力及其演化规律，根据这些规律研究固体和结构的破坏 (刚度、强度、疲劳、断裂以及稳定性等) 根据研究对象的不同：弹性力学，塑性力学，断裂力学，冲击 力学; 材料力学，理论力学等 根据采用的方法：实验，理论和计算4

固体力学的任务(续) 重点：建立固体在外部作用下的变形和应力 以及演化规律的数学模型(控制方程) 例如： 应力~外力之间的关系：平衡方程(运动方程) 应力~应变之间的关系：本构方程 研究变形的机理，变形的诱因(外部作用)、应 力和应变的定量关系： 弹性问题：Hooke定律 热弹性问题：热膨胀规律 塑性问题：屈服条件；强化准则；流动准则 断裂问题：起裂条件；扩展规律5

变形的描述以及几何关系 主要研究变形的描述方式(应变，位移，转角等) 建立变形与位移之间的定量关系 应变与位移之间的定量关系例如：小变形条件下：

1 ui u j ij = 2 x j xi

1 ui u j 1 um um = 有限变形条件下： ij 2 x j xi 2 x j xi 边界条件： 位移边界：应力边界：

ui ( x) ui

x u

ij ( x)n j ti

x 6

求解方法以及对应的控制方程 (1)力法未知量：以应力作为基

本未知量 控制方程：平衡方程；相容方程(变形协调方程)

(2)位移法：未知量：以位移作为基本未知量 控制方程：位移表示的平衡方程； 边界条件

对应原理，变分原理

研究：微分方程的积分形式， 泛函变分与基本方程的对应 建立各种问题所对应的变分原理 任务：国体力学是建立固体变形规律所必须 满足的规律以及数学模型，为各种求解策略 提供理论基础。

计算固体力学的任务和研究内容 任务：以固体力学的基本理论为基础，研究利用计算机 科学与技术求解固体力学中各类问题的数值分析 理论、方法、建模、软件实现；

研究内容：(1)研究固体力学中各类问题的数值计算方法，基 本原理； (2)采用数值模拟技术，分析固体的变形演化规律、 破坏规律、应力分布规律，揭示新的力学现象， 包括材料性能揭示；工程中的力学问题等。 (3)工程问题的模型化、可视化、虚拟现实9

结构分析问题 各种工程结构常见的结构元件：(1)杆、梁、柱(长>>宽和高) (2)板(中厚板)、壳(厚<<长和宽) (3)三维体 (4)薄壁结构(飞机机翼与机身等) (5)以上结构类型的复合体

结构分析问题包括：

(1)强度问题(应力) (2)刚度问题(变形) (3)稳定性问题 (4)振动问题10

计算固体力学的主要方法 有限元法(位移协调元,杂交元，应力元，拟协调元) 边界元法 无网格法(mesh-free method): Non-structural finite difference(Orkisz,2001); Element-free Galerkin(Belytschko,1994) Smooth particle hydrodynamic(Gingold,1997) Partition of Unity((Melenk,1996) Finite Point method(Onate,1996) Meshless finite element(Onate,2003) Finite sphere(Bathe,2001) Natural element (Belytschko,1998) 扩展的有限元法(x-FEM) 等几何法(isogeometric method) 变分法 加权残数法11

有限元法的发展历史近似求解偏微分方程的数值方法： Lord Rayleigh and Ritz , Galerkin 采用试函数(trial functions) 对偏微分方程的解进行近似 Courant 引入子域内分片连续试函数( piecewise-continuous functions) 的概念，标志着有限元方法的起始 1960s. Clough在平面应力分析中引入 finite element 的名称 1960s-1970s. 板弯曲、壳弯曲、压力容器、三维弹性问题、流 动、热传导等采用有限元方法求解 ；美国空间计划支持 Nastran的开发 1970s. 开发了ANSYS , ALGOR , COSMOS/M , SAP, NONSAP and ABAQUS etc. 目前. FEM系统可在微机上解决大规模结构分析问题12

计算力学发展展望 计算力学研究采用计算机和相应计算方法求解力学问题、认识力学现象的方 法、理论、软件实现和工程应用。计算力学是力学学科和计

算机科学技术交 叉而形成的力学分支，是计算科学和工程的核心学科。

计算力学起始于有限元法 . 有限元法的诞生可追溯到 50 年代中期 Martin, Clough,Turner(1956),Argyris(1955) 等的 …… 此处隐藏：2551字，全部文档内容请下载后查看。喜欢就下载吧 ……