PFC3d手册命令中文详解

PFC3D学习记录5－COMMANDREFERENCE通用规定和特征命令语法默认情况下，命令行忽略大小写，然而可以通过SETcase命令激活区别大小写情况。斜体字代表数字，以i,j,mn为开头的字母代表整数，其他情况为实数命令词，关键字和数值可用空格，以及（），＝分隔。在分号（;）后为注释一行最多可包含80个字符。而&表示续行Range语法定义一定限制范围的物体可运用Range的命令：Change/initial/propertyClumpFix/FreeGroupJsetModelPlot(运用于不同输出项目)PrintRange(创建一个命名的范围)一个Range定义一系列物体。由许多范围元素构成。如果使用多个范围元素，最终物体几何为不同元素集合的交集。但是可通过关键词any改变选择特性。如果使用关键词any选项，任何包含在范围内的单元可被认为最终物体集合的一部分。即求合集除此之外，通过在定义范围单元后接关键词not确定所有不在指定范围的单元通常还可以用定义group方式定义复杂的range。内置范围单元：Annuluscenter(x,y,z)radius=r1,r2球形空心区域中心在（x,y,z）内径r1,外径r2Coloriliu球体或墙体在间隔（il,iu）内颜色索引。若无iu则iu=il球体颜色索引列表通过property命令建立。墙体颜色索引列表必须用FISH内置的w_color设置Cylinderend1(x1,y1,z1)end2(x2,y2,z2)radius=r圆柱范围。圆柱中心轴过end1(x1,y1,z1)end2(x2,y2,z2)半径为r.Fishfname调用用户定义的范围单元。通过FISH函数fname执行。而此函数中必须包含一个range_element指定范围单元函数。此时fname就会调用所有在范围内的各个物体。物体的地址将通过fc_arg(0)传入函数，而坐标x-,y-,和z-分量将通过fc_arg(1),fc_arg(2),fc_arg(3)分别传入函数。物体类型由FISH内置的pointer_type指定。如果物体不在指定范围内，则fname的值必须在函数内设置为零，如果未指定，默认值为１。groupgname定义属于group名字gname的颗粒球体，（粒子组通过group命令定义）idiliuid编号在il，iu范围内的物体，如果iu没有指定，则默认il=iujsetiliu节理编号在在il，iu范围内的接触，如果iu没有指定，则默认il=iu节理通过JSET命令创建namername创建一个命名范围rname，通过在命令行中范围元素集合定义，通常一个范围用range命令创建。planekeyword无限平面包围区域Keyword:Originx,y,z平面上一点（x,y,z）dddd方向角x-y平面上方向角，从y轴正向顺时针方向测量dipdip俯角从x-y平面，向z轴负向。Distanced在距离d范围内区域Above平面上方below平面下方normalnx,ny,nz平面的单位法向矢量radiusrlrutolerancet半径在（rl,ru）范围内球体。ru未指定，则范围（rl–t,ru–t）t未指定，则t=rl×10-6.注意误差范围只有在未指定ru条件下指定。rname命名范围rnamespherecenter(x,y,z)radius=r球形区域，中心（x,y,z）半径rx=xl,xu,tolerance=ty=yl,yu,tolerance=tz=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可在任何分析阶段给出，但是必须在调用时指定计算模式。CppudmC++用户定义接触模型CppuwcC++用户编写代码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控制表的创建和修改。成对数值索引命令集合建立球体Ballradius=rdensity=**hertzid=**x=**y=**z=**Radius:球体半径density:密度（m/V）Hertz:指定Hertz接触模型否则为线性接触模型Id:编号，唯一性，默认值大于当前编号最大值的整数x,y,z：x-,y-z-坐标，默认为0.0用ball建立球体允许重叠，球体属性由property命令设置Build控制执行AC/DC操作Buildassemblyxa,ya,za通过复制粒子和包含在pbrick中的力（必须已经存在）建立粒子集合，而初始粒子由Callfname调用data文件，默认为pfc3d.dat,而若扩展名未指定，等于*.dat文件调用允许多级调用，但不允许相互迭代。Change同property命令一样创建团块Clumpid=**full=on/offpermanentreleaseupdate_cycle=nrange…创建一个团块，若id等于已有团块编号，则为修改团块属性Range为要建立团块的范围，范围内原来属于其他团块的粒子转化到当前团块，必须指定range防止所有粒子变成一个团块Printclump输出团块信息，plotclump输出团块1、Fullon在计算团块运动时，使用全部惯性向量力矩，这会耗费更多计算时间。因为cross-couplingterm使用，必须进行矩阵转换。当旋转速度较大和陀螺效应影响大时使用。Fulloff计算团块运动时，只有使用diagonalinertia，而忽略旋转速度乘积。只有3D计算受full关键词影响。当fullon时，应设置update_cycle=12、Permanent通过阻止构成团块粒子间生成接触使得构成团块的粒子在随后模拟过程保持接触。这样不需要存储接触信息，并且避免经常性更新接触信息及其对应的矢量。当在模拟过程中团块有个形成和分解过程，就不能指定这个关键词然而，指定颗粒间永久接触有利于节约内存和计算时间。当一个团块有其permanent属性设置，其设置为永久的然而永久团块可被删除且团块内球体可被释放。Release关键词用于从团块提取球体，释放时，球体将不存在接触3、release范围内所有属于同一团块将被释放。如果id未指定，所有范围内属于任何团块的粒子将被释放。团块内所有粒子被释放时，团块本身不被删除，形成“空团块”,这种情况在计算过程中忽略，而在粒子重新加入之后，被激活。Deleteclump删除团块4、update_cyclen控制数据更新速度，默认值n=10。N越大，计算时间小。可被更新速度有Momentofinertiatensor惯性张量力矩（3D）contactcoordinates接触坐标contactnormalvector接触法向矢量接触剪力矢量(3D)平行粘结力矩（3D）平行粘结剪力矢量(3D)注意团块阻尼系数为构成团块粒子的阻尼系数的算术平均值当模型中存在团块时，在本版本中就不能使用密度放大系数Fix/free命令用于团块内任何粒子时，此工况用于团块重心。Printclump可看到各个团块的fix/free状况。Configkeywordkeyword指定可选计算模式Keyword:CppudmC++用户定义模型CppuwcC++用户编写代码Fluid流体分析Thermal热分析Continue重启由pause中止的data文件过程Cycle/Stepn执行n时间步对于热分析和热力分析计算为热计算步，否则为力学步计算过程输入Esc，中止计算过程。计算过程输入空格，继续中止当前cycle.当前计算过程由setmechanical和setthermal命令指定。对于热力耦合计算，每个热步后跟着最多n个力学计算子步，命令：thermalsetsubstep.Dampkeywordkeyword…range…设置局部阻尼的阻尼属性（球体）和粘性阻尼的阻尼属性（接触）局部阻尼：施加一个与各个球体自身不平衡力成比例大小的阻尼力。粘性阻尼：在每个接触施加一个法向和切向粘壶，这些粘壶与原有接触作用平行，从而生成与两接触实体（球－球或球－墙）相对速度差值成比力的粘滞阻力两种阻尼可以同时存在当未指定阻尼系数时，局部阻尼发生作用，且所有球体阻尼系数为0.7，而粘性阻尼不发生作用。Keyword1、default指定新建球的阻尼参数defaultlocal=v局部阻尼系数默认v=0.7viscousnormal=v法向临界阻尼系数，默认v=0notensionoff/on无拉力标记on–无张力默认offshear=v切向临界阻尼系数，默认v=02、local=v范围内所有球体局部阻尼系数为v3、viscousnormal=vnotension=off/onshear=v设定范围内粘滞阻尼系数DefinefunctionnameEnd定义fish函数所有define…end内的语句被编译并以压缩方式存储以便后续过程使用Deletekeyword…删除球体，团块，时程关系，墙体Keyword:1、Ballsid删除id号球体range…删除质心在范围内球体未指定任何参数时，模型内所有球体被删除2、Clumpid删除编号为id的团块，不能使用range参数。团块内所有球体被释放，但是球体未删除。4、histories删除所有时程记录。时程记录内容可以用historyreset命令删除。Wallsid1id2删除编号为id1,id2的墙体Id1toid2id3toid4删除编号在id1,id2间墙体Fixkeyword…range固定球体指定自由度速度值。即赋值速度常量。速度属性可以由property命令设置。关键词：x，y,zx-,y-,z-方向速度xspin,yspin,zspin:x-,y-,z-方向角速度例子，固定球体（id=5）x-方向速度Fixxrangeid=5Propertyxvel=1.5rangeid=5Fluidkeywordkeyword…流体分析关键词:1、modelxl,xh,yl,yh,zl,zhsiz