FLAC3d中非线性边界条件的施加

pqt01

各位大牛，在FLAC3D里怎么施加非线性边界条件（含根号），比如

SXX=11+ 0.033*H + 0.509*H^(1/2) +0.000415*H^2     (H为垂直深度)

程序代码怎么编写呢？
许多书籍里只提到线性边界条件的施加，如apply sxx -2.892e7 grad 0 0 0.024e6 range x -0.1 0.1
对于非线性边界，如何设置呢？谢谢！

pqt01

怎么没有回复？

wqq8889

好像真没有这样的功能，我看了一下命令参考，apply的关键词里没有直接实现非线性边界的关键词。不知道新的版本里有没有。可能你得尝试用fish实现这个功能。

pqt01

def stress_boundary
    pnt = gp_head
    loop while pnt # null
    H=gp_zpos(pnt)
    ax=11+ 0.033*H + 0.509*H^(1/2) +0.000415*H^2
         command
         apply sxx ax ran x ym1_a ym1_b
         apply sxx ax ran x ym2_a ym2_b
         endcommand
    endloop
end
stress_boundary
fish语句实现的是对每个节点赋值，遍历整个模型。边界荷载是在边界小范围内施加的。fish语句还是有些困难

lovemm

太厉害了，学习了！

benjackxu

pqt01 发表于 2012-11-28 14:03
def stress_boundary
    pnt = gp_head
    loop while pnt # null

加个if语句限定边界范围就可以实现了

pqt01

按楼上的意见，代码如下：
def stress_boundaryxx
    pnt = gp_head
   
        loop while pnt # null
        xx=gp_xpos(pnt)
        if xx>ym1_a
          if xx<ym1_b
                H=gp_zpos(pnt)
                     ax=11+ 0.033*H + 0.509*H^(1/2) +0.000415*H^2
        endif
        endif
        command
              apply sxx ax ran x ym1_a ym1_b z H
            endcommand
        pnt=gp_next(pnt)
    endloop
end
stress_boundaryxx
还是有问题，主要的问题在下面的ran范围上
            command
              apply sxx ax ran x ym1_a ym1_b z H
            endcommand
如何将某深度处的应力施加到同一深度处的单元上？

rhgbgrimm

pqt01 发表于 2012-12-19 21:45
按楼上的意见，代码如下：
def stress_boundaryxx
    pnt = gp_head

回复楼主，将z的范围调整一下，就可以了
command
     apply sxx ax ran x ym1_a ym1_b z H Hmax（Hmax模型最高高度）
endcommand

rhgbgrimm

pqt01 发表于 2012-11-28 14:03
def stress_boundary
    pnt = gp_head
    loop while pnt # null

高手啊，学习了，谢谢