第7章 内力和应力.pptVIP

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