Harbin Institute of Technology材料力学电算大作业课程名称：材料力学设计题目：组合截面几何性质计算作者院系：作者班级：作者姓名：作者学号：指导教师：完成时间：一、软件主要功能X4，X5，X6分别是n1个圆形截面，n2个圆环形截面，n3个矩形截面的形心位置X与面积的乘积Y4，Y5，Y6分别是n1个圆形截面，n2个圆环形截面，n3个矩形截面的形心位置Y与面积的乘积Xc,Yc是总截面的形心坐标Ix1,Ix2,Ix3分别是n1个圆形截面，n2个圆环形截面，n3个矩形截面对通过形心且与x轴平行的轴的惯性矩Iy1,Iy2,Iy3分别是n1个圆形截面，n2个圆环形截面，n3个矩形截面对通过形心且与y轴平行的轴的惯性矩Ixy1,Ixy2,Ixy3分别是n1个圆形截面，n2个圆环形截面，n3个矩形截面对通过形心且与x,y轴平行的两轴的惯性积a是通过形心的主轴与x轴的夹角Imax,Imin分别是截面对形心主轴的主惯性矩软件截图：二、程序源代码Dim n1 As DoubleDim d1(10) As DoubleDim X1(10) As DoubleDim Y1(10) As DoubleDim n2 As DoubleDim d2(10) As DoubleDim d3(10) As DoubleDim X2(10) As DoubleDim Y2(10) As DoubleDim n3 As DoubleDim h(10) As DoubleDim d(10) As DoubleDim X3(10) As DoubleDim Y3(10) As DoubleDim S1 As Double, S2 As Double, S3 As DoubleDim X4 As Double, Y4 As Double, X5 As Double, Y5 As Double, X6 As Double, Y6 As DoubleDim Xc As Double, Yc As DoubleDim Ix1 As Double, Iy1 As Double, Ix2 As Double, Iy2 As Double, Ix3 As Double, Iy3 As Double, Imax As Double, Imin As DoubleDim Ixy1 As Double, Ixy2 As Double, Ixy3 As DoubleDim a As DoublePrivate Sub Text1_Change()n1 = Val(Text1.Text)For i = 1 To n1d1(i) = Val(InputBox("输入第" & (i) & "个圆的直径"))X1(i) = Val(InputBox("输入第" & (i) & "个圆的圆心的x坐标值")) Y1(i) = Val(InputBox("输入第" & (i) & "个圆的圆心的y坐标值")) Next iFor i = 1 To n1S1 = S1 + 3.14159 * d1(i) * d1(i) / 4X4 = X4 + X1(i) * 3.14159 * d1(i) * d1(i) / 4Y4 = Y4 + Y1(i) * 3.14159 * d1(i) * d1(i) / 4Next iEnd SubPrivate Sub Text2_Change()n2 = Val(Text2.Text)For i = 1 To n2d2(i) = Val(InputBox("输入第" & (i) & "个圆环的外径"))d3(i) = Val(InputBox("输入第" & (i) & "个圆环的内径"))X2(i) = Val(InputBox("输入第" & (i) & "个圆的圆心的x坐标值")) Y2(i) = Val(InputBox("输入第" & (i) & "个圆的圆心的y坐标值")) Next iFor i = 1 To n2S2 = S2 + 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4X5 = X5 + X2(i) * 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4Y5 = Y5 + Y2(i) * 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4Next iEnd SubPrivate Sub Text3_Change()n3 = Val(Text3.Text)For i = 1 To n3h(i) = Val(InputBox("输入第" & (i) & "个矩形的高"))d(i) = Val(InputBox("输入第" & (i) & "个矩形的宽"))X3(i) = Val(InputBox("输入第" & (i) & "个矩形形心的x坐标值")) Y3(i) = Val(InputBox("输入第" & (i) & "个矩形形心的y坐标值")) Next iFor i = 1 To n3S3 = S3 + h(i) * d(i)X6 = X6 + X3(i) * h(i) * d(i)Y6 = Y6 + Y3(i) * h(i) * d(i)Next iEnd SubPrivate Sub Command1_Click()Xc = (X4 + X5 + X6) / (S1 + S2 + S3)Yc = (Y4 + Y5 + Y6) / (S1 + S2 + S3)Text4 = XcText5 = YcEnd SubPrivate Sub Command2_Click()Picture1.BackColor = vbWhitePicture1.Scale (-Max - 10, Max + 10)-(Max + 10, -Max - 10)Picture1.Line (0, Max + 10)-(0, -Max - 10)Picture1.Line (Max + 10, 0)-(-Max - 10, 0)Picture1.FillStyle = 0For i = 1 To n1 Step 1Picture1.FillColor = vbBluePicture1.Circle (X1(i), Y1(i)), d1(i) / 2, vbRedNext iFor i = 1 To n2 Step 1Picture1.FillColor = vbGreenPicture1.Circle (X2(i), Y2(i)), d2(i) / 2, vbRedPicture1.FillColor = vbWhitePicture1.Circle (X2(i), Y2(i)), d3(i) / 2, vbRedNext iFor i = 1 To n3 Step 1Picture1.FillColor = vbRedPicture1.Line (X3(i) - d(i) / 2, Y3(i) - h(i) / 2)-(X3(i) + d(i) / 2, Y3(i) + h(i) / 2), vbRed, BFNext iEnd SubPrivate Sub Command3_Click()For i = 1 To n1Ix1 = Ix1 + 3.14159 * d1(i) * d1(i) * d1(i) * d1(i) / 64 + (Y1(i) -Yc) * (Y1(i) - Yc) * 3.14159 * d1(i) * d1(i) / 4Iy1 = Iy1 + 3.14159 * d1(i) * d1(i) * d1(i) * d1(i) / 64 + (X1(i) - Xc) * (X1(i) - Xc) * 3.14159 * d1(i) * d1(i) / 4Ixy1 = Ixy1 + (X1(i) - Xc) * (Y1(i) - Yc) * 3.14159 * d1(i) * d1(i) / 4Next iFor i = 1 To n2Ix2 = Ix2 + 3.14159 * (d2(i) * d2(i) * d2(i) * d2(i) - d3(i) * d3(i) * d3(i) * d3(i)) / 64 + (Y2(i) - Yc) * (Y2(i) - Yc) * 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4Iy2 = Iy2 + 3.14159 * (d2(i) * d2(i) * d2(i) * d2(i) - d3(i) * d3(i) * d3(i) * d3(i)) / 64 + (X2(i) - Xc) * (X2(i) - Xc) * 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4Ixy2 = Ixy2 + (X2(i) - Xc) * (Y2(i) - Yc) * 3.14159 * (d2(i) * d2(i) - d3(i) * d3(i)) / 4Next iFor i = 1 To n3Ix3 = Ix3 + h(i) * h(i) * h(i) * d(i) / 12 + (Y3(i) - Yc) * (Y3(i) - Yc) * h(i) * d(i)Iy3 = Iy3 + d(i) * d(i) * d(i) * h(i) / 12 + (X3(i) - Xc) * (X3(i) - Xc) * h(i) * d(i)Ixy3 = Ixy3 + (X3(i) - Xc) * (Y3(i) - Yc) * h(i) * d(i)Next iImax = (Ix1 + Ix2 + Ix3 + Iy1 + Iy2 + Iy3) / 2 + Sqr((Ix1 + Ix2 + Ix3 - Iy1 - Iy2 - Iy3) * (Ix1 + Ix2 + Ix3 - Iy1 - Iy2 - Iy3) / 4 + (Ixy1 + Ixy2 + Ixy3) * (Ixy1 + Ixy2 + Ixy3))Imin = (Ix1 + Ix2 + Ix3 + Iy1 + Iy2 + Iy3) / 2 - Sqr((Ix1 + Ix2 + Ix3 - Iy1 - Iy2 - Iy3) * (Ix1 + Ix2 + Ix3 - Iy1 - Iy2 - Iy3) / 4 + (Ixy1 + Ixy2 + Ixy3) * (Ixy1 + Ixy2 + Ixy3))Text6 = ImaxText7 = IminIf Ix <> Iy Thena = Math.Atn(-2 * (Ixy1 + Ixy2 + Ixy3) / (Ix1 + Ix2 + Ix3 - Iy1 - Iy2 - Iy3)) / 2Picture1.Line (Xc, Yc)-(Xc + 20, Yc + 20 * Tan(a)), vbGreenPicture1.Line (Xc, Yc)-(Xc - 20 * Tan(a), Yc + 20), vbGreenElsePicture1.Line (Xc, Yc)-(Xc, Yc + 20), vbGreenPicture1.Line (Xc, Yc)-(Xc + 20, Yc), vbGreenEnd IfEnd Sub程序见/share/link?shareid=734076929&uk=100682977 三、例题的计算结果数据已知一工字型钢如图所示上下两矩形一样，宽是1，高是0.4，中间的矩形宽是0.4，高是2 求形心及主惯性矩解：Xc=0.5 Yc=1.4Ix1= 0.58 Iy1=0.0333Ix2= 0.267 Iy2=0.0107Ix3= 0.58 Iy3=0.0333Ixy1= Ixy2= Ixy3=0夹角a=0°Ix总= Ix1+ Ix2+ Ix3=1.427Iy总= Iy1+ Iy2+ Iy3=0.0773Imax=1.427Imin=0.0773计算结果如图：。