离散元软件入门指导.pptVIP

JSETid1dip50,2number5spacing50,3上述命令是生成5条间距为单位50的节理面，第2至第5条节理面生成的过程中，倾角允许偏差为2度，间距值允许偏差为单位3。赋予节理面属性时，同样采用PROPERTY命令，如：propertyn_bond=0s_bond=0rangejset1或：propertyfric0.25rangejset1,2第126页，共194页，星期日，2025年，2月5日9.加载PFC2D荷载中的类型一般分为两类：主动荷载和被动荷载。被动荷载如重力，采用SETGRAVITY设置。主动荷载则包括施加颗粒体速度和外力。颗粒体可以直接施加速度和外力，但外力却不能直接施加在墙体上，因为运动方程对于墙体是无效的，但可以通过其他等效的方法来模拟墙体受集中力或应力。第127页，共194页，星期日，2025年，2月5日加载的操作主要从这四方面来执行：1).墙体的控制；2).颗粒体外力的控制；3).颗粒体速度控制；4).混合受力的控制。上述几方面，其具体编程方法与前面所介绍的边界条件与初始条件类似，在这里就不详细介绍，具体的现实问题或工程所需的特殊的加载方式，需采用FISH来进行编写和操作，同时需从简单的模型来验证其实际控制效果。第128页，共194页，星期日，2025年，2月5日10.求解PFC2D中所谓的求解即为执行许多个计算循环步(使用CYCLE和STEP)，同时观察计算后的模型响应。对于静力问题，模型系统要么计算达到平衡，要么破坏坍塌。对于动力问题，就显得较复杂些，特别需要注意如何消除个别已脱离模型颗粒体。虽然程序会对这些已脱离模型的颗粒体持续计算，但他们的变化对整个模型的影响会越来越小，因为只有当颗粒体处于计算的约束区域内才会产生其应有的效应。一种简易的消除逃逸颗粒体的方法是每隔几千步计算就在一定范围内执行DELETEBALL命令。第129页，共194页，星期日，2025年，2月5日1).静力求解对于求解命令CYCLE或STEP，用户也可以用SOLVE命令来进行静力求解。该命令是给定某个不平衡力准则，当达到这个标准时，SOLVE命令即停止。同时可配合一些关键词共同执行。如Average，允许指定平均不平衡力与平均接触力的比值；maximum则允许指定最大不平衡力与最大接触力的比值，默认的值为0.01(1%)。solveaverage=0.005maximum=1e-10第130页，共194页，星期日，2025年，2月5日SOLVEratio命令可以设置模型计算不平衡力的比率大小，默认值为1×10-5，ratio的值可通过SETsolve_ratio命令来设置，一般该命令在热力学求解模式下使用较多。SolveClockt(以分钟计)Solvetimet(以秒计)……对于静力计算的时间步，是与颗粒体的大小和接触刚度有关的，如果模型中颗粒体的大小和接触刚度变化幅度较大，时间步将不统一，为了能较快的收敛平衡，可以设置setdtdscale，该命令使颗粒体所有的自由度均采用相等的时间相应，这样有效的加快了模型的收敛。第131页，共194页，星期日，2025年，2月5日2).动力求解动力求解和静力求解的差异主要是动力求解存在一个现实的阻尼值。但一般动力求解之前总是要进行静力求解，其目的是为了获得初始应力平衡分配的计算模型。在此基础上来研究其进一步的动力上的力学响应。对于阻尼值的设定，我们需采用实际值，它的获取通过试验室测定或现场试验结果，程序默认的值为0.7。采用PROPERTY命令执行。第132页，共194页，星期日，2025年，2月5日defstartupold_time=timeomega=freq*2.0*pienddefshakewhile_steppingreal_time=time-old_timeifreal_timetlengthvel=ampl*sin(omega*real_time)elsevel=0.0end_ifbp=ball_headloopwhilebp#nullifb_xfix(bp)=1