应力状态分析2016课.docVIP

应力状态分析 问题的提出 题目7: 应力状态分析 输入:1.平面应力状态输入：（ax、oy、Txy ）;某截面方位角a 2 ?空间应力状态输入（ax、oy、az、Txy、Tyz、tzx ） 输出： 5.输出主应力（al、a2、＜?3 ） 6.最大切应力[Ttnax=Tl3= （ al-a3 ） /2] 7. 如为平面应力状态则需输入方位角a斜截面上的应力aa、Ta及主方向角aO。 8.画出应力圆示意图 该程序可以解决的问题： 平面应力状态，已知ax , oy , Txy和角度的情况下计算不同任意截面的应力分量，还可以计算平面应力状态主应力的大小和方向 并能画出应力圆； 空间应力状态，已知ax、ay、az、Txy、Tyz、tzx的情况下计算主应力和最大切应力并能画出应力圆。 数学模型 1、平面应力状态任一斜截面上既有正应力又有切应力，公式为： 5Mi2a+x cosSz 2、平面应力状态主应力大小及方向 1、公式 圆心： 主应力计算公式 主平面方位计算公式 三向应力计算公式: 3、最大切应力：最大切应力发生在与主平面夹45度角的平 面方位。 其公式 特征方程： a-/ia+M-Zi=0 三个不变量yxzxXXI(J2yz 三个不变量 yx zx X XI (J 2 yz 程序概况 2菜单 选择是平面应力状态或者空间应力状态。 如果选择平面应力状态，输入ax、ay、Txy、方位角a并单击计算画应力图，得到主应力及应力圆，输出主应力(a2、a3) 最大切应力[1^=下13= ( ara3) /2] (3 )如果选择空间应力状态，输入ax、ay、az、Tyz、tzx )并单击计算画应力图，得到主症力及应力圆，输出主应力(a2、 a3)最大切应力 Trnox 以下为一^普通情况下的计算结果 以下为A、a2、＜j3相同情况下的计算结果，以及画出的点圆 以下为计算平面应力情况下的计算结果，以及画出的应力圆 以下为q、a2、(j3全部为0情况下的计算结果，以及画出的圆 3程序介绍 程序设计难点 (1)平面: 在计算过程中对三个主皮力进行排序(事后反思其实可以先算出三个主皮力再排序); 当画点圆的时候发现绘图程序自动画出的点圆看不见，于是在程序中对点圆进行了判断， 在需要画点圆的地方画了一个半径为0.5的圆，使其看起来像是一个点圆。 If r = 0 Then Picture2.Circle (O, 0), 0.06, vbRed Else Picture2.Circle (O, 0), r, vbRed End If 添加(ax, Txy) (oy，yyx) (aa, Ta)到圆心的连线，并添加了标注，使绘制的应力圆更加易读 Picture2.Line (0, 0)-(X / k，-tl / k), vbBlue Picture2.Print (ax, Txy) -APicture2.Line (O, 0)-(Y / k，tl / k), vbBlue Picture2.Print (oy，Tyx) -A Picture2.Line (O, 0)-(out(4) / k, -out(5) / k)，vbGreen Picture2.Print (aa，to) ⑵空间： 针对盛金公式计算平面应力状态会出现问题的情况，把平面的情况单独提出来计算 If al 2 + a4 ~ 2 + a5 2 = 0 Then If 1 / 2 * (a2 + a3) + Sqr(((a2 - a3) / 2) A 2 + a6 △ 2) 0 Then X(2)= Format$(l / 2 * (a2 + a3) + Sqr(((a2 - a3) / 2) 2 + a6 2), 0.0) If 1 / 2 * (a2 + a3) - Sqr(((a2 - a3) / 2) 2 + a6 2) 0 Then X(l) = Format$(0) X(0)= Format$(l / 2 * (a2 + a3) - Sqr(((a2 - a3) / 2) 一 2 W 2)，0.0) Else X(0) = Format$(0) X(l)= Format$(l / 2 * (a2 + a3) - Sqr(((a2 ? a3) / 2广 2 + a6 2)，0.0) End If Else X(2) = Format$(0) X(l)= Format$(l / 2 * (a2 + a3) + Sqr(((a2 - a3) / 2) A 2 + a6 A 2), 0.0) X(0)= Format$(l / 2 * (a2 + a3) - Sqr(((a2 - a3) / 2) A 2 ? a6 A 2)，0.0) End If Elself a2/2 + a