2有限元法的基本概念与求解方法.ppt

总刚的奇异性 对任意的vi ，ui上式恒等于0，因此其系数分别等于0 即： 由对称性可知： 上式表示了总刚中各行元素之和均为零，即K对应的行列式的各 行线性相关，根据行列式的性质：“若行列式的各行线性相关，或某 一行是其余各行的线性组合，则行列式=0” 有： ?总刚是奇异矩阵 位移边界条件的处理 处理位移边界条件的常用方法 由于K和R中的各元素均已按照一定的顺序分别存储在相应的数 组中，在对K及R处理时，应尽量不打乱原有的存储顺序，并希望处 理的元素越少越好。常用的方法有3种。 1 .降阶法：降低结构刚度方程阶次的方法 若结构刚度方程为：K ? =R 在节点位移中，令?A——未知位移， ?B已知位移，利用矩阵分块： 位移边界条件的处理 处理位移边界条件的常用方法 式中RB为未知载荷，并按第一行展开： KAA ?A + KAB ?B =RA 令 则： （阶次比K ? =RA低） 若?B为零位移，则上式变为： 上式相当于在原结构刚度方程中，将与零位移约束对应的行与列 划去得到，由该式可以解出未知位移 ?A 采用计算机解题时，降阶会打乱原K及R的存储顺序，且需重新 安排KAA ，RA及?A ?在程序设计中一般不采用 位移边界条件的处理 处理位移边界条件的常用方法 2 .对角置一法： 对结构刚度方程 位移边界条件的处理 处理位移边界条件的常用方法 已知位移边界条件： ，将其引入刚度方程。 为了不改变列数，处理第i列诸元素。 位移边界条件的处理 处理位移边界条件的常用方法 为了改变行数，对第i行处理，使其体现 位移边界条件的处理 处理位移边界条件的常用方法 对于平面问题，只要进行三次即可使k成为非奇异阵，故可求出 待解位移。 位移边界条件的处理 处理位移边界条件的常用方法 3 .对角元乘大数法 若： 对第i行的主对角元Kii乘以一个大数，如1020，并将对应的Ri改为： 其它各行元素均保持不变。这样将第i行展开得到： 同除1020得： ? 位移边界条件的处理 处理位移边界条件的常用方法 若 则在Kii处乘以大数，Ri处置零即可 这种方法应用的最为普遍。 位移边界条件的处理 基本公式 由结构刚度方程解出?后可得到用? e单元内应变、应力与节点 位移间的关系。就可得出单元中任一点处的应变与应力： 当不考虑温度影响时： 或者 ?=B?e （B几何矩阵） 应力计算 载荷移置的普遍公式 推导：（１）根据单元内位移与节点位移关系　　 f=[u, v]Ｔ=N ?e ? f*= N ?*e ?P 　　　（２）根据能量等效原则： ?*e ＴRe =f*Ｔ?P ? ?*e ＴRe =(N ? ?*e ) Ｔ?P = ?*e Ｔ? NＴ?P ? Re = NＴ?P 这就是集中力Ｐ的移置公式，式中Ｎ为单元的形函数矩阵。 载荷移置与等效节点载荷 载荷移置的普遍公式 ２.体力g的移置公式 设：单元ijm上作用有体力　　g=[gx, gy ]Ｔ 推导（１）将单元体tdxdy上的体积力gdxdy当作集中力，应用集中 力的移置公式　Re = NＴ?P 有微元体上 　　　　　　　dRe = NＴ? g tdxdy 　　（２）积分在整个单元上有 Re = ??? NＴ? g tdxdy = t??? NＴ? g dxdy 载荷移置与等效节点载荷 载荷移置的普遍公式 3.表面力q的移置公式 设单元ijm的jm边上作用有表面力q=[qx，qy]Ｔ 可将微元面积上tds上的面力qtds当作集中力 则　Re = ??sjmdRe = ?sjm qt ds = t?sjm q ds 上述公式适用于任何单元及任意坐标方向。 载荷移置与等效节点载荷 载荷移置举例 　　以单元自重（或作用在单元形心处的集中力）为例。 　　设一个均质等厚的三角形单元ijm，其厚度为t，面积为?，材料 比重为?，则单元的自重为：W=t????,且其作用在单元形心c处 载荷移置与等效节点载荷 载荷移置举例 思路： 　　（１）欲求哪个节点在哪个方向上的载荷分量，就在该方向加 一单位虚位移，其他自由度为０。 　　（２）利用