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

lfs03020315

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

lovemm

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

abblcn

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

jiangjwhl

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

jiangjwhl

你的那个小命令流大概思路差不多啊，你可以用f_accum = f_accum + in_nstr(p_in)命令试试看