有限元课程设计..doc

数学与统计学院信息与计算科学系 《有限元课程设计》 姓 名： 学 号： 专 业： 班 级： 指导老师： 时 间： 成绩评定 指 导 教 师 评 语 与 成 绩 学号 姓名 成绩 指导教师： 时 间： 设计题目： （1）平面应力问题上的单元节点位移求法 （2）双曲方程满足Dirichlet的有限元方法解答过程 摘要: 随着计算机技术的飞速发展，编制高效的程序求解各种偏微分方程问题已经成为可能，偏微分方程数值解法也成为科学和工程计算中的重要分支。然而，对于广大应用工作者来说，从偏微分方程模型出发，使用有限元法或有限差分法求解都要经过诸多步骤，耗费很大的工作量，才能得到数值解。进一步将结果可视化也是一个迫切需要解决的问题。目前，MATLAB?中的PDE toolbox已实现了二维偏微分方程定解问题的准确高效求解，人们只需使用图形用户界面或编写M文件，即可显示解的图形或输出解的数值。MATLAB为求解偏微分方程带来了极大方便，并在此基础上可以解决更为复杂的问题，偏微分方程在科学和工程上有着广泛的应用。有限元法是一种重要的偏微分方程数值解法。编程实现从偏微分方程到有限元求解全过程需要很好的理论基础和编程技巧,难度较高,PDE toolbox实现了偏微分方程的有限元解法。 关键词：MATLAB，PDE toolbox，有限元 问题背景： 如下图所示的结构受集中力P作用，试按平面应力问题分析，用3节点三角形单元，求图中结构2，4节点位移，取E为常量，μ =1/6，t=1。 2.详细推导出 双曲方程满足Dirichlet的有限元方法解答过程，并求解 正方形区域上的波动方程 ， 初始条件 ，， 边界条件 ，Dirichlet条件，，Neumann条件。 详细解答 第一题： 第一步：对单元①，i,j,m对应节点3,2,1 A=1/2 几何矩阵B为： 弹性矩阵D为： 利用公式： 求得： 从而求得： 单元贡献矩阵 注意：这里的上标代表单元的号码。 第二步：对于单元②，i,j,m对应节点3,4,2. A=1/2 几何矩阵B为： 利用公式： 求得： 从而求得： 单元贡献矩阵： 第三步：求得总刚度矩阵： 有了总体刚度矩阵后，再形成载荷矩阵，即可得整体刚度方程，约束处理后就可求解节点位移。 载荷矩阵为： 位移列阵为： 形成整体平衡方程： 处理后得到： 解以上；联立方程得 第二题 正方形区域上的波动方程 ， 初始条件 ，， 边界条件 ，dirichlet条件，，Neumann条件。 推导： 有限元法的基本步骤： [1] 把微分方程定解问题转化为变分形式 [2] 选定单元的形状，对求解区域做剖分 [3] 构造基函数或单元形状函数 [4] 形成有限元方程 [5] 求解有限元方程 [6] 收敛性及误差分析 考虑热传导方程的初边值问题 (1) (2) (3) 其中，G是光滑边界的平面有界域。 [1] 把微分方程定解问题转化为变分形式 设对固定的t0关于(x,y)属于以乘(1)两端并积分，得 利用Green公式和边值条件(3)，得 引进双线性形式和内积： (6) 则初边问题(1)-(3)的变分形式为：求，满足 [2] 对区域G作网格剖分 在中取一n维子空间, 从而Galerkin方程为 (9) (10) [3] 构造