哈工大材料力学上机大作业.docx

1、哈工大材料力学上机大作业材料力学上机实验作业 组合截面几何性质计算 指导教师：王秋生 班级：学号：姓名：题目3：组合截面几何性质计算输入:1.各几何图形形状（圆形、矩形）及个数2.截面参数输入（形状、尺寸及位置）即圆截面-d，圆心位置;圆环截面-D,d,圆心位置;矩形截面h,b,形心位置输出：1. 截面构形图2. 形心位置坐标（图示）3. 在形心坐标系下画出主轴，计算主惯性矩程序说明：本程序用于计算材料截面几何量，是使用VB编程语言汇编而成。进入程序后，可选择三种模式： 1.常用模式(此模式用于计算常用正多边形的几何量) 2.直线模式，此模式下可利用线段绘制多边形(包括镂空多边形及多个不相交的

2、多边形) 3.曲线模式，此模式下可用鼠标任意绘制平面图形(同样包括镂空图形及多个不相交图形）一.绘图模式 1.常用模式 可鼠标绘制正多边形，或利用键盘精确输入正多边形。 2.直线模式 可用鼠标绘制任意多边形，或利用键盘精确输入点坐标。输入完成倒数第二个点后，单击淡绿色圆圈(起点)，在弹出的对话框选择是，完成第一个图形。可继续绘制下一图形或开始计算。 3.曲线模式 可用鼠标绘制任意曲线以构成图形。最后一点，单击淡绿色圆圈(起点)，在弹出的对话框选择是，完成第一个图形。可继续绘制下一图形或开始计算。二.计算 当绘制完某一图形后，点击“开始计算”按钮，即可计算出面积、静矩、形心坐标、惯性矩、惯性半径

3、、惯性积、主惯性轴角度、主惯性矩、形心主惯性轴角度、形心主惯性矩等16个几何量。欲重新开始，可点击“重置”按钮。点击任意按钮后屏幕下方都会有相应的提示。三.外部图形 首先将所需图形大小修改为500470，格式为jpg，命名为“截面”，并将其保存在程序根目录下，可载入外部图形以精确计算其各个量。源程序：注：以下为form1的代码Private X As Single, Y As Single, Xp As Single, Yp As SinglePrivate N As Integer, Num As Integer, M As Integer, missN(10000) As Integer,

4、 j As Integer, k As IntegerPrivate stepX As Single, stepY As Single, Step As SinglePrivate P(1 To 30000, 1 To 2) As Single, Q1 As Single, Q2 As SinglePrivate minX As Single, maxX As Single, minY As Single, maxY As SinglePrivate runX As Single, runY As SinglePrivate X1 As Single, Y1 As Single, X2 As

5、Single, Y2 As SinglePrivate Lx As Single, Ly As SinglePrivate Lpoint(200) As Single, Lpointmemory As SinglePrivate Nn As Integer, Rr As SinglePrivate Xx As Single, Yy As Single, T As Integer, Ang As SinglePrivate Sx As Single, Sy As Single, Ax As Single, Ay As SinglePrivate ShapeCenterX As Single, S

6、hapeCenterY As SinglePrivate Ix As Single, Iy As Single, IxC As Single, IyC As Single, IxyC As SinglePrivate Rx As Single, Ry As SinglePrivate Ixy As Single, Iyx As Single, LIxy As Single, LIyx As SinglePrivate Angle As SinglePrivate Imin As Single, Imax As SinglePrivate AngleC As SinglePrivate Icmi

7、n As Single, Icmax As SinglePrivate Can As IntegerPrivate Msg As IntegerPrivate MoveN As Long, MoveM As Long, Decider As Integer, Shape As IntegerPrivate A As Single, B As Single, C As Single, E As Single, F As SinglePrivate Usual As IntegerPrivate Sub PaintNnRr()M = 0: N = 1For i = 0 To Nn - 1 Text

8、1.Text = Xx + Rr * Cos(Ang * 3.1415926535898 / 180 + 2 * 3.1415926535898 / Nn * i) Text2.Text = Yy + Rr * Sin(Ang * 3.1415926535898 / 180 + 2 * 3.1415926535898 / Nn * i) Command1_ClickNextEnd SubPrivate Sub Shorten(A As Single, B As Single, C As Single)i = 0Do While A = 10 A = A / 10 i = i + 1LoopB

9、= Format(A, 0.00)C = iEnd SubPrivate Sub Calculate()ShapeCenterY = CLng(Sx / Ax)ShapeCenterX = CLng(Sy / Ay)If Ix = Iy And Ixy 0 Then Angle = -3.1415926535898 / 4ElseIf Ix = Iy And Ixy 0 Then AngleC = -3.1415926535898 / 4ElseIf IxC = IyC And IxyC 0 Then AngleC = 3.1415926535898 / 4ElseIf IxC = IyC A

10、nd IxyC = 0 Then anflec = 0Else AngleC = 0.5 * Atn(-2 * IxyC / (IxC - IyC)End IfEnd SubPrivate Sub Output()Picture1.Line (Sy / Ay + 1000, 6000 - Sx / Ax)-(Sy / Ay + 1000, 6000), RGB(255, 121, 121)Picture1.Line (Sy / Ay + 1000, 6000 - Sx / Ax)-(1000, 6000 - Sx / Ax), RGB(255, 121, 121)Label6.Left = S

11、y / Ay + 100 + 1000: Label6.Top = 5775 + 20Label6.Caption = X = & ShapeCenterXLabel6.Visible = TrueLabel7.Left = 1080: Label7.Top = 6000 - (Sx / Ax - 20)Label7.Caption = Y = & ShapeCenterYLabel7.Visible = TrueText3.Text = ShapeCenterX: Text4.Text = ShapeCenterYCall Shorten(Ax, E, F)Text5.Text = E &

12、& 10 & & FCall Shorten(Sx, E, F)Text6.Text = E & & 10 & & FCall Shorten(Sy, E, F)Text7.Text = E & & 10 & & FCall Shorten(Ix, E, F)Text8.Text = E & & 10 & & FCall Shorten(Iy, E, F)Text9.Text = E & & 10 & & FText10.Text = CLng(Rx * 10) / 10: Text11.Text = CLng(Ry * 10) / 10Call Shorten(Ixy, E, F)Text1

13、2.Text = E & & 10 & & FText13.Text = Format(Angle * 180 / 3.1415926535898, 0.0)Call Shorten(Imin, E, F)Text14.Text = E & & 10 & & FCall Shorten(Imax, E, F)Text15.Text = E & & 10 & & FText16.Text = Format(AngleC * 180 / 3.1415926535898, 0.0)Call Shorten(Icmin, E, F)Text17.Text = E & & 10 & & FText18.

14、Text = IcmaxCall Shorten(Icmax, E, F)Text18.Text = E & & 10 & & FEnd SubPrivate Sub Sort()For i = 1 To j - 2 For k = i + 1 To j - 1 If Lpoint(k) Lpoint(i) Then Lpointmemory = Lpoint(i) Lpoint(i) = Lpoint(k) Lpoint(k) = Lpointmemory End If NextNextEnd SubPrivate Sub CloseLine()Picture1.Circle (Q1 + 1

15、000, 6000 - Q2), 90, RGB(255, 255, 255)Label1.Caption = 输入第 & N + 1 & 点坐标:P(N, 1) = Q1: P(N, 2) = Q2Picture1.Line (P(N - 1, 1) + 1000, 6000 - P(N - 1, 2)-(Q1 + 1000, 6000 - Q2), RGB(0, 0, 255)missN(Num) = NM = 0N = N + 1Num = Num + 1Decider = 0Label21.Caption = 点击计算“截面的几何量”按钮开始计算,或单击下一点以绘制下一图形的起点End

16、 SubPrivate Sub GetstepXstepY()minY = P(1, 2)For i = 2 To N - 1 If P(i, 2) maxY Then maxY = P(i, 2) Else End IfNextstepY = (maxY - minY) / StepminX = P(1, 1)For i = 2 To N - 1 If P(i, 1) maxX Then maxX = P(i, 1) Else End IfNextstepX = (maxX - minX) / StepEnd SubPrivate Sub GetLxLIxy()j = 1: Lx = 0:

17、LIxy = 0For i = 1 To N - 1 Can = 0 For k = 1 To Num - 1 If i = missN(k) Then Can = 1 Else End If Next If Can = 0 Then If (P(i, 2) - runY) * (P(i + 1, 2) - runY) = 0 Then runY = runY + stepY / 1000 Else End If If (P(i, 2) - runY) * (P(i + 1, 2) - runY) 0 Then X1 = P(i, 1): Y1 = P(i, 2): X2 = P(i + 1,

18、 1): Y2 = P(i + 1, 2) Lpoint(j) = X1 + (X2 - X1) * (runY - Y1) / (Y2 - Y1) j = j + 1 Else End If Else End IfNextCan = 0Call SortFor i = 1 To j - 1 Step 2 Lx = Lx + Lpoint(i + 1) - Lpoint(i) LIxy = LIxy + 0.5 * (Lpoint(i + 1) 2 - Lpoint(i) 2)NextEnd SubPrivate Sub PaintX()For i = 1 To j - 2 Step 2 Pi

19、cture1.Line (Lpoint(i) + 1000, 6000 - runY)-(Lpoint(i + 1) + 1000, 6000 - (runY + stepY), RGB(255, 113, 116), BFNextEnd SubPrivate Sub GetAxSxIx()Ax = Ax + Lx * stepYSx = Sx + (runY + 0.5 * stepY) * (Lx * stepY)Ix = Ix + (runY + 0.5 * stepY) 2 * (Lx * stepY)Ixy = Ixy + (runY + 0.5 * stepY) * (LIxy *

20、 stepY)End SubPrivate Sub GetLyLIyx()j = 1: Ly = 0: LIyx = 0For i = 1 To N - 1 Can = 0 For k = 1 To Num - 1 If i = missN(k) Then Can = 1 Else End If Next If Can = 0 Then If (P(i, 1) - runX) * (P(i + 1, 1) - runX) = 0 Then runX = runX + stepX / 1000 Else End If If (P(i, 1) - runX) * (P(i + 1, 1) - ru

21、nX) 0 Then X1 = P(i, 1): Y1 = P(i, 2): X2 = P(i + 1, 1): Y2 = P(i + 1, 2) Lpoint(j) = Y1 + (Y2 - Y1) * (runX - X1) / (X2 - X1) j = j + 1 Else End If Else End IfNextCall SortFor i = 1 To j - 1 Step 2 Ly = Ly + Lpoint(i + 1) - Lpoint(i) LIyx = LIyx + 0.5 * (Lpoint(i + 1) 2 - Lpoint(i) 2)NextEnd SubPri

22、vate Sub PaintY()For i = 1 To j - 2 Step 2 Picture1.Line (runX + 1000, 6000 - Lpoint(i)-(runX + stepX + 1000, 6000 - Lpoint(i + 1), RGB(118, 89, 255), BFNextEnd SubPrivate Sub GetAySyIy()Ay = Ay + Ly * stepXSy = Sy + (runX + 0.5 * stepX) * (Ly * stepX)Iy = Iy + (runX + 0.5 * stepX) 2 * (Ly * stepX)I

23、yx = Iyx + (runX + 0.5 * stepX) * (LIyx * stepX)End SubPrivate Sub Command1_Click() 画线Label1.Caption = 输入第 & N + 1 & 点坐标:Text1.Text = CLng(Text1.Text * 10) / 10: Text2.Text = CLng(Text2.Text * 10) / 10Xp = Text1.Text: Yp = Text2.TextX = Xp: Y = YpP(N, 1) = X: P(N, 2) = YIf M = 0 Then Q1 = X: Q2 = Y

24、Picture1.PSet (Q1 + 1000, 6000 - Q2), RGB(0, 0, 255) Picture1.Circle (Q1 + 1000, 6000 - Q2), 90, RGB(180, 255, 180) If Shape = 0 Then Decider = 0 Else Decider = 1 End IfElseIf M 0 Then Picture1.Line (P(N - 1, 1) + 1000, 6000 - P(N - 1, 2)-(P(N, 1) + 1000, 6000 - P(N, 2), RGB(0, 0, 255)ElseEnd IfM =

25、1N = N + 1Command2.Enabled = FalseCommand4.Enabled = FalseCommand5.Enabled = FalseCommand9.Enabled = FalseLabel21.Caption = 点击以输入下一个点的坐标 .或点击当前图形起点(淡绿色圆圈内)以封闭图形End SubPrivate Sub Command10_Click()Nn = Text19.TextRr = Text20.TextXx = Text22.TextYy = Text23.TextAng = Text21.TextCall PaintNnRrCall Clos

26、eLineCommand2_ClickCommand10.Enabled = FalseEnd SubPrivate Sub Command2_Click()N = N - 1Call GetstepXstepYrunX = minX: runY = minYAx = 0: Sx = 0: Ix = 0Do While runY maxY Call GetLxLIxy Lx为与x轴平行的L Call PaintX Call GetAxSxIx runY = runY + stepYLoopAy = 0: Sy = 0: Iy = 0Do While runX maxX Call GetLyLI

27、yx Call PaintY Call GetAySyIy runX = runX + stepXLoopCall CalculateCall OutputCommand1.Enabled = FalseCommand2.Enabled = FalseN = N + 1Label21.Caption = 单击“重置程序”按钮以初始化程序Form1.Width = 9750End SubPrivate Sub Command3_Click()Form_ActivateEnd SubPrivate Sub Command4_Click()Decider = 0Shape = 0Label21.Caption = 进入直线模式If T = 1 Then Frame6.Visible = False T = 0ElseEnd IfCommand1.Enabled = TrueCommand4.Enabled = FalseCommand5.Enabled = FalseCommand9.Enabled = FalseEnd SubPrivate Sub Command5_Click()Decider = 0Shape = 1Label21.Caption = 进入曲线模式If T = 1 Then Frame6.Visible = False T = 0ElseEnd IfCommand