当前位置:首页 >> 机械/仪表 >>

PFC3D中文讲解第三天


先对 plot 补充点说明: 1、plot create ***命令可用于创建多个图像显示窗口,但每次只能看一种结果,可以通 过菜单栏里的 window 选项卡来切换窗口,查看想看的结果 2、如只想显示小球,则使用命令 plot ball,如只想显示墙面,则使用命令 plot wall,如 二者都想显示则 plot ball wall, 如运行完毕后想显示别的结果, 如接触力则

plot cforce contact ( force 的简写) 。总而言之,plot 后面加你想显示的东东即可,用户需灵活使用。 此外,程序中的注释用英文的分号; ,后面加注释的内容即可。

第二天的讲解中小球都是自由状态,无外力的作用,因此有必要讨论施加外力的情况。 (1) 对小球施加重力作用

;fname: footing.DAT (tutorial example for PFC3D) new ; clear program state to begin new problem set random ; reset random-number generator title 'Tutorial Example' wall id=1 face ( 0, 0, 0) (10, 0, 0) (10, 0 -5) ( 0, 0,-5) 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) wall id=5 face ( 0, 0,-5) (10, 0,-5) (10, 2,-5) ( 0, 2,-5) wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0) gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0
; plot create Contact_Force plot set title text 'Model at equilibrium' plot add ball yellow plot add wall white plot add cforce black plot add axes brown plot show ;

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 ; prop density 2000 kn 1e8 ks 1e8 ini rad mul 1.47 ; set hist_rep=5 hist ball zvel 3 1 0 hist diagnostic muf set dt dscale

;
set grav 0 0 -9.81 prop fric 1.0 ; cyc 1000

此程序与第二天的程序几乎相同,唯一的不同是多加了两句程序: set grav 0 0 -9.81 prop fric 1.0 第一句是设定重力加速度为 9.81m/s2,方向为 Z 轴负方向,grav 是重力 gravity 的简写; 第二句是设定球与球之间的摩擦系数,fric 是 friction 的简写,在该程序中设定摩擦系 数为 1.0。运行结果如下图所示:

(2) 将上顶面换做一小块下降的板,压缩箱子里的小球 程序如下: ;fname: footing.DAT (tutorial example for PFC3D) new ; clear program state to begin new problem set random ; reset random-number generator title 'Tutorial Example' wall id=1 face ( 0, 0, 0) (10, 0, 0) (10, 0 -5) ( 0, 0,-5) 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) wall id=5 face ( 0, 0,-5) (10, 0,-5) (10, 2,-5) ( 0, 2,-5) wall id=6 face ( 0, 0, 0) ( 0, 2, 0) (10, 2, 0) (10, 0, 0) gen id=1,750 rad 0.15,0.20 x=0,10 y=0,2 z=-5,0 ; plot create Contact_Force plot set title text 'Model at equilibrium' plot add ball yellow plot add wall blue plot show ;

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 ; prop density 2000 kn 1e8 ks 1e8 ini rad mul 1.47 ; set hist_rep=5 hist ball zvel 3 1 0 hist diagnostic muf set dt dscale ; set grav 0 0 -9.81 prop fric 1.0 ; delete wall 6 wall id 10 face (0,0,0) (0,2,0) (2,2,0) (2,0,0) wall id 10 ks 1e8 kn 1e8 fric 1.0 wall id 10 zvel=-0.5e-3 property n_bond=5e5 s_bond=5e5 hist wall zforce id 10 cyc 6000

运行结果如图所示

与只施加重力的情况相比,这个程序多加了如下几句代码: delete wall 6

wall id 10 face (0,0,0) (0,2,0) (2,2,0) (2,0,0) wall id 10 ks 1e8 kn 1e8 fric 1.0 wall id 10 zvel=-0.5e-3 property n_bond=5e5 s_bond=5e5 hist wall zforce id 10 其中,delete wall 6 是删除编号为 6 的面,即箱子的顶面。 wall id 10 face (0,0,0) (0,2,0) (2,2,0) (2,0,0) wall id 10 ks 1e8 kn 1e8 fric 1.0 wall id 10 zvel=-0.5e-3 是建立一个 ID 为 10 的墙面,设定面的法向刚度、切向刚度、摩擦系数和 Z 方向 的速度,这个速度的单位是 m/step,因为开启了 density scaling(程序 set dt dscale 的作用) ,如果是完全的动力学分析的话,单位将会是 m/s(这个地方可解释为: 静力学只考虑系统达到平衡稳态的情况, 故单位与时间无关; 动力学关注的是每一 时刻的状态,故单位与时间有关) property n_bond=5e5 s_bond=5e5 是为所有的接触面定义法向和切向强度, 以此来模拟 初始的接触力(应该是类似于预紧力的作用,值得商榷) ,单位为 N hist wall zforce id 10 用于监测下降的板的 Z 方向的力的变化情况, 计算完毕后如想 查看,只需输入 plot hist 3 即可。结果如图所示

这个力的曲线图与英文帮助手册上的图明显不同, 但不知道是什么原因出的错, 有待 思考。下面附上英文帮助手册上的图。


相关文章:
PFC3D中文讲解第三天
PFC3D中文讲解第三天_机械/仪表_工程科技_专业资料。先对 plot 补充点说明: 1、plot create ***命令可用于创建多个图像显示窗口,但每次只能看一种结果,可以通 ...
PFC3D中文讲解第二天
;fname: footing.DAT (tutorial example for PFC3D) new ; clear program ...PFC3D中文讲解第一天 PFC3D中文讲解第三天 PFC3D中文讲解第四天 PFC3D中文讲解...
PFC3D中文讲解
(5)通过例子学习几种常见命令 例二:PFC3D 目录下的 Guide\Start\footing.dat...PFC3D中文讲解第三天 暂无评价 4页 1下载券 PFC3D中文讲解全五天 15页 免费...
PFC3D中文讲解第五天
PFC3D中文讲解第五天_机械/仪表_工程科技_专业资料。所学所想(...PFC3D中文讲解第一天 PFC3D中文讲解第二天 PFC3D中文讲解第三天 PFC3D中文讲解...
PFC3d手册命令中文详解
PFC3d手册命令中文详解_冶金/矿山/地质_工程科技_专业资料。PFC3D手册中的所有命令(不包含fish命令部分)的中文翻译版本PFC3D 学习记录 5- COMMAND REFERENCE 通用规...
PFC3d手册命令中文详解
PFC3d手册命令中文详解_计算机软件及应用_IT/计算机_专业资料。PFC3D 学习记录 ...(1)球体属性 1)color 颜色索引(非零整数,0 对应指定球体颜色列表第一项) ...
每日一词:中考英语易错词汇讲解(第三天)
每日一词:中考英语易错词汇讲解(第三天)_英语_初中教育_教育专区。每日一词:中考英语易错词汇讲解(第三天) 2014-10-22 学而思上海中考指导 principal adj. 主要的...
PFC3D_手册中文_fishtank
PFC3D_手册中文_fishtank_数学_自然科学_专业资料。pfc翻译资料AUGMENTED FISHTANK 执行固体模型,需要特别支持算法扩展设置,固体以粒子集合代表, 粒子接触中设置粘结。...
PFC3D 程序解析
PFC3D 程序解析 PFC3D 是基于命令驱动模式(COMMAND-DRIVEN FORMAT)的软件,各种命令控制着 程序的运行,这部分内容将介绍 PFC3D 软件的内嵌命令。 本说明文件译自 ...
_PFC3D命令集
例如,命令 CYCL 10000 意为执行 10000 步循环计算,如果用户在第 5000 步 的时候按下了<Esc>键, PFC3D 将在此时停止计算并返回命令交互模式; 则如 果在 5000...
更多相关标签:
pfc3d中文手册 | pfc3d5.0中文手册 | pfc电路讲解 | pfc电路讲解视频 | pfc2d中文手册 | pfc5.0用户手册中文版 | pfc2d中文教程 | pfc中文手册 |