数值模拟技术3xin.ppt

注意:(1）上式中[K]e是平面问题中计算单元刚度矩阵的一般形式；（2）对常应变三角形单元，[K]可简化为： 将[B]、[D]代入（5）式中，可得到单刚矩阵的具体表达式，[K]e可写成分块矩阵形式： 2.单元刚度矩阵的物理意义 将（4）式写成下面形式： 单刚矩阵的2*2阶子矩阵[Krs]e表示节点s（s=i,j,m)产生单位位移时，在节点r（r=i,j,m)上所需施加的节点力的大小。 将（8）式中各列阵进一步展开，则有： 其中各个元素的物理意义为： 3.单元刚度矩阵的特性 （1）对称性： （2）奇异性： 物理意义：若仅知节点力，不知约束条件，则单元除产生本身的形变外，还可产生任意的刚体位移。这部分刚体位移分量按给定的节点力是无法唯一确定的。 （3）分块性。 五、荷载向节点的移置—单元荷载列阵 移置原因：有限元讨论的是节点的受力、平衡和位移。 移置原则：静力等效原则。 静力等效原则： （1）对于弹性体，原荷载与节点荷载在任何虚位移上的虚功都相等。 （2）原荷载与节点荷载在任一轴上的投影之和及对任一轴的力矩之和都相等，具相同的主矢量和主矩。 荷载分类：（1）集中力；(2)体积力：重力、地震力 (3)面力：构造应力、自重应力、风力。 1、集中力的移置 如图，将集中荷载移置到节点上后形成的节点荷载列阵用{R}e表示： 假设(1)M点发生的虚位移： (2)各节点与之相应产生虚位移为 由静力等效原则： 两边同乘 则上式变为： 2、体力的移置 设单元内分布体力的密度为 在单元内取一小的微分体，体积为tdxdy,t为厚度 将{r}tdxdy当作一集中力P，由（10）有： 3、面力的移置 图示，设在单元ij边上作用分 布面力 在ij边上取微 段tds，作用于此段上的荷载为 将其作为集中荷载，则可由(10)式 积分得： 4、线性位移模式下的荷载移置 对于本章讨论的三节点三角形单元线性位移模式情况，以上的通用公式可进一步简化。 （1）重力的移置:单元为均质等厚的三角形单元，W为作用于三角形形心的集中荷载。 首先求W移置到节点i的垂直荷载Pyi.如图3-14(a). 设j,m不动，i有垂直向下的单位虚位移。由于bc/bi=1/3，由三角形法则可知c的虚位移为1/3，由静力等效原则 同理可得： 其次，求重力向节点移置 的水平荷载，设节点i产生单位水平虚位移，如图3-14(b),此时形心c所产生的虚位移 cd为水平向，与W垂直，故W所作的虚位移功为0，即： （2）边界上集中荷载的移置 如图三角形一边受水平向集中 荷载p。 首先,求向i的水平移置荷载Rxi 设i产生单位虚位移，j,m不动。 由相似三角形原理，集中荷载 作用点产生水平位移为Li/L,由虚功原理： (3)分布荷载的移置 首先，求合力，然后，进行移置。 如图，ij边受三角形分布荷载作用，其合力为： ，其作用点位置为： 将其代入（12）式，则节点 荷载列阵为： 对于任意荷载可先分解为X， Y方向，然后利用上述结果，由 叠加原理求得节点荷载列阵。 第二章 平面线弹性问题的有限单元法 第一节 有限单元法的基本思路 第二节 位移模式及单元分析 第三节 整体分析—整体刚度矩阵及其修正 第四节 小结及算例 第一节 有限单元法的基本思路 平面线弹性问题的有限单元法是针对平面问题，依据弹性力学的基本方程所构造的有限单元法。（Finite Element Method_FEM) 有限单元法的基础：结构力学 对于岩土体这样的连续介质，怎样用结构力学的分析方法来解决它的弹性力学问题？ 其步骤如下： 1、将连续体离散化：即将连续的弹性体变换成一个离散的结构物（构件）。这个结构物（构件）是由有限多个有限大小的构件，在有限个连续点上联系起来而构成的。（如图） 其中：这有限大小的构件就称为单元； 有限个连