受冲击荷载作用下土体的动态响应.docxVIP

受冲击荷载作用下土体的动态响应 1 钻头冲击碎岩动态动力学分析 影响方法是土壤样品的一定结构和输出，其中一些可切割土层深度，以及钻头对柱状土芯的扰动。根据钻头冲击碎岩特点，钻头高速冲击土，而后又离开土层，是一个高速动力过程，可以采用有限元中瞬时动态动力学分析来模拟分析。 有限元分析的关键是建立合适的计算模型。把土当作一种弹塑性材料来看，则钻头冲击土体的过程可以抽象为一个半无限体在圆环面积上受均布力作用的问题。 2 分析条件 2.1 钻体 钻头与土体接触部分的尺寸，外径D=129 mm，内径d=91mm。 2.2 单次激发能量和影响时间 根据室外传感器测冲击力实验可得，作用在钻头上的冲击力为4.92 kN，接触时间为0.029 s。据此分析土体内部应力应变状况。 2.3 boussi现行修正 把土当作一种弹塑性材料来看，则钻头冲击土体的过程可以抽象为一个半无限体在圆环面积上受均布力作用的问题，则土体变形分析体可以是一段圆柱体，只需确定其径向尺寸和轴向尺寸即可。为了确定土体的受力影响范围，笔者从受冲击土体的轴向位移着手，利用著名的Boussinesq解答, 推导出了土体表面上某点和土体样品轴线上某点的轴向位移即沉陷, 并由此作了解答。具体推导过程如下 设圆环形区域的内外半径分别为r z方向的位移分量ω可由Boussinesq位移解答的第二式确定： 积分可得 对于式 (3) ，设v=0.27，经换算最终可得如下结果： 当z=0时， 当z=12r 由此可得，在深为6倍的钻头外半径处的沉陷为地表沉陷的1/10，可忽略。 对于钻头外围土体表面上任意一点，由半无限体表面任意一点的法向位移积分并变换，最终得到如下公式： 对于式 (6) ，当r=r 当r=2r 当r=3r 由以上结果可知，水平向3倍钻头外径处土体的沉陷为钻头周围的1/5，因为水平方向不是我们关心的重点，故范围可定到此处。 因此，受冲击影响的土体范围可按如下式子确定： 水平方向：宽度=6r 垂直方向：高度=12 r 2.3 土体的泊松比和抗剪强度 土体的弹性变形模量：E=5.6 MPa 土体的泊松比：ν=0.27 土的重度：γ=19.6 kN/m 土体的抗剪强度：τ=0.23 MPa 土体的内摩擦角：φ=13o 2.4 土体分析体结构 利用2.3所述土的力学性能参数，经对比得出采用能考虑土体塑性的轴对称4节点平面单元对整个土体进行自动粗剖分，该模型结构简单对称，可采用Mapping方式进行，通过不断试计算，并细化应力较大的部位，剖分结果为12 800个单元，13 041个节点。网格划分后的土体有限元模型见图1 2.5 土体接触分析 根据实验，由动量定律可算得钻头冲击土体的接触时间为0.029 s。为分析方便，钻头对土体的冲击可用一个脉冲力来模拟，变化曲线如图2所示。其中，F 2.6 ．应用广泛 ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元分析软件，可广泛用于土木工程、机械制造等科学研究。软件主要包括三个部分：前处理模块、求解模块（即分析计算模块）和后处理模块 3 土体应力分析 有限元分析计算土体节点密塞斯应力分布状况和土体的轴向位移分布状况。根据以上建立的几何模型，确定土体的边界条件，使钻头在冲击土体的过程中，只有在圆柱体范围内的土才受冲击力的影响，在其内有应力的重分布、质点位移。 有代表性的时间子步上土体节点密塞斯应力分布状况和土体的轴向位移分布状况见表1。密塞斯应力分布状况和轴向位移分布状况见图3和图4。 4 地层冲击试验结果 在时间后处理器中，可以对分析土体中某一个感兴趣的点进行时间追踪，从而得出在冲击过程中应力、应变和时间变化曲线；应力、应变和位移关系。由于钻头几何模型和加载的轴对称性，我们只取模型的右半边进行跟踪分析。由分析土体的网格剖分图可知，土样中心所在点节点号为2 019，记为M2019。为了研究冲击钻进对土样的扰动影响，沿分析土体的上端面取一路径，长度为钻头的外径，记为Surface。 Surface路径上的密塞斯应力分布情况如图5所示。由图可知，密塞斯应力在土样部位几乎都保持在应力最小值的水平，不到钻头底唇面部位应力值的1/10，而且分布情况几乎没变；而在土中钻头内外壁和土层接触处的应力几乎是突变的。可以看出，冲击对土样的扰动只是在一个很小的范围内的。 Surface路径的轴向位移分布情况如图6所示。该路径的长度是钻头的外径，其中点就是土样的中心。由图可知，分析土体的轴向位移都是负的，即为轴向压缩，其压缩量以钻头底唇面处为最大；土样中心部位为最小，而且土样中心部位的位移相对比较平缓。 对比分析发现，在同一时刻，土体中轴向和径向位移的值绝大部分都相差两个数量级。这意味着在冲击条件下冲击功主要用于轴向压缩土体，而径向的挤压作用则相对很小。