双轴代码转换

孤亭远见

已于 2023-08-06 21:13:13 修改

阅读量777

点赞数 1

分类专栏： PFC 文章标签：开发语言

于 2023-05-12 09:16:22 首次发布

本文链接：https://blog.csdn.net/weixin_52546971/article/details/130634514

版权

PFC 专栏收录该内容

30 篇文章

订阅专栏

双轴代码转换

二维转三维

5.0转6.0

双轴代码转换

二维转三维

添加z方向数值

domain extent [-width*2.0] [width*2.0] [-width*2.0] [width*2.0] ...
        [-height*2.0] [height*2.0]
set random 10001
；比较标准的话可能设置40000个颗粒

更新墙体指针，并重新定义y和z方向的尺寸/应力

def wall_ini
    wpup=wall.find(2)
    wpdown=wall.find(1)
    wpleft=wall.find(3)
    wpright=wall.find(4)
    wpfront=wall.find(5)
    wpback=wall.find(6)
end
@wall_ini

def computer_chicun
    wlx=wall.pos.x(wpright)-wall.pos.x(wpleft)
    wly=wall.pos.y(wpback)-wall.pos.y(wpfront)
    wlz=wall.pos.z(wpup)-wall.pos.z(wpdown)
end

伺服系数除以面积

gx=1.0*servo_factor*(wly*wlz)/(zongKNX*global.timestep) 
gy=1.0*servo_factor*(wlx*wlz)/(zongKNY*global.timestep)

除此之外，因为2D的y变为了现在3D的z，一些对应的细节也需要更改。

[strainRate=1e-2]
wall attribute zvel [strainRate*wlz] range id 1
wall attribute zvel [-strainRate*wlz] range id 2

5.0转6.0

new 转为 model new

restore sample 转为 model restore "sample"

关于模型的操作在前面需要加上model

字符串需要加上双引号
model cycle 2000 calm 50
ball property "fric" 0.5
model solve
model save "sample"
cmat需要在前面加上contact
contact cmat add 1 model linearpbond method deform "emod" 1e8 "kratio" 1.5 ...
                            pb_deform "emod" 1e8 "kratio" 1.5 ...
                            property "pb_coh" 1e6 "pb_ten" 1e6 "pb_fa" 40 "fric" 0.5 ...
                            range contact type "ball-ball"
contact cmat apply
fish call back -1.0 @computer_strain 改为 fish callback add @computer_strain -1.0

在PFC6.0中，yvel 被改为了velocity-y

→找不到对应字符的情况应该到手册中查找