PFC3d手册命令中文详解.doc
《PFC3d手册命令中文详解.doc》由会员分享,可在线阅读,更多相关《PFC3d手册命令中文详解.doc(54页珍藏版)》请在咨信网上搜索。
1、PFC3D 学习记录5 COMMAND REFERENCE通用规定和特征命令语法默认情况下,命令行忽略大小写,然而可以通过SET case命令激活区别大小写情况。斜体字代表数字,以i, j, m n 为开头的字母代表整数,其他情况为实数命令词,关键字和数值可用空格,以及(),分隔。在分号(;) 后为注释一行最多可包含80个字符。而&表示续行Range 语法定义一定限制范围的物体可运用Range的命令:Change/ initial/ propertyClumpFix/ FreeGroupJsetModelPlot ( 运用于不同输出项目)PrintRange( 创建一个命名的范围)一个Rang
2、e 定义一系列物体。由许多范围元素构成。如果使用多个范围元素,最终物体几何为不同元素集合的交集。但是可通过关键词any改变选择特性。如果使用关键词any选项,任何包含在范围内的单元可被认为最终物体集合的一部分。即求合集除此之外,通过在定义范围单元后接关键词not确定所有不在指定范围的单元通常还可以用定义group 方式定义复杂的range。内置范围单元:Annulus center ( x, y, z) radius = r1, r2球形空心区域中心在(x, y, z ) 内径r1, 外径r2Color il 球体或墙体在间隔(il, iu )内颜色索引。若无iu 则iu = il球体颜色索引
3、列表通过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
4、( 1 ), fc_arg( 2 ), fc_arg( 3 )分别传入函数。物体类型由FISH内置的pointer_type指定。如果物体不在指定范围内,则fname的值必须在函数内设置为零,如果未指定,默认值为。group gname定义属于group名字gname的颗粒球体,(粒子组通过group命令定义)id il id编号在il,iu范围内的物体,如果iu没有指定,则默认il = iujset il 节理编号在在il,iu范围内的接触,如果iu没有指定,则默认il = iu节理通过JSET命令创建name rname创建一个命名范围rname,通过在命令行中范围元素集合定义,通常一个范
5、围用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 半径在(rl, ru )范围内球体。ru未指定,则范围(rl t, ru t )t 未指定,则t = rl 10-6. 注意误差范围只有在未指定ru 条件下指定。rname 命名范围rnamesp
6、here center( x, y, z ) radius = r 球形区域,中心(x, y, z ) 半径rx = xl, , y = yl, , z = zl, , 指定x, y, z 坐标在指定范围内的物体。其他同radius 条重点理解fish定义range部分。在线帮助Command ? 适用command的各个关键词?/ help 显示所有可用命令命令分类程序控制命令这些命令支持批处理执行模式,允许改变程序状态。其与用程序建立实际数值模型无关。Cycle & Step 执行指定时间步Solve执行计算,直至满足收敛条件Build 控制AC/DC 操作 call 调用data文件,执
7、行命令,批处理模式 return 从批处理模式中返回到控制台。或在多级调用中返回上一级调用文件continue- pause重新/ 中止调用data文件 Fluid设置流体分析条件 / MPI 控制MPI 操作/ PARALLEL控制并行计算程序操作NEW 重启新模型Quit & Stop / Save Restore Set 改变全局程序设置/ Thermal 设置热分析/ Title 设置模型名称,这些将保存在随后的输出图并记录在save文件中/ Ucode注册用户编译的C+ 类特殊计算模式可选计算模式,通过config命令调用。Config 可在任何分析阶段给出,但是必须在调用时指定计算
8、模式。Cppudm C+用户定义接触模型Cppuwc C+用户编写代码fluid 流体分析thermal 热分析模型属性部分模型建立和修改Ball 建立新球体Clump 建立新团块,或修改新团块属性Delete 删除球体,团块,墙体和时程记录Generate 生成球体集合Jset 生成一系列节理Wall 生成墙体或修改墙体属性Fix 对球体设置固定速度Free 将球体固定速度清零Model 对指定接触指定用户定义接触模型Property 修改已有球体属性(材料属性和施加外力和速度),修改属性,粘结属性,接触属性模型监控History 指定执行过程变量(预定或用户定义fish变量Measure
9、建立一个测量球或改变测量球属性Movie 对屏幕输出进行控制抓图,从而实现动画输出Plot 控制模型几何和反应属性输出Print 生成输出变量和模型信息Trace 追踪模型中能量及功的输出命令工具Define 定义fish函数Group 定义和修改组集合Macro 命令宏Range 定义范围Table 控制表的创建和修改。成对数值索引命令集合建立球体Ball radius = r Radius :球体半径density : 密度(m/V)Hertz: 指定Hertz接触模型否则为线性接触模型Id: 编号,唯一性,默认值大于当前编号最大值的整数x, y, z :x-, y- z-坐标,默认为0.
10、0用ball 建立球体允许重叠,球体属性由property命令设置Build 控制执行AC/DC 操作Build assembly xa, ya, 通过复制粒子和包含在pbrick中的力(必须已经存在)建立粒子集合,而初始粒子由Call 调用data文件,默认为pfc3d.dat, 而若扩展名未指定,等于*.dat文件调用允许多级调用,但不允许相互迭代。Change 同property命令一样创建团块Clump full = on / off permanent release update_cycle = n创建一个团块,若id等于已有团块编号,则为修改团块属性Range 为要建立团块的范围
11、,范围内原来属于其他团块的粒子转化到当前团块,必须指定range防止所有粒子变成一个团块Print clump 输出团块信息,plot clump 输出团块1、Full on 在计算团块运动时,使用全部惯性向量力矩,这会耗费更多计算时间。因为cross-coupling term使用,必须进行矩阵转换。当旋转速度较大和陀螺效应影响大时使用。Full off 计算团块运动时,只有使用diagonal inertia,而忽略旋转速度乘积。只有3D计算受full关键词影响。当full on 时,应设置update_cycle = 12、Permanent 通过阻止构成团块粒子间生成接触使得构成团块的
12、粒子在随后模拟过程保持接触。这样不需要存储接触信息,并且避免经常性更新接触信息及其对应的矢量。当在模拟过程中团块有个形成和分解过程,就不能指定这个关键词然而,指定颗粒间永久接触有利于节约内存和计算时间。当一个团块有其permanent属性设置,其设置为永久的然而永久团块可被删除且团块内球体可被释放。Release关键词用于从团块提取球体,释放时,球体将不存在接触3、release 范围内所有属于同一团块将被释放。如果id未指定,所有范围内属于任何团块的粒子将被释放。团块内所有粒子被释放时,团块本身不被删除,形成“空团块”, 这种情况在计算过程中忽略,而在粒子重新加入之后,被激活。Delete
13、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状况。Con
14、fig keyword 指定可选计算模式Keyword:Cppudm C+用户定义模型Cppuwc C+用户编写代码Fluid 流体分析Thermal 热分析Continue 重启由pause中止的data文件过程Cycle/Step n执行n 时间步对于热分析和热力分析计算为热计算步,否则为力学步计算过程输入Esc,中止计算过程。计算过程输入空格,继续中止当前cycle.当前计算过程由set mechanical 和set thermal命令指定。对于热力耦合计算,每个热步后跟着最多n个力学计算子步,命令:thermal set substep.Damp keyword 设置局部阻尼的阻尼属
15、性(球体)和粘性阻尼的阻尼属性(接触) 局部阻尼:施加一个与各个球体自身不平衡力成比例大小的阻尼力。粘性阻尼:在每个接触施加一个法向和切向粘壶,这些粘壶与原有接触作用平行,从而生成与两接触实体(球球或球墙)相对速度差值成比力的粘滞阻力两种阻尼可以同时存在当未指定阻尼系数时,局部阻尼发生作用,且所有球体阻尼系数为0.7,而粘性阻尼不发生作用。Keyword1、 default 指定新建球的阻尼参数default local = v 局部阻尼系数默认v= 0.7 viscous normal = v 法向临界阻尼系数,默认v = 0 notension off/ on 无拉力标记on 无张力默认o
16、ff shear = v 切向临界阻尼系数,默认v = 02、 local = v 范围内所有球体局部阻尼系数为v3、 viscous normal = v notension = off/on shear = v设定范围内粘滞阻尼系数Define functionnameEnd定义fish函数所有define end 内的语句被编译并以压缩方式存储以便后续过程使用Delete keyword删除球体,团块,时程关系,墙体Keyword:1、Balls 删除id号球体 删除质心在范围内球体未指定任何参数时,模型内所有球体被删除2、Clump id删除编号为id的团块,不能使用range参数。团
17、块内所有球体被释放,但是球体未删除。4、 histories 删除所有时程记录。时程记录内容可以用history reset命令删除。Walls id1 删除编号为id1,id2的墙体Id1 to id2 删除编号在id1,id2间墙体Fix keyword 固定球体指定自由度速度值。即赋值速度常量。速度属性可以由property命令设置。关键词:x,y, z x-,y-,z-方向速度xspin, yspin, zspin : x-, y-, z- 方向角速度例子, 固定球体(id = 5 )x-方向速度Fix x range id = 5Property xvel = 1.5 range i
18、d = 5Fluid keyword 流体分析关键词:1、 model xl, xh, yl, yh, zl, zh size nx, ny, nz设置x-, y-, z-方向模型尺寸和单元数目。在xl, xh yl, yh zl, zh 范围内建立总数为nx* ny * nz数目流体单元,每个流体单元x-, y-,z-方向尺寸为(xh-xl)/ nx, (yh yl ) / ny, ( zh zl ) / nz,注意模型总单元数为(nx + 2)*( ny + 2) * (nz+2), 边界上附加一个单元为边界条件2、 property density d 流体密度d Viscosity 流
19、体粘滞系数v3、boundary nonslip keyword 在keyword指定位置设置无滑动边界Slip keyword 在keyword指定位置设置滑动边界 Pressure pv keyword 在keyword指定位置设置压力为pv Velocity xv, yv, zv keyword 在keyword指定位置指定apparent速度为xv, yv, zv注:apparent速度定义:对于进入网格的流体将flow rate除以单元长度(或单元面积)边界位置指定:xl : x 方向下边界xu: x 方向上边界yl : y 方向下边界yu: y 方向上边界zl : z 方向下边界z
20、u: z 方向上边界at xp,yp,zp 在点(xp, yp, zp )4、set keyword 控制流体计算循环参数buoterm off( 默认 )/ on on: 打开浮力作为外部体力(流体计划中 grterm off(默认)/on off: 在Navier-Strokes方程中忽略重力on: 在流体计算中考虑重力影响 set gravity visterm off(默认)/on off: 忽略粘滞影响下面参数控制SIMPLE项目,注意如果迭代不收敛,pfc3D发出一个警告消息而不是结束计算 Crit cr 指定SIMPLE 计划收敛标准,默认10-6 Intval intv 如果在
21、tdel auto 时间步模式下指定。则流体计算将在每intv计算步调用,而流体计算步有模型状况指定Iteration it 指定SIMPLE scheme迭代限制,默认it = 5000 P_relax pr 在压力调整计算中指定松弛因子默认pr = 0.8Por_relax por 在孔隙率中指定松弛因子por = 1.0 Tdel auto 流体分析计算中时间步为自动计算 Tdel tdel 指定流体时间步tdel V_relax vr 指定速度计算的松弛因子, 默认vr= 0.5Free keyword 消除指定范围内指定自由度的固定速度约束关键词:x/y/z/ xspin/yspin
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- PFC3d 手册 命令 中文 详解
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【精***】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【精***】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。