平面单元-等效结点荷载计算.ppt

北京大学力学与工程科学系 陈斌 1. 位移函数 1. 位移函数 1. 位移函数 1. 位移函数 1. 位移函数－例题 1. 位移函数－例题 2. 关于节点等效力的一点说明 3. 总体刚度矩阵的形成与特点—整体分析的一般步骤 3. 总体刚度矩阵的形成与特点—整体分析的一般步骤 3. 总体刚度矩阵的形成与特点—整体分析的一般步骤 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的形成 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的形成 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的特点 1. 总体刚度矩阵的形成与特点—整体刚度矩阵的特点 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的特点 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的特点 3. 总体刚度矩阵的形成与特点—整体刚度矩阵的特点 4. 边界条件的处理 无约束结构的整体刚度矩阵是奇异的，即整体平衡方程的解不唯一。位移约束常分为：节点固定和给定节点位移两种约束。 由于引入位移约束条件通常在整体刚阵及节点载荷形成后进行（也有在此之前进行的，如直接删除法），即此时[K]、{R}中的元素均已按一定顺序分别储存于相应的数组，故引入位移约束时，要求尽量不要打乱[K]、{R}的储存顺序。 引入约束的方法常有： 1）直接带入法 2）对角元素置1法 3）大数法 4）直接删除法（降阶法） 5）罚单元法 4. 边界条件的处理 1）直接带入法（降阶法） 4. 边界条件的处理 2）对角元素置1法 处理ai=b形式的边条件。 4. 边界条件的处理 4. 边界条件的处理 3）大数法（适合于计算机处理） α1 4. 边界条件的处理 4）直接删除法 只能处理ai=0的情况。在单刚集成总刚时，对应与自由度ai的元素不进入 总刚，即在建立节点自由度与方程号之间的对照表时，把对应自由度的ID设成0。 这种方法降低了总刚的阶数。 4. 边界条件的处理 5）罚单元法：适于处理自由度耦合的约束 5. 面积坐标 记 L1=A1/A L2=A2/A L3=A3/A P点对应的面积坐标(L1, L2, L3) 面积坐标相互不完全独立： L1+L2+L3=1 显然 Li(xi,yi)=δij (i,j=1,2,3) 5. 面积坐标 5. 面积坐标 5. 面积坐标 5. 面积坐标 5. 面积坐标 6. 六节点三角形单元 6. 六节点三角形单元 6. 六节点三角形单元—例题 6. 六节点三角形单元—例题 4. 其它单元 完全二次多项式。 一次项保证了完备性。单元边界为二次变化，完全￡边界节点决定，保证了协调性 1(1,0,0) 2(0,1,0) 3(0,0,1) 4(1/2,1/2,0) 5(0,1/2,1/2) 6(1/2,0,1/2) x y 3 (x3,y3) 1 (x1,y1) 2 (x2,y2) 边长L q 4 5 6 例题. 求图所示二次三角形平面单元在142边作用有均布侧压q(线荷载密度)时的节点等效力，假设节点坐标已知。 * x y 1 2 3 u1 v1 u2 v2 u3 v3 (x1, y1) (x2, y2) (x3, y3) 六个节点位移只能确定六个多项式的系数，所以取这样的位移函数。该位移函数，将单元内部任一点的位移设定为坐标的线性函数，该位移模式很简单。其中α1-6为广义坐标或待定系数，可据节点1、2、3的位移值和坐标值求出。 为2A第1行各个元素的代数余子式 插值函数矩阵或形函数矩阵 x y 1 (a,0) 2 (0,a) 3 (0,0) 例题：图示等腰三角形单元，求其插值函数矩阵[N]。 连续弹性体离散为单元组合体时，需把弹性体承受的任意分布的载荷都向节点转移，而成为节点等效载荷（或节点等效力）。如果弹性体承受的载荷全都是集中力，则将所有集中力的作用点取为节点，就不存在转移的问题，集中力就是节点等效载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在节点上。因此，必须进行载荷转移。如果集中力的作用点未被取为节点，该集中力也要向节点转移。 将载荷转移到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。前面推导时使用的能量泛函пp对{ae}进行变分之后产生的{δae}实际上就是虚位移，以上公式可以适用于任意复杂的荷载情况。 如果单元为线性单元（如，本章的三节点三角形单元），则可以采用直接的静力等效法和虚功等效法。 图示结构的网格共有四个单元和六个节点。在节点1、4、6共有四个支杆支承。结构的载荷已经转换为节点载荷。 整体分析的四个步骤： 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度