怎么在模型底面添加抛物线型的边界条件？

苍穹之上

假设路基在填筑完成后，路基底面的位移是抛物线型的，用什么命令才能给路基底面各个节点或者单元施加位移边界条件？用来反算路基在抛物线型沉降条件下，受力情况。初步考虑是根据抛物线计算处位移值的大小，再通过fish施加到各个节点或单元上。
def aa
    p_z=z_near( 0,0,0)
    z_szz(p_z)=-1000
    z_sxx(p_z)=-1000
    z_syy(p_z)=-1000
end
aa
def bb
    p_z=z_near( 0,0,0)  
    gp_zdis(p_z)=-0.01
    gp_xdis(p_z)=0
    gp_ydis(p_z)=0
end
bb

尝试了一下，发现单元的力可以加上，施加位移的话提示“user-defined scalar  cannot  have  arguments”
希望各位大神能给点思路，看看这种边界条件怎么加上去。

dongxx0221

你要搞清基本概念。位移边界实际是加在节点上的，不是单元面上，虽然通常命令写加在面上。
你按照网格尺寸和抛物线坐标，近似确定各节点位置，fix gp zdisp x y z, 就会固定 x y z 最近的gp

苍穹之上

dongxx0221 发表于 2018-1-21 11:39
你要搞清基本概念。位移边界实际是加在节点上的，不是单元面上，虽然通常命令写加在面上。
你按照网格尺寸 ...

好的，多谢大神的帖子，等我按你说的试一下，再来请教

苍穹之上

dongxx0221 发表于 2018-1-21 11:39
你要搞清基本概念。位移边界实际是加在节点上的，不是单元面上，虽然通常命令写加在面上。
你按照网格尺寸 ...

fix gp zdis x y z ,这条命令，flac3d不识别，提示Unrecognized parameter 1(gp)，是不是版本有问题还是，，，还得多多指教下哈