本书以一维杆单元为例,系统地阐述了有限单元法的基本原理、数值方法、程序实现和固体力学领域各类问题中的应用。
全书共13章。前6章为有限单元法的理论基础,包括直接刚度法,一维杆的“强”形式与“弱”形式,单元和插值函数的构造,加权余量法与虚功原理建立有限元格式,变分原理建立有限元格式。后7章为专题部分,包括线性静态有限元分析,线性动态有限元分析,几何非线性有限元分析,材料非线性有限元分析,复合材料多尺度分析,结构灵敏度分析,桁架结构有限元教学软件EFESTS。本书通过一维杆单元详尽地展示了有限单元法的细节,使读者更容易地学习有限元理论,这是作者的基本出发点,也是本书的特色。
本书适用于力学专业本科生、硕士生教学,可作为机械、土木、航空航天、材料工程和交通运输等理工科专业的研究生教材,也可作为相关专业高等学校教师和工程技术人员的参考书。
在当当、京东等网站可以搜索并购买本书。
有限单元法是工程分析中应用最广泛的微分方程数值解法,不仅被普遍地列为工科专业本科生和研究生的学位课程,而且也是相关工程技术人员和教师继续学习的重要内容。自20世纪50年代有限单元法问世以来,相关专著和教材层出不穷,涵盖的基本内容有:有限元一般原理和表达格式、单元和插值函数、等参元和数值积分、线性问题及其解法、非线性问题及其解法。由于有限元是一种数值方法,对自由度超过3个的问题,就不容易在纸面上展开详细推导,所以这些专著和教材基本上都重点关注有限元表达格式的推导,而对有限元理论结合具体单元的推导,以及有限元理论对具体案例的应用过程,细节展示不够。因此给大多数读者的感觉是有限元理论过于抽象和复杂,即使其本身是非常具体的数值方法。作者学习有限单元法以来,也有类似的感受,相关书籍很多,但是容易上手、容易看懂、适合自学的并不多,尤其是非线性有限元,需要很多书籍结合着阅读才能有所理解。
鉴于这些原因,作者选择最简单的一维轴力杆单元来详细地展示有限单元法的基本原理与求解方法。本书首先通过直接刚度法感性地展示有限元的基本过程;接着针对杆单元的“强”形式,构造了各种“弱”形式以及各种变分形式;其次采用杆单元给出了线性静态分析、线性动态分析、几何非线性分析、材料非线性分析、复合材料多尺度分析的有限元格式、解法以及案例的详细步骤;然后推导了有限元方程的导数,也即灵敏度分析;最后给出了作者开发的开源EFESTS软件和部分Matlab代码。此外,结构灵敏度分析本应是属于结构优化的内容,但是对有限元方程求其导数是个很自然的过程,所以将结构灵敏度分析列在本书中也是一个新的尝试。EFESTS是个完全面向有限元教学的软件,具有友好的图形用户界面、快速的二维桁架结构有限元建模、详细的有限元分析过程与数据导出、求解过程与结果的可视化,比如单元刚度矩阵的组装过程就是动态可视化的。对于有限元的每种分析方法,书中都通过一个小规模的桁架结构,详细地展示了求解的每一步骤,有助于读者深入理解其基本理论。总之,通过杆单元这只“五脏俱全”的“小麻雀”,较为具体地展示了有限单元法的细节,使读者能够容易地学习有限单元法。同时,建议读者结合其他较为抽象的有限元书籍,互为补充,一起学习有限单元法,避免杆单元“一叶障目”带来的不足。
前言
符号表
第1章 绪论
1.1 有限单元法的发展史
1.2 有限元软件的发展史
1.3 本书的写作目的
1.4 关于学习本书的建议
第2章 直接刚度法
2.1 一维杆单元
2.2 结构平衡方程的组装
2.3 边界条件施加与方程求解
2.4 二维杆单元
2.5 三维杆单元
第3章 一维杆的“强”形式与“弱”形式
3.1 一维杆的微分方程“强”形式
3.2 一维杆的积分方程“弱”形式
3.3 近似函数的连续性
3.4 “强”“弱”形式的等价性
3.5 算例
第4章 单元和插值函数的构造
4.1 近似函数的插值
4.2 2节点线性单元
4.3 3节点2次单元
4.4 直接构造形函数
4.4.1 总体坐标内的插值函数
4.4.2 自然坐标内的插值函数
第5章 加权余量法与虚功原理建立有限元格式
5.1 伽辽金加权余量法建立有限元格式
5.2 虚位移原理建立有限元格式
5.2.1 一维杆弹性力学基本方程
5.2.2 平衡方程的等效积分“弱”形式━━虚位移原理
5.2.3 虚位移原理推导平衡方程和面力边界条件
第6章 变分原理建立有限元格式
6.1 最小势能原理
6.2 最小势能原理推导杆单元平衡方程和面力边界条件
6.3 H-W变分原理
6.4 H-R变分原理
6.5 加权余量法与能量解法的源流关系
第7章 线性静态有限元分析
7.1 数值积分方法
7.2 单元刚度矩阵与载荷列阵的解析积分和数值积分
7.3 等参数单元
7.3.1 等参变换
7.3.2 等参变换与换元积分的关系
7.4 组装总体刚度矩阵与载荷列阵
7.5 矩阵压缩存储
7.6 优化结构的最小半带宽
7.6.1 基本概念
7.6.2 节点编号步骤
7.7 引入位移边界条件
7.8 求解线性方程组
7.9 杆单元节点位移、轴力、轴向应变、结构约束反力的求解
7.10 算例
7.10.1 2杆桁架解析算例
7.10.2 10杆桁架数值算例
第8章 线性动态有限元分析
8.1 基本方程
8.1.1 弹性动力学的基本方程
8.1.2 杆件离散化
8.1.3 有限元格式
8.2 质量矩阵
8.2.1 协调质量矩阵
8.2.2 集中质量矩阵
8.3 直接积分法
8.3.1 中心差分法
8.3.2 Wilson-θ方法
8.3.3 Newmark方法
8.3.4 数值算例
8.4 振型叠加法
8.4.1 求解固有频率和固有振型
8.4.2 求解模态坐标下的动力响应
8.4.3 2杆桁架模态分析的解析算例
8.5 大型特征值问题的解法
8.5.1 子空间迭代法
8.5.2 化广义特征值问题为标准特征值问题
8.5.3 Jacobi旋转法求解标准特征值问题
8.5.4 自由模态问题的移频
8.5.5 子空间迭代法的数值算例
第9章 几何非线性有限元分析
9.1 应变与应力的度量
9.2 完全拉格朗日格式
9.2.1 最小势能原理推导平衡方程
9.2.2 切线刚度矩阵
9.3 非线性方程组的解法
9.4 解析算例
9.4.1 求解单元切线刚度矩阵
9.4.2 组装切线刚度矩阵
9.4.3 求解内部力并组装
9.4.4 求解平衡方程
9.5 数值算例
9.5.1 2杆桁架结构
9.5.2 19杆桁架结构
第10章 材料非线性有限元分析
10.1 各向同性强化模型与随动强化模型
10.2 单向应力状态的弹塑性本构方程
10.2.1 屈服准则
10.2.2 强化法则
10.2.3 流动法则
10.2.4 一致性准则
10.3 弹塑性状态的决定
10.4 弹塑性分析的增量有限元格式
10.5 弹塑性有限元分析的基本流程
10.6 几种常见的等向强化模型
10.6.1 线性强化模型
10.6.2 二次强化模型
10.6.3 指数强化模型
10.6.4 Ramberg-Osgood强化模型
10.7 算例
10.7.1 3杆桁架解析算例
10.7.2 3杆桁架数值算例
10.7.3 10杆桁架循环加载算例
10.7.4 桁架桥梁结构算例
第11章 复合材料多尺度分析
11.1 摄动法
11.2 均匀化方法
11.3 算例
第12章 结构灵敏度分析
12.1 静态位移灵敏度分析
12.1.1 差分法
12.1.2 解析法
12.1.3 伴随变量法
12.1.4 半解析法
12.1.5 复数法
12.1.6 算例
12.1.7 各种算法的综合比较
12.2 应力灵敏度分析
12.3 模态灵敏度分析
12.4 应变能灵敏度分析
12.5 单元层次灵敏度分析
12.5.1 尺寸变量灵敏度分析
12.5.2 形状变量灵敏度分析
12.5.3 材料拓扑灵敏度分析
12.6 灵敏度分析算例
12.6.1 静态位移灵敏度分析
12.6.2 频率灵敏度分析
第13章 桁架结构有限元教学软件EFESTS
13.1 软件介绍
13.2 前处理
13.3 求解器
13.3.1 线性静态分析
13.3.2 几何非线性静态分析
13.3.3 材料非线性静态分析
13.3.4 模态分析
13.3.5 动力响应分析
13.4 位移灵敏度分析
13.4.1 解析法
13.4.2 伴随变量法
13.5 应力灵敏度分析
13.6 模态灵敏度分析
参考文献
附录1 算例数据
附录2 桁架结构有限元分析Matlab代码与实例
左文杰,1981年出生,吉林大学机械与航空航天工程学院教授,博士生导师。2001~2010年在吉林大学获得学士、硕士、博士学位;2014~2015年,美国密歇根大学机械工程系国家公派访问学者;2017年晋升吉林大学教授。从事结构优化设计研究,提出了完整的车身结构正向设计方法,解决了结构薄壁化设计难题。开发了具有自主知识产权的车身结构正向设计CAE软件CarFrame,得到了车企的认可与应用。主持国家自然科学基金2项,省部级与车企项目6项。在AIAA、IJVD等期刊发表SCI论文25篇、EI论文15篇;获发明专利4项、软件著作权10项;行业学术会议特邀报告6次;获省部级奖励2项;Google Scholar被引800余次,ESI高被引论文1篇。
汽车正向设计与轻量化团队
左文杰
zuowenjie@jlu.edu.cn
吉林省长春市人民大街5988号
吉林大学机械学院力学系