岩土工程数值方法有限单元法_第1页
岩土工程数值方法有限单元法_第2页
岩土工程数值方法有限单元法_第3页
岩土工程数值方法有限单元法_第4页
岩土工程数值方法有限单元法_第5页
已阅读5页,还剩11页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、7.2.2 单元矩阵的变换 概述 在有限元分析中,为建立求解方程,需要进行各个单元体积内和面积内的积分,它们的一般形式可表示为而G和g中还常包含着场函数对于总体坐标x,y,z的导数。 由于在现在的情况下,场函数是用自然坐标表述的,又因为在自然坐标内的积分限是规格化的,因此我们希望能在自然坐标内并按规格化的数值积分方法进行上述积分。为此需要建立两个坐标系内导数、体积微元、面积微元之间的变换关系。 具体分为以下三种情况: 1. 导数之间的变换; 2. 体积微元、面积微元的变换; 3. 面积(或体积)坐标与笛卡尔坐标之间的变换。 导数之间的变换 按照通常的偏微分规则,函数 对 的偏导数可表示成: 对

2、于其它两个坐标 ,可写出类似的表达式,将它们集合成矩 阵形式,则有 上式中 称为Jacobi矩阵,可记作 ,利用7.2.2式, 可以显式地表示为自然坐标的函数 这样一来, 对于x,y,z的偏导数可用自然坐标显式地表示为其中 的逆矩阵,它可按下式计算得到 的行列式,称为Jacobi先烈式, 的伴随矩阵,它的元素 的元素 的代数余子式。 体积微元、面积微元的变换从图7.2可以看到 在笛卡尔坐标系内所形成的体积微元是设:定义 的混合积,且有其几何意思是以 为棱的平行立面体的体积,另外: 而 其中i,j,k是笛卡尔坐标x,y,z方向的单位向量。将7.2.10式代入7.2.9式,得到关于面积微元,例如在

3、 常数(c)的面上其它面上的 可以通过轮换 而得到。 称 矢量积, 的方向按右手法则确定。向量积的模 表示向量 为边所构成的平行四边形的面积。 在有了以上几种坐标变换的关系式以后,积分7.2.3和7.2.4式最终可以变换到自然坐标系的规则化域内进行,它们可分别表示成 其中对于二维情况,以上各式将相应蜕化,这里Jacobi矩阵是二个坐标之间的偏导数关系是 在笛卡尔坐标内形成的面积微元是 的曲线上, 在笛卡尔坐标内的线段微元长度是 面积(或体积)坐标与笛卡尔坐标 之间的变换 以上关于 等的公式原则上对于任何坐标和笛卡尔坐标之间的变换都是适用的,但是当自然坐标是面积或体积坐标时要注意两点: (1).面积或体积坐标都不是完全独立的,分别存在关系式; (2). 7.2.13,7.2.14等式的积分限应根据体积坐标和面积坐标特点,作必要的改变,这样一来,上述各式将成为。 (1).面积或体积坐标都不是完全独立的,分别存在关系式: 因此可以重新定义新的自然坐标,例如对于三维情况,可令 这样一来,7.2.5-7.2.12式形式上都保持不变, 也仍保持它的原来形式,只是它对 的导数应作如下替换 对于二维情况,则因为可令所以 (2). 7.2.13,7.2.14等式的积分限应根据体积坐标和面积坐标特点,作必要的改变,这样一来,上述各式将成为: 7.2.24式用于 的表面,类似

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论