pb模型

模型简介

linear contact bond model ：模拟胶结物和被胶结物分界不是很清楚

linear parallel bond model ：岩石（线性弹簧+）

控制变形:pb_kn、pb_ks

控制强度：pb_ten、pb_coh、pb_fa

胶结模型可以承受一定拉力。

图3：平行结合力和力矩的力-位移规律：（a）法向力与平行结合面间隙；（b）剪切力与相对剪切位移；（c）扭曲力矩与相对扭曲旋转；（d）弯曲力矩与相对弯曲旋转。

法向验证

[ct=contact.find("ball-ball",1)] 查找接触

pb模型在纯压状态永远不会破坏。

model.clean命令是一个全局命令，可以在模型创建过程中的任何时候给出，以创建模型中的所有接触，初始化所有件属性，确保数据结构存在以进行空间搜索，并根据几何形状更新接触活动状态。

pb_stage 在右侧plot窗口中选择查看破坏状态。

new ;纯压，8s应力值为8Mpa
;纯拉，2s应力值为1Mpa
domain extent -100 100 ball create position 0 -10 radius 10
ball create position 0 10 radius 10ball fix vel spin range id 1
ball fix vel spin range id 2
ball attribute vel 0 1 range id 2 ；可以多次调试竖向速度为1或者-1，横向速度为1或者-1ball attribute density 2.3e3 damp 0.7;力fbond=pb_kn*disp*A---------bond（对于bp胶结模型）其中A表示面积
;力flinear=kn*disp-------linear（对于线性模型）[pb_kn=1e7/20.0]
；pb_kn=kn/A=1e7/20.0
cmat default model linearpbond property ...kn 0 pb_kn @pb_kn pb_ks @pb_kn ... pb_ten 1e6 pb_fa 45 pb_coh 1e6
;pb_kn和pb_ks控制变形，pb_ten 1e6 pb_fa 45 pb_coh 1e6 控制强度
clean
contact method bond gap 1
set timestep fix 1e-3
[ct=contact.find("ball-ball",1)]
def jiance   ；进行监测whilestepping  ；进行回调faxiang_force=contact.force.normal(ct)  ；法向力faxiang_yingli=faxiang_force/20.0     ；法向应力qiexiang_force=contact.force.shear(ct)qiexiang_yingli=qiexiang_force/20.0qiexiang_fengzhi=1e6+math.tan(math.pi*0.25)*faxiang_yingli
；   计算切向力的峰值，来给出包络线c+σ*tanφ
;    lin_stress=comp.x(contact.prop(ct,"lin_force"))/20.0
;    pb_stress=contact.prop(ct,"pb_sigma")
end
history id 1 @faxiang_yingli
history id 2 @qiexiang_yingli
history id 3 @qiexiang_fengzhi
solve time 8
save faxiang1

全量计算与增量计算

全量：累积之前的记录

增量：不考虑之前的记录

线性部分可以有全量和增量，而胶结部分只有增量。

def jiance   ；进行监测whilestepping  ；进行回调faxiang_force=contact.force.normal(ct)  ；法向力faxiang_yingli=faxiang_force/20.0     ；法向应力qiexiang_force=contact.force.shear(ct)qiexiang_yingli=qiexiang_force/20.0qiexiang_fengzhi=1e6+math.tan(math.pi*0.25)*faxiang_yinglilin_stress=comp.x(contact.prop(ct,"lin_force"))/20.0pb_stress=contact.prop(ct,"pb_sigma")；comp.x返回矢量的x方向的数值，应力等于力除以面积；线性力不是从0开始计算的；胶结力是从0开始计算的；pb_sigma指的是Normal stress at bond periphery [stress]
end

切向验证

计算切向力和切向应力，并绘制强度包络线。

def jiance   ；进行监测whilestepping  ；进行回调faxiang_force=contact.force.normal(ct)  ；法向力faxiang_yingli=faxiang_force/20.0     ；法向应力qiexiang_force=contact.force.shear(ct)qiexiang_yingli=qiexiang_force/20.0qiexiang_fengzhi=1e6+math.tan(math.pi*0.25)*faxiang_yingli；计算切向力的峰值，来给出包络线c+σ*tanφend

加胶结方式

首先是无胶结生成颗粒，然后再给颗粒添加胶结。

无胶结生成颗粒

new
domain extent -100 100 ball create position 0 -9 radius 10
ball create position 0 10 radius 10
；生成颗粒有一定的重叠量wall generate box -10 10 -19 20    ；给一个墙体约束ball attribute density 2.3e3 damp 0.7
cmat default model linear property kn 1e7
；生成线性模型cycle 1
solve
save sample
；松散体内部达到平衡，此时模型内部存在力

加胶结

边界和颗粒之间无胶结，但是颗粒和颗粒之间有胶结。随着颗粒运动，线性弹簧由压迫状态变为松弛；而胶结弹簧一开始就是松弛的，由于可以运动，弹簧被拉伸，产生法向力，当大于一定程度时，将会产生破坏。

负号表示拉伸，正号表示压缩。

把时步调快能够增快计算速度。

cmat对新生成的接触赋值，如果已经有旧的cmat，则需要使用cmat apply。

restore sample
[pb_kn=1e9/20.0]
cmat default type ball-facet model linear property ...kn 1e8 ；设置墙体刚度
cmat default type ball-ball model linearpbond property ...kn 1e9 pb_kn @pb_kn pb_ks @pb_kn ... pb_ten 1e6 pb_fa 45 pb_coh 1e6 ；kn和ks控制变形；ten，fa，coh控制强度
cmat apply
；更新模型后，将会把原先的力进行重置为0，而实际计算时会根据kn*A来重新计算
clean
set timestep fix 1e-5   ；减小时步来查看问题
[ct=contact.find("ball-ball",1)] ；去看contact的标号，查找后为1
def jiancewhilestepping    lin_stress=comp.x(contact.prop(ct,"lin_force"))/20.0 ；线性应力，应该有原始积累，全量计算pb_stress=contact.prop(ct,"pb_sigma")                ；胶结应力，增量计算
end
history id 1 @lin_stress
history id 2 @pb_stress
cycle 1
solvecontact method bond gap 1cycle 1
solve
save problem

避免胶结破坏的方法

① 将接触力清零

calm
contact property lin_force 0.0 0.0
ball attribute contactforce multiply 0.0 contactmoment multiply 0.0
;弊端是线性力和胶结力增加的模式不太符合物理实际

② 在改完刚度（cmat apply）之后，再加一个胶结——刚度变化带来的不平衡力计算平衡之后，再添加胶结。

history id 1 @lin_stress
history id 2 @pb_stress
cycle 1
solvecontact method bond gap 1
；之前改变刚度带来的不平衡力已经抵消
；此时再添加胶结，形成bondcycle 1
solve
save problem

PFC教程_pb模型相关推荐

MySQL什么是关系_MySQL教程-关系模型
MySQL教程-关系模型 Date 2021-1-22 Author Haihan Gao 什么是关系模型关系模型本质上就是多个存储数据的二维表表的每一行称为记录记录是一个逻辑意义上的数据表的每 ...
php.ini开启命名空间,Zend Framework教程之模型Model基本规则和使用方法
本文实例讲述了Zend Framework教程之模型Model基本规则和使用方法.分享给大家供大家参考,具体如下: 这里讲讲Zend中的model.其实Zend中的Model处理是相当简单的. 这主要 ...
实践教程 | Pytorch 模型的保存与迁移
实践教程 | Pytorch 模型的保存与迁移在本篇文章中,笔者首先介绍了模型复用的几种典型场景:然后介绍了如何查看Pytorch模型中的相关参数信息:接着介绍了如何载入模型.如何进行追加训练以及进 ...
数据库系统教程——关系模型
数据库系统教程--关系模型这里是理论,建议直接跳转到白话!!! 2.1 关系模型的基本概念基本术语关系的定义与性质关系代数并差交笛卡尔积投影操作选择连接除法以下是白话建议直接 ...
飞控教程基于模型的开发（across）
摘自:https://mp.weixin.qq.com/s/SESLziI-qk_UYLX56Upl9A 飞控教程&基于模型的开发原创 across超越者 across说 7月28日 ACR ...
视图添加字段_Odoo开发教程2-为模型添加字段
开启开发者模式后,我们可通过菜单 Settings > Technical > Database Structure > Models来查看模型的定义.这时搜索 res.partne ...
[翻译：ASP.NET MVC 教程]理解模型、视图和控制器
本篇教程为你提供了ASP.NET MVC的模型.视图和控制器的高级概述.换句话说,即本文向你解释了在ASP.NET MVC中"M"."V"和"C&qu ...
Word2Vec教程-Skip-Gram模型
上海站 | 高性能计算之GPU CUDA培训 4月13-15日三天密集式学习快速带你晋级阅读全文 > 正文共1930个字,6张图,预计阅读时间5分钟. 原文:Word2Vec Tutor ...
二、MMsegmentation 配置教程+训练教程+模型测试（服务器）
时间: 2022年4月8日内容:训练 MM Segmentation 中的deeplabv3 深度神经网络如果想要系统的学习,可以参考官方文档: https://mmsegmentation.re ...

PFC教程_pb模型