有限元基础-讲-第2讲-之2.pptVIP

? Copyright 2000 PTC. Visit Visit ? Copyright 2000 Parametric Technology Corporation 板材单向拉伸试验右图所示。试验的目的是确定载荷与位移量的关系曲线，从而得到材料的应力应变关系。 选取一个典型的单元来进行单元的分析。首先要建立用单元结点位移表示单元内部任一点位移的关系式。设单元相关的结点号为i、j、m。 如果用u和v表示单元内部任一点的位移，它们都应该是x、y的函数，可以用对结点位移插值的方法来确定。这里选择最简便的线性函数，即： 确定α1,β1,α2,β2,α3,β3 已知各结点坐标为（xi，yi），（xj，yj），（xm，ym）将它 们代入 Thanks for your listening 二、平面问题 一般性概念 1 关于单元类型与位移模式 图示4结点四边形单元的位移模式应取： 通过该单元的四个结点，定出式中8个系数α1、β1 、… α4、β4 。而不单独取x2或y2作二次项，以保证几何各向同性。 二、平面问题 一般性概念 1 关于单元类型与位移模式 巴斯卡三角形 一般地说，用多项式位移模式处理平面问题应按照巴斯卡三角形中的对称原则来选择其项数，以保证x和y的项数（也包括幂次数）对等。 二、平面问题 一般性概念 1 关于单元类型与位移模式 例如： 平面8结点单元可选1、x、y、x2、y2、xy、x2y、xy2这8项，或1、x、y、x2、y2、xy、x3、y3； 9结点的单元则可选1、x、y、x2、y2、xy、x3、y3、x2y2这9项，其余情况依次类推 在空间问题中，同样地选择合理的项数，保证其位移模式对x、y、z三个坐标的对称性。 巴斯卡三角形 二、平面问题 一般性概念 2 整体刚度矩阵的性质 整体刚度矩阵 二、平面问题 一般性概念 2 整体刚度矩阵的性质 性质一：[K]是对称阵 由于刚度矩阵[K]的对称性，使实际计算时只需求出主对角线及其以上（或以下）的一半元素，节省了近一半的计算工作量。 性质二：主元恒正 这是因为结点上的力和位移方向一致的关系。 性质三：[K]是稀疏阵 有利于节约计算机的存贮单元，并缩短计算的时间。 二、平面问题 一般性概念 2 整体刚度矩阵的性质 性质四：[K]是奇异阵，排除刚体位移后[K]是正定阵 因为在{R}力的作用下，结构处于平衡状态，对于平面问题应满足三个平衡方程。反映在[K]中就存在三个线性相关的行列式，因而它是奇异阵，不存在逆阵。 在排除刚体位移后，[K]将是正定阵，这是因为： 式中：n—单元数 [K]e为单元的刚度矩阵 二、平面问题 一般性概念 2 整体刚度矩阵的性质 性质四：[K]是奇异阵，排除刚体位移后[K]是正定阵 两边同乘 代 入 二、平面问题 一般性概念 2 整体刚度矩阵的性质 性质四：[K]是奇异阵，排除刚体位移后[K]是正定阵 由于弹性矩阵[D]是正定的，而且t和△都是正的，所以只有在每个单元中都有{ε}=0时才有 ，否则它大于零。也就是说，排除了刚体位移{ε}=0之后，若{q}≠0，则二次型 恒大于零，于是[K]必定为正定阵。 二、平面问题 作 业 试分析方板对角拉伸过程中各结点的位移（分2个单元） 已知： E=2×105MPa=2×105N/mm2 μ=0.3 t=1mm P=10kN (1,0) ① ② (-1,0) (0,-1) (0,1) ? Copyright 2000 PTC. Visit 下次内容 2.1 有限元基本方法 2.2 平面问题 2.3 轴对称问题 2.4 空间问题 2.5 等参元和板壳元 * ? Copyright 2000 PTC. Visit 第二节 2.1 有限元基本方法 2.2 平面问题 2.3 轴对称问题 2.4 空间问题 2.5 等参元和板壳元 二、平面问题 实例—板材单向拉伸 1 结构离散化 利用对称性，取试件的四分之一为分析对象，坐标系可选 在试件的对称中心处。可采用三角形单元，也可以采用矩形单元。如图所示，将试件有效部分划分成16个三角形单元，通过设置在其上的15个结点连线分割而成。最后将单元和结点分别依次编号，即完成结构离散过程。 二、平面问题 实例—板材单向拉伸 2