5土压力程序设计
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
22
土木学院报告
１.５画土压力分布图 Picture1.Cls A = A + Val(Text7(i)) Picture1.Scale (0, 0)-(Val(Text10(Val(Text11) - 1)) + 2, Picture1.Line -(Val(Text10(i)), A) A + 1) Picture1.Line (0, A)-(Val(Text10(i)), A) A=0 For i = 0 To Val(Text11) - 1 Picture1.Print Val(Text10(i)) If i = 0 Then Picture1.CurrentX = Val(Text10(i)) A = A + Val(Text7(i)) Picture1.CurrentY = A Picture1.Line (0, 0)-(Val(Text10(i)), A)‟斜线 End If‟梯形情况 Picture1.Line (0, A)-(Val(Text10(i)), A)„横线 Next Picture1.Print Val(Text10(i))
创建地层
15
土木学院报告
１.２
创建地层
Private Sub Text11_Change() If Val(Text11) > 0 Then '加载土层数 ElseIf Val(Text11) < Text7.Count Then If Val(Text11) > Text7.Count Then For i = Text7.UBound To Val(Text11) Step 1 For i = Text7.UBound + 1 To Val(Text11) - 1 Load > 0 Then If i Text7(i): Load Text8(i): Load Text10(i) Text7(i).Top =Text7(i): 1).Top + 300 Unload Text7(i - Unload Text8(i): Text8(i).Top = Text8(i - 1).Top + 300 Unload Text10(i) Text10(i).Top = Text10(i - 1).Top + 300 End If Text7(i).Visible = True Next Text8(i).Visible = True End If End IfText10(i).Visible = True Next
18
土木学院报告
１.４计算合力与作用点 '计算合力大小 E0 = 0 For i = 0 To Val(Text11) - 1 If i = 0 Then E0 = E0 + 0.5 * Val(Text7(i)) * Val(Text10(i)) Else E0 = E0 + 0.5 * Val(Text7(i)) * (Val(Text10(i)) + Val(Text10(i - 1))) End If Next Text9(1) = Format(E0, "0.000kN/m")
5
土木学院报告
1.4 演示数据 Private Sub Command3_Click() Text2 = 4 Text2_Change Text1(0) = "0,100" Text1(1) = "50,90" Text1(2) = "50,50" Text1(3) = "0,0" End Sub
21
土木学院报告
１.４.２求合力作用点 x = x / Val(Text9(1)) A=0 For i = 0 To Val(Text11) - 1 A = A + Val(Text7(i)) Next x=A-x Text9(2) = Format(x, “0.000m”) 求解出合力作用点 距底端的距离
End Sub
16
土木学院报告
１.３计算每层的静止土压力
Private Sub Command7_Click() Dim E0 As Double, P0() As Double，H() As Double, y() As Double ReDim H(Val(Text11) - 1), y(Val(Text11) - 1), P0(Val(Text11) - 1) Dim i As Integer For i = 0 To Val(Text11) - 1 H(i) = Val(Text7(i))： y(i) = Val(Text8(i)) If i > 0 Then P0(i) = P0(i - 1) + K0 * H(i) * y(i) Else P0(i) = K0 * H(i) * y(i) End If Text10(i) = Format(P0(i), "0.000kPa") '输出每层的土压力 Next Text9(0) = P0(i - 1)
17
土木学院报告
１.４计算合力与作用点 (1)合力大小的求解方法 合力大小就是左图中的面 积。 A.i=0时，第一层是一个 三角形， B.i>0时，是梯形，采用 梯形的面积计算公式 E0 = E0 + 0.5 * Val(Text7(i)) * (Val(Text10(i)) + Val(Text10(i - 1)))
4
土木学院报告
1.3 加载多边形的角点数
Dim i% ElseIf Text1.Count > Val(Text2) Then If Val(Text2.Text) > Text1.Count And Val(Text2) <= 20 For i = Text1.Count - 1 To Val(Text2) Step -1 Then Unload Label1(i) For i = Text1.Count To Val(Text2) - 1 Unload Text1(i) Load Label1(i) ：Load Text1(i)'加载多边形 Next Label1(i).Caption = i + 1 ElseIf Val(Text2) < 3 Then Label1(i).Top = Label1(i - 1).Top + 450 MsgBox "边数小于3，构不成多边形", vbCritical, "多 边形输入出错" Text1(i).Top = Text1(i - 1).Top + 450 EndText1(i) = "" If Label1(i).Visible = True Text1(i).Visible = True Next i
土木学院报告
1
土木学院报告
2
土木学院报告
1.1设计界面 有时路基是多边形的，要 用平均断面法或平均距离 法计算土石方量，就应用 到多边形面积求解。
3
土木学院报告
1.2 计算原理 Ｓ＝０.５*∑(x(i) + x(i - 1)) * (y(i) - y(i - 1)) 最后一点与最初一点相连。 按逆时针或顺时针输入，注意要能构成多边形。
超固结粘性土： k OCR (1 sin f ' ) 0
11
土木学院报告
静止土压力系数取值 《公路桥涵设计通用规范》给出了静止土压力系 数的参考值为： 砾石、卵石：0.20 砂土：0.25 亚砂土：0.35 亚粘性土：0.45 粘土：0.55
12
土木学院报告
１.１
静止土压力系数求解
13
土木学院报告
１.１
静止土压力系数求解
基本方法： Private Sub Option2_Click(Index As Integer) 。。。。。。。 当某个单选按扭被选中时，就用该方 Case 4： K0 = 0.2 Case 5： K0 = 0.25 法来计 Private1 To Text1.ubound For i = Sub Command1_Click() Dim AA() Sum + (x(i)x()x(i - Double, y()y(i - 1)) Sum = As String, + As 1)) * (y(i) - As Double，Sum#,i%'把输入数据转存到数组之中 Next For ij= 0 To Text1.ubound If Text1(i).Text Like Dim As Double "*,*" Then j = Text1.ubound Sum = Sum + (x(0) + x(j)) * ",") - y(j)) AA = Split(Text1(i).Text, (y(0) Sum = Sum / 2 ReDim Preserve x(i):ReDim Preserve y(i) Text3 = Format(Abs(Sum), "0.###") x(i) = AA(0): y(i) = AA(1) End If End Sub Next
20
土木学院报告
１.４.２求合力作用点
Dim x As Single, A As Single x = 0 '力矩 A=0 For i = 0 To Val(Text11) - 1 If i = 0 Then x = x + 0.5 * Val(Text10(i)) * Val(Text7(i)) * 2 / 3 * Val(Text7(i)) A = A + Val(Text7(i)) Else x = x + Val(Text10(i - 1)) * Val(Text7(i)) * (0.5 * Val(Text7(i)) + A) '矩形产生的力矩 x = x + 0.5 * ((Val(Text10(i)) - Val(Text10(i - 1))) * Val(Text7(i)) * (2 / 3 * Val(Text7(i)) + A)) '三角形 产生的力矩 A = A + Val(Text7(i)) End If Next i
7
土木学院报告
运行界面
8
土木学院报告
9
土木学院报告
静止土压力理论 静止土压力设计步骤： １.静止土压力系数确定求解 ２.创建土层数，输入相关参数 ３.计算每层的静止土压力，设Ｋ０不变 ４.计算合力，计算合力作用点。 ５.绘制静止土压力分布图。
10
土木学院报告
3.1 静止土压力系数取值 对于侧限应力状态： sv z p0=sh sh =K0sv sh =u/(1-u)sv sv = u/(1-u)gz K0=n/(1-n) 由于土的 n 很难确定，K0常用经验公式计算 对于砂土、正常固结粘土：K0≈1-sinf‟
Case 6： K0 = 0.35 Case 7： K0 = 0.45 Case 8： K0 = 0.55 End Select If K0 = 0 Then Text12 = "出错" Else Text12 = K0 End If End Sub
14
土木学院报告
１.２
19
土木学院报告
１.４.２求合力作用点
方法： ①可对顶点取力矩或对底端点求力 矩。 本程序对顶点取矩。 ②梯形分解成矩形和三角形分别求 力矩。 x = x + Val(Text10(i - 1)) * Val(Text7(i)) * (0.5 * Val(Text7(i)) + A) '矩形产生的 力矩 x = x + 0.5 * ((Val(Text10(i)) - Val(Text10(i - 1))) * Val(Text7(i)) * (2 / 3 * Val(Text7(i)) + A)) '三角形 产生的力矩
土木学院报告
１.５画土压力分布图 Picture1.Cls A = A + Val(Text7(i)) Picture1.Scale (0, 0)-(Val(Text10(Val(Text11) - 1)) + 2, Picture1.Line -(Val(Text10(i)), A) A + 1) Picture1.Line (0, A)-(Val(Text10(i)), A) A=0 For i = 0 To Val(Text11) - 1 Picture1.Print Val(Text10(i)) If i = 0 Then Picture1.CurrentX = Val(Text10(i)) A = A + Val(Text7(i)) Picture1.CurrentY = A Picture1.Line (0, 0)-(Val(Text10(i)), A)‟斜线 End If‟梯形情况 Picture1.Line (0, A)-(Val(Text10(i)), A)„横线 Next Picture1.Print Val(Text10(i))
创建地层
15
土木学院报告
１.２
创建地层
Private Sub Text11_Change() If Val(Text11) > 0 Then '加载土层数 ElseIf Val(Text11) < Text7.Count Then If Val(Text11) > Text7.Count Then For i = Text7.UBound To Val(Text11) Step 1 For i = Text7.UBound + 1 To Val(Text11) - 1 Load > 0 Then If i Text7(i): Load Text8(i): Load Text10(i) Text7(i).Top =Text7(i): 1).Top + 300 Unload Text7(i - Unload Text8(i): Text8(i).Top = Text8(i - 1).Top + 300 Unload Text10(i) Text10(i).Top = Text10(i - 1).Top + 300 End If Text7(i).Visible = True Next Text8(i).Visible = True End If End IfText10(i).Visible = True Next
18
土木学院报告
１.４计算合力与作用点 '计算合力大小 E0 = 0 For i = 0 To Val(Text11) - 1 If i = 0 Then E0 = E0 + 0.5 * Val(Text7(i)) * Val(Text10(i)) Else E0 = E0 + 0.5 * Val(Text7(i)) * (Val(Text10(i)) + Val(Text10(i - 1))) End If Next Text9(1) = Format(E0, "0.000kN/m")
5
土木学院报告
1.4 演示数据 Private Sub Command3_Click() Text2 = 4 Text2_Change Text1(0) = "0,100" Text1(1) = "50,90" Text1(2) = "50,50" Text1(3) = "0,0" End Sub
21
土木学院报告
１.４.２求合力作用点 x = x / Val(Text9(1)) A=0 For i = 0 To Val(Text11) - 1 A = A + Val(Text7(i)) Next x=A-x Text9(2) = Format(x, “0.000m”) 求解出合力作用点 距底端的距离
End Sub
16
土木学院报告
１.３计算每层的静止土压力
Private Sub Command7_Click() Dim E0 As Double, P0() As Double，H() As Double, y() As Double ReDim H(Val(Text11) - 1), y(Val(Text11) - 1), P0(Val(Text11) - 1) Dim i As Integer For i = 0 To Val(Text11) - 1 H(i) = Val(Text7(i))： y(i) = Val(Text8(i)) If i > 0 Then P0(i) = P0(i - 1) + K0 * H(i) * y(i) Else P0(i) = K0 * H(i) * y(i) End If Text10(i) = Format(P0(i), "0.000kPa") '输出每层的土压力 Next Text9(0) = P0(i - 1)
17
土木学院报告
１.４计算合力与作用点 (1)合力大小的求解方法 合力大小就是左图中的面 积。 A.i=0时，第一层是一个 三角形， B.i>0时，是梯形，采用 梯形的面积计算公式 E0 = E0 + 0.5 * Val(Text7(i)) * (Val(Text10(i)) + Val(Text10(i - 1)))
4
土木学院报告
1.3 加载多边形的角点数
Dim i% ElseIf Text1.Count > Val(Text2) Then If Val(Text2.Text) > Text1.Count And Val(Text2) <= 20 For i = Text1.Count - 1 To Val(Text2) Step -1 Then Unload Label1(i) For i = Text1.Count To Val(Text2) - 1 Unload Text1(i) Load Label1(i) ：Load Text1(i)'加载多边形 Next Label1(i).Caption = i + 1 ElseIf Val(Text2) < 3 Then Label1(i).Top = Label1(i - 1).Top + 450 MsgBox "边数小于3，构不成多边形", vbCritical, "多 边形输入出错" Text1(i).Top = Text1(i - 1).Top + 450 EndText1(i) = "" If Label1(i).Visible = True Text1(i).Visible = True Next i
土木学院报告
1
土木学院报告
2
土木学院报告
1.1设计界面 有时路基是多边形的，要 用平均断面法或平均距离 法计算土石方量，就应用 到多边形面积求解。
3
土木学院报告
1.2 计算原理 Ｓ＝０.５*∑(x(i) + x(i - 1)) * (y(i) - y(i - 1)) 最后一点与最初一点相连。 按逆时针或顺时针输入，注意要能构成多边形。
超固结粘性土： k OCR (1 sin f ' ) 0
11
土木学院报告
静止土压力系数取值 《公路桥涵设计通用规范》给出了静止土压力系 数的参考值为： 砾石、卵石：0.20 砂土：0.25 亚砂土：0.35 亚粘性土：0.45 粘土：0.55
12
土木学院报告
１.１
静止土压力系数求解
13
土木学院报告
１.１
静止土压力系数求解
基本方法： Private Sub Option2_Click(Index As Integer) 。。。。。。。 当某个单选按扭被选中时，就用该方 Case 4： K0 = 0.2 Case 5： K0 = 0.25 法来计 Private1 To Text1.ubound For i = Sub Command1_Click() Dim AA() Sum + (x(i)x()x(i - Double, y()y(i - 1)) Sum = As String, + As 1)) * (y(i) - As Double，Sum#,i%'把输入数据转存到数组之中 Next For ij= 0 To Text1.ubound If Text1(i).Text Like Dim As Double "*,*" Then j = Text1.ubound Sum = Sum + (x(0) + x(j)) * ",") - y(j)) AA = Split(Text1(i).Text, (y(0) Sum = Sum / 2 ReDim Preserve x(i):ReDim Preserve y(i) Text3 = Format(Abs(Sum), "0.###") x(i) = AA(0): y(i) = AA(1) End If End Sub Next
20
土木学院报告
１.４.２求合力作用点
Dim x As Single, A As Single x = 0 '力矩 A=0 For i = 0 To Val(Text11) - 1 If i = 0 Then x = x + 0.5 * Val(Text10(i)) * Val(Text7(i)) * 2 / 3 * Val(Text7(i)) A = A + Val(Text7(i)) Else x = x + Val(Text10(i - 1)) * Val(Text7(i)) * (0.5 * Val(Text7(i)) + A) '矩形产生的力矩 x = x + 0.5 * ((Val(Text10(i)) - Val(Text10(i - 1))) * Val(Text7(i)) * (2 / 3 * Val(Text7(i)) + A)) '三角形 产生的力矩 A = A + Val(Text7(i)) End If Next i
7
土木学院报告
运行界面
8
土木学院报告
9
土木学院报告
静止土压力理论 静止土压力设计步骤： １.静止土压力系数确定求解 ２.创建土层数，输入相关参数 ３.计算每层的静止土压力，设Ｋ０不变 ４.计算合力，计算合力作用点。 ５.绘制静止土压力分布图。
10
土木学院报告
3.1 静止土压力系数取值 对于侧限应力状态： sv z p0=sh sh =K0sv sh =u/(1-u)sv sv = u/(1-u)gz K0=n/(1-n) 由于土的 n 很难确定，K0常用经验公式计算 对于砂土、正常固结粘土：K0≈1-sinf‟
Case 6： K0 = 0.35 Case 7： K0 = 0.45 Case 8： K0 = 0.55 End Select If K0 = 0 Then Text12 = "出错" Else Text12 = K0 End If End Sub
14
土木学院报告
１.２
19
土木学院报告
１.４.２求合力作用点
方法： ①可对顶点取力矩或对底端点求力 矩。 本程序对顶点取矩。 ②梯形分解成矩形和三角形分别求 力矩。 x = x + Val(Text10(i - 1)) * Val(Text7(i)) * (0.5 * Val(Text7(i)) + A) '矩形产生的 力矩 x = x + 0.5 * ((Val(Text10(i)) - Val(Text10(i - 1))) * Val(Text7(i)) * (2 / 3 * Val(Text7(i)) + A)) '三角形 产生的力矩