击实配水程序设计.ppt

击实试验原理v①以塑限为中心，相差2%含水量制备５个试 样； v v现在某土样的天然含水率为15%v已知试样的塑限为16%v击实采用重型击实试验，每个筒需要试样约 6000gv试求配水方案v解：根据塑限为16%，则配水方案应该为：v12%，14%，16%，18%，20%v则加水时应加多少水？v当目标含水率小于天然含水率时，怎么办？第一种情况：当目标含水率高于天然含水率情况v公式推导：v则需要添加的水量为：v夏天时要考虑蒸发，需多加2～3%的水量v界面实现，为了防止出现太多的控件，可以 有表格实现，只要记住表格相应位置的数据 ，则直接引用即可，没有什么难度v以本题为例：v为了方便计算，一般提供一组演示数据：v'“当目标含水率高于天然含水率情况“vSpreadsheet1.Rows(1).Font.Color = vbBluevSpreadsheet1.Cells(1, 1) = “当目标含水率高于天然含水率情况“vSpreadsheet1.Cells(2, 1) = “试样的天然含水率w0(%)“vSpreadsheet1.Cells(3, 1) = “试样的目标含水率w1(%)“vSpreadsheet1.Cells(4, 1) = “天然试样的质量m0(g)“vSpreadsheet1.Cells(5, 1) = “考虑蒸发需增加水百分比(%)“vSpreadsheet1.Cells(2, 2) = 15vSpreadsheet1.Cells(3, 2) = 20vSpreadsheet1.Cells(4, 2) = 6000vSpreadsheet1.Cells(5, 2) = 2计算v‘第一种情况，把表格给变量vDim W0#, W1#, M0#, AddB#, Mw#vW0 = Spreadsheet1.Cells(2, 2) / 100vW1 = Spreadsheet1.Cells(3, 2) / 100vM0 = Spreadsheet1.Cells(4, 2)vAddB = Spreadsheet1.Cells(5, 2) / 100变量参与运算，根据公式写出程序则行vIf W1 = WW2 And WW2 >= WW1 ThenvMM0 = MM * (1 + WW0) * (WW2 - WW1) / (1 + WW2) / (WW0 - WW1)输出vSpreadsheet1.Cells(8, 4) = “需要天然含水率试样 (g):“vSpreadsheet1.Cells(8, 5) = Format(MM0, “####0.000“)vSpreadsheet1.Cells(9, 4) = “需要风干试样(g):“vSpreadsheet1.Cells(9, 5) = Format(MM - MM0, “####0.000“)vElsev MsgBox “根据你提供的含水率，无法得出目标 含水率“, vbCriticalvEnd If土木学院报告土木学院报告 ２.１朗肯土压力计算Ø①考虑朗肯土压力计算公式适用条件 Ø②分析输入的数据 Ø③朗肯主动土压力每层计算Ø难点，当是粘性土时，存在拉裂现象 Ø④分析计算合力与合力作用点Ø难点：粘性土，负压力情况下怎样计算 Ø⑤朗肯被动土压力计算，并计算合力与作用点。

Ø⑥绘制主动土压力与被动土压力分布图土木学院报告 ２.２朗肯土压力适用条件Ø适用于墙背光滑，墙背垂直，填土面光滑 Ø输入参数的确定： Ø主要有墙后填土荷载，每层土厚度，容重，c，内 摩擦角及墙后填土分层数 Ø输出参数： 每层上土压力，每层下土压力 Ø合力与合力作用点土木学院报告 ２.３土层数的确定Ø此处代码略，已讲Ｎ次了只介绍步骤： Ø①当目前的土层数小于目标土层数时： Ø加载load 对象(i) Ø定位 对象.top Ø显示控件 对象.visible=true Ø②当目前的土层数大于目标土层数时： Ø卸载相应控件 unload 对象(i)土木学院报告 2.4 演示数据 ØPrivate Sub Command2_Click() ØText22 = 3 ØText22_Change ØText29 = “0kPa“ ’ ØText23(0) = 1:Text23(1) = 2:Text23(2) = 3’高度 ØText24(0) = 16.5：Text24(1) = 18：Text24(2) = 19‘容重 ØText25(0) = 0：Text25(1) = 20：Text25(2) = 10‘粘聚力 ØText26(0) = 30：Text26(1) = 0：Text26(2) = 20‘fai ØEnd Sub土木学院报告 ２.５每层的主动土压力求解 Ø计算方法： Ø定义一个动态数组 ，接收相关输入的 数据，参与计算Ø Dim E0 As Double, Pa() As Double Ø Dim H() As Double, y() As Double, C() As Double, fai() As Double, Ka() As Double, H1() Ø Dim Z As Integer Ø If Val(Text22) > 1 Then‘多组数据的情况下 Ø ReDim H(Val(Text22) - 1), y(Val(Text22) - 1), Pa(Val(Text22) - 1), Ka(Val(Text22) - 1), C(Val(Text22) - 1), fai(Val(Text22) - 1), H1(Val(Text22) - 1) Ø Z = Val(Text22) Ø Else‘只有一层土 Ø ReDim H(0), y(0), Pa(0), Ka(0), C(0), fai(0), H1(0) Ø Z = 1 Ø End If ’定义相关变量土木学院报告 ２.６用变量接收相应数据 ØDim i As Integer ØFor i = 0 To Z - 1 Ø H(i) = Val(Text23(i)) Ø H1(i) = Val(Text23(i)) Ø y(i) = Val(Text24(i)) Ø C(i) = Val(Text25(i)) Ø fai(i) = Val(Text26(i)) Ø Ka(i) = Tan((45 - fai(i) / 2) * 3.14159265 / 180) * Tan((45 - fai(i) / 2) * 3.14159265 / 180) Ø ØNext土木学院报告 ２.７每层土压力的计算Ø难点１，当一层出现负值时，要计算开裂深度， 此处暂不考虑该层上下都是负数的情况。

ØH1(i) = H(i) - ((2 * C(i) * Sqr(Ka(i)) - Ka(i) * G)) / Ka(i) / y(i) ‘计算开裂剩余深度 Ø难点２：每一层土都要作为下一层土的荷载 ØG = G + y(0) * H(0) + Val(Text29) ØG = G + y(i) * H(i)土木学院报告 ２.７每层土压力的计算Ø Text27(0) = Ka(0) * Val(Text29) - 2 * C(0) * Sqr(Ka(0)) '深度为0 Ø Text30(0) = Ka(0) * Val(Text29) + Ka(0) * y(0) * H(0) - 2 * C(0) * Sqr(Ka(0)) '深度为H1，H1上 Ø Dim G As Double Ø G = G + y(0) * H(0) + Val(Text29) Ø For i = 1 To Z - 1 Ø Text27(i) = G * Ka(i) - 2 * C(i) * Sqr(Ka(i)) 'H1下 Ø Text30(i) = Ka(i) * G + Ka(i) * y(i) * H(i) - 2 * C(i) * Sqr(Ka(i)) 'H2上 Ø If Val(Text27(i)) 0 Then’梯形情况，分两部分 Ø Liju = Liju + Val(Text27(i)) * H(i) * (H2 + H(i) / 2) + 0.5 * H(i) * (Val(Text30(i)) - Val(Text27(i))) * (H2 + H(i) / 3) Ø Else‘三角形情况 Ø Liju = Liju + 0.5 * H1(i) * (Val(Text30(i))) * (H2 + 1 / 3 * H1(i)) Ø End If Ø Next Ø Text28(1).Text = Format(Liju / S, “距底部0.000m“)土木学院报告 ２.10画土压力分布图Ø 难点分解 Ø 确定最大值与最小值，不然画图时可能画 到图片框外面去。

Ø H2 = 0 '计算总高度 Ø Dim MinPa As Double, MaxPa As Double Ø MinPa = Val(Text27(0)): MaxPa = Val(Text30(0)) Ø For i = 0 To Z - 1 Ø H2 = H2 + H(i) Ø If Val(Text27(i)) > MaxPa Then MaxPa = Val(Text27(i)) Ø If Val(Text27(i)) MaxPa Then MaxPa = Val(Text30(i)) Ø If Val(Text30(i)) 0 Then Picture3.Line (0, H2)- (Val(Text30(i - 1)), H2) ØNext土木学院报告 2.10 ３把计算值标注到图上 Ø '打印计算值 Ø H2 = 0 Ø For i = 0 To Z - 1 Ø H2 = H2 + H(i) Ø Picture3.CurrentX = Val(Text27(i)): Picture3.CurrentY = H2 - H(i) Ø Picture3.Print Format(Val(Text27(i)), “0.0“) Ø Picture3.CurrentX = Val(Text30(i)): Picture3.CurrentY = H2 Ø Picture3.Print Format(Val(Text30(i)), “0.0“) Ø Next土木学院报告 运行界面土木学院报告土木学院报告 ２.１１Ø难点分析： Ø其基本代码与求解主动土压力相差不大，大部分 可以直接使用。

Ø其比主动土压力要简单，因为被动土压力不可能 出现负值情况，都按梯形求解Ø考试内容：土木学院报告 总体设计界面土木学院报告 运行界面土木学院报告土木学院报告 ３.１基本计算原理土木学院报告 ３.１基本计算原理土木学院报告 ３.１基本计算原理土木学院报告 ３.１基本计算原理土木学院报告 ３.２程序设计难点Ø考虑输入的参数，此方法是用于路堤计算，所以 要考虑上部荷载的情况：土木学院报告 ３.３ 设计界面土木学院报告 ３.４ 运行界面给定安全系数，求剩余下滑力土木学院报告Ø设计计算： Ø注意： Ø当不是最后一块时，剩余下滑力为负时，说明没 有剩余下滑力，该力不能传递给下一块可令其 为０. Ø最后一块剩余的下滑力为负时，说明安全，否则 是不安全的土木学院报告 ３.５ 运行界面求解安全系数土木学院报告 ３.６ 求解安全系数算法Ø设安全系数从０.１增加到２０，增幅０.００１ Ø当最后一块剩余下滑力正好为０时，或接近于０ 时，则此时假设的安全系数就是滑坡的安全系数 。