结构静力分析有限元软件设计与开发PDF电子书下载

第1章 绪论 1

1.1 有限单元法的产生及发展 1

1.1.1 概述 1

1.1.2 有限单元法的前期积累 1

1.1.3 有限单元法的思想火花 4

1.1.4 有限单元法的正式诞生 6

1.1.5 有限单元法的蓬勃发展 8

1.2 通用有限元软件的产生及发展 9

1.2.1 IBM 704计算机的推出 9

1.2.2 高级语言的问世 10

1.2.3 通用有限元软件的诞生 11

1.2.4 有限元软件产业的繁荣 12

1.3 中国的有限单元法及有限元软件 14

1.3.1 早期的艰难探索 14

1.3.2 改革开放后的进展 15

1.4 有限元软件开发的重要论著 16

参考文献 17

第2章 杆件结构有限单元法理论 20

2.1 基本概念及基本符号 20

2.1.1 有限单元法中的一些基本概念 20

2.1.2 单元杆端内力和杆端位移的表示方法 23

2.1.3 坐标变换矩阵 24

2.2 单元分析 26

2.2.1 单元分析之Ⅰ——局部坐标系中的单元刚度矩阵 26

2.2.2 单元分析之Ⅱ——整体坐标系中的单元刚度矩阵 28

2.2.3 单元分析之Ⅲ——单元荷载的等效节点荷载向量 31

2.3 整体分析 33

2.3.1 整体分析之Ⅰ——结构的整体刚度矩阵 33

2.3.2 整体分析之Ⅱ——结构的综合节点荷载向量 33

2.3.3 整体分析之Ⅲ——结构刚度方程的求解 33

2.4 单元再分析 35

2.4.1 单元再分析之Ⅰ——单元的最终杆端内力 35

2.4.2 单元再分析之Ⅱ——单元的内力变化方程 36

2.4.3 单元再分析之Ⅲ——单元的变形曲线方程 36

2.5 三维杆单元 37

2.5.1 三维杆单元的杆端内力向量和杆端位移向量 37

2.5.2 三维杆单元的坐标变换矩阵 38

2.5.3 三维杆单元局部坐标系中的单刚 38

2.5.4 三维杆单元整体坐标系中的单刚 39

2.6 三维梁单元 40

2.6.1 三维梁单元的杆端内力向量和杆端位移向量 40

2.6.2 三维梁单元的坐标变换矩阵 41

2.6.3 三维梁单元的第三个节点 42

2.6.4 三维梁单元局部坐标系中的单刚 44

2.6.5 三维梁单元整体坐标系中的单刚 45

2.7 杆件结构有限单元法计算步骤总结 46

参考文献 47

第3章 杆件结构内力分析程序 48

3.1 FRAME2D简介 48

3.1.1 程序总体结构及界面 48

3.1.2 有限元模型数据 50

3.2 FRAME2D的计算模块入口程序 57

3.3 读入有限元模型数据子程序 58

3.3.1 存储模型数据的重要全程变量和数组 58

3.3.2 读入有限元模型数据文件的子程序 59

3.4 计算前预处理子程序 61

3.5 编写节点总码子程序 63

3.6 确定总刚主对角线上元素地址子程序 68

3.7 数值计算模块总控子程序 71

3.8 计算总刚子程序 72

3.9 形成节点荷载向量子程序 74

3.1 0求解线性方程组子程序 79

3.1 1利用节点位移求单元杆端内力 81

3.1 2单元调度程序 82

3.1 3两节点杆单元模块 83

3.1 3.1 杆单元主控程序 83

3.1 3.2 杆单元的计算单刚子程序 84

3.1 3.3 强制节点位移在杆单元中引起的等效节点荷载 87

3.1 3.4 杆单元上单元荷载产生的等效节点荷载 88

3.1 3.5 计算杆单元局部坐标系中杆端内力子程序 91

3.1 4两节点梁单元模块 92

3.1 4.1 梁单元的主控程序 92

3.1 4.2 梁单元的计算单刚子程序 93

3.1 4.3 强制节点位移在梁单元中引起的等效节点荷载 97

3.1 4.4 计算梁单元上单元荷载的等效节点荷载子程序 98

3.1 4.5 计算梁单元杆端内力子程序 101

3.1 5算例 102

3.1 5.1 算例1 102

3.1 5.2 算例2 110

参考文献 113

第4章 连续体结构有限单元法理论 114

4.1 二维弹性力学问题有限元分析基础 114

4.1.1 二维问题的一些基本概念 114

4.1.2 二维问题的一些基本方程 116

4.1.3 二维连续体的势能及有关定理 119

4.1.4 平面单元积分形式的单刚 121

4.1.5 平面单元的等效节点荷载 121

4.1.6 二维连续体的整体刚度方程 122

4.2 三节点三角形单元 123

4.2.1 三节点三角形单元节点力向量和节点位移向量 123

4.2.2 三节点三角形单元的形函数 123

4.2.3 面积坐标 124

4.2.4 三节点三角形单元的应变矩阵 124

4.2.5 三节点三角形单元的应力矩阵 124

4.2.6 三节点三角形单元单刚的显式表达式 125

4.2.7 三节点三角形单元的等效节点荷载 125

4.2.8 局部节点应力平滑技术 126

4.2.9 单元内的应力分布 126

4.2.1 0应力等值线算法 127

4.3 四节点四边形单元 127

4.3.1 四节点四边形等参单元的形函数 128

4.3.2 等参单元的一些基本变换 129

4.3.3 四节点四边形等参单元的单刚 131

4.3.4 高斯积分 131

4.3.5 四节点四边形等参单元的等效节点荷载 133

4.3.6 单元应力及最佳应力点 134

4.4 三维弹性力学问题有限元分析基础 134

4.5 四节点四面体单元 136

4.5.1 基于直角坐标系的四节点四面体单元形函数 136

4.5.2 四节点四面体单元的一些基本变换 137

4.5.3 四节点四面体单元的应变矩阵 139

4.5.4 四节点四面体单元的应力矩阵 139

4.5.5 四节点四面体单元的显式单刚 139

4.5.6 四节点四面体单元的等效节点荷载 140

4.6 八节点六面体单元 140

4.6.1 八节点六面体单元的形函数 140

4.6.2 八节点六面体单元的应变矩阵 141

4.6.3 八节点六面体单元的雅可比矩阵 142

4.6.4 八节点六面体单元的单刚 142

4.7 全域应力平滑技术 142

参考文献 143

第5章 连续体结构应力分析程序 144

5.1 SOLID2D的计算模块入口程序 144

5.2 SOLID2D的计算模块主控程序 145

5.3 SOLID2D的计算总刚子程序 146

5.4 SOLID2D的单元调度程序 147

5.5 三节点三角形单元模块 148

5.5.1 三节点三角形单元主控模块 148

5.5.2 计算三节点三角形单元单刚子程序 149

5.5.3 计算三节点三角形单元应力子程序 151

5.5.4 计算三节点三角形单元光滑矩阵和光滑向量子程序 155

5.6 四节点四边形单元模块 157

5.6.1 四节点四边形单元主控模块 157

5.6.2 计算四节点四边形单元单刚子程序 158

5.6.3 计算四节点四边形单元高斯点处应力子程序 162

5.6.4 计算四节点四边形单元光滑矩阵和光滑向量子程序 165

5.7 计算六节点三角形单元单刚子程序 170

5.8 计算十节点三角形单元单刚子程序 174

5.9 计算八节点四边形单元单刚子程序 180

5.1 0 SOLID3D程序 185

5.1 0.1 SOLID3D简介 185

5.1 0.2 计算八节点六面体单元单刚子程序 186

5.1 0.3 计算八节点六面体单元应力子程序 192

5.1 1算例 199

5.1 1.1 算例1 199

5.1 1.2 算例2 201

5.1 1.3 算例3 202

参考文献 205

第6章 温度场问题有限单元法理论 206

6.1 传热学基本理论 206

6.1.1 傅里叶定律 206

6.1.2 二维问题导热微分方程 206

6.1.3 三维问题导热微分方程 207

6.1.4 初始条件和边界条件 208

6.2 二维温度场有限元分析理论 209

6.2.1 二维温度场的泛函 209

6.2.2 二维温度场问题有限元方程的推导 210

6.3 三节点三角形热单元 213

6.3.1 平面求解区域的离散 213

6.3.2 三节点三角形热单元的形函数 213

6.3.3 三节点三角形热单元的导热刚度矩阵 214

6.3.4 三节点三角形热单元的蓄热刚度矩阵 214

6.3.5 三节点三角形热单元的边界贡献矩阵 215

6.3.6 三节点三角形热单元的温度荷载向量 217

6.4 四节点四边形热单元 218

6.4.1 四节点四边形热单元的形函数 218

6.4.2 四节点四边形热单元的导热单刚 219

6.4.3 四节点四边形热单元的蓄热单刚 219

6.4.4 四节点四边形热单元的边界贡献矩阵 219

6.4.5 四节点四边形热单元的温度荷载向量 222

6.5 六节点三角形热单元 223

6.5.1 六节点三角形热单元的形函数 223

6.5.2 六节点三角形热单元的导热单刚 224

6.5.3 六节点三角形热单元的蓄热单刚 225

6.5.4 六节点三角形热单元的边界贡献矩阵 225

6.5.5 六节点三角形热单元的温度荷载向量 227

6.6 二维热阻单元 227

6.6.1 二维热阻单元导热单刚的积分公式 227

6.6.2 四节点矩形热阻单元导热单刚的显式 229

6.6.3 六节点矩形热阻单元导热单刚的显式 230

6.7 三维温度场有限元分析理论 230

6.7.1 三维温度场的泛函 230

6.7.2 三维温度场问题有限元方程的推导 231

6.8 八节点六面体热单元 233

6.8.1 八节点六面体热单元的形函数及其对自然坐标的导数 233

6.8.2 八节点六面体热单元的雅可比矩阵及B矩阵 234

6.8.3 三维曲面面积微元的变换公式 235

6.8.4 八节点六面体热单元的导热单刚 236

6.8.5 八节点六面体热单元的蓄热单刚 237

6.8.6 八节点六面体热单元的边界贡献矩阵 237

6.8.7 八节点六面体热单元的温度荷载 237

6.9 三维热阻单元 237

6.9.1 三维热阻单元导热单刚的积分公式 237

6.9.2 特殊三维热阻单元导热单刚的显式 239

参考文献 241

第7章 温度场分析程序 242

7.1 二维温度场分析程序TEMP2D 242

7.1.1 TEMP2D简介 242

7.1.2 TEMP2D的模型数据 242

7.2 TEMP2D的计算模块入口程序 243

7.3 向后差分法数值计算模块总控子程序 244

7.3.1 向后差分算法 244

7.3.2 向后差分算法的总控子程序 245

7.4 计算导热总刚子程序 248

7.5 计算蓄热总刚子程序 249

7.6 形成总体温度荷载向量子程序 250

7.7 三节点三角形热单元模块 251

7.7.1 三节点三角形热单元主控子程序 251

7.7.2 三节点三角形热单元计算导热单刚子程序 252

7.7.3 三节点三角形热单元计算蓄热单刚子程序 255

7.7.4 三节点三角形热单元计算荷载向量子程序 256

7.8 四节点四边形热单元模块 258

7.8.1 四节点四边形热单元计算导热单刚子程序 258

7.8.2 四节点四边形热单元计算蓄热单刚子程序 261

7.9 四节点矩形热阻单元模块 263

7.10 三维温度场分析程序TEMP3D 265

7.10.1 计算八节点六面体热单元导热单刚子程序 265

7.10.2 计算八节点六面体热单元蓄热单刚子程序 273

7.10.3 计算八节点六面体热单元温度荷载向量子程序 276

7.11 算例 279

7.11.1 算例1 279

7.11.2 算例2 281

7.11.3 算例3 282

参考文献 287

第8章 前处理技术 288

8.1 二维网格自动剖分的映射方法 288

8.2 二维网格自动剖分子程序 291

8.3 确定绘图区坐标系统 296

8.3.1 确定整体坐标系的方法 296

8.3.2 建立整体坐标系的子程序 298

8.4 二维网格的显示、填充、放缩及移动 300

8.4.1 二维有限元网格的显示 300

8.4.2 二维有限元网格的填充 301

8.4.3 二维有限元网格的放缩及移动 302

8.5 交互编辑材料号和边界号 303

8.5.1 编辑单元材料号 303

8.5.2 编辑温度场问题的边界号 304

8.5.3 编辑力学问题的边界号 305

8.5.4 绘制矩形选择框子程序 305

8.6 交互编辑荷载 309

8.6.1 编辑节点荷载 309

8.6.2 编辑单元荷载 309

8.7 三维网格的坐标变换 310

8.8 三维网格的消隐计算 312

8.8.1 一些基本概念 313

8.8.2 自遮挡隐藏面的确定 314

8.8.3 最终可见线段的确定 315

8.9 深度优先级排序表的确定 317

参考文献 318

第9章 后处理技术 320

9.1 画云图方法概述 320

9.2 颜色与场量的关系 321

9.2.1 颜色模型 321

9.2.2 颜色与场量的函数关系 322

9.3 扫描母元法原理 323

9.3.1 扫描母元法理论基础 323

9.3.2 牛顿－拉夫逊法求自然坐标的过程 325

9.3.3 画单元云图的过程 326

9.4 画数据场云图的主要源程序 327

9.5 用扫描母元法画三角形单元云图 331

9.5.1 画三角形单元云图的原理 331

9.5.2 画三角形单元云图的过程 333

9.6 三维数据场可见表面上云图 336

9.7 六面体单元任意剖面上云图 337

9.7.1 六面体等参单元插值公式 337

9.7.2 坐标变换及单元包围盒 337

9.7.3 相关单元的判断及点到平面的距离 338

9.7.4 确定像素的自然坐标（r,s,t）的方法 339

9.7.5 画单元任意剖面上云图的步骤 339

9.7.6 应用实例 340

参考文献 341

第10章 有限元理论与软件发展展望 342

参考文献 344