当前位置：文档之家 › PFC3d手册命令中文详解95406

PFC3d手册命令中文详解95406

  • 格式：doc
  • 大小：765.00 KB
  • 文档页数：54

PFC3D 学习记录5－COMMAND REFERENCE通用规定和特征命令语法默认情况下，命令行忽略大小写，然而可以通过SET case命令激活区别大小写情况。

斜体字代表数字，以i, j, m n 为开头的字母代表整数，其他情况为实数命令词，关键字和数值可用空格，以及（），＝分隔。

在分号（;）后为注释一行最多可包含80个字符。

而&表示续行Range 语法定义一定限制范围的物体可运用Range的命令：Change/ initial/ property Clump Fix/ Free Group JsetModel Plot ( 运用于不同输出项目)PrintRange( 创建一个命名的范围)一个Range 定义一系列物体。

由许多范围元素构成。

如果使用多个范围元素，最终物体几何为不同元素集合的交集。

但是可通过关键词any改变选择特性。

如果使用关键词any选项，任何包含在范围内的单元可被认为最终物体集合的一部分。

即求合集除此之外，通过在定义范围单元后接关键词not确定所有不在指定范围的单元通常还可以用定义group 方式定义复杂的range。

内置范围单元：Annulus center ( x, y, z) radius = r1, r2球形空心区域中心在（x, y, z ）内径r1, 外径r2Color il < iu >球体或墙体在间隔（il, iu ）内颜色索引。

若无iu 则iu = il球体颜色索引列表通过property命令建立。

墙体颜色索引列表必须用FISH 内置的w_color设置Cylinder end1( x1, y1, z1 ) end2( x2, y2, z2 ) radius = r圆柱范围。

圆柱中心轴过end1( x1, y1, z1 ) end2( x2, y2, z2 ) 半径为r.Fish fname调用用户定义的范围单元。

通过FISH函数fname 执行。

而此函数中必须包含一个range_element指定范围单元函数。

此时fname就会调用所有在范围内的各个物体。

物体的地址将通过fc_arg( 0 )传入函数，而坐标x-, y-, 和z-分量将通过fc_arg( 1 ), fc_arg( 2 ), fc_arg( 3 )分别传入函数。

物体类型由FISH内置的pointer_type指定。

如果物体不在指定范围内，则fname的值必须在函数内设置为零，如果未指定，默认值为１。

group gname定义属于group名字gname的颗粒球体，（粒子组通过group命令定义）id il < iu >id编号在il，iu范围内的物体，如果iu没有指定，则默认il = iujset il < iu >节理编号在在il，iu范围内的接触，如果iu没有指定，则默认il = iu节理通过JSET命令创建创建一个命名范围rname，通过在命令行中范围元素集合定义，通常一个范围用range命令创建。

plane keyword无限平面包围区域Keyword:Origin x, y, z 平面上一点（x, y, z ）dd dd 方向角x-y平面上方向角，从y轴正向顺时针方向测量dip dip 俯角从x-y平面，向z轴负向。

Distance d 在距离d范围内区域Above 平面上方below 平面下方normal nx, ny, nz 平面的单位法向矢量radius rl < ru > <tolerance t >半径在（rl, ru ）范围内球体。

ru未指定，则范围（rl – t, ru – t ）t 未指定，则t = rl ×10-6. 注意误差范围只有在未指定ru 条件下指定。

rname 命名范围rnamesphere center( x, y, z ) radius = r球形区域，中心（x, y, z ）半径rx = xl, <xu >, <tolerance = t >y = yl, < yu >, < tolerance = t >z = zl, < zu >, < tolerance = t >指定x, y, z 坐标在指定范围内的物体。

其他同radius 条重点理解fish定义range部分。

在线帮助Command ? 适用command的各个关键词？/ help 显示所有可用命令程序控制命令这些命令支持批处理执行模式，允许改变程序状态。

其与用程序建立实际数值模型无关。

Cycle & Step 执行指定时间步Solve执行计算，直至满足收敛条件Build 控制AC/DC 操作call 调用data文件，执行命令，批处理模式return 从批处理模式中返回到控制台。

或在多级调用中返回上一级调用文件continue- pause重新/ 中止调用data文件Fluid设置流体分析条件/ MPI 控制MPI 操作/ PARALLEL控制并行计算程序操作NEW 重启新模型Quit & Stop / Save – RestoreSet 改变全局程序设置/ Thermal设置热分析/ Title 设置模型名称，这些将保存在随后的输出图并记录在save文件中/ Ucode注册用户编译的C++ 类特殊计算模式可选计算模式，通过config命令调用。

Config 可在任何分析阶段给出，但是必须在调用时指定计算模式。

Cppudm C++用户定义接触模型Cppuwc C++用户编写代码fluid 流体分析thermal 热分析模型属性部分模型建立和修改Ball 建立新球体Clump 建立新团块，或修改新团块属性Delete 删除球体，团块，墙体和时程记录Generate 生成球体集合Jset 生成一系列节理Wall 生成墙体或修改墙体属性Fix 对球体设置固定速度Free 将球体固定速度清零Model 对指定接触指定用户定义接触模型Property 修改已有球体属性（材料属性和施加外力和速度），修改属性，粘结属性，接触属性模型监控History 指定执行过程变量（预定或用户定义fish变量Measure 建立一个测量球或改变测量球属性Movie 对屏幕输出进行控制抓图，从而实现动画输出Plot 控制模型几何和反应属性输出Print 生成输出变量和模型信息Trace 追踪模型中能量及功的输出命令工具Define 定义fish函数Group 定义和修改组集合Macro 命令宏Range 定义范围Table 控制表的创建和修改。

成对数值索引命令集合建立球体Ball radius = r <density = ** hertz id = ** x = ** y = ** z = **>Radius :球体半径density : 密度（m/V）Hertz: 指定Hertz接触模型否则为线性接触模型Id: 编号，唯一性，默认值大于当前编号最大值的整数x, y, z ：x-, y- z-坐标，默认为0.0用ball 建立球体允许重叠，球体属性由property命令设置Build 控制执行AC/DC 操作Build assembly xa, ya, < za >通过复制粒子和包含在pbrick中的力（必须已经存在）建立粒子集合，而初始粒子由Call <fname> 调用data文件，默认为pfc3d.dat, 而若扩展名未指定，等于*.dat文件调用允许多级调用，但不允许相互迭代。

Change 同property命令一样创建团块Clump <id = **> full = on / off permanent release update_cycle = n<range …>创建一个团块，若id等于已有团块编号，则为修改团块属性Range 为要建立团块的范围，范围内原来属于其他团块的粒子转化到当前团块，必须指定range防止所有粒子变成一个团块Print clump 输出团块信息，plot clump 输出团块1、Full on 在计算团块运动时，使用全部惯性向量力矩，这会耗费更多计算时间。

因为cross-coupling term使用，必须进行矩阵转换。

当旋转速度较大和陀螺效应影响大时使用。

Full off 计算团块运动时，只有使用diagonal inertia，而忽略旋转速度乘积。

只有3D计算受full关键词影响。

当full on 时，应设置update_cycle = 12、Permanent 通过阻止构成团块粒子间生成接触使得构成团块的粒子在随后模拟过程保持接触。

这样不需要存储接触信息，并且避免经常性更新接触信息及其对应的矢量。

当在模拟过程中团块有个形成和分解过程，就不能指定这个关键词然而，指定颗粒间永久接触有利于节约内存和计算时间。

当一个团块有其permanent属性设置，其设置为永久的然而永久团块可被删除且团块内球体可被释放。

Release关键词用于从团块提取球体，释放时，球体将不存在接触3、release 范围内所有属于同一团块将被释放。

如果id未指定，所有范围内属于任何团块的粒子将被释放。

团块内所有粒子被释放时，团块本身不被删除，形成“空团块”, 这种情况在计算过程中忽略，而在粒子重新加入之后，被激活。

Delete clump 删除团块4、update_cycle n控制数据更新速度，默认值n= 10。

N越大，计算时间小。

可被更新速度有Moment of inertia tensor 惯性张量力矩（3D）contact coordinates 接触坐标contact normal vector 接触法向矢量接触剪力矢量(3D) 平行粘结力矩（3D）平行粘结剪力矢量(3D )注意团块阻尼系数为构成团块粒子的阻尼系数的算术平均值当模型中存在团块时，在本版本中就不能使用密度放大系数Fix/free命令用于团块内任何粒子时，此工况用于团块重心。

Print clump 可看到各个团块的fix/free状况。

Config keyword < keyword >指定可选计算模式Keyword:Cppudm C++用户定义模型Cppuwc C++用户编写代码Fluid 流体分析Thermal 热分析Continue 重启由pause中止的data文件过程Cycle/Step n执行n 时间步对于热分析和热力分析计算为热计算步，否则为力学步计算过程输入Esc，中止计算过程。

计算过程输入空格，继续中止当前cycle.当前计算过程由set mechanical 和set thermal命令指定。

合集下载

最新版PFC Command 命令大全

最新版PFC Command 命令大全

PFC3D命令说明（COMMON COMMAND REFERENCE）PFC3D是基于命令驱动模式（COMMAND-DRIVEN FORMA T）的软件，各种命令控制着程序的运行，这部分内容将介绍PFC3D软件的内嵌命令。

本说明文件译自PFC3D软件2.0版使用手册中的COMMAND REFERENCE部分，并补充了一些手册中没有的命令（如设置粘性阻尼、生成圆柱、螺旋壁面等，这些命令在更高版本的手册中有说明）。

命令说明的顺序没有采用原手册中按字母排序的方式，而是根据创建PFC3D模型解决实际物理问题的一般过程，对相关命令加以说明。

PFC3D手册中COMMON COMMAND REFERENCE只对每个命令的格式和基本功能做了简单介绍，本说明文件对每个命令做了更详细的解释，为保持文件的可读性，对命令的详细解释都以附录形式给出。

本文件介绍了PFC3D的基本功能，对初学者有较大帮助，但要进行高级应用，还需清楚了解DEM 和PFC3D相关功能的基本原理以及软件的结构等。

由于只关注使用PFC3D解决颗粒流的问题，涉及颗粒流的命令介绍的比较详细，而用于岩土工程等其他领域的命令，由于关注很少且专业知识缺乏，只是不求甚解。

0.关于命令格式的说明：PFC3D中每个命令基本包含3部分：命令名、必需关键词和可选关键词。

本说明文件中，命令名和必需关键词写在第一行，尖括号<Keyword…>表示该命令具有可选关键词，罗列于命令名的下一行次级列；带参数的关键词，参数紧随关键词之后用加粗斜体表示。

名字较长的“命令”和“关键词”，PFC3D允许不写出全部字母，只需简写前面的一些字母，最短可简写名（即PFC3D能识别的最短简写“命令”和“关键词”名）用加粗字体表示。

例如命令：BALL rad r<keyword…>…hert z …id idx xy yz z说明：命令名为BALL的命令，具有必需关键词rad（参数r）和hert z、id、x、y、z等可选关键词。

PFC3D中文讲解-推荐下载

PFC3D中文讲解-推荐下载

4、定义墙面的法向刚度（normal stiffness）和切向刚度（shear stiffness） wall id=1 kn=1e8 ks=1e8 wall id=2 kn=1e8 ks=1e8 wall id=3 kn=1e8 ks=1e8 wall id=4 kn=1e8 ks=1e8
wall id=5 kn=1e8 ks=1e8 wall id=6 kn=1e8 ks=1e8
wall id=2 face ( 0, 0, 0) ( 0, 0,-5) ( 0, 2,-5) ( 0, 2, 0)
wall id=3 face (10, 0,-5) (10, 0, 0) (10, 2, 0) (10, 2,-5)
wall id=4 face (10, 2,-5) (10, 2, 0) ( 0, 2, 0) ( 0, 2,-5)
每个命令的作用），以下是程序代码。
;fname: footing.DAT (tutorial example for PFC3D)
new
; clear program state to begin new problem
set random ; reset random-number generator
其中 kn 表示法向刚度，ks 表示切向刚度，单位为 N/m 5、定义球的密度，法向刚度和切向刚度 prop density 2000 kn 1e8 ks 1e8 使用命令 property，可简写为 prop，定义密度 density 为 2000kg/m3，法向刚 度和切向刚度均为 108N/m 6、改变小球半径 ini rad mul 1.47 ini 应该是 initial 的简写，用于初始化小球半径（有待考证） rad 表示半径 mul 是 multiply 的简写，是倍数的意思。 1.47 是倍数因子，表示现在半径是原来半径的 1.47 倍。 7、设定监测变量和采样间隔时间 set hist_rep=5 hist ball zvel 3 1 0 hist diagnostic muf hist 是 history 历史记录是简写，set hist_rep=5 表示每隔 5 个运算步记录一 次变量的值（默认的是 10 个运算步记录一次）；本程序设定了两个监测对象， 一个是 hist ball zvel 3 1 0，表示的是中心坐标最接近 x=3，y=1，z=0 的球的 z 方向的速度被监测和记录（在这个例子中，282 号球最接近这个坐标，因此作 为监测球）；另一个是 hist diagnostic muf，是用来监测系统的不平衡力，主要用 于静力学分析，当这个值趋于很小时，表示系统达到稳定状态。 8、设定求解选项 set dt dscale 因为这是一个静态分析，因此我们开启 density-scaling 选项（简写为 dscale）来优化求解效率，dt 应该是这个选项中的某个模式（有待考证）。 9、设定求解步数 cycle 3000 当写入这句程序后，计算就开始了，计算的步数，时间步，总时间，不平 衡力平均值和最大不平衡力这些信息均显示在屏幕上，并且每隔 5 个计算步更 新一次（这是 set hist_rep=5 的作用），图像动画每隔 20 个计算步更新一次（默 认），这个更新时间可用命令 set pinterval 来改变。 另外，阻尼已经在运动方程中被默认设定，用于吸收振动能量，更快地达 到平衡状态。 10、查看结果 求解完毕后，用户可输入命令 plot his 1 来查看历史记录中第一个变量的数 值曲线。
PFC3D建模流程和代码解释

PFC3D建模流程和代码解释

本示例表示生成球体在一定的空间内进行运动的部分代码以及代码表示意义。

新建一个PFC3D模型必要部分生成的整体日志。

PFC3D建模需要的部分和代码解释*********************************************** Logging started at 周三八月30 16:50:55 2017* By pfc3d Version 5.00 Release 025*** Job Title: Pandect example**********************************************pfc3d>pfc3d>; 3-设置模型名称，可有可无。

pfc3d>title 'Pandect example'pfc3d>pfc3d>; 4-设定计算区域(必要条件)pfc3d>; Set the domain extentpfc3d>domain extent -10.0 10.0 -20 20 ; domain当中有两个关键词是condition和extent，condition指的是边界条件periodic ;destroy ;stop ;reflect；condition默认值是stop。

extent指的是模型区域pfc3d>pfc3d>; 5-指定随机种子（若不指定，种子随机）pfc3d>set random 10001;默认值是一万，且数量级是相同的数量级，也就是该随机数值10000不能过大，不能过小。

设置随机种子以后，生成的颗粒的半径和位置都在范围内随机。

pfc3d>pfc3d>; 6-生成及修改模型组件(必要条件)，有三种命令是来生成颗粒的，generate，create，distribute。

pfc3d>; Generate 30 balls in a box，wall就是一个生成颗粒的容器。

PFC3D操作指导

PFC3D操作指导

PFC3D操作指导与flac3d的2.0和2.1版本相比，flac3d3.0有哪些新功能？答：在计算方面：1.所有的计算和数据均采用双精度浮点制。

2.其运行速度较2.1约快10-15%功能方面：1.动力模块中增加了hysteretic阻尼在动力学计算中，增加了一个新的阻尼特性：超圆柱阻尼。

有了这种形式的阻尼，数值模拟产生的应变可以基于模型模量和阻尼函数的共同作用。

这允许用户将等效线性计算方法的结果与完全非线性非简化本构模型的结果进行比较。

另外，对一些其它形式的阻尼，如rayleigh阻尼，进行了较大的折减，这样在使用hysteretic阻尼时可以有效地节省计算时间。

2.岩石本构模型中增加了Hook Brown模型加入了霍克-布朗屈服准则，使用户对岩石材料计算结果较为合理。

它的屈服面是非线性的，同时是考虑最大、最小主应力的关系的基础上产生的。

该模型综合了塑性流准则，其特性随着不同的应力水平，呈一变化的函数。

3.热/流体水平对流模型flac3d之前已经能提供非线性固体与多孔介质渗流的流固耦合模拟，机械地耦合流体和固体。

而3.0版本的新功能增加了温度可受流体密度影响和流体中温度的水平对流。

4.3 dshop生成的网格导入3dshop是一种能力强大的固体建摸和六面体网格的软件包，也是itasa的产品。

flac3d原始的网格仍可用，但是用基于windows操作系统的3dshop建复杂模型更为简单方便。

3dshop生成的网格能被flac直接读取。

5.动画命令：FLAC3D现在可以生成AVI和DCX动画，这是以前版本中不可用的功能。

6.记录颗粒轨迹：颗粒的路径能被记录和显示此外，flac3d3 0还包括以下新功能（在以前的版本中无法实现）：1）网络版本2）应用于混凝土加工模拟的水合作用模型。

Q：如何显示变形等高线命令？答：ploskemagf10其中10为放大系数问：如何查看各个时段不平衡力的具体数值？答：his用于记录计算结果，包括位移、应力和其他命令hisunbalhisgp(zone)zdisrange(000)或者id=？导出数据命令hiswitenvsmbegin time step end time step filename Hisn表示记录的IDM，并表示推导不平衡力hisunbal特定值的时间步长step100000orsolve他的write1vsstep开始于1000file123。