PFC3D概述简介.docVIP

PFC3D概述 简介（INTRODUCTION） PFC2D/3D软件是美国ITASCA公司出品的商用DEM软件，适用于研究颗粒集合体（Particles assemblies）的破裂和破裂发展问题，以及大位移的颗粒流问题，是一款用于模拟复杂固体力学问题和颗粒流问题的有效工具。 PFC3D研究的基本对象是颗粒（Particles）和颗粒间的接触（contact），它能直接模拟球形颗粒间的运动和相互作用的物理问题；可以通过“连接”（Attach）两个或多个小颗粒来创建任意形状的大颗粒，“连接”而成的“组合颗粒”（group of particles）可以作为独立的颗粒体研究。PFC3D可以模拟固体的破裂问题。通过“粘结”（bond）相邻颗粒得到的颗粒集合体可作为具有弹性属性的“固体”（solid），当颗粒间的“粘结”（bond）逐渐破坏时，该固体即产生“破裂”（fracture）。 …… PFC3D的基本特征 颗粒（Particle） PFC3D可模拟任意大小球形颗粒的集合体的动力学行为，可以自动生成统计学上特定分布形式的颗粒集合，颗粒的半径可以均匀分布（即颗粒半径都相等），也可遵循高斯分布。任何状态下，颗粒的半径都能以任意大小的系数进行膨胀，并由此改变空隙率（porosity，或填充率）的大小。通过膨胀颗粒半径，不用通过反复试验即可得到指定空隙率的致密填充状态。 PFC3D对单个颗粒（particles）或接触（contact）赋予物性参数，而不是“型数”（type numbers）。因此，颗粒物性以及半径可按照指定的方式连续且有梯度地变化。绘图时的颜色也被作为一种物性赋予颗粒，因此用户可以使用更多样化的分类方案（marking schemes）【注：我们的研究中就使用多种颜色的分类方案，详见后续说明】。 接触（Contact） 颗粒间的接触是DEM中的另一个主要研究对象，PFC3D软件中提供了多种用于计算颗粒间接触列的接触模型，主要有： 线性弹簧—阻尼器模型（linear springs）和简化Hertz-Mindlin模型（simplified Hetz-Mindlin law）； 库伦滑移模型（Coulomb sliding）； 两种可选择的“粘结”模型：接触粘结模型和平行粘结模型； 任何用户自定义的接触模型。 接触也被赋予物性参数，并可以读取和修改。如使用粘性阻尼时的粘性系数是赋予接触点出的能量耗散系数。 壁面（Wall） PFC3D建立的模型中，其边界条件由壁面组成。与颗粒、接触一样，壁面也被赋予物性参数，如法向和切向刚度系数、摩擦系数、颜色编号等。 PFC3D可生产多种形式的壁面，主要包括： 有限平面，一般是由有限个点连接而成的任意多边形，有限平面在使用时应注意遵从右手法则； 无限平面，指定平面的法向方向以及平面上的一个点，即可生成一个无限大的平面。无限平面； 圆柱面、圆锥面、圆台面，这些面可以通过使用同一个命令生产，这些回旋面作为一个整体而存在，并不是若干个分割小平面拼接而成； 其他形式的复杂面，通过FISH语言允许用户编写形式更复杂的边界面； 每个壁面有其特有的接触特性。可以指定壁面速度，并可以监测作用于每个壁面的合力和合力矩。颗粒和壁面可以在模拟过程中的任意时刻创建或删除。 时间步长（Timestep） PFC3D默认自动确定计算的时间步长。时间步长自动确定，并包含了由于Hertz接触模型产生的刚度系数变化效应(The timestep calculation is automatic and includes the effects of changing stiffnesses due to the Hertz contact model )。另外，模拟过程中，时间步长也随每个颗粒周围的接触数量以及瞬时刚度值而变化。 PFC3D中也可以使用固定时间步长。PFC3D自动计算得到的时间步长对于特定问题并一定是最佳（？是如何计算的？为何不是最佳？），自动时间步长为了保证数值稳定性，往往偏保守，致使模拟的计算时间很长。对于特定问题，可以使用比最大自动时间步长更大的合适的固定时间步长值。固定时间步长的选取是经验性的，需要先验证其正确性。 阻尼（Damping） PFC3D使用阻尼模型耗散颗粒的动能，使颗粒在合力的迭代步数内达到稳定的运动状态（静止[stationary]或恒稳运动[steady motion]）。PFC3D有3种阻尼模型可选择：粘性阻尼、局部非粘性阻尼和组合阻尼。 局部非粘性阻尼（Local non-viscous Damping）被作为一种属性赋予颗粒，局部非粘性阻尼具有以下优点： 稳态运动下的体力为0（只在加速运动时有阻尼作用）？？； 阻尼常数是无量纲量； 因为局部阻尼与频率无关，使用相