平面高阶单元.ppt

  1. 1、本文档共19页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
3.1 四结点矩形单元 (1) 位移模式选择:两个,每个四项, (2) 形函数计算 在式(3-1)中代入结点位移和结点坐标后,可解出待定系数。将这些系数再代入式(3-1),可得形函数: (3) 单元应变 (4) 单元应力: (5) 单元刚度矩阵 把[B]、[D]代入上式,整理后可得: (6)等效结点荷载 单元的体积力和表面力引起的结点力仍可用式(1-36)和(1-37)进行计算。由于位移分量在x为常数及y为常数的直线上是线性变化的,因此,载荷向结点的分配也符合静力等效的原则。 (7) 整体平衡方程 根据各单元的刚度矩阵[k]、等效结点力列阵,按对号入座的方式叠加组装整体刚度矩阵和结点荷载列阵,从而得到整体平衡方程: 3.2 六结点三角形单元 1、位移模式 在三角形单元i、j、m的各边中点上,各增设一个结点,使每个单元具有6个结点,则得到图示的六结点三角形单元。这种单元具有12个自由度,可以采用完全二次多项式的位移模式: 收敛条件:所取位移模式反映了单元的刚体位移和常应变。另外,在单元边界上位移分量是按抛物线变化的,而每条公共边界上有3个公共结点,正好可以保证相邻两单元位移的连续性。因此,上述位移模式满足收敛的完备协调性条件。 根据假定的位移模式,按照前面的方法和过程,可以确定形函数[N]、应变矩阵[B]、应力矩阵[S]、单元刚度矩阵[k]、等效结点荷载列阵{f},然后组装总刚和荷载列阵。但是,过程非常复杂。采用采用面积坐标可以大大简化计算。 2、面积坐标 三角形单元ijm中,任一点P(x,y)的位置,可以用如下的三个比值来确定: (1)面积坐标与直角坐标的关系 三角形Pjm的面积为 (2)面积坐标的导数公式 根据面积坐标与直角坐标的关系,由复合函数的求导公式,有 (3)面积坐标的积分公式 在三角形单元上进行积分时,有 在三角形某一边(设边ij,边长l 为)上进行积分时,有 3. 六结点三角形单元计算公式 对于六结点三角形单元,形函数可用面积坐标表示为 可以利用形函数的性质检验式(3-29)的正确性 。 (在i点等于1,在其他结点等于0 ,见下) 形函数检验 单元位移、应变 单元应力、单元刚度、等效结点荷载 单元重力引起的等效结点力: 它表示各边中点承担单元重力的1/3。 它表示边中点承担2/3,载荷集度大的角结点承担1/3。 六结点三角形单元中的应变、应力不为常量,因此可以应用于应力梯度较大的地方,精度较高。 但其计算复杂。 * 上述的阶数为二次型,比三角形单元高,故可以更好地反映物体的应力、应变状态,得到更高的精度。它反映了单元的刚体位移和常应变,在单元的边界上位移是按线性分布的,因此,相邻单元在公共边上的位移是连续的。这样,位移模式满足了解答收敛性的充分必要条件。 第三章 平面高阶单元 [B]是ξ、η的函数,即是x,y的函数。因此单元中的应变不再是常数。 引入位移约束条件,解上述线性方程组可得结点位移,进而可求各单元应力。 四结点矩阵单元采用较高阶的位移模式,具有比三结点三角形单元较高的计算精度。但矩形单元也有缺点,一是不能适应斜线及曲线边界,二是不便于采用大小不同的单元。 面积坐标为 将上式中的3个式子分别乘以xi、xj、xm 并利用ai、bi、ci得 同理 等效结点荷载:由于位移模式是非线性的,因此体积力和表面力引起的结点力向量不能采用静力等效原理进行分配,而应采用相应公式进行计算。 单元体积力引起的结点力: 单元表面力引起的结点力 设在ij边上受有x方向的均匀分布力ps,则 它表示边中点分担4/6,两角点各承担载荷的1/6。 设在ij边上受有x方向的三角形分布力ps,其集度在i点为ps ,在j点为0。则 *

文档评论(0)

139****6621 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档