第7章节 内力和应力.pptVIP

* * 第16章 内力和应力计算 调用第15章解方程的子程序，可得到总体坐标系下结构的节点位移?。根据这一位移向量，可计算出每个单元上的节点力或单元内的应力。 §16.1 计算公式 16.1.1 单元节点位移 解平衡方程得到的节点位移，是总体坐标系下定义的整体位移。对每个单元而言，要从中取出单元位移，才能计算节点力。对于杆单元和梁单元还必须转换到局部坐标系下，这样得到的节点力才具有明确的物理意义。 1．总体坐标系下的单元节点位移 总体坐标系下单元节点位移向量?e与整个结构的节点位移向量? 的对应关系，可通过IS数组联系。根据第14章IS数组的定义，单元内部编号为I的节点位移，对应的结构整体位移编号为IS(I)，所以 I = 1 , 2 , ? , NFD (16.1) 2．局部坐标下的单元节点位移 已知单元总体坐标系下的节点位移 以后，可利用坐标变换求出局部坐标系下单元的节点位移 (16.2) 式中R为坐标变换矩阵。 16.1.2 单元节点力 根据结构设计的需要，对于杆单元和梁单元一般是以内力作为承载能力的依据，所以需要计算和输出单元节点力。节点力是在单元局部坐标系下定义的，它们分别对应于轴力、剪力和弯矩等内力。 1．杆单元 杆单元节点力与节点位移之间的关系，由局部坐标系下的单元刚度矩阵联系 (16.3) 2．梁单元 梁单元上若有非节点荷载作用，节点力还需在式(16.3)的基础上叠加各非节点荷载引起的固端反力 ，即 (16.4) 16.1.3 单元应力 1．杆单元 对于桁架杆单元，通常计算输出单元（杆段）轴力N和正应力? 的数值。由图2.1关于杆单元局部坐标系下节点力的定义，可知轴力N等于j节点的节点力，即 (16.5) 所以，轴向正应力为 (16.6) 2．平面应力单元 平面应力单元的应力向量?通过弹性矩阵D和应变向量相联系，而应变向量又和应变矩阵B、节点位移向量?e联系在一起，即 (16.7) 式中S=DB称为应力矩阵。 对于三角形单元，应变和应力在单元内为常量，一般认为它是三角形形心处的应力代表值。对于三角形某一顶点上的应力可由相关单元应力的算术平均值来表示。对于八节点等参单元，应力在单元内是变化的，所以需要计算指定点的应力。通常计算单元内高斯积分点的应力，也可以计算八个节点的应力，随需要而定。 强度计算的依据不是式(16.7)所示的正应力和剪应力，而是由此产生的主应力或等效应力。这是二向应力状态，由材料力学理论可得主应力和主方向分别为 (16.8) 和 (16.9) 是较大主应力 与x轴的夹角，由莫尔应力圆分析可知：当 ? 时， ? 90?；当 ? 时， ?90?； 当 = 时， ?0、2 = - 90?， ? 0、2 = 90? 。 将应力 、 和0按代数值的大小排序，得到第一、第二和第三主应力 、 和 。第一主应力为最大主应力、第三主应力为最小主应力，第二主应力可由应力的第一不变量计算出来，即： (16.10) 由三个主应力可以计算相应的Tresca等效应力（第三强度理论）和Mises等效应力（第四强度理论）： (16.11) §16.2 程序设计 16.2.1 变量说明 NL2 ： 局部坐标系下单元自由度（桁架杆单元=2 ；梁单元=6） PPE(NL2) ：局部坐标系下单元节点力 UPE(NL2)：局部坐标系下单元节点位移 UE(NFD) ：总体坐标系下单元节点位移 S(3 , 6) ：平面应力三角形单元的应力矩阵S = DB S(3 , 16) ：平面八节点等参单元应力矩阵S = DB SS(8)、TT(8)： 平面八节点等参单元标准单元的八个节点上的 值。 ST(3) ：平面应力单元的应力向量 16.2.2 程序段 1．梁杆单元 DO 800 IE = 1 , NE 对单元循环 CALL KEP(NN , NA , IE , MEA , AEU , AAI , X , Y , AKEP) 计算局部单元刚度 CALL CR(NN , IE , MEA , X , Y , R) 计算坐标变换矩阵R（空间桁架有z坐标） CALL FIS( IE , MEA , NF , ND , NFD ,