三维弹性应变能计算的fish程序

lessse

有哪位老兄做过三维弹性应变能计算，望指教啊！

[ 本帖最后由 lessse 于 2008-1-10 17:17 编辑 ]

wolf2005

根据弹性力学推导出的弹性应变能公式比较简单，很多参考书上都有推导过程，下面给出公式和FISH代码及结果。
res elasticroadway.sav
config zextra 1
def elasenergy
array ai(6)
p_z=zone_head
loop while p_z#null
dum=z_fsi(p_z,ai)
sxx=abs(z_sxx(p_z))
syy=abs(z_syy(p_z))
szz=abs(z_szz(p_z))
sxy=abs(z_sxy(p_z))
syz=abs(z_syz(p_z))
sxz=abs(z_sxy(p_z))
exx=abs(ai(1))
eyy=abs(ai(2))
ezz=abs(ai(3))
exy=abs(ai(4))
eyz=abs(ai(5))
exz=abs(ai(6))
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz)
z_extra(p_z,1)=ee
p_z=z_next(p_z)
endloop
end
elasenergy
save newelasenergy.sav

huixing_183

看起来很难

zzy518

文章发表在，2007‘河南采矿与安全国际会议(ISMSST)。文章名字为Strain energy density method applications to forecast rock burst and design parameters of pressure relief holes。由science press正式出版

jonavy

原帖由 zzy518 于 2008-1-5 20:50 发表
文章发表在，2007‘河南采矿与安全国际会议(ISMSST)。文章名字为Strain energy density method applications to forecast rock burst and design parameters of pressure relief holes。由science press正式出版

你能把你的大作贴出来共享下么，那个没下载到。另外你考虑弹塑性状态时的应变能求解办法了么？（开挖区不规则情况）

chunhuo420

只要把应变能的公式给出，然后对应输入相应的应力应变，应该可以得到的。

zzy518

声明：我在模拟时plot的是弹性应变能，模型屈服准则是摩尔库伦谈属性准则。这边文章我会近期贴给大家。请关注。

zzy518

忘了说一点，我在做模拟时plot了弹性应变能密度。再在空间体积上积分才是弹性应便变

zzy518

想plot弹性应变能密度，应变应该只用弹性应变，不包括塑性应变。请大家想清楚。而摩尔库伦是理想的弹塑性，在受压状态达到屈服点时应力不变。正利用了这一点得到了弹性应变能密度

zzy518

论文贴给大家，本想的检索结果出来给大家，文章好像被ISTP检索了，检索证明明天才能拿到。

zzy518

呵呵，请大家批评指正

lessse

这是弹性应变能密度吧

[ 本帖最后由 lessse 于 2008-1-9 10:07 编辑 ]

lessse

我觉得不应该加上abs，能量标量在计算时带符号相乘即可：另外要得到弹性应变能用z_volume与密度相乘即可

[ 本帖最后由 lessse 于 2008-1-8 20:58 编辑 ]

zzy518

原帖由 wolf2005 于 2008-1-8 16:15 发表
根据弹性力学推导出的弹性应变能公式比较简单，很多参考书上都有推导过程，下面给出公式和FISH代码及结果。
res elasticroadway.sav
config zextra 1
def elasenergy
array ai(6)
p_z=zone_head
loop while p_ ...

呵呵，那样你的应变应该包括塑性应变了，你所绘出的就不是弹性应变能密度了

lessse

原帖由 zzy518 于 2008-1-10 16:45 发表




呵呵，那样你的应变应该包括塑性应变了，你所绘出的就不是弹性应变能密度了

也对啊，包括了塑性应变呢，以上列出来的公式岂不是不能用了，zzy518兄不介意把你的程序贴出来给大家分享啊！

[ 本帖最后由 lessse 于 2008-1-10 17:20 编辑 ]

jonavy

原帖由 wolf2005 于 2008-1-8 16:15 发表
根据弹性力学推导出的弹性应变能公式比较简单，很多参考书上都有推导过程，下面给出公式和FISH代码及结果。
res elasticroadway.sav
config zextra 1
def elasenergy
array ai(6)
p_z=zone_head
loop while p_ ...

1. z_fsi 给出的是应变增量，这个是包括了剪应变和体应变了么？
2. 既然z_fsi  是应变增量那么我们应该用累积吧吧 strzz=strzz+ai(3)
3.兄台能否给出把6个应变输出到txt文件里fish语言。

希望兄台继续

jonavy

原帖由 zzy518 于 2008-1-10 16:45 发表




呵呵，那样你的应变应该包括塑性应变了，你所绘出的就不是弹性应变能密度了

兄台，因为弹性应变能是对弹性体而言的，我们在使用弹塑性的前提下，应变能降低了，这部分损耗一部分是开挖释放了，另一部分是转化成塑性功了吧，不知这个是不是对，那么你怎么认识这部分塑性应变的，你又是怎么只提取弹性应变的？弹性应变能和应变能密度也请你确切解释一下。对此很感兴趣，能把你的论文中的程序贴出来讲更好。先谢了！

wolf2005

计算时是用的弹性模型（model elastic），当然就只有弹性应变，没有塑性应变啊
如果采用的是莫尔库仑这样的弹塑性模型，才会出现塑性应变，塑性应变能的定义是应力在塑性应变增量方向上做的功，反映在应力应变曲线上就是应力峰值后的应力应变曲线围成的面积，数学公式应该用积分表示。

jonavy

望兄台wolf2005回复16楼的困惑，先谢谢了

[ 本帖最后由 jonavy 于 2008-1-11 12:42 编辑 ]

wolf2005

1.z_fsr是全应变，包括各个方向上的应变分量（xx,yy,zz,xy,yz,xz），FLAC3D中的剪应变和体应变的数值分别是应变偏量第二不变量的平方根和应变第一不变量。剪应变和体应变可以根据应变分量求得。
2.个人认为在小变形情况下，应变增量是个累计值，也就是应变值。
3.可在论坛里搜搜，有很多这方面的帖子。