基于切应力条件的广义协调等参元.pdfVIP

2005拄 6月 郑州大学学报(工学版) Jun．2005 of V01．26No．2 第26卷第2期 JournalZhengzhouUniversity(EngineeringScience) 文章编号：1671—6833(2005)02—0092—03 基于切应力条件的广义协调等参元 秦力一，许德刚，周爱民 (郑州大学工程力学系，河南郑州450002) 摘要：非协调位移模式可有效改善计算精度，但对任意非规则网格常不能满足收敛性条件．根据线性 边界力状态下的广义协调条件，对任意不规则四边形单元构造出一种广义协调等参元．该单元的导出考 虑了线性切应力条件，使单元间位移协调条件在加权残数意义上得以满足．进一步给出了相应的广义协 调应变形函数矩阵，模式紧凑而不易出现奇异性．该等参元对任意不规则四边形网格能通过分片检验， 当单元为平行四边形时蜕化为Q6元．算例表明，该类型等参元精度较高，对不规则网格剖分能保持良好 的数值性态． 关键词：广义协调；等参元；线性切应力 242．21 中图分类号：TB21；0 文献标识码：A 式中：Ⅱ。为协调元相应的能量泛函；T为单元边 0 引言 界力函数；F。和Pi分别为给定的体积力和面力． 为改进平面四结点等参元(Q4元)性态，Wil—用加权残值法可使单元间位移协调条件在某种积 son等引进非协调位移模式(Q6元)，与内部自由分意义上得以满足： 度相关联，对弯曲问题的计算精度显著改善，其位 4)T7(M一五)ds：0(2) 移模式由两部分组成： J孙。 {u}={u}+{五}(1) 为保证单元收敛并保证计算精度，考虑如下 其中协调位移为双线性插值，非协调位移模拟弯 线性边界力状态： 曲引起的形状变化．对非平行四边网格，06元不 能通过分片检验．文献[1]和文献[2]中，又分别根 据常应力分片检验的强形式和广义协调条件导出 式中：凡。，11,，表示单元边界外法线的方向余弦；常 非协调元NQ6和广义协调等参元GC—Q6． 应力项(^，兀，f5)以保证单元通过小片检验，线性 本文作者仍应用常应力和线性应力下的广义 项为提高计算精度，f6叩和厶善二项相应于线性切 协调条件，导出另外一种广义协调等参元．在单元 应力状态；其余各项(A)略去．将式(3及式(1)代 导出过程避开形心位移条件，而补充以边界线性 人式(2)后，注意到协调位移满足单元边界位移条 切应力广义协调条件．该非协调模式形式紧凑且 件M一“=0(在3A，域)，可得 不易_出现奇异性． 中：儿dS=0 J aA 1一线性边界应力下的广义协调条件 9)a^inyridS=o 记M和面为单元内部和边界位移变量．非协 调元不保证单元边界上处处满足条件M一¨=0， 9)。^u4nxridS_0 其相应的能量泛函应修正为[3] 4)锄。}d5=0 (4)