一套可行的PFC力控制实现循环加载方法(附同步视频教程+离散元学习资料)
- 作者优秀
- 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
- 内容充实
- 工程案例
- 平台推荐
- 主编推荐/内容稀缺/全网独家
导读:
我认识Lobby老师(仿真秀优秀讲师)是通过他连续在仿真秀平台发布原创文章,并且深受用户喜爱(持续有用户订阅),这让我很惊奇!于是我翻看了他发布的所有文章,毫不犹豫发起对他约稿。
后来进一步获悉:他毕业于上海同济大学,一直在使用离散元技术搞科研和解决复杂工程问题,尽管工作十分饱和,工作之余在仿真秀平台分享他这些年积累的离散元经验和学习心得,希望能对理工科学子和研发工程师成长有些许帮助,如有不当,欢迎批评指正,共同进步!一、基本原理
PFC中的计算元素包括ball、clump、wall,其中ball和clump是参与力学计算的,也就是说可以由接触力来计算其速度位移。但wall是不可以进行力学计算的,只能使用位移进行控制,我们当然可以通过伺服原理来实现位移和力的相互转化,但是使用wall直接进行力加载始终有些不顺畅,原因在于力的传递需要时间,而伺服原理根本而言给的是试样内部的力,而不是边界的力,这个和现实的加载又有些出入。并且利用伺服实现力的加载,对伺服系数敏感度较高,有可能会出现力的波动情况。clump是PFC内置的可以模拟不规则形状的单元,其力学原理和ball一样,只是形状不是圆形而已。这里便采用矩形的clump来模拟加载板,这样和实际也可以对的上。当然模拟也可以高出现实,在于我们可以设置加载板间是可以没有力的,这样便可以实现加载板的互相穿透,摆脱了现实真三轴试验的仪器限制。最终的加载板如图所示,加载板中的pebble数量一共800个,相对于模拟的上万颗粒数而言不算很多了,所以对计算效率的影响不会很大。![](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1606363999118_ece4zt.jpg?imageView2/0)
二、clump实现预压和围压
Clump直接可以施加外力,所以根据尺寸和应力计算出力就可以直接施加在加载板上面了。但是有时候如果应力变化比较大的话,对于散体材料肯定会出现很大的扰动,所以可以运用分级加载来实现应力缓慢的变化。这里采用二次函数实现,如果是指数函数的话可能还有点偏大。图2 分级加载曲线
这里可以根据当前应力wsxx00和目标应力mubiaoxx来计算系数a:a= (wsxx00-mubiaoxx)/ 10^2txx=mubiaoxx ((math.abs(n-10))^10)*a三、等偏应力幅的循环加载
首先我们一下加载路径,如图6所示,对试样的纵向施加余弦应力振动。则纵向应力的值为:tyy=txx fengzhi*0.5-fengzhi*0.5*math.cos(2*math.pi*pinlv*mech.age)
其中:tyy为竖向应力,txx为横向应力,fengzhi为偏应力峰值,pinlv为加载频率,即1s内有多少个余弦波,mech.age为当前时间。我们监测加载过程中tyy与wsyy的变化,如图7所示,蓝色线为tyy,红色线为wsyy。其中tyy是我们预设的应力,也就是施加给加载板的力,wsyy为加载板和试样的接触力。这里给的围压为100kPa,偏应力峰值为50kPa,则按照假设,纵向应力应该在100-150kPa间作余弦波动。从图中可以看出,给加载板施加的力与加载板监测到的力之间,有一定的滞后,在数值上也不是完全一致,材料的动力特性以及加载频率影响比较大,图8和图9为频率为2以及50的曲线,可以看出加载频率越小,滞后性越小。 图8 频率为2的应力曲线
图9 频率为50的应力曲线
图10为频率为10 的加载频率的横向应力与竖向应力的变化,可以看到由于加载的滞后性以及加载速度的问题,横向应力wsxx并不是一直维持在100kPa,而是有一定的波动,这个波动随着频率的减小也会不怎么明显。这也是一样的道理,wsxx为试样内部力,而txx为仪器给的力,所以txx可以维持不变,而wsxx由于惯性肯定会有所变化。图10 横向与纵向的应力曲线
试样循环加载应力应变曲线如图11所示,注意这里的应力用的是wsyy,因为wsyy有一定的滞后性,所以采用tyy的话效果会好很多,tyy对应的也是现实中仪器给的应力值,而wsyy对应的是试样内部的应力值。由图可以看出,频率越小,这种应力的滞后性越不明显,频率比较大的时候,不光应力会出现滞后性,应力的值也会有所区别。在做不排水的循环试验时,我们更关注其累计变形。图12为频率为10的轴向应变随时间的关系,即粉色的曲线,首先曲线增大,随后在某一个区间内往复运动。循环加载中的变形可以分为两个部分,一部分为不可恢复的塑性变形,即蓝色的区域,一部分为可恢复的弹性变形,即粉色曲线的区域。从图中可以看出累计塑性变形先增大,最后稳定在一个数值。总的变形先是大部分转化为塑性变形,之后都是转化为弹性变形。![](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1606363998105_9sibrg.jpg?imageView2/0)
四、参数分析
对试样进行2Hz、10Hz、50Hz的加载,整理后可得下图的曲线对比,图13为tyy曲线的对比,可以看到频率越小,第一次加载产生的应变越大,其最终的应变也是越大的。10Hz表现出的滞回圈最大,这是因为频率比较小的时候,试样的变形比较彻底,而频率比较大时,试样来不及大的变形。图14为wsyy曲线,这里也可以看出为什么大频率下产生的应变会小,这是因为试样内部的应力没有达到仪器所给的应力,而低频率下应力基本差不多,所以最后的变形值也不会有太大的区别。相比于wsyy曲线,tyy曲线更加圆润一点,但是wsyy曲线更加能够体现试样内部的应力应变情况。图14 wsyy曲线
图15为不同加载频率的轴向变形曲线,可以看出2Hz下的曲线能够更快的达到一个平稳阶段,这个和应力曲线也是相对应的。虽然2Hz下的累计变形最大,但是累计塑性变形却最小,而50Hz的加载频率可以发生更大的累计塑性变形。在做单元试验的时候,我们通常会采用damp去耗散系统中的动能,使其能够更快的平衡。这个对于准静态条件来说是完全可以成立的,但是对于动力试验我们需要仔细去考虑damp的取值。一般来讲damp模拟的是介质对颗粒的粘滞效应,也就是说damp的值决定了你的材料是放在空气中,还是放在水里,或者是放在黏液中。这里采用三种不同的damp值对其动力特性进行研究,即取damp为0.2、0.5、0.7。图16和17分别为试样的tyy应力曲线以及wsyy应力曲线。可以看出,阻尼系数越小,其产生的累计变形越大。这个也是可以理解的,颗粒变形的速度会受到阻尼的影响,阻尼越大,颗粒速度减小的越快,于是其产生的变形越小。![](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1606363996671_hd7u3u.jpg?imageView2/0)
图16 tyy应力曲线
图17为不同阻尼系数下的累计变形曲线,可以看出不同的阻尼系数其累计塑性变形不会有很大的区别,但是可恢复的弹性变形随着阻尼系数的减小而增大。CSR为偏应力与正应力的比值,这里做了三组试验,分别为CSR=0.1、0.2、0.5。图18为三种CSR的应力应变曲线,可以看出CSR越大,其产生的应变越大。图19也能看出一样的效果。![](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1606363996872_sky340.jpg?imageView2/0)
五、PFC模拟仿真实践与技巧
本文使用clump模拟加载板,直接对加载板施加力来实现对试样的加载。提出了一套可行的模拟方法,并进行等偏应力的循环加载方法。对模拟进行参数分析,分析了加载频率、阻尼系数和CSR对模拟的敏感性,也证明了模型的正确性。为此笔者在仿真秀平台同步发布了视频教程,感兴趣的朋友可以关注我的专栏。点击下图即可查看视频教程学习福利:曾几何时,我在仿真秀平台发布了一篇原创文章《我的PFC岩土颗粒流离散元分析攻略(附赠学习资料)》,给离散元技术学习者赠送了一份资料,如下图:请在本文附件下载即可(云盘地址失效请在文章下方留言即可)。![](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1606363996819_zmy9gx.jpg?imageView2/0)
声明:原创文章,本文首发仿真秀App,部分图片源自网络,如有不当请联系我们,欢迎分享,禁止私自转载,转载请联系我们。获赞 9713粉丝 21108文章 3396课程 215