（原创）逐个取得模型中interface面上的法向力和切向力！

tg0215

;file
;Int_force_setout.dat
;Used for FLAC2D only
;#########################################################################
; Get the normal froce and shear shear of interfaces.Set out the normal and

; shear force by plot.In the plot , X-axes is the X-cordinate of the ;

interface and  Y-axes is the absolute value of the force.
;#########################################################################
;-------------- written by tg0215 on 2006-9-4 -----------------------------

set echo off
cal int.fin
set echo on

; example model
g 10 11
m e
prop dens 2000 sh 1e8 bu 2e8
m n j=6
ini x add 2.5  y add -0.9 j=7,12
int 4 aside from 1,6 to 11,6 bside from 1,7 to 11,7
int 4 kn 5e8 ks 2.5e8 fric 10
set grav 10,5  large
fix x y j=1
step 1000
;
def save_data
  
  ip=int_pnt
  loop while ip # 0
    id_int=imem(ip+$kicid)  
    A_side_pointer=imem(ip+$kicapt)
    B_side_pointer=imem(ip+$kicbpt)
    n=1
    loop while A_side_pointer # 0
      jj=imem(A_side_pointer+$kidj)
      ii=imem(A_side_pointer+$kidi)
      xtable(id_int,n)=x(ii,jj)
      xtable(id_int+10,n)=x(ii,jj)
      ytable(id_int,n)=abs(fmem(A_side_pointer+$kidfn))
      ytable(id_int+10,n)=abs(fmem(A_side_pointer+$kidfs))
      A_side_pointer=imem(A_side_pointer)
      n=n+1
    endloop
    ip=imem(ip)
  endloop
end
save_data

def set_out
  ip=int_pnt
  loop while ip # 0
    id_int=imem(ip+$kicid)
    ppp=id_int+10
    command
      pl hold table id_int both
      pl hold table ppp both
    endcommand
    ip=imem(ip)
  endloop
end
set_out

tg0215

模型及interface所在位置。

tg0215

法向力－normal force

tg0215

切向力－shear force

tg0215

两点说明
1，interface的两个面，只取了aside面的normal force和shear force；
2，无法取liner和beam单元的interface上的力。

yts1

过来支持你！

zuo-0000

顶一下！

phlow

先学习一下flac2d，不错的资料，哈哈

huixing_183

有人做3D的吗

异乡客

谢谢版主讲解，谢谢了。。。

yisihai

牛人啊,敬仰

深夜来访电话

只是用于2d，先是很高兴，后是很失望

深夜来访电话

不过还是可以学习研究一下的