tg0215 发表于 2006-9-4 16:35:10

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

本帖最后由 cloud328 于 2010-6-8 13:29 编辑

;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 andY-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.5y 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,5large
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 发表于 2006-9-4 16:38:32

模型及interface所在位置。

模型及interface所在位置。

tg0215 发表于 2006-9-4 16:39:23

取出的interface法向力－normal force

法向力－normal force

tg0215 发表于 2006-9-4 16:40:00

取出的interface切向力－shear force

切向力－shear force

tg0215 发表于 2006-9-4 17:27:36

两点说明

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

yts1 发表于 2006-9-5 14:40:46

过来支持你！

zuo-0000 发表于 2007-12-21 10:19:44

:) 顶一下！

phlow 发表于 2008-6-29 09:28:46

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

huixing_183 发表于 2008-7-3 11:26:22

厉害

有人做3D的吗

异乡客 发表于 2008-9-5 09:01:07

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

yisihai 发表于 2009-10-26 17:42:35

牛人啊,敬仰

深夜来访电话 发表于 2010-6-8 12:10:45

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

深夜来访电话 发表于 2010-6-8 12:11:40

不过还是可以学习研究一下的
页: [1]
查看完整版本: （原创）逐个取得模型中interface面上的法向力和切向力！