应力状态分析-2016最新课.docVIP

应力状态分析 问题的提出 题目7: 应力状态分析 输入:1.平面应力状态输入：（σx、σy、τxy）；某截面方位角α 2.空间应力状态输入（σx、σy、σz、τxy、τyz、τzx） 输出： 5. 输出主应力（σ1、σ2、σ3） 6. 最大切应力[τmax=τ13=（σ1-σ3）/2] 7. 如为平面应力状态则需输入方位角α斜截面上的应力σα、τα及主方向角α0。 8. 画出应力圆示意图 该程序可以解决的问题： 平面应力状态，已知σx，σy，τxy和角度的情况下计算不同任意截面的应力分量，还可以计算平面应力状态主应力的大小和方向并能画出应力圆； 空间应力状态，已知σx、σy、σz、τxy、τyz、τzx的情况下计算主应力和最大切应力并能画出应力圆。 数学模型 1、平面应力状态任一斜截面上既有正应力又有切应力，公式为： 2、平面应力状态主应力大小及方向 ??? １、公式 ??? 主应力计算公式 ??? ??? 主平面方位计算公式 ??? 3、最大切应力：最大切应力发生在与主平面夹45度角的平面方位。 ??? 其公式 ??? 4、应力圆 应力圆原理： 圆心： 半径： 5、三向应力状态 三向应力计算公式： 特征方程： 三个不变量  程序概况 2菜单 选择是平面应力状态或者空间应力状态。 （2）如果选择平面应力状态， 输入σx、σy、τxy 、方位角α并单击“计算”“画应力图”，得到主应力及应力圆，输出主应力（σ1、σ2、σ3）最大切应力[τmax=τ13=（σ1-σ3）/2] （3）如果选择空间应力状态， 输入σx、σy、σz、τxy、τyz、τzx）并单击“计算”“画应力图”，得到主应力及应力圆，输出主应力（σ1、σ2、σ3）最大切应力τmax 以下为一个普通情况下的计算结果 以下为σ1、σ2、σ3相同情况下的计算结果，以及画出的点圆 以下为计算平面应力情况下的计算结果，以及画出的应力圆 以下为σ1、σ2、σ3全部为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 添加(σx, τxy) (σy, τyx) (σa, τa)到圆心的连线，并添加了标注，使绘制的应力圆更加易读 Picture2.Line (O, 0)-(X / k, -t1 / k), vbBlue Picture2.Print (σx, τxy) Picture2.Line (O, 0)-(Y / k, t1 / k), vbBlue Picture2.Print (σy, τyx) Picture2.Line (O, 0)-(out(4) / k, -out(5) / k), vbGreen Picture2.Print (σa, τa) ⑵空间： 针对盛金公式计算平面应力状态会出现问题的情况，把平面的情况单独提出来计算 If a1 ^ 2 + a4 ^ 2 + a5 ^ 2 = 0 Then If 1 / 2 * (a2 + a3) + Sqr(((a2 - a3) / 2) ^ 2 + a6 ^ 2) 0 Then X(2) = Format$(1 / 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(1) = Format$(0) X(0) = Format$(1 / 2 * (a2 + a3) - Sqr(((a2 - a3) / 2