lfs03020315 发表于 2009-10-6 22:40:46

如何得到接触面切向应力和法向应力的合力?

最近在研究抗滑桩与边坡土之间的相互作用,在桩土之间建立了四个接触面。计算后可以显示接触面的切向应力和法向应力,但我想得到的是合力的在竖向的分布,即土体在桩正面对桩的推力和桩侧面对桩的摩擦力。
       我的想法:通过吧桩沿竖向分成几个区,把每个区内的接触面节点上的法向应力(或切向应力)与这个节点特有的面积相乘,来得到这个区内法向应力的合力(或切向应力的合力),我查看了彭文斌的FLAC3D一书p71的接触面变量,接触面节点法向应力的变量为in_nstr(p_in),切向应力的变量为in_sstr(p_in,dof),节点Node的特有面积变量为in_area(p_in)。最后把每个区的法向应力(或切向应力)合力连成线就是桩所受的推力和桩侧面对桩的摩擦力。
       但现在我的问题是不能用fish实现上述想法,主要是:我找不到接触面接触面节点的头指针,不能实现循环。p71上有分界面第一个节点的的地址,但我还是实现不了循环,因为运行时指针错误。
       请高手们指点我的fish循环程序(见附件,这只是实现我的想法的第一小步,就仅仅是想实现接触面节点的遍历),急切盼望早日解决此问题。另附我的接触面上剪应力分布图。

lovemm 发表于 2010-2-25 17:07:23

试试壳体、衬砌等支护体吧!彭文斌P222的例子可以参考一下!

abblcn 发表于 2010-5-26 15:37:38

这个问题一直困扰着。。。。

jiangjwhl 发表于 2011-5-18 15:13:08

楼主的点子很好啊,我也在思考呢,想把整个面的法向或者剪切等效应力搞出来

jiangjwhl 发表于 2011-5-18 15:57:47

你的那个小命令流大概思路差不多啊,你可以用f_accum = f_accum + in_nstr(p_in)命令试试看
页: [1]
查看完整版本: 如何得到接触面切向应力和法向应力的合力?