材力上机报告.docVIP

材料力学上机报告 应力状态分析 院 系： 能源科学与工程学院 班 级： 设 计 者： 学 号： 指导教师： 刘一志 王春香 软件要求： Visual Basic 6.0 设计时间： 作业要求 一.实验问题提出： 本程序使用Microsoft Visual Basic编写，由输入平面应力状态，输入空间应力状态，实现了平面应力状态分析，空间应力状态分析，可输出主应力，最大剪应力，输出应力圆图形，如为平面应力状态则输出方位角斜截面上的应力及主方向角。能进行简单计算平面应力状态和空间应力状态各应力及角度的大小，并利用平面应力圆来求解主应力和切应力。 所用到的相关公式： 1，平面应力状态：σx σy τxy α σmax = σmin = τmax = τmin = （主平面的方位角） tan 2 = （主切平面的方位角） 2，空间应力状态：σx σy σz τxy τyz τzx I1 =σx + σy + σz I2 = + I3 = I1 + I2I3 = 0 (主应力的三次方程) tmax=1/2*(σ1-σ3) 二，VB程序源代码 Private Sub Comman() Dim X As Single, Y As Single, W As Single, Z As Single, A As Single, B As Single （定义变量） End Sub Private Sub Command1_Click() Dim X As Single, Y As Single, W As Single, Z As Single, A As Single, B As Single X = Val(Text1.Text) Y = Val(Text2.Text) W = Val(Text3.Text) （给变量赋值） Z = Val(Text4.Text) A = (X + Y) / 2 + Sqr(((X - Y) / 2) ^ 2 + W ^ 2) B = (X + Y) / 2 - Sqr(((X - Y) / 2) ^ 2 + W ^ 2) （计算主应力的大小） If A O Then Text11.Text = A Else Text12.Text = A Text11.Text = 0 End If If B 0 Then Text12.Text = B （进行比较大小，并输出答案） Text13.Text = 0 Else Text13.Text = B Text12.Text = 0 End If Text14.Text = Sqr(((X - Y) / 2) ^ 2 + W ^ 2) Text15.Text = (X + Y) / 2 + 0.5 * (X - Y) * Cos(2 * Z * 3.1415 / 180) + W * Sin(2 * Z * 3.1415 / 180) Text16.Text = -0.5 * (X - Y) * Sin(2 * Z * 3.14159 / 180) + W * Cos(2 * Z * 3.1415 / 180) If X = Y Then （计算主切应 Text17.Text = 45 力及角度） Else Text17.Text = 0.5 * Atn(2 * W / (X - Y)) * 180 / 3.1415 End If Text18.Text = Text17.Text + 90 End Sub Private Sub Command2_Click() Dim X(6) As Double Dim Y(3) As Double Dim J As Integer (定义变量) Dim A As Double, B As Double, C As Double, D As Double Dim Z(3) As Double, W As Double, T As Double X(0) = Val(Text5.Text) X(1) = Val(Text6.