弹性力学与有限单元法试卷1

1、 至学年第学期出卷教师： 王建伟 适应班级：土木10级考试方式：开卷报告形式 本试卷分数考试占总评成绩的 80%姓名： 悦敢超 学号：201048040603第一题：以平面应力弹性力学问题为例，说明弹性力学的所研究问题的数学模型，并简单导出弹性力学位移法与应力法的数学模型。（20分）第二题：采用半逆解法求解下面薄壁梁（其参数为厚度t ，高度h ，弹性模量E ，泊松比等）的变形后的应力、应变和位移？（20分） 第三题：以平面应力弹性力学问题为例说明最小位能原理（能量法泛函极值）对问题的描述完全等价于第一题中的位移法描述（微分形式）。（20分）第四题：说明有限单元法的基本思想与基本步骤。要求：必须

2、准确说明采用有限元法解题的每一步（不必计算），如你认为有必要，可以结合例子说明。（20分）第五题（上网查资料）谈一谈你对有限元法的认识，说明有限单元法在工程上的应用，并举出两个以上的你有兴趣的工程实例。说明今后有限元法的发展趋势与方向。（20分）报告要求：1、报告中需图形表述的可以借用讲义中的图形，如果无法找到合适图形，可手绘于标准打印纸上，注意标注清晰。除图形外所有文字与公式必须打印。如不满足上述要求，无成绩处理。2、考勤不符合要求者，无成绩处理。3、交报告时将此试卷作为报告首页（填写姓名与学号），其它报告内容按顺序装订后提交。一、解：平面应力弹性问题的数学模型，就是8个方程应力边界条件和位

3、移边界条件，目标是求解8个未知（场）函数。位移法 二、解：（1）应力函数 （2）求应力几何方程物理方程（3）求应变 四、解：有限单元法的思想就是先对位能泛函做离散近似，然后在变分求得位移解。位移有限元法最重要就是要建立Ka=P这个有限元方程（本质为2n 个线性代入方程，后引入Su 条件求解即可），因而K和P的结成是关键。 具体解题步骤如下： 1单元剖分 把连续弹性体分割成许多个有限大小的单元，并为单元和节点编号。 2单元特征分析 以节点位移e 为基本未知量，设选一个单元位移函数，之后：（1）用节点位移表示单元位移，f=N e。（2）通过几何方程用节点位移表示单元应变，=B e。 （3）通过物理

4、方程用节点位移表示单元应力，=G e。（4）通过虚功方程用节点位移表示节点力，Fe=K ee，得出单元刚度矩阵。3总体结构合成（1）分析整理各单元刚度矩阵，通过节点的平衡方程形成节点载荷列阵、合成总体刚度矩阵，建立以节点位移为未知量的、以总体刚度矩阵为系数的线性代数方程：K =F。 （2）对线性代数方程组进行边界条件处理，求解节点位移。进而由=G e可求得单元应力。五、解：有限元法（finite element method）是一种高效能、常用的求解微分方程的计算方法。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值

5、问题有着紧密的联系）。其基本思想是由解给定的泊松方程化为求解泛函的极值问题。基本原理是将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达，从而使一个连续的无限自由度问题变成离散的有限自由度问题。岩土工程应用岩 土 工 程 较 为 复 杂， 这 一 方 面 是 因 为 岩 土 材 料 力 学 特 性 复 杂， 另 一方 面 也 因 为 地 质 结 构 复 杂， 且 具 隐 蔽 性， 无 法 完 全 勘 测， 从 而 导 致 岩 土 工 程 对 经 验 的 依 赖 性 非 常 强。

6、 尽 管 如 此， 有 限 元 法 在 岩 土 工 程 中 也 得 到 了 普 遍 的 应 用， 发 挥 的 作 用 越 来 越 大。 有 限 元 法 较 适 合 分 析 连 续 介 质 的 力 学 问 题， 而 岩 体 类 材 料 属 不 连 续 介 质， 其 中 存 在 断 层 和 众 多 的 结 构 面。 对 此 类 问 题， 人 们 提 出 了 一 些 特 殊 单 元 来 模 拟 它 们， 如 薄 层 单 元、Goodman 单 元 等， 较 好 地 解 决 了 这 一 问 题， 使 有 限 元 法 可 以 成 功 地 应 用 到 岩 土 工 程 中， 为 岩 土 工 程 的 设 计

7、和 施 工 提 供 了 很 好 的 理 论 分 析 方 法 和 计 算 手 段。再 如 在 地 下 工 程 中， 通 过 线 性 或 非 线 性 有 限 元 法 对 围 岩 进 行 结 构 分 析， 可 以 预 测 地 下 洞 室 开 挖 后， 围 岩 的 变 形、 应 力 重 分 布 情 况 及 洞 室 的 稳 定 性， 为 地 下 洞 室 的 可 行 性 研 究 提 供 理 论 根 据。 有 限 元 法 与 大 变 形、 流 变 和 损 伤 等 理 论 相 结 合， 还 可 研 究 洞 室 的 长 期 稳 定 性、 破 坏 规 律 和 过 程， 为 经 济、 科 学 地 设 计 有 效 支

8、 护 形 式 和 分 析 加 固 效 果 提 供 可 靠 参 数。 在 边 坡、 土 坝 和 地 基 等 工 程 的 研 究 过 程 中 也 较 多 使 用 了 有 限 单 元 法， 用 来 分 析 在 荷 载 作 用 下 的 变 形、 稳 定 和 加 固 等 问 题。 如 在 边 坡 的 开 挖 过 程 中， 稳 定 性 致 关 重 要。 过 去 由 于 受 计 算 手 段 的 限 制， 常 采 用 刚 体 极 限 平 衡 法， 往 往 过 于 保 守。 采 用 非 线 性 有 限 元 法 后， 可 以 分 析 整 个 边坡 的 变 形 状 态， 确 定 危 险 的 区 域 或 岩 块， 找

9、 出 可 能 的 滑 动 面， 计 算 相 应 的 稳 定 安 全 系 数， 还 可 分 析 加 固 效 果 和 施 工 程 序 对 稳 定 性 的 影 响。 由 于 岩 土 类 材 料 力 学 参 数 的 随 机 性 和 地 质 结 构 的 不 确 定 性， 人 们 在 有 限 元 法 中 引 进 了 概 率 理 论， 产 生 了 随 机 有 限 元 法， 其 结 果 更 具 参 考 价 值。有 限 元 法 在 求 解 具 体 问 题 时， 首 先 要 形 成 所 谓 的 网 格 信 息。 对 于 连 续 介 质 和 比 较 规 则 的 区 域， 目 前 的 有 限 元 法 软 件 都 能

10、较 方 便 地 自 动 生 成 网 格 和 坐 标 信 息。 对 于 岩 体 这 类 不 连 续 介 质 材 料， 其 众 多 的 结 构 面 对 网 格 的 生 成 带 来 了 很 大 的 麻 烦。 而 这 些 地 质 结 构 面 对 岩 体 的 力 学 性 能 有 直 接 的 影 响， 在 利 用 有 限 元 法 进 行 分 析 时 必 须 加 以 模 拟。 如 何 根 据 地 质 信 息， 由 计 算 机 自 动 形 成 网 格， 这 方 面 工 作 还 值 得 研 究。 为 避 免 上 述 困 难， 已 经 有 些 学 者 提 出 了 一 些 新 的 数 值 分 析 方 法， 如 块

11、体 单 元 法、DDA 和 无 单 元 网 格 法 等。机械工程应用有限元技术的出现为机械工程结构的设计、制造提供了强有力的工具，它可以解决许多以往手工计算根本无法解决的问题，为企业带来巨大的经济效益和社会效益。在现代机械工业中要设计生产出性能优越、可靠的机械产品，不应用计算及进行辅助设计分析是根本无法实现的，因此目前各生产设计部门都非常重视在设计制造过程中采用先进的计算机技术。近年来，国内外许多学者对机械零部件的有限元分析进行了大量的研究，归纳起来主要是以下几个方面：a 静力学分析。当作用在结构上的载荷不随时间变化或随时间的变化十分缓慢，应进行静力学分析。这是对机械结构受力后的应力、应变和变

12、形的分析，是有限元法在机械工程中最基本、最常用的分析类型。 b 动力学分析。机械零部件在工作时不仅受到静载荷作用，当外界有与其固有频率相近的激励时，还会引起共振，严重破坏结构从而引起失效。故零部件在结构设计时，对复杂结构，在满足静态刚度要求条件下，要检验动态刚度。C 热应力分析。这类分析用于研究结构的工作温度不等于安装温度时或工作时结构内部存在温度分布时，结构内部的温度应力。d 接触分析。接触分析用于分析两个结构物发生接触时的接触面状态、法向力等。由干机械结构中结构与结构间力的传递均是通过接触来实现的，所以有限元法在机械结构中的应用很多都是接触分析。 这是一种非线性分析， 以前受计算能力的制约

13、， 接触分析应用的较少。 e 屈 曲分析。 这是一种几何非线性分析，用于确定结构开始变得不稳定时的临界载荷 和屈曲模态形状，例如压杆稳定性问题。 有限元法最初应用在求解结构的平面问题， 发展至今已由二维问题扩展到三 维问题、板壳问题，由单一物理场的求解扩展到多物理场的耦合，由静力学问题 扩展到动力学问题、稳定性问题，由结构力学扩展到流体力学、电磁学、传热学 等学科，由线性问题扩展到非线性问题，由弹性材料扩展到弹塑性、塑性、黏弹 性、黏塑性和复合材料，从航空技术领域扩展到航天、土木建筑、机械制造、水 利工程、 造船、 电子技术及原子能等， 其应用的深度和广度都得到了极大的拓展。 有限元法的发展过程是与计算机技术的发展紧密相联的。只有计算机技术高度 发展以后，有限元法才得到广泛的应用。一个复杂的问题的求解，过去用小型机 花费几天才能得到结果， 现在用PC机几个小时就能完成同样的工作。商业有限元 软件也由只能在大中型计算机