击实配水程序设计

《击实配水程序设计》由会员分享，可在线阅读，更多相关《击实配水程序设计（48页珍藏版）》请在金锄头文库上搜索。

1、击实试验原理v以塑限为中心，相差2%含水量制备个试 样； v v现在某土样的天然含水率为15%v已知试样的塑限为16%v击实采用重型击实试验，每个筒需要试样约 6000g。v试求配水方案。v解：根据塑限为16%，则配水方案应该为：v12%，14%，16%，18%，20%v则加水时应加多少水？v当目标含水率小于天然含水率时，怎么办？第一种情况：当目标含水率高于天然含水率情况v公式推导：v则需要添加的水量为：v夏天时要考虑蒸发，需多加23%的水量v界面实现，为了防止出现太多的控件，可以 有表格实现，只要记住表格相应位置的数据 ，则直接引用即可，没有什么难度。v以本题为例：v为了方便计算，一般提供一

2、组演示数据：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.Cell

3、s(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变量参与运算，根据公式写

4、出程序则行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

5、“根据你提供的含水率，无法得出目标 含水率“, vbCriticalvEnd If土木学院报告土木学院报告 .朗肯土压力计算考虑朗肯土压力计算公式适用条件 分析输入的数据 朗肯主动土压力每层计算难点，当是粘性土时，存在拉裂现象 分析计算合力与合力作用点难点：粘性土，负压力情况下怎样计算 朗肯被动土压力计算，并计算合力与作用点。 绘制主动土压力与被动土压力分布图。土木学院报告 .朗肯土压力适用条件适用于墙背光滑，墙背垂直，填土面光滑。 输入参数的确定： 主要有墙后填土荷载，每层土厚度，容重，c，内 摩擦角及墙后填土分层数。 输出参数： 每层上土压力，每层下土压力 合力与合力作用点。土木学院报告

6、.土层数的确定此处代码略，已讲次了。只介绍步骤： 当目前的土层数小于目标土层数时： 加载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容重

7、 Text25(0) = 0：Text25(1) = 20：Text25(2) = 10粘聚力 Text26(0) = 30：Text26(1) = 0：Text26(2) = 20fai 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(Tex

8、t22) 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 Intege

9、r 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 *

10、 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

11、 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) - V

12、al(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(

13、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.Cur

14、rentX = 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土木学院报告 运行界面土木学院报告土木学院报告 .难点分析： 其基本代码与求解主动土压力相差不大，大部分 可以直接使用。 其比主动土压力要简单，因为被动土压力不可能 出现负值情况，都按梯形求解。

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