资源描述：
《材料力学上机大作.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、材料力学电算大作业題目名称:：应力状起分析班号：1108103学号：姓名:指导老师：张桂莲完庆时间：2013.06.07哈尔滨工业大学题耳7应力状态分析输入：1、平面应力状态输入：（5、勾、Txy）；某截面方位角（X2^空间应力状态输入：心、勺、6、Txy＞珈、Tzx）输出：1、输出主应力（巧、立、a3）；2、最大剪应力（Tmax=t13=^）；3、如为平面应力状态则需输出方位角cc斜截面上的应力％、s及主方位角；4、画出应力圆示意图。一、数学模型1、平面应力状态任一斜截面上既有正应力又有切应力，公式为:o严丑空+^空Zx-T、sin2x??t(i=662a+[cos2cz2、平
2、面应力状态主应力大小及方向主应力计算公式主平面方位计算公式3、最大切应力：最大切应力发生在与主平面夹45度角的平面方位。其公式6一64、应力圆应力圆原理:/2匕一如圆心:6+6,a=2半径：5、三向应力状态三向应力计算公式：特征方程：a-ha+La~L=°三个不变量心+中rr,yyx56二、程序流程图1、平面应力计算打开界面选择平面应力计算出现数据输入对话框►输入。x,0y,Txy和角度►点击计算可计算出任意截面的应力分量以及三个主应力和最大切应力记录数据►点击清除可进行重复计算。2、空间应力状态计算打开界面选择空间应力状态计算►出现数据输入对话框►输入ox,oy,oz,TXV
3、,Tyz,TTZ点击计算得i出结果可计算出三个主应力及最大切应力记录数据一►点击清除可以重复计算。三、程序说明本程序分为平面应力状态与空间应力状态，点击相应主菜单进行相应的计算。平面应力状态：本程序能进行多次重复输入与计算根据提示输入数据进行计算输入X,Y,XY方向上的应力值及X'转动的角度值，点击"计算〃按扭，进行数值计算。可画出与数据相对应的应力圆空间应力状态：本程序能进行多次重复输入与计算根据捉示输入数据进行计算输入X,Y,Z,XY,YZ,ZY方向上的应力值，点击〃计算〃按扭,可以计算出三个主应力及最大切应力。四、程序源代码运用VisualBasic语言编程：Private
4、Subclear3_Click(IndexAsInteger)DimjAsIntegerText1.TextText2.TextText3.Text=Clear=Clear=ClearText11.Text=ClearPiel.CisFori=0To6Step1Text10(i).Text=ClearNextiEndSubPrivateSubCombo1Click()Piel.CisTfCombo1.ListTndex=1ThenFrame2.Visible=FalseFramel.Visible=TrueFrame3.Visible=FalseElseFramel.Visibl
5、e=FalseFrame2.Visible=TrueFrame3.Visible二FalseEndIfEndSubPrivateSubCommandlClick0Textl0(4).Visible=FalseText10(5).Visible=FalseTextlO(6).Visible=FalseDima(6)AsDoubleDimi(3)AsDoubleDimjAsIntegerDimx(3)AsDouble,wAsDouble,tAsDoublea(0)=Vai(Text4.Text)a(1)=Vai(Text5.Text)a(2)=Vai(Text6.Text)a(3)=
6、Vai(Text7.Text)a(4)=Vai(Text8.Text)a(5)=Vai(Text9.Text)i(0)=-1*(a(0)+a(l)+a⑵)’lli(1)=a(0)*a(l)-a(3)*a(3)+a(1)*a(2)+a(0)*a(2)-a(4)*a(4)-a(5)*a(5)'12i(2)=-1*(a(0)*a(l)*a(2)+2*a(3)*a(4)*a(5)-a(l)*a(5)*a(5)-a(0)*a(4)*a(4)-a(2)*a(3)*a(3))'13'以下解一元三次方程。3+1(0)*。2+1(1)*。+1(2)二0'abcd'一元三次方程aX^3+bX^2+c
7、X+d=0,(a,b,c,dWR,HaHO)。'重根判别式:'A二b"2—3ac；'B=bc—9ad；'C=c"2—3bd,'总判别式：'A=B^2-4ACo'当A=B=O时，盛金公式①(WhenA二B二0,Shengjin'sFormulaQ:'X1二X2二X3二-b/(3a)=一c/b=一3d/c。八'当厶二B2—4AC二0吋，盛金公式③(When△二B2—4AC=0,Shcngjin'sFormula③):'Xl=-b/a+K；''X2二X3二一K/2,'其中K=B/A,(
