《平面单元等效结点荷载计算.ppt》由会员分享，可在线阅读，更多相关《平面单元等效结点荷载计算.ppt（35页珍藏版）》请在三一文库上搜索。
1、2019年4月,1,1. 位移函数,六个节点位移只能确定六个多项式的系数，所以取这样的位移函数。该位移函数，将单元内部任一点的位移设定为坐标的线性函数，该位移模式很简单。其中1-6为广义坐标或待定系数，可据节点1、2、3的位移值和坐标值求出。,2019年4月,2,1. 位移函数,2019年4月,3,1. 位移函数,为2A第1行各个元素的代数余子式,2019年4月,4,1. 位移函数,插值函数矩阵或形函数矩阵,2019年4月,5,1. 位移函数例题,例题：图示等腰三角形单元，求其插值函数矩阵N。,2019年4月,6,1. 位移函数例题,2019年4月,7,2. 关于节点等效力的一点说明,连续弹性
2、体离散为单元组合体时，需把弹性体承受的任意分布的载荷都向节点转移，而成为节点等效载荷（或节点等效力）。如果弹性体承受的载荷全都是集中力，则将所有集中力的作用点取为节点，就不存在转移的问题，集中力就是节点等效载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在节点上。因此，必须进行载荷转移。如果集中力的作用点未被取为节点，该集中力也要向节点转移。 将载荷转移到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。前面推导时使用的能量泛函p对ae进行变分之后产生的ae实际上就是虚位移，以上公式可以适用于任意复杂的荷载情况。 如果单元为线性单元（如，本章的三节
3、点三角形单元），则可以采用直接的静力等效法和虚功等效法。,2019年4月,8,3. 总体刚度矩阵的形成与特点 整体分析的一般步骤,图示结构的网格共有四个单元和六个节点。在节点1、4、6共有四个支杆支承。结构的载荷已经转换为节点载荷。 整体分析的四个步骤： 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度矩阵； 3、解方程组，求节点位移； 4、根据节点位移求出应力。,2019年4月,9,3. 总体刚度矩阵的形成与特点 整体分析的一般步骤,1、建立整体刚度矩阵 上图中的结构有六个节点，共有12个节点位移分量(自由度)和12个节点力分量，它们之间的关系为：,总体刚度方程中的自由度与节点位移之间的
4、对应关系,2019年4月,10,3. 总体刚度矩阵的形成与特点 整体分析的一般步骤,2、根据支承条件修改整体刚度矩阵。 建立整体刚度矩阵时，每个节点的位移当作未知量看待，没有考虑具体的支承情况，因此进行整体分析时还要针对支承条件加以处理。 在上图的结构中，支承条件共有四个，即在节点1、4、6的四个支杆处相应位移已知为零：u1=u4=v4=v6=0 建立节点平衡方程时，应根据上述边界条件进行处理。 3、解方程组，求出节点位移。 通常采用消元法和迭代法两种方法。 4、根据节点位移求出应力。,2019年4月,11,3. 总体刚度矩阵的形成与特点 整体刚度矩阵的形成,1、总刚形成的物理背景： 刚度矩阵
5、中的元素，即由节点作单位位移时引起的节点力。在单刚Ke中，Kije表示第j个位移(自由度)给一单位位移，其它位移为零时，单元在第i位移方向上引起的节点力；类似，在整体刚阵中，Ki,j表示第j个自由度给一单位位移，其它自由度为零时，整体结构在第i个自由度上引起的节点力（即所有与第i、j个自由度相关的单元在第i个自由度上引起的节点力之和）。 如上图结构，计算K3,5时（第3和5个自由度分别对应第2和3号节点的u,即x向位移），与节点2和3相关的单元有单元和，当节点3发生x向单位位移时，相关单元和同时在节点2的x向引起节点力，将这两个力相加，就得出 K3,5 = K511 + K153 。由此看出，
6、总刚的刚度系数是相关单刚的刚度系数的集成。,2019年4月,12,2、刚度矩阵的集成规则： 1）在整体离散结构变形后，应保证各单元在节点处仍然协调地相互连接，即在该节点处所有单元在该节点上有相同位移。 2）整体离散结构各节点应满足平衡条件。即环绕每个节点的所有单元作用其上的节点力之和应等于作用于该节点上的节点载荷Ri。,3. 总体刚度矩阵的形成与特点 整体刚度矩阵的形成,2019年4月,13,1、对称性。由Kij的物理意义和互易定理可以很容易得到此结论。 利用对称性可以只存贮矩阵的上三角部分，节省近一半的存贮容量。 2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。,3. 总体刚度
7、矩阵的形成与特点 整体刚度矩阵的特点,节点1只与周围的两个节点(2、3)用三角形单元相连，它们是1的相关节点。在矩阵K中，第1行的非零元素只有6个(对应于相关节点的x,y向自由度)。,问题的规模越大，矩阵中的非零元素所占的比例就越小,2019年4月,14,3、带形分布规律。 右图中，矩阵K的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中(包括对角线元素在内)，每行具有的元素个数叫做半带宽，用d表示。半带宽的一般计算公式是： 半带宽 d = (相邻结点码的最大差值 + 1) * 2 左图中相邻节点码的最大差值为4，故d=(4+1)*2=10 利用带形矩阵的特点并利用对称
8、性，可只存贮上半带的元素，叫半带存贮。若每行都取不同的半带宽则称作轮廓线存储。,1. 总体刚度矩阵的形成与特点 整体刚度矩阵的特点,2019年4月,15,3. 总体刚度矩阵的形成与特点 整体刚度矩阵的特点,图(a)中的矩阵K为n行n列矩阵，半带宽为d。半带存贮时从K中取出上半带元素，按图(b)中的矩阵K的排列方式进行存贮，即将上半部斜带换成竖带。存贮量n*d，存贮量与K中元素总数之比为d/n，d值越小，则存贮量约省。,矩阵K 矩阵K* 对角线 第1列 r行 r行 r列 45度斜线 r行s列 r行s-r+1列元素 元素,2019年4月,16,3. 总体刚度矩阵的形成与特点 整体刚度矩阵的特点,同
9、一网格中，如果采用不同的节点编码，则相应的半带宽d也可能不同。如图，是同一网格的三种节点编码，相邻节点码的最大差值分别为4、6、8，半带宽分别为10、14、18。因此，应当采用合理的节点编码方式，以便得到最小的半带宽，从而节省存贮容量。,2019年4月,17,3. 总体刚度矩阵的形成与特点 整体刚度矩阵的特点,4、Kii0 5、带入边条件之前，总刚各行（列）元素之和等于0 6、带入边条件之前，总刚奇异,2019年4月,18,4. 边界条件的处理,无约束结构的整体刚度矩阵是奇异的，即整体平衡方程的解不唯一。位移约束常分为：节点固定和给定节点位移两种约束。 由于引入位移约束条件通常在整体刚阵及节点
10、载荷形成后进行（也有在此之前进行的，如直接删除法），即此时K、R中的元素均已按一定顺序分别储存于相应的数组，故引入位移约束时，要求尽量不要打乱K、R的储存顺序。 引入约束的方法常有： 1）直接带入法 2）对角元素置1法 3）大数法 4）直接删除法（降阶法） 5）罚单元法,2019年4月,19,4. 边界条件的处理,1）直接带入法（降阶法）,改变了原方程的顺序，只适用于一些简单的问题,2019年4月,20,4. 边界条件的处理,2）对角元素置1法 处理ai=b形式的边条件。,2019年4月,21,4. 边界条件的处理,图示结构，对边界支承条件处理后，整体刚度矩阵修改为：,2019年4月,22,4
11、. 边界条件的处理,3）大数法（适合于计算机处理） 1,2019年4月,23,4. 边界条件的处理,4）直接删除法 只能处理ai=0的情况。在单刚集成总刚时，对应与自由度ai的元素不进入 总刚，即在建立节点自由度与方程号之间的对照表时，把对应自由度的ID设成0。 这种方法降低了总刚的阶数。,2019年4月,24,4. 边界条件的处理,5）罚单元法：适于处理自由度耦合的约束,2019年4月,25,5. 面积坐标,记 L1=A1/A L2=A2/A L3=A3/A P点对应的面积坐标(L1, L2, L3) 面积坐标相互不完全独立： L1+L2+L3=1 显然 Li(xi,yi)=ij (i,j=
12、1,2,3),三角形的高次单元如果仍然直角坐标系来定义插值函数Ni，其公式将变得很复杂，若采用面积坐标则很简单,1,2,3下标轮转,2019年4月,26,5. 面积坐标,L1、L2、L3实际上就是推导三节点三角形平面单元时的N1、N2、N3 面积坐标与直角坐标之间的变换关系：,2019年4月,27,5. 面积坐标,2019年4月,28,5. 面积坐标,利用此式可以重新计算线性三角形单元侧边均布压力的节点等效力,2019年4月,29,5. 面积坐标,例：q为线荷载密度，利用面积坐标计算节点1、2的等效节点力。,2019年4月,30,5. 面积坐标,所以三角形线荷载在1,2号节点的节点等效力分别为：,2019年4月,31,完全二次多项式。 一次项保证了完备性。单元边界为二次变化，完全边界节点决定，保证了协调性,6. 六节点三角形单元,2019年4月,32,6. 六节点三角形单元,2019年4月,33,6. 六节点三角形单元 例题,例题. 求图所示二次三角形平面单元在142边作用有均布侧压q(线荷载密度)时的节点等效力，假设节点坐标已知。,2019年4月,34,6. 六节点三角形单元 例题,P1:P4:P2=1:4:1,2019年4月,35,4. 其它单元,(1). 10节点三次三角形单元,完全三次多项式。 插值函数,