企业官网建设流程全解析

PFC模拟技术在岩土工程地基处理中的应用，最能体现颗粒流价值的方向之一就是振动密实法。我在某模拟项目X里要评估振冲密实方案对一处松散砂层地基的加固效果，设计文件上写着“影响半径约1m”，可这个数字是怎么来的、能不能通过调整振动参数扩大，现场经验给不了定量解释。于是我把振冲器建模、土体颗粒模型、振动输入三个部分全部搬进PFC，打算在颗粒尺度上完整演算一遍“振冲头插进去→高频振动→砂层被挤密”的过程。真正动手之后才发现，土体颗粒模型反而是相对成熟的部分，最容易被低估、也最影响成败的，是振冲器本身怎么建——它既不是墙，也不是普通颗粒，而是需要准确描述几何形态、质量分布和运动输入的刚性结构。这篇内容是我用Clump技术构建振冲器/振动杆的完整复盘，涉及建模选型、参数标定、结果解读和调试经验，适合正在做PFC离散元模拟、或者想用颗粒流方法评估地基处理方案的工程师和研究生参考。

1. 振动密实方案论证中，传统手段给不了的那“半句话”

1.1 经验公式和连续介质方法为什么不够用

振动密实法（行业内也叫振冲密实法）加固砂土，本质上是靠振冲器产生的高频振动，让砂粒摆脱粒间锁定状态，重新排列到更密实的骨架位置上。设计阶段最常用的依据是两条线：一条是经验关系，比如根据振冲器功率、地层类别和桩间距来估算加固影响半径，另一条是现场检测，通过标准贯入、圆锥动力触探或者载荷板试验来验证完工后的加固效果。这两条线各有短板：经验关系只能给一个大致范围，没法区分具体土性差异、振动参数差异和施工工艺差异；现场检测当然最可靠，但它发生在施工之后，如果发现效果不达标，返工成本已经落地了。

所以在设计阶段，很多人想把数值模拟当作一条新的论证路径。但常见的连续介质方法，无论有限元还是有限差分，遇到振动密实问题都有两个绕不开的坎。第一，振冲器在土体中高频扰动时，周围土体会发生大位移、大变形，网格畸变、单元交错、接触面反复开合，收敛性极差。第二，连续介质土体本构模型大多是从静力或准静态试验中拟合出来的，很难把振动频率、振幅、加速度与颗粒重排、孔隙率变化直接挂钩。即便强行把振冲荷载折算成边界条件，也只能得到一个“荷载作用下的宏观变形场”，模拟不了土体颗粒被挤得重新排列这一本质过程。

1.2 颗粒流模拟补上的那一块拼图

PFC的做法，是把土体看成由大量离散刚性圆球（ball）组成的系统，每个颗粒的运动都遵守牛顿第二定律，颗粒之间的相互作用由接触模型表达。这样一来，振动密实就不再只是一个“荷载-变形”问题，而是真真切切的颗粒位置重排、接触形成与破坏、力链网络转移的演化过程。你可以在模拟中直接统计某个深度、某个半径范围内颗粒孔隙率的实时变化，也可以观察接触数量的增减，这些微观观测量在连续介质方法里根本拿不到。

更关键的是，PFC可以把土体和振冲机具放在同一个物理框架下处理。振冲器不再被简化成一个固定边界或等效荷载，而是有质量、有形状、有表面摩擦特性的刚性结构，土体颗粒在振动下产生的力学响应反过来也会作用在振冲器上。这种“机具-土体耦合”的模拟，对地基处理类问题意义很大——毕竟振冲器发出的振动能量到底有多少传进土体、以什么方式传递，取决于它和土体的接触状态，而不是我们单方面给的边界位移。

当然，PFC也不是拿来就能用的万能工具。它对颗粒数量、接触参数、阻尼设置极其敏感，颗粒模型本身如果标定不好，结果就只是一道漂亮的伪命题。所以我这篇内容的重点，是把实际操作中“把振冲器建模做对”的经验梳理出来，尤其是用Clump构建振冲器/振动杆背后的逻辑、流程和常见坑。

2. 不是一个铁疙瘩那么简单：振冲器建模选型分析

2.1 振冲器的运动特征决定了建模思路

先回到工程现场想一想，振冲器在地基处理时到底处于什么运动状态。振冲器内部有偏心块，电机驱动偏心块高速旋转，产生周期性激振力；整个振冲器连同上部振动杆被吊车吊着，在自重和激振力共同作用下逐步沉入土层。激振力是往复的，有竖直分量也有水平分量，施工过程中还要反复上拔、下插，因此振冲器和土体之间的接触界面一直处于高应力、高频交替的加载状态。

这种运动特征对模拟提出的直接要求是：振冲器不能是静止的“占位空间”。它必须在振动输入下于土体内部做高频往复运动，同时外表面几何必须足够准确，尤其前端锥体的形状。锥角大一点还是小一点，直接决定颗粒是被排向两侧还是被压向下方，这种几何差异在连续介质建模里通常被一笔带过，但在PFC的颗粒级模拟里根本无法回避。

2.2 四种可选表示方式：墙、bonded-ball、刚性块与Clump的取舍

我实际比较过PFC里可以模拟振冲器的几类对象，差异相当明显。

wall的实现速度确实最快，用三角网格片拼出圆柱表面就行。但墙本质上没有质量和惯性，你让它动它就动，土体反力对它的影响无法体现，做高频振动模拟时墙面与颗粒的接触行为也容易产生不自然的约束，除非只是演示步骤，否则我不推荐。

bonded-ball方案是把一批ball用平行粘结键连成一个整体，优点是接触计算与土体颗粒完全统一。但它的代价是内部接触大量参与迭代，计算量明显偏大；更麻烦的是，在30Hz高频强振动下平行键反复承受交变应力，很容易超过强度发生破裂。虽然可以人为把键强度调得极高来避免破坏，但那种“强制不裂”的刚体，动态响应会变得很别扭。

rigid block适合描述可以转动、滑移的刚体，但它对表面形状的还原依赖网格划分质量，振冲头这类锥台过渡几何建起来并不顺手，还要额外处理与ball之间的接触参数匹配。

最终我在模拟项目X里选了Clump。Clump由多个pebble刚性拼成，pebble之间没有相对位移和转动，整体作为一个刚体运动，接触计算只发生在对外表面与外部颗粒之间。它既能用不同半径的pebble逼近复杂几何形状，又不会引入内部接触的额外迭代，计算效率和稳定性都符合高频振动场景的需求。

2.3 Clump的底层原理与关键概念

给刚接触PFC的朋友补一句底层原理。Clump可以理解成“用若干圆球颗粒粘成的一个刚体糖块”，这些圆球叫pebble。pebble之间的重叠不会产生接触力，生成后位置关系保持不变，整个clump只有一个重心和一套惯性张量，运动由一个刚体方程控制。

这里的核心优势在于，你可以把振冲器这种金属机具近似为完全刚性体来处理。金属结构在振动过程中的自变形相对于土体颗粒的位移确实可以忽略，直接采用刚性假设，既符合物理，又能节省大量计算量。

使用Clump时有几个参数必须理解清楚。pebble半径决定表面分辨率和计算效率，半径太大表面就会呈锯齿状；pebble重叠度会影响体积和质量统计；整体密度则决定clump总质量。振冲器的总质量和转动惯量会影响它与土体的耦合响应——如果我只关心影响半径倒还好，但如果要分析振动器在自重下的下沉速度、上拔时的阻力，质量分布就必须跟真实机具匹配，不能随手给个大密度了事。

3. 基于Clump构建振冲器和振动杆的实操流程

3.1 先把振冲器拆成三段几何

一个大尺寸clump一次成型不是不行，但效率和可调性都很差。实际操作中我习惯把振冲器拆成三个几何段，分别布置pebble，最后合并到同一个clump模板里。

• 振冲头（前端）：多为截锥形或半球形，是直接挤压土体的部分，这一段pebble半径需要细分，尤其锥形母线的曲率要尽量平滑，否则颗粒会在端头附近出现不自然的卡滞。
• 振动器本体（中段）：圆柱形，外径最大，侧壁是振动能量向四周传递的主要界面，这一段pebble可以适当放大，但不能出现明显棱边。
• 振动杆（上段）：细长圆柱，起传力作用。因为杆径较细，pebble半径要比中段更小，否则表面粗糙度会被成倍放大，影响细杆周围颗粒的排布。

3.2 用pebble填充几何轮廓并生成Clump模板

三段几何轮廓确定后，就可以用脚本批量生成pebble。这里给一段我在项目中用的思路，PFC新版本支持Python for PFC，用脚本控制比鼠标点击高效得多，也方便参数扫描时快速重建。

# 示意代码：按几何尺寸生成pebble列表并建立clump模板 import math pebbles = [] # 振冲头：截锥段，锥角约60度，从尖端向本体渐扩 z_start, z_end = -0.35, -0.05 for z in np.arange(z_start, z_end, 0.015): r = 0.10 * (z - z_start) / (z_end - z_start) pebbles.append([0.0, 0.0, z, max(r, 0.015)]) # 振动器本体：圆柱段，半径0.10m for z in np.arange(-0.05, 0.45, 0.02): pebbles.append([0.0, 0.0, z, 0.10]) # 振动杆：细圆柱段，半径0.03m for z in np.arange(0.45, 2.0, 0.01): pebbles.append([0.0, 0.0, z, 0.03]) # 创建clump模板，后续实例化时再赋物理属性 clump_template_create(name='vibrator') for p in pebbles: clump_template_add_pebble('vibrator', p)

这段代码的思路，是把几何体每一个高度截面近似为一个圆，圆的中心作为pebble位置，半径取该截面处的半径。截锥段通过线性插值从尖端的小半径逐渐放大到本体半径，保证表面连续；中段和细杆段用固定半径循环填充，相邻pebble在轴向有重叠也没关系，重叠加不会在clump模板中产生接触力，只影响体积和质量计算。

模板生成后一定要做一次几何质检。我的做法是取几个代表性高度，用PFC的几何接口输出clump表面点坐标，跟理论半径逐点校验，误差控制在小颗粒粒径量级以内。锥面如果出现局部半径波动偏大，就把该处的pebble间距缩小。这一步偷懒的话，后面土体颗粒在振冲头附近的排布就会出现假象。

3.3 物理属性与振动输入：质量、速度时程和激振力

clump模板只定义了形状，实例化之后还要赋物理属性。重点是整体密度，要让总质量与真实振冲器质量一致。比如某个型号振冲器总质量大约2t，轮廓体积算出来0.18立方米，那么密度直接按质量除以体积来设，不需要去纠结钢材真实密度，因为内部空腔和偏心块已经被体积平均化了。质量对动态问题影响很直接：质量太小，高频振动下clump容易被土体弹开；质量太大，惯性主导一切，土体来不及响应，影响半径会偏小。

振动输入方面，我优先用速度时程驱动。原因很简单：位移驱动在PFC的运动积分过程中容易在每一步产生速度突变，速度驱动则平滑得多。现场实测的振冲器加速度时程，可以先积分成速度时程，再按时间采样点加载到clump上。

如果采用简化的正弦振动假设，速度与位移的关系是相位超前90度，也就是v_z(t) = 2πf·A·cos(2πft)。在PFC里大致可以这样处理：

while step_condition: time = current_time() v_z = 2 * math.pi * freq * amp * math.cos(2 * math.pi * freq * time) clump.set_vel(0, 0, v_z) step(1)

freq和amp分别取施工振动频率和振幅，比如30Hz、3mm振幅时峰值速度约0.565m/s。如果要模拟振冲器的下沉过程，还需要叠加一个竖向贯入速度，但必须控制好叠加方式，别把它变成“反复打桩”。

如果要用激振力驱动，就更贴近现场。激振力大小由偏心块的偏心质量和转速决定，Fe = me·e·ω²，这个力可以加在clump重心上，配合自重和重力场让系统自行振动起来。这种方式的好处是振幅不由外部限制，可能更真实地反映变幅特性；代价是数值收敛更难，参数稍微不对就会出现发散。

3.4 clump与土体接触的三个关键设置

第一，clump与ball之间的接触模型建议用线性接触（linear contact），有效模量设为土体颗粒间接触模量的5到10倍，摩擦系数按钢-砂界面取值，一般0.3~0.5。接触间隙设为零，避免初始嵌入产生巨大的排斥力。

第二，放置clump之前，要把它所在位置的土体颗粒先删除，再把clump“预埋”进去。如果不删直接压入，初始接触力会大得离谱，系统平衡阶段就会把颗粒弹出几米远，后面就全乱套了。

第三，高频振动会让clump周围颗粒的接触频繁开合，如果发现颗粒不断穿越clump表面，多半是接触模型的有效模量太低或者时间步太大。穿透不光是动画难看，更会让接触力累积失真，必须尽早处理。

4. 土体颗粒模型与接触参数标定：决定成败的地基

4.1 接触模型选型要结合土性

振冲密实处理的对象以砂土、砾石土为主，我在PFC里通常用线性接触模型描述颗粒间的基本力学响应。它的物理图像很清楚：法向力正比于颗粒重叠量，切向力的上限由摩擦系数决定。再加上滚动阻力，就能近似模拟真实砂粒棱角的咬合效应。

如果对象是黏性土或含细粒较多的粉土，线性接触模型会明显偏弱，因为缺乏黏聚和抗拉强度，这时要考虑平行键模型（parallel-bond）。但加了bond以后模型标定更繁琐，而且振动反复作用下键会断裂，你得判断断裂后的堆积体是否还等于原来的土体。所以，只要目标地层是砂土，就别画蛇添足。

我个人的选型判断标准是：孔隙比大、颗粒互锁强的纯净砂土，linear加rolling resistance足够；饱和砂土要引入CFD耦合或对结果做排水条件修正；黏性土建议慎重使用干颗粒PFC模拟的结果。

4.2 从三轴试验到PFC参数：一套可复用的标定流程

参数标定不能靠感觉。我的流程分四步。

第一步，建立数值三轴试样，颗粒级配和现场模型保持一致，试样尺寸参考室内三轴标准，比如直径50mm、高100mm，生成时目标孔隙率设为现场松散状态，例如0.42。

第二步，让试样在目标围压下固结平衡，再施加轴向压缩，记录偏应力-轴向应变曲线和体应变曲线。

第三步，把数值曲线和室内试验曲线对照。先调有效模量匹配初始切线刚度，再调摩擦系数匹配峰值强度，最后用刚度比和阻尼来匹配剪胀特征和峰后软化。每个参数都有明确的宏观对应，不要同时往多个方向乱调。

第四步，如果数值强度始终低于室内结果，通常是因为圆球颗粒滚动自由度太大，需要引入滚动阻力，或者把摩擦系数略微提高，来等效模拟真实砂粒的棱角咬合。我用滚动阻力系数0.2~0.4、摩擦系数0.5~0.7的组合，通常能得到和标准砂三轴曲线比较接近的结果。

标定完成不是终点。颗粒模型参数往往带有压力依赖性，振冲模拟涉及的深度范围较大时，我会在浅层、中层、深层分别标定一组，确保不同深度处的刚度响应都合理。

4.3 阻尼设置：振动模拟里的隐形调节阀

PFC里的能量耗散主要靠接触粘性阻尼和局部阻尼。静力模拟常用局部阻尼0.7来快速“刹车”，但这个值放到振动模拟里，会把振动能量耗散掉大半，颗粒几乎动不起来。我在做振动密实时局部阻尼通常设在0.15~0.3之间。数值很讲究：太大会让颗粒重排偏慢、影响半径偏小；太小会让系统长时间振荡，计算效率很差。

一个实用校核方法是：先把阻尼设到极小，跑一小段看颗粒速度的量级和波动特征，再逐步加大阻尼，找到一个“系统能稳定下来但颗粒位移量没有明显减少”的临界值，就用那个值。接触粘性阻尼一般只给法向很小的值，切向给更小，因为振动时切向摩擦本身的耗散已经足够，再加粘性阻尼容易让响应发“僵”。

4.4 初始条件与边界：为什么先平衡那么重要

地基模型生成后的第一件事，不是开振冲器，而是让土体在重力和边界约束下充分平衡。很多新手跳过这一步，直接加振动，结果会看到颗粒在振动初期疯狂飞溅——那不是土体被振松，而是初始接触力还没稳定，系统内部应力在释放。

边界设置方面，底面用固定墙或粘滞边界，侧面尽量远离振冲器。模型平面尺寸至少取预期影响半径的3倍，深度至少在被加固深度以下再加2倍影响半径。模型太小，边界反射波会明显改变影响半径曲线，数据可信度大打折扣。

5. 数值结果如何回答工程设计的三个核心问题

5.1 加固深度够不够：密实度随深度的演化

密实度变化是振冲模拟最直观的输出。在PFC里用测量球（measurement sphere）沿深度布置一排，统计每个球覆盖范围内颗粒的孔隙率，记录振动前后的差值Δn。

我在模拟项目X里的一组结果是：初始孔隙率0.43的松散砂土，振冲头在1.0m深度持续振动10秒后，振冲头所在位置的孔隙率降到0.37左右，往上0.5m处降到0.39，往下0.5m处降到0.41。这种梯度分布说明加固效果集中在振冲头活动段附近，符合现场经验。如果发现孔隙率整体下降缓慢但能量输入很高，通常不是土体“不吃力”，而是刚度或阻尼把能量以数值形式耗散掉了，得回头查参数。

5.2 影响半径有多大：水平向密实度剖面

平面上的影响半径，可以沿中心向外的不同半径处放测量球，统计孔隙率降幅，得到一条径向衰减曲线。按孔隙率降幅5%作为有效边界，就能画出一个等效影响半径。

这个虚拟值与工程经验的对照非常有用。同样在模拟项目X里，设定激振力300kN、频率35Hz，模拟得到的影响半径约1.0~1.2m，跟同一参数下现场试桩后检测出的加固范围非常接近。这里要说清楚：模拟不能替代试桩，但它能在试桩之前把试验方案筛得更准，减少无效试验的浪费。

5.3 振冲参数怎么优化：频率、振幅与持续时间的敏感性

有了稳定的PFC模型，参数敏感性分析就变得非常便宜。把地基模型复制成若干份，分别跑不同频率、振幅、振动时长，然后对比孔隙率、影响半径、有效加固时间。

我在实际对比中发现：频率从25Hz提到50Hz，影响半径的增长不如振幅从3mm提到6mm来得明显；持续振动超过15秒后，孔隙率下降速率明显变缓，继续延长时间意义不大。这些结论如果靠现场试验验证，至少要打几十根试验桩；而用PFC跑，一组脚本就能看到趋势。虽然绝对数值不能直接当验收依据，但可以大幅缩小现场试验的搜索空间。

这种参数扫描的关键前提，是基准工况已经经过标定和验证。标定参数小幅偏移，可能会改变不同工况的优劣排序，所以敏感性分析之前务必先把基准工况校核扎实。

6. 高频振动数值模拟的典型问题和调试方法

6.1 clump刚度与颗粒穿透

高频振动下最容易出现的现象，是振冲器周围的颗粒逐渐“钻进”clump表面。颗粒坐标不断穿越边界，基本可以断定是clump与ball接触的有效模量不足，颗粒在接触力作用下过度嵌入。解决方法是逐步提高有效模量，同时注意时间步也会变小，计算时间变长。

如果穿透只发生在某个局部，比如振动杆上段，那多半是pebble半径和土体颗粒半径的比例失衡。杆子太细、pebble太粗时，表面呈锯齿状，小颗粒容易卡进齿缝。处理办法是局部加密小半径pebble，让杆面更平滑。

6.2 计算时间太长怎么办

高频振动模拟的计算量很头痛。30Hz振动10秒，每个周期大约0.033秒，而PFC的临界时间步通常只有10的负6次方秒量级，一个周期就要跑上万步。我总结出三个对策，按优先级排序。

第一，缩短模拟时长。观察密实度变化率，只要孔隙率在连续几个周期内下降不足0.5%，就可以提前结束，没必要跑完真实振冲时长。第二，在不改变模型尺寸的前提下适当增大土体颗粒半径，减少颗粒总数，但必须做颗粒尺寸敏感性验证，确认影响半径和密实度结果仍在可接受范围内。第三，局部质量缩放，只对远离振冲器的边界区域颗粒放大质量，可以提高全局时间步，又不影响核心区域的动态响应。

6.3 边界反射波的识别与消除

如果影响半径曲线上出现孔隙率降幅在某处不降反升，或者曲线有周期性起伏，多半是边界反射波造成的。把振动结束后不同时刻的颗粒速度分布画出来，如果看到明显的波前从边界返回中心，就说明反射已经干扰结果了。

处理办法是在边界设置粘滞阻尼墙，或者在边界外侧加一圈吸收层，厚度按主频和土体波速估算。加吸收层之后，影响半径曲线会明显更平滑。这个细节虽然不起眼，但直接影响影响半径的判读精度。

6.4 饱和砂土与排水条件问题，必须提前想清楚

PFC纯颗粒模拟默认颗粒之间没有流体介质，相当于干砂、排水条件极好的工况。但真实饱和砂土地基里，高频振动会引发超静孔隙水压力，有效应力降低，颗粒更易移动；排水条件好时，超孔压快速消散，挤密效果明显。干颗粒模拟和饱和工况的结果可能差出一个数量级。

所以如果工程对象是饱和砂土，至少要做到两点：一是解读纯PFC结果时明确标注“排水极限工况”的适用范围，二是条件允许时叠加CFD-DEM耦合，把流固耦合纳入考虑。这是PFC在振动密实领域最值得深入研究的方向，也最容易被低估。

最后再分享一个我个人习惯的验证小技巧：每次跑完一组振动密实模拟，我都会把频率、振幅、clump刚度、局部阻尼、影响半径这几个关键量记在同一个表格里，立即做横向对比。这样当别人问“你那个模拟到底靠不靠谱”时，你可以拿出参数和曲线告诉他：这个数据跟某深度标贯结果吻合，那个结果提示频率需要下调。数值模拟本质上是在给工程论证提供辅助依据，自己先把流程和数据管理得足够严谨，它才能真正成为顺手的工具。