Flac3D怎样给边界施加面荷载（荷载是关于坐标的正弦函数）

66847789

flac3D的grad功能可以很容易实现，给边界施加三角形的梯度荷载。但如果是面荷载，而且荷载是关于坐标的正弦函数或者是二次函数，应该如何做？请高手指点算法。

xwjwy7833

编写fish函数

66847789

原帖由 xwjwy7833 于 2007-12-29 14:32 发表
编写fish函数

这个问题问题困扰了我比较久，能不能把过程说一下？
通过哪几个变量？
比如说把szz_=sin(x1)的应力加到z=5的顶面
我想了几个办法，但感觉都不对，希望大家多多指点
1、利用ad1=gp_near(x,y,z),x1=gp_zforce(ad)，用循环把顶面所有的节点加上zforce=sin(x1)的节点力，用节点力来近似面力，这个肯定会有误差。
2、用x1=z_xcen(p_z),zstress=sin(x1),然后z_szz(p_z)=sin(x1),用loop把没个顶点单元的应力更改？这个我不知道对不对。
我现在都糊涂了，思路也没有章法，请大家指点错误啊。
谢谢大家。

niebaisong

不是很会啊