abaqus地应力平衡计算不过去

1、“地应力平衡”的含义、目的、作用

我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。 真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当

*initial conditions,type=stress,input=0.csv

mdb.models['Model-1'].setValues(noPartsInputFile=ON)

（将网格划分精密一些，地应力平衡精度更高一些，2010-7-12注）

注意：导入的inp、csv或其他文件里不得有空行，否则会出现element 0。允许有空格。 --------------------------------------------------------------------

方法一：

桩土地应力平衡不容易通过，可先直接指定应力： *initial conditions,type=stress,geostatic Set-pile-soil,0,0,-510000,30,0.6,0.6

将桩土摩擦系数设为0，容易平衡，输出含有应力s11、s22等的应力文件；

导入应力文件，将摩擦系数提高为一个较小的数值（低于正常的摩擦系数），再计算，输出应力文件； 再提高摩擦系数，导入应力文件，计算，重复上述操作，直到摩擦系数达到正常值。

上述如果不行，可以干脆将桩与的接触改为tie，2010-8-18的模型按照上面的步骤操作不行，后来tie就

*initial conditions,type=stress,input=0.csv

mdb.models['Model-1'].setValues(noPartsInputFile=ON)

（将网格划分精密一些，地应力平衡精度更高一些，2010-7-12注）

注意：导入的inp、csv或其他文件里不得有空行，否则会出现element 0。允许有空格。 --------------------------------------------------------------------

方法一：

桩土地应力平衡不容易通过，可先直接指定应力： *initial conditions,type=stress,geostatic Set-pile-soil,0,0,-510000,30,0.6,0.6

将桩土摩擦系数设为0，容易平衡，输出含有应力s11、s22等的应力文件；

导入应力文件，将摩擦系数提高为一个较小的数值（低于正常的摩擦系数），再计算，输出应力文件； 再提高摩擦系数，导入应力文件，计算，重复上述操作，直到摩擦系数达到正常值。

上述如果不行，可以干脆将桩与的接触改为tie，2010-8-18的模型按照上面的步骤操作不行，后来tie就

运动会的时候，有30个学生加5个家长的集体跳绳比赛项目。在此之前，班里展开了训练和选拔。

第一次我碰巧跳过去了，第二次绊住了，一个趔趄，差点摔倒，心里便有了阴影。每每跳绳的时候总觉得自己跳不过去，这样一暗示，果然就跳不过去，唉，烦人！

校运会跳长绳时我没有报名，可李焕桂却把我的名字写上去了。我和珊珊那时是跳不过去的人，我就去找老师把名字画掉，老师也同意了。参加跳长绳的人，要六点到架空层下去练习，我觉得他们好辛苦啊，我就去加油。同宿舍的杨莎、张晓丽和我一样，都属于笨人一族”，都是不会跳的，可她们担任了甩绳的工作。

我班不会跳长绳的笨人”太多，情况紧急，人数还是不够啊。杨莎、张晓丽到底还算聪明人，边甩绳边看着别人跳，竟然掌握了跳绳的窍门，再一跳，居然都跳过去了。我和珊珊只好光荣地接过甩绳的艰巨任务，可是，手无缚鸡之力的我们甩绳也不够作文力气，没办法，换了个男生，我只好再次退出。

比赛情况，正如想象中那么差，第一名的班级3分钟跳了近300下，可我们班，说不出口啊，3分钟仅仅跳了101下。田老师安慰我们，重在参与，继续加油。

我以为校运会比完集体跳绳就不用大课间的时候跳长绳了，可事实不是这样的，大课间我们又跳起了长绳。这次，田老师给我们

地应力计算公式

（一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生 2003年6月） 主应力计算

根据泊松比?、地层孔隙压力贡献系数V、孔隙压力P0及密度测井值?b可以计算三个主应力值：

????H???A???v?VP0??VP0

?1???????h???B???v?VP0??VP0

?1???H?v???b?dh

0相关系数计算：

应用密度声波全波测井资料的纵波、横波时差（?tp、?ts）及测井的泥质含量Vsh可以计算泊松比?、地层孔隙压力贡献系数V、岩石弹性模量E及岩石抗拉强度ST。 ① 泊松比

??0.5?ts2??t2p2(?t??t)2s2p

2?b?ts2(3?ts2?4?tp)② 地层孔隙压力贡献系数 V?1? 222?m(?tms??tmp)22?b3?ts?4?tp③ 岩石弹性模量 E??t2s??t??t22s2p

④ 岩石抗拉强度 ST?a??b?(3?ts?4?tp)?[b?E?(1?Vsh)?c?E?Vsh]

注：?,?m,?tms,?tmp分别为密度测井值，地层骨架密度，横波时差和纵波时差值。a,b,c为地区试验常数。 其它参数

不同地区岩石抗

3.1 地应力与地应力测量方法简介

地应力，又称原岩应力，也称岩体初始应力或绝对应力，是在漫长的地质年代里，由于地质构造运动等原因产生的。在一定时间和一定地区内，地壳中的应力状态是各种起源应力的总和。主要由重力应力、构造应力、孔隙压力、热应力和残余应力等耦合而成,重力应力和构造应力是地应力的主要来源。地应力的形成主要与地球的各种动力运动过程有关，其中包括：板块边界受压、地幔热对流、地球内应力、地心引力、地球旋转、岩浆侵入和地壳非均匀扩容等。另外，温度不均、水压梯度、地表剥蚀或其他物理化学变化等也可引起相应的应力场。而重力作用和构造运动是引起地应力的主要原因，其中尤以水平方向的构造运动对地应力的形成影响最大。

地应力测量，就是确定拟开挖岩体及其周围区域的未受扰动的三维应力状态，这种测量通常是通过多个点的量测来完成的。地应力测量是确定工程岩体力学属性、进行围岩稳定性分析、实现岩土工程开挖设计和决策科学化的前提。地应力对矿山开采、地下工程和能源开发等生产实践均起着至关重要的作用，所以地应力研究是当前国际采矿界上的一个前沿性课题，近几十年来，世界上许多国家均开展了地应力的测量及应用研究工作,取得了众多的成果。

随着矿区开采现代化进程的不断提高和开采深

ABAQUS建模如何施加预应力

本文参考了百度文库中的文章：

/link?url=dt_VLOGCUf8hUo7A9THhyv7BuSHry71EbLVtBtkWpoiYtkmLxbfk0Io63jsygs6vWbFU7x22HHFv8pIGgPMYkv1lyFXWbgPJqvCodSioUqa

关键字格式：

“*initial conditions, type=stress, input=bb.dat”

上面的关键字，即绿色部分，全部插于*STEP语句之前（如下图），两语句之间不能有空格。施加预应力场只是initial conditions关键字的一个应用，详见abaqus6.8帮助文档，《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。

实例：点焊所产生的焊点中存在着残余应力，本文就是教大家如何完成焊点残余应力的模拟。 原理说明：先在模型上施加一个任意载荷（记为状态1），可得出此载荷作用下模型上的等效节点载荷，然后通过keywords让这个等效节点载荷作用于模型上，使它与之前施加在模型上的载荷相平衡，便得到了一个位移为0的初始状态（记为状态2），该状态下，模型中含有的应力场与状态1相

有意义

第26卷 第5期

岩石力学与工程学报 Vol.26 No.5

2007年5月 Chinese Journal of Rock Mechanics and Engineering May，2007

深部矿井地应力测量方法研究与应用

康红普，林 健，张 晓

(煤炭科学研究总院 北京开采研究所，北京 100013)

摘要：目前地应力测量方法有很多种，可分为力学法、地球物理法及地质构造信息法等。在分析各种测量方法特点的基础上，介绍适用于深部煤矿巷道的地应力快速测量方法——小孔径水压致裂地应力测量法，包括小孔径水压致裂地应力测量原理、测量装置的组成部分和技术特征。在我国煤矿典型的深部矿区——新汶矿区井下巷道中，完成9个测点的地应力测量。其中5个测点的埋深超过1 000 m，最深达1 220 m。测量结果表明，7个测点的最大水平主应力大于垂直主应力；最大水平主应力高达42.1 MPa；最大水平主应力与垂直主应力的比值为1.004～1.550；新汶矿区超过1 000 m深井地应力以水平应力为主，最大水平主应力方向主要集中在N3°E～N43.5

2000年 第24卷 石油大学学报(自然科学版) Vol.24 No.5 第5期 JournaloftheUniversityofPetroleum,China Oct.2000

文章编号:1000-5870(2000)05-0001-03

声发射地应力测量中凯塞点的确定

薛亚东, 高德利

(石油大学石油工程系,北京102200)

摘要:准确判断凯塞点是岩石声发射地应力测量中的重点和难点,目前还没有一种有效的确定方法。对砂岩声发射试验中声发射信号的重标区间(R/S)进行了分析,发现lg(R/S)与时间间隔对数值具有很好的线性相关性,相关系数均大于0.97,这表明声发射信号曲线具有分形特征。研究还发现,不同载荷阶段声发射信号曲线的赫斯特指数(H)可反映声发射信号强弱变化的趋势。当H>0.5时,声发射信号趋于加强;H<0.5时,声发射信号趋于减弱。根据这一规律,通过计算声发射信号曲线的赫斯特指数(或分形维数),便可确定声发射试验中的凯塞点。实践证明这是一种确定凯塞点的有效方法。

关键词:地应力;测量;声发射;分形;凯塞效应中图分类号:TU453 文献标识码:A

引 言

岩石同许多金属材

2000年 第24卷 石油大学学报(自然科学版) Vol.24 No.5 第5期 JournaloftheUniversityofPetroleum,China Oct.2000

文章编号:1000-5870(2000)05-0001-03

声发射地应力测量中凯塞点的确定

薛亚东, 高德利

(石油大学石油工程系,北京102200)

摘要:准确判断凯塞点是岩石声发射地应力测量中的重点和难点,目前还没有一种有效的确定方法。对砂岩声发射试验中声发射信号的重标区间(R/S)进行了分析,发现lg(R/S)与时间间隔对数值具有很好的线性相关性,相关系数均大于0.97,这表明声发射信号曲线具有分形特征。研究还发现,不同载荷阶段声发射信号曲线的赫斯特指数(H)可反映声发射信号强弱变化的趋势。当H>0.5时,声发射信号趋于加强;H<0.5时,声发射信号趋于减弱。根据这一规律,通过计算声发射信号曲线的赫斯特指数(或分形维数),便可确定声发射试验中的凯塞点。实践证明这是一种确定凯塞点的有效方法。

关键词:地应力;测量;声发射;分形;凯塞效应中图分类号:TU453 文献标识码:A

引 言

岩石同许多金属材