离散元软件PFC2D学习问题总结
关于颗粒流软件PFC的离散元数值模拟参数标定
关于颗粒流软件PFC的离散元数值模拟参数标定摘要：简要介绍了运用离散元理论及颗粒流软件PFC进行岩土类工程问题数值模拟的研究思路和方法，简述了运用PFC软件进行数值模拟时通常会遇到的关键性问题，通过结合笔者在研究过程中的实践经验，针对数值模拟的首要步骤即参数标定这一步，给出了相应的解决方案及研究结论。
关键词：颗粒流软件PFC；离散元；数值模拟；参数标定1工程背景随着工程规模的日益扩大，作为建筑材料以及地基介质的岩土体一直以来都是人们研究的热点。
目前，解决岩土类工程问题的主要研究手段有三大类，分别是理论分析、室内实验和数值模拟。
其中，数值模拟又因其便捷、高效以及可控性成为当今最为常用的研究方法。
常用的数值模拟方法主要基于有限元分析理论或离散元分析理论，运用计算机软件对实际的工程问题进行模拟，从而发现问题、得到内在规律并作出研究结论，对于指导工程实践具有重大意义。
2研究方法针对不连续介质问题以及大变形问题[1]，近年来新兴的离散元理论逐步代替传统的有限元理论被人们采用。
其主要理论依据是牛顿第二定律和力与位移的关系，求解思路为：将求解空间离散为分离的单元阵，相邻两个单元之间用合理的元件相连接；基本变量是单元间的相对位移，通过位移和力的关系可以分别得到单元之间的法向、切向作用力；求单元上受到的合力及合力矩，由牛顿运动定律计算出单元体的加速度；最后，积分得到单元的速度和位移。
通过以上计算过程可以得到任意时刻所有单元体的位移、速度和加速度。
在基于离散元分析理论开发出来的计算机软件中，最为广泛使用的是美国ITASCA公司开发的UDEC（univeralditinctelementcode）、3DEC（3-dimenionalditinctelementcode）块体离散元程序，和PFC2D （particleflowcodein2dimenion）、PFC3D（particleflowcodein3dimenion）软件。
PFC二维模拟的缺点
PFC二维模拟的缺点PFC二维模拟的缺点在分析PFC2D模拟的结果时，需要考虑PFC2D的这些局限性，对结果做出合理的解释。
与三维模拟具有三个力的分量和三个弯矩分量不同，二维模拟只有两个力的分量和一个弯矩，其余几个分量并未考虑在平衡方程中。
默认情况下，生成的BALLS是圆心位于同一平面上的球体；当SET DISK ON命令打开后，BALLS是具有单位厚度的圆柱体，该命令影响到球的质量和内在特性的计算。
应力和应变PFC2D的应力和应变不是存在于颗粒组的任何地方，而是离散的，因此需要设置测量圆对其进行平均化计算。
平均化计算中，不考虑命令SET DISK，默认所有颗粒为单位厚度，这种算法在SET DISK OFF 时并不精确。
（算法见“PFC2D学习笔记之理论与背景”）PFC2D的运动方程中没有考虑平面外的力分量和应力应变，而是采用力-位移关系，为确定平面应变状态所必需的平面外约束在这里没有考虑。
填密和孔隙比孔隙比用来表示颗粒填密的程度，但是二维和三维孔隙比的计算没有明显关系，这使得PFC2D在模拟三维效应明显的材料时很难让人满意。
在三维下，均匀球体在矩形区域最紧密状态（无侵入）的孔隙比为0.2595；在二维下，则为0.0931，这说明三维可以容纳更多的空隙。
在真实材料中，是达不到这种最紧密状态的，因为在填密过程中颗粒之间会形成“锁合”，而三维情况下形成“锁合”的机率也比二维要多得多。
因此，三维条件下，孔隙比会更大。
三维条件下，小颗粒会从大颗粒的空隙中漏出，而在二维条件下则不可能发生这种渗漏。
因此，对于三维敏感的问题，有必要先使用PFC3D做少量的计算，得到某些物理响应的相关参数，然后再利用这些参数进行PFC2D的大量的二维模拟。
质量特性每个颗粒的质量特性取决于嗦采用的颗粒是球体还是圆盘，这会影响到运动计算和重力。
质量特性受到影响，但并不意味着其他模型特性，如应力和空隙比也必然受到影响。
pfc知识问答
与FLAC3D 的2.0和2.1版本相比，FLAC3D 3.0有哪些新功能？答：计算方面：1.所有的计算和数据均采用双精度浮点制。
2.其运行速度较2.1约快10-15%功能方面：1. 动力模块中增加了hysteretic阻尼对于动力计算，加入了一个新的阻尼特性：hysteretic阻尼。
采用这种形式的阻尼，数值计算模拟产生的应变可基于模型模量和阻尼函数的共同作用。
这就可以使用户将等价线性计算方法的结果与完全非线性无折减本构模型的计算结果进行比较。
另外，对一些其它形式的阻尼，如Rayleigh阻尼，进行了较大的折减，这样在使用hysteretic 阻尼时可以有效地节省计算时间。
2. 在岩石本构模型中增加了Hook-Brown（霍克-布朗）模型加入了霍克-布朗屈服准则，使用户对岩石材料计算结果较为合理。
它的屈服面是非线性的，同时是考虑最大、最小主应力的关系的基础上产生的。
该模型综合了塑性流准则，其特性随着不同的应力水平，呈一变化的函数。
3.热/流体水平对流模型FLAC3D之前已经能提供非线性固体与多孔介质渗流的流固耦合模拟，机械地耦合流体和固体。
而3.0版本的新功能增加了温度可受流体密度影响和流体中温度的水平对流。
4. 3Dshop生成的网格导入3DSHOP是一种能力强大的固体建摸和六面体网格的软件包，也是ITASA的产品。
FLAC3D 原始的网格仍可用，但是用基于WINDOWS操作系统的3DSHOP建复杂模型更为简单方便。
3DSHOP生成的网格能被FLAC直接读取。
5. 动画命令：FLAC3D现在能产生AVI和DCX动画，这是以前的版本没有的功能。
6. 记录颗粒轨迹：颗粒的路径能被记录和显示此外，FLAC3D 3.0也包含下面的新特点（这也是以前版本所不能实现的）：1）网络版2）应用于混凝土加工模拟的水合作用模型。
问：如何显示变形轮廓线的命令？答：plo ske magf 10 其中10为放大系数问：如何查看各个时段不平衡力的具体数值？答：采用his来记录计算，包括位移应力等命令his unbalhis gp(zone) zdis range (0 0 0) 或者id=？导出数据命令his write n vs m begin 时步 end 时步 file filename.hisn表示纪录的id m表示时步要导出不平衡力的具体数值his unbalstep 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值问：initial 与 apply 有何区别？答：initial初始化命令，如初始化计算体的应力状态等；apply边界条件限制命令，如施加边界的力、位移等约束等。
PFC数值软件基础教学总结
块体逻辑支持附属粒子组或块体的创建，促进了
程序的推广普及。块体内粒子可以任意程度的重叠， 作为刚性体具有可变形边界的每一个块体，可作为一 般形状的超级粒子。通过指定墙的速度、混合的粒子 速度、施加外力和重力来给系统加载。“扩展的FISH 库”提供了在集合体内设置指定应力场或施加应力边 界条件的函数。时步计算是自动的，包括因为Hertz接 触模型刚度变化的影响。模拟过程中，根据每个粒子 周围接触数目和瞬间刚度值，时步也在变化。基于估 计的粒子数，单元映射策略采用最佳的单元数目，自 动调整单元的外部尺寸来适应粒子缺失和指定的新对 象。单元映射方案支持接触探测算法以保证求解时间 随粒子数目线性增加，而不是二次方增加。
粒子生成器根据粒子的指定分布规律自动概 率地生成。粒子半径按均匀分布或按高斯分布规 律分布。
初始孔隙度一般比较高，但通过控制粒子半 径的扩大可以获得密度压实。在任何阶段任何因 素都可以改变半径。所以不需反复试验就可以获 得指定孔隙度的压实状态。
属性与各个粒子或接触有关，而不是与“类型号” 有关。
因此，可以指定属性和半径的连续变化梯度。 “节理生成器”用来修改沿指定轨迹线的接触特性。 假定这些线叠加在颗粒集合体上。用这种方法，模型 可以被成组的弱面，如岩石节理切割。
PFC2D (Particle Follow Code 2 Dimension)即二 维颗粒流程序，是通过离散单元方法来模拟圆形颗 粒介质的运动及其相互作用。最初，这种方法是研 究颗粒介质特性的一种工具，它采用数值方法将物 体分为有代表性的数百个颗粒单元，期望利用这种 局部的模拟结果来研究边值间题连续计算的本构模 型。以下两种因素促使PFC2D方法产生变革与发 展:(1)通过现场实验来得到颗粒介质本构模型相当 困难:(2)随着微机功能的逐步增强，用颗粒模型模 拟整个问题成为可能，一些本构特性可以在模型中 自动形成。因此，PFC2D便成为用来模拟固体力学和 颗粒流问题的一种有效手段。
关于颗粒流软件PFC的离散元数值模拟参数标定
关于颗粒流软件PFC的离散元数值模拟参数标定作者：孙婧何佩珊齐梦菊来源：《山东工业技术》2016年第10期摘要：简要介绍了运用离散元理论及颗粒流软件PFC进行岩土类工程问题数值模拟的研究思路和方法，简述了运用PFC软件进行数值模拟时通常会遇到的关键性问题，通过结合笔者在研究过程中的实践经验，针对数值模拟的首要步骤即参数标定这一步，给出了相应的解决方案及研究结论。
关键词：颗粒流软件PFC；离散元；数值模拟；参数标定DOI：10.16640/ki.37-1222/t.2016.10.0391 工程背景随着工程规模的日益扩大，作为建筑材料以及地基介质的岩土体一直以来都是人们研究的热点。
目前，解决岩土类工程问题的主要研究手段有三大类，分别是理论分析、室内实验和数值模拟。
其中，数值模拟又因其便捷、高效以及可控性成为当今最为常用的研究方法。
常用的数值模拟方法主要基于有限元分析理论或离散元分析理论，运用计算机软件对实际的工程问题进行模拟，从而发现问题、得到内在规律并作出研究结论，对于指导工程实践具有重大意义。
2 研究方法针对不连续介质问题以及大变形问题[1]，近年来新兴的离散元理论逐步代替传统的有限元理论被人们采用。
其主要理论依据是牛顿第二定律和力与位移的关系，求解思路为：将求解空间离散为分离的单元阵，相邻两个单元之间用合理的元件相连接；基本变量是单元间的相对位移，通过位移和力的关系可以分别得到单元之间的法向、切向作用力；求单元上受到的合力及合力矩，由牛顿运动定律计算出单元体的加速度；最后，积分得到单元的速度和位移。
通过以上计算过程可以得到任意时刻所有单元体的位移、速度和加速度。
在基于离散元分析理论开发出来的计算机软件中，最为广泛使用的是美国ITASCA公司开发的UDEC（universal distinct element code）、3DEC（3-dimensional distinct element code）块体离散元程序，和PFC2D（particle flow code in 2 dimensions）、PFC3D（particle flow code in 3 dimensions）软件。
PFC2D学习笔记之流固耦合
因此必须采用一些手段消除这些宏观摩擦。
流体公式
公式体现了可变形流体和颗粒间的全耦合。在使用公式时，需要选择局部流体变量。 在 PFC2D 程序中，压力和体积与颗粒间的孔隙相关，流动与接触相关；在 PFC3D 中，则是 压力、体积和流动都与颗粒间的孔隙相关。
流固耦合
流体与颗粒的相互作用方式
流体流动的细节模拟及其与非连续介质的耦合，需要的计算量很大，但是在很多情况下 只需要得到近似模拟即可，因此下面介绍各种复杂问题及其近似方法。 1、 均匀压力或静水压力
当颗粒只是受到重力左右，以休止角浸没在流体中，由于静水梯度的影响，应当使用颗 粒的浮容重。 2、流体中的分散颗粒 当颗粒在流体中的运动彼此独立时，即流体中的颗粒很分散，应当将粘滞力施加在颗粒 上，该粘滞力是颗粒和流体相对速度的函数。 3、 颗粒组中的分散流体 当流体体积占颗粒体积的很小一部分时，流体在颗粒的接触处受到表面张力作用。此时， 需要定义一个特殊的接触法则，以考虑流体的影响。该法则应包括粘滞分量和摩擦分量， 其中摩擦分量依赖于接触处颗粒的分离程度和流体体积。如果有两种流体相，则需要综 合考虑第 1 和 3 条. 4、 小压力梯度的饱和颗粒组 当压力在平均颗粒半径的距离内变化很小时，可以使用基于多个颗粒的孔隙比和渗透 性，当作连续流体来计算。拖曳力由压力梯度计算。由颗粒材料程序计算得到平均渗透性， 并将其赋给连续性材料计算程序，然后从连续性材料计算程序得到流体平均流速向量，再计 算得到颗粒体力。该弱耦合方式可以得到表面侵蚀、管道流和沟槽流的机理。但是，颗粒组 平均化参数计算时，所取的颗粒组范围需要自己决定。 5、 大压力梯度的饱和、粘聚性颗粒组 考虑孔隙的“虚”和“实”，可以将该问题分成两种情况讨论，两种情况下都具有大压 力梯度和粘聚性。 a) 流体在名义裂缝中的流动
PFC2D培训课件
学习中的难点及解决方法
遇到概念和技术难题
软件操作不熟练
当遇到PFC2D概念和技术难题时，可以通 过查阅相关资料、请教老师或与同学讨论等 方法解决。
当软件操作不熟练时，可以ห้องสมุดไป่ตู้过观看教程、 练习实例或参加培训课程等方式进行提高。
无法实现预期功能
学习缺乏兴趣和动力
在实现预期功能时，可能会遇到一些问题， 可以尝试调整设计思路、寻求设计辅助或者 咨询专业人士帮助。
易用性
图形界面
PFC2D提供直观的图形界面，用户可通过简单的拖拽操作轻 松建立模型，设置计算参数和运行模拟。
文档和教程
PFC2D提供详细的文档和教程，帮助用户快速上手并掌握使 用技巧。
可扩展性
可扩展性架构
PFC2D采用模块化、可扩展性架构，使得添加新功能或扩展现有功能变得简 单和灵活。
开放API
准备学习资料
03
收集和整理PFC2D学习所需的资料，如教材、教程、视频和在
线课程等。
学习资料与计划
制定学习计划
根据个人学习需求和学习时间，制定详细的学习计划，包括学习 内容和时间安排。
选择合适的学习资料
根据学习需求和计划，选择适合的学习资料，如教材、教程、视 频和在线课程等。
建立学习目标
在学习过程中，建立明确的学习目标，以便更好地把握学习方向 和提高学习效率。
PFC2D通过将颗粒视为离散的刚体，利用物理模型和算法来 模拟颗粒之间的相互作用和运动，可以应用于各种工业领域 。
PFC2D发展历程
1
PFC2D的前身是Particle Flow Code，由美国 国家实验室开发，最初用于核工业领域。
2
随着技术的发展和应用领域的扩展，PFC2D逐 渐被应用于采矿、冶金、化工、制药等工业领 域。
PFC2D学习笔记之颗粒生成
使用命令BALL和GENERATE可以生成颗粒。
BALL命令只能单个生成的颗粒，且可以相互重叠；GENERA TE命令可以生成颗粒组，但自动生成的一组颗粒间没有重叠。
生成的颗粒组有两种，一种是规则的排列，用于模拟结构部件，这种阵列的初始接触力可以精确的预测；另一种是不规则的排列，用于模拟实体或内部结构无序的粒状材料，要预估其初始接触力量值是不可能的。
对于第二种颗粒组，尽管颗粒排列是任意的，但某些特性还是需要在颗粒生成后添加，比如弱面或各向异性强度等。
规则排列通过用户自定义FISH函数，在循环体内使用BALL命令生成。
示例如下loop row (1,n_row)loop col (1,n_col)commandball id=idc x=xc y=yc rad=rcend_commandidc = idc + 1xc = xc + r2end_loopyc = yc + yincxc = x0 + radius * (row - (row/2) * 2)end_loop不规则排列一般，边界都是由墙定义好的，要在给定空间内生成颗粒，既要保证孔隙比符合要求，又要使组合达到平衡。
显然，孔隙比不能任意的小。
无法将圆颗粒在紧密堆积时将他们设置在任意紧密的区域（低孔隙比）。
这里讨论两种方法来以给定的半径和孔隙比生成颗粒，这两种方法都是先固定边界，即半径扩大法和颗粒排斥法。
还有一种方法是移动边界法，可以得到精确的孔隙比，但这种方法有三个缺点：边界条件改变了；达到平衡的时间很慢；最终分布很不均匀。
半径扩大法和颗粒排斥法中，颗粒半径按照某种概率分布生成，如均匀分布或高斯分布，也可以通过FISH自定义其分布形式。
BALL和GENERA TE的具体用法，参考《PFC2D学习笔记之命令篇》。
半径扩大法其基本FISH函数语句如下：mult = 1.6 ; initial radius multiplication factorn0 = 1.0 - (1.0 - poros) / multˆ2r0 = sqrt(height*width*(1.0 - n0)/(pi*num))rlo = 2.0 * r0 / (1.0 + rat)rhi = rat * rlo... ...gen id=1,num rad=rlo,rhi x=0,width y=0,heightprop dens=1000 ks=s_stiff kn=n_stiff... ...sum = 0.0bp = ball_headloop while bp # nullsum = sum + pi * b_rad(bp)ˆ2bp = b_next(bp)end_looppmeas = 1.0 - sum / (width * height)mult = sqrt((1.0 - poros) / (1.0 - pmeas))ini rad mul multcycle 1000prop fric 0.2cycle 2000注意，为了使颗粒组在比较快的时间内达到平衡，必须设置在CYCLE前使用FRIC设置摩擦系数，用于消散颗粒能量，但这种干预方式阻碍了颗粒的自由运动，对其最终组合状态有影响。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
macro ’wall_fric’ ’0.5’ wall id 4 wall_props node pt1 pt2 ;bottom
new
gen id 1 400 x 0 5 y 0 5 rad 0.05 0.05 ini rad mul 1.4 ；定义一组范围的球体，名为Tunnel group Tunnel range x 1.5 3.5 y 0 5 plot create the_view plot add ball yellow plot show
pause
；删除先前定义的组
delete ball range group Tunnel
；只产生一个球，半径是必需的，如果ID未指定，则自动赋一个比现在ID大的数，如果x与y未指 定，则在坐标原点，指定参数hertz后，则质点将服从Hertz-Mindlin contact而不是默认的Linear law ，指定该法则后，则弹性体必需给出其属性，通过PROPERTY给出。产生的该质点可以与其它质 点重叠，而不像GENERATE命令，当执行循环后，重叠的质点将产生排斥力，除非给定边界限制 。
PFC2D学习初期总结[geotechman]
欢迎您： HBT2008 隐身 | 我的主页 | 控制面板 | 搜索 | 风格 | 展区 | 插件 | 退出 geotechman → 数值模拟 → PFC2D → PFC2D学习初期总结
左栏 右栏 我的收件箱 (1)
主题：PFC2D学习初期总结
总数 19 1 2 下一页 1 << 上一主题 下一主题 >> 共有2492人关注过本帖 树形 打印 收藏 修复 推荐
PFC2D学习初期总结[geotechman]
add ball lblue add wall black add axes brown show ；建立视点Footing2 create Footing2
add ball lblue
add wall black add axes brown set size -1.0 1.0 -2.0 2.0 set title bottom
hist 1 hist 2
；将显示当前的几个视点，如在plot 模式下。 print list ；如果在命令模式下，则如下： plot print list ；在命令模式下，如下的命令将显示所有的历史信息 print hist ；提供重力，加速度为9.81m/s2 set grav 0,-9.81 ；提供球之间的磨擦力,磨擦系数为1.0 property fric 1.0 cycle 1000
PFC2D学习初期总结[geotechman]
macro pt3 ’1.0 1.0’ macro pt4 ’0.0 1.0’ macro wall_props ’kn wall_n_stiff ks wall_s_stiff fric wall_fric’ wall id 1 wall_props node pt2 pt3 ;side wall id 2 wall_props node pt3 pt4 ;top
内存分配 当载入PFC程序后，程序将动态分配内存，随着模型的构建，内存相应地增大。Print memory 将显 示当前内存大小。 内存分配 当载入PFC程序后，程序将动态分配内存，随着模型的构建，内存相应地增大。Print memory 将显 示当前内存大小。 图形控制 在图形模式下，按键M为放大图形，缩小为Shift+M；箭头来移动；Ctrl+G为在彩图与黑白图形间 转换；Ctrl+C显示照相机对话框；Ctrl+L为cutting plane；Ctrl+R为恢复图形显示为默认值；Ctrl+Z 为用一矩形框来放大图形局部；F9为重画。 实例 new ；将球体看成一定厚度的盘子；Note that only mass properties are affected by the SET disk command set disk on ；下面为建立四个直线边壁，分别赋予1~4的编号。每一直线壁以两点连线建立。 wall id=1 nodes ( 0, -5) (10, -5) wall id=2 nodes (10, -5) (10, 0) wall id=3 nodes (10, 0) ( 0, 0) wall id=4 nodes ( 0, 0) ( 0, -5) ；下面为产生500个球，每个球的编号被赋予1~500，产生球的地方在0<x<10,-5<y<-0.5的范围内，
;球为白色但有影子 plot add ball white shade on
一． 粒子创建 两种创建粒子的命令：BALL and GENERATE BALL 创建单一的粒子，并指定其位置，其不受已存在粒子的制约，也就是说其创建的粒子
可以与已存在的粒子叠加。 GENERATE 创建一系列的粒子，其不可重叠。 有两种粒子集：规则的和不规则的。
property n_bond=5e5 s_bond=5e5 ;记录ID号为10的wall在Y方向的力 hist wall yforce id 10
cycle 2000 plot current 0 plot hist 3 ；拷贝Contact_ Force 到Foot_Load下，并定义标题为：Loaded footing plot copy Contact_ Force Foot_Load both current 4 set title text ’Loaded footing’ show ；显示所有的球 print ball ；显示所有的墙 print wall ；显示模型的信息 print info
wjwhappy
加好友 发短信
等级：岩土一级 帖子：32 积分：0 威望：0 精华：1 注册：8/24/2005 7:33:42 AM
个性首页 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC
小 大 1楼
PFC2D学习初期总结
Post By：3/11/2006 9:39:55 PM
最近学了一段时间PFC2D，由于做课题要用到它，时间又很紧，只是粗粗的学了一下，现将我的学习笔记抄录如下，比 较乱一些。以期对比我更初级的学习者有所帮助。
;查看当前视图的详细信息
plot print view
;查看当前视图的具体模型
plot print item
;改变当前视图为另一视图
plot current Footing
;显示任何已存在的视图
plot show <viewname>
；以下命令将壁3删除，并建立壁10，因密度缩放比例是开着的，所以其速度单位为m/step，如 果是动态分析的话，则速度单位为length/time delete wall 3 wall id 10 nodes (2,0) (0,0) wall id 10 ks 1e8 kn 1e8 fric 1.0 wall id 10 yvel = -0.5e-3 ；在接触面上提供法向强度normal strength和切向强度shear strength，即内聚力cohesion，单位为： N.当超过法向或切向强度后，bond就会失稳破坏，球之间的摩擦力即发生效果。
ball rad v id i x v y v <hertz> ；建立球的最大与最小半径相同；如果半径不同的话，则产生的半径是最小的或是最大的，随机分布的。如果指定gaus s后，则按高斯分布，即按(R1+R2)/2与 (R2-R1)/2确定半径。 gen big rad=.6,.6 x=0,10 y=0,10 gen small rad=.3,.3 x=0,10 y=0,10 ；将大球定义Index 为0，小球为1 change c_index 0 range big change c_index 1 range small ；定义的wall的光面是向上的，故而在定义球时的no_shadow即这光面上方。 wall id 1 nodes -10 -10 10 10 gen no_shadow gauss x 0 6 y 0 6 rad .2 .22 id=1,50 del wall 1 ；定义的wall的光面是向下的，故而在定义球时的no_shadow即这光面上方。 wall id 1 nodes 10 10 -10 -10 gen no_shadow gauss x 0 6 y 0 6 rad .3 .32 id=51 75
；给球赋属性，密度为1000kg/m3,刚度都为108N/m
prop density 1000 kn 1e8 ks 1e8
；给球的半径乘以一个系数，从而减小球与球间的空隙
prop rad mul 1.51
；每5步监测一下变量的变化，从而将变化的数据存入一数组，默认为每10步；两个变量被监测， 一个是圆点在（3，0）附近的球在Y方向的速度，一个是平均不不衡力mean unbalanced force(muf) 的诊断
set hist_rep = 5
history ball yvel (3, 0) history diagnostic muf ；因为这是一个静态分析，用密度缩放参数(density-scaling option)使问题最优化 set dt dscale cycle 4000 ；视图Footing中，图像每20步进行一次更新，更新间隔可以通过SET pinterval 命令改变 ；以下4条命令进入到视点0中，显示刚才设定的历史1的图形，在该例中是圆点在（3，0）附近的 球在Y方向的速度和hist 2的图形 plot current 0 ；当前视点
； 定义宏 macro wall_s_stiff ’1e6’ macro wall_n_stiff ’1e6’ macro wall_fric ’0.0’ macro pt1 ’0.0 0.0’ macro pt2 ’1.0 0.0’
/old/dispbbs.asp?boardid=49&Id=3606[2009-4-26 23:31:23]