第八章平面应力问题的有限单元法详解.pptx

船舶结构力学;§8-1弹性体的应力、位移与应变;2012-12-29;2012-12-29;2012-12-29;2012-12-29;2012-12-29;2012-12-29;2012-12-29; 1.平面应力问题; 因为板面上(z=± t/2 )不受力，所以有应力分量;; 在船体结构中，很多问题可以简化为平面应力问题处理。例如甲板开口、舷侧、横梁开孔和肘板的强度问题等等。;; 首先以通过中心C，并平行于z轴的直线为矩轴，列力矩平衡方程?Mc=0。; 其次，以x轴为投影轴，列出力投影的平衡方程?Fx=0：; 其次，以y轴为投影轴，列出力投影的平衡方程?Fy=0：; 由此，得到平面应力问题的平衡方程：;; 线段PA的正应变：; 剪应变?xy：; 剪应变?xy：; 由式（8-11）可见，当板内各点的位移分量u、v为已知函数时，就可确定各点的应变分量?x、 ?y、 ?xy；反之，假定三个应变分量函数，那么按式（8-11）的前两个式子就可以求出位移函数u、v，若用此两个位移分量函数代入第三个方程式求?xy，就会与假定的?xy不同。这样就出现了矛盾。这一矛盾是因为板内任一点的应变分量之间有相互联系所造成的，如果在假定各应变分量函数时不反映出这种联系，那就将使变形不连续，即板变形将发生空隙或裂缝。从数学上讲，式（8-11）的应变分量是三个，而位移分量只有两个，因此三个应变分量不能相互独立，而必然有一定的联系，这个关系叫应变协调方程或应变相容条件。现推导如下：; 应变相容方程或叫应变协调方程：; （3）物理方程式; 或写成：; 用矩阵形式表示为：; （4）边界条件; （4）边界条件; （4）边界条件; 在物体边界上取直角三角形的微块PAB，它的斜面AB与物体边界面重合，如图6-4所示。用N表示边界面的外法线法向，其方向余弦为：; 设该边界面AB 的长度为ds，则截面PA和PB的长度分别为lds和mds，另设微块的厚度为1。由微块的平衡条件?Fx=0，得;; 1.解题方法简介; 和：;;;;;2012-12-29;2012-12-29;2012-12-29;2012-12-29;3、计算等效节点力;§ 8-4 三角形单元的位移函数与刚度矩阵; 2.位移函数;式中 a1, a2... a6, 是待定常数，它们可以由单元的边界条件，即节点的位移值来确定。为此，只要将节点的坐标值代入式（a)，就得到节点的位移值：;其中，节点的坐标值是已知的，令;;2012-12-29;;;4. 单元 刚度矩阵;;;;; 结构的平衡条件可用所有结点的平衡条件表示。假定i 结点为结构中的任一公共结点，则该结点平衡条件为：; 每个结点由两个平衡方程组成，若结构共有n个结点，则有2n个平衡方程。整个结构的平衡条件上面的公式求和得到，即：;进而可得到;作用在弹性体上的载荷不外乎是：* 体力（自重、惯性力）* 面力* 集中力三种。 用有限元法解题时，既然全部问题都归结到节点来处理，那么，当单元上作用有外载荷时，也应把它们移置到节点上来，成为节点载荷。 这种移置必须按照静力等效原则进行; 所谓“静力等效”，系指原载荷与移置后的节点载荷，在弹性体的任何虚位移过程中的虚功相等。 当插值函数已经确定时，这种移置的结果是唯一的。 在取线性位移插值时，符合刚体静力等效原则，即：载荷与节点载荷在任一轴上投影之和相等，对任一轴的力矩之和也相等，也就是说，原载荷与节点载荷将具有相同的主矢和主矩。;(1)体力的移置;首先，我们求移置到i节点上的垂直节点载荷 。为了便于计算虚功，假想该单元在节点i处沿y方向产生一个单位虚位移，而其它两点不动；这相当于在j点及m点安置了铰支座，在i点安置了水平连杆支座，如右图。 由于我们在三角形单元中采用的位移插值函数是线性的，所以任一条直线上各点位移都呈线性变化。; 现在m点及j点的位移都等于0，所以在 边上各点位移都等于0； 线上各点的垂直位移也按线性变化，在b点等于0，在i点为1。因此c点的垂直位移将为1/3。 按静力等效原则，体力载荷W的虚功应等于 的虚功，即有： 或 Yei=-W/3负号表示Yei的方向与图上所画的方向相反。 ;用同样的方法可以得到 Yej=-W/3 , Yem=-W/3;下面来求移置到节点i