空间轴对称问题.pptxVIP

会计学; 由于各分量都对称于z轴，与 无关，因此弹性体内各点只可能存在着径向位移u和轴向位移w（此时，u, w只是r, z的函数），而环向位移v = 0。则必有剪应变 ，否则这种轴对称问题的弹性体将不能保持轴对称状态而发生歪扭，这在实际中是不可能出现的。;几何方程为; 根据轴对称问题的特点，我们只须考查轴对称结构的任意子午面上各单元的各应力分量和各结点的位 移分量。在轴对称结构的任意子午面上 （即rz平面）任取一个三角形单元i,j,m， 基本未知量仍然取结点位移，单元的结 点位移可用列阵表示为 ; 仿照平面问题，取线性位移模式;而;(b)式也可写成矩阵形式;其中; 由此可见，单元中的应变分量 都是常量，但是环向正应变 不是常量，它与 中的r有关。;其中; 显然，只有应力分量 在单元中为常量外，其余三个正应力在单元中都不是常量。在实用上，为了简化计算和消除对称轴上由于r = 0所引起的麻烦，常把各个单元中的r及z近似地当作常量，并且分别等于各单元形心的坐标，即; 现在，再运用虚功原理求导轴对称结构上任意单元的刚度矩阵;将上式代入(a)式，并注意到;上式右边与单元结点位移列阵 相乘的矩阵便是单元刚度矩阵; 由于在轴对称问题的矩阵 中出现坐标r、z，所以(4-11)式的积分运算比平面问题要复杂得多。现在仍取单元形心的坐标; 对于整体刚度矩阵，如果弹性体被划分为 个单元和 n 个结点，于是就可得到 个型如(d)式的方程组。与平面问题的情况完全相类似的处理，把各单元的 、 、 等都加以扩大到整个结构的自由度的维数，然后叠加得到; 整体刚度矩阵;;其中子矩阵为; (g)式右边的载荷列阵展开的形式为;式中 为集中载荷 作用点的径向坐标。将上节(b)式代入上式并考虑到 ，上式可以化成; 采用平面问题中相同的符号： 集中力的等效结点力;于是(c)式可以改写成; 1. 体积力 (1) 自重。在此情况下 ；其中 为密度。于是单元的自重移置到结点i,j,m上的等效结点力???;这样就得到; (2) 离心力。在此情况下 ，其中 为角速度。于是单元的离心力移置到结点i,j,m上的等效结点力;;注意到在ij边上 ，于是积分; 经过类似推导可得移置到结点j和m上的等效结点力; 两种特殊情况 (1) ; (2) ; 3．考虑温度改变的影响，由于应力应变关系包含温度应变，于是公式(e)还应加上热载荷;上式 中的两个积分通常可用数值积分法求积，或简单地取近似表达式 ; 在§4-3中我们导得单元刚度矩阵 ，它可以划分为3个2×2的子矩阵，即;其中 是象§4-3中那样，用单元形心坐标 代入 后得到的；;而 是它的变化部分。由(4-7)式可知;把(a)式代入(4-11)式，有;故有;;上式使用了缩写符号;其中; 利用(4-33)式进行计算时，有两种特殊情况会出现奇异性。一种情况是当单元的某个结点i位于对称轴上（即ri=0），于是(4-33)式中包含对数的项; 结合以上两种情况，可以将(4-34)式改写成如下的形式以消除计算中的奇异性，即; 对于有两个结点在对称轴上的单元，例如结点i和j，则有ri=rj=0，此时除了子矩阵 外，其它子矩阵的积分均发散。但是，由于此时ui=uj=0，正好可将刚度矩阵中的第2i-1行和2i-1列以及2j-1行和2j-1列划去。因此，在计算程序中，可以统一应用公式(4-35)使(4-33)式的右端获得一个有限值。这些有限值出现在将要划去的行和列上，因此，并不影响计算的正确性。; 实际上，一些例子表明，精确刚度矩阵与近似刚度矩阵，对于主对角元，修正影响甚微，修正项在精确刚度矩阵的相应元素中所占的比例均在5%以下；对于非主对角元，一些元素的修正项所占的比例可达33%。尽管如此，对于位移与应力的计算，其影响并