有限单元法报告1.0.docVIP

《有限单元法》报告 姓名 张仁磊 班级 土建研1401 学号 第一题 以平面应力弹性力学问题为例，说明弹性力学的所研究问题的数学模型，并简单导出弹性力学位移法与应力法的数学模型。（10分） 第二题 以平面问题为例，推导应力边界条件。（10分） 第三题 采用逆解法求解下面薄壁梁（其参数为厚度t，高度h，弹性模量E，泊松比等）的变形后的应力、应变和位移？（20分） 第四题 以平面应力弹性力学问题为例说明最小位能原理（能量法－泛函极值）对问题的描述完全等价于第一题中的位移法描述（微分形式）。（20分） 第五题 通过例子详细说明有限单元法的理论和具体计算步骤。（20分） 第六题（上网查资料） 谈一谈你对有限元法的认识，说明有限单元法在工程上的应用，并举出两个以上的你有兴趣的工程实例。说明今后有限元法的发展趋势与方向。（20分） 报告要求： 报告中需图形表述的可以借用讲义中的图形，如果无法找到合适图形，可手绘于标准打印纸上，注意标注清晰。除图形外所有文字与公式必须打印。如不满足上述要求，无成绩处理。 考勤不符合要求者，无成绩处理。 交报告时将此试卷作为报告首页（填写姓名与学号），其它报告内容按顺序装订后提交。 第一题 解：（1）微元平衡微分方程（2个） 微元受力分析如右图 由切应力互等定理或→τxy=τyx 即 化简得 同理由 得 故可得两个平衡方程 物理方程（3个） 平面应力弹性力学问题应力应变关系式： （3）几何方程（3个） 同理 又因为 同理得 所以 即平面的几何方程为： 平面力学问题的几何条件 ?应力边界条件 由切应力互等定理或→τxy=τyx 同理由 平面弹性力学应力边界条件为： ?位移边界条件： 位移法的数学模型 思路：?变换物理方程，用应变来表示应力得方程?； ?将几何方程代入?得用位移来表示应力的方程?； ?将方程?代入平衡方程和应力边界条件即可得到位移法的数学模型。 （6）应力解法的数学模型 思路：应力未知函数有3个，因此应力解法的数学模型需要3个以应力表示的微分方程。已经有两个平衡方程差一个，即相容方程。 相容方程的推导： 故应力解法的数学模型为： 第二题 解：应力边界条件 由切应力互等定理或→τxy=τyx 同理由 平面弹性力学应力边界条件为： 第三题 解：（1）确定应力函数 前面我们求过类似模型，可知在如上坐标系下应力函数为 应力求解 应变求解 将应力函数代入几何方程即得位移场函数 （4）位移求解 应变代入位移函数 将上面求得的u，v代入中得 确定待定系数 位移约束条件： 放松边界条件 第第四题 解：（1）位能指标： W为外力做的功 为约束力做功；为主动力做功；为主动表面力做功；为主动力做功。 于是位能表达式 将上式中的U1用位移u,v表示可得宗量为u(x,y),v(x,y)的两宗量泛函 虚位移原理 给宗量以微小变化容易求的，但不好求。 对于变换成应变复合宗量层面得 最小位移原理 对于上式的前半部分由格林公式可知 其中L为区域Ω的边界，方向为逆时针方向，并且∠2=∠α，∠β+∠1=180°所以l=cosα=cos∠2,m=cosβ=-cos∠1 将上述结果代入中可得 整理可得 整理最后可得 由变分法中的引理可知上式等价于 故位移解法(微分型)和位移能量法(积分型)是等价的。 第五题 解：有限单元法分析的基本步骤：在工程或物理问题的数学模型（基本变量、基本方程、求解域和边界条件等）确定以后，有限单元法作为对其进行分析的数值计算方法的基本步骤如下： 位能的矩阵表示及离散化 位能是一个积分形式的能量泛函指标，我们对它进行离散(如果不离散，原问题仍然不可解)，离散的方式有很多种，最基本的是三角形离散。 (2) 单元分析 位移插值函数 通过插值的方式实现用节点位移表示单元内任意一点的位移。做位移插值，目标是实现用单元节点位移来表示单元内任意一点的位移。使单元所有的变形信息（应力、应变和位移）都决定于节点位移或者说由节点位移来表示。 建立单元位能的位移表示 3）对单元位能的表达式做进一步的处理 利用结构位移转化矩阵，建立局部单元节点位移同总体节点位移之间的