原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
3.4 二维有限元解 这里以典型的二维Dirichler问题为例，论述二维有限元法的基本步骤。 [问题的提法] 求二维Ω域内位函数u(x,y)的解 其边界条件是 Dirichlet条件 对应的能量泛函 和变分提法 3.4.1 区域划分 首先，将域Ω共剖分为若干个三角形单元，其中结点数p分成两类： 引入未知函数ut矩阵 (5-18-97) 十分清楚，上式中uI是域内函数，而uII是边界函数。这时，能量泛函可写成 (5-18-98) 典型单元称之为e，共有ne个三角形单元。其中 3.4.2 典型单元e分析 设典型单元e结点为i、j、m时对应的结点参数为ui、uj、um，于是典型单元的位函数矩阵有 (5-18-101) 而i、j、m单元内任一点(x,y) 的位函数u(x,y)可表示为 (5-18-102) 进一步有 (5-18-103) (5-18-104) 十分清楚，位函数的导数只作用在N函数上。因此，在单元e 中 (5-18-105) 其中 (5-18-106) 另一方面 (5-18-107) 根据单元能量泛函定义式，把其中的u(x,y)用p(x,y)逼近，有 (5-18-108) 若进一步引入 则可得到 的简洁矩阵表式 (5-18-111) 3.4.3 和 的计算公式 上述引入的k(e)计算比较简单，因为被积矩阵s(e)是常数矩阵。显见 (5-18-112) 关于r(e)矩阵计算可以由ri (e)分量来举例说明，即 (5-18-113) 式中，ρ表示电荷密度，在所求的问题中是已知的，ρ在单元e中的计算可以有以下几种近似情况。 情况1：常数元近似 设 (5-18-114) 这时 (5-18-115) 由附录可知 (5-18-116) 最后得到 (5-18-117) 情况2：线性元近似 这种情况下，认为ρ在单元e中近似作线性变化，也即 (5-18-118) 把式(5-18-118)代入ri (e)定义式(5-18-113)可得 (5-18-119) 同样由附录可知 于是得到 (5-18-122) 于是电场储能亦分为两部分 (5-18-28) 其中 (5-18-29) 则写出 (5-18-30) 于是电场总储能ω可表示为 (5-18-31) 在问题中μ0已知，故ω是μ1的一元函数。于是ωmin要求 (5-18-32) 容易解出 (5-18-33) 这是位函数分布 (5-18-34) 明显看出，这种情况下位函数呈直线分布，且可统一写出 (5-18-35) 特别应指出，式(5-18-35)正是此问题的精确解。 3.2 一维有限元法 通过研究[例3]的推广来讨论一维有限元法，把x0和xN分为N段，如图5-18-4所示。 图5-18-4 电位u(x)的分段展开函数 设两端边界条件已知 (5-1
文档评论（0）