ansys地应力平衡

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就

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

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

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

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

地应力计算公式

（一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生 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为地区试验常数。 其它参数

不同地区岩石抗

SX:X-Component of

stress;SY:Y-Component of stress;SZ:Z-Component of

stress--X,Y,Z轴方向应力。

SXY:XY Shear stress;SYZ:YZ Shear

stress;SXZ:XZ Shear stress--X,Y,Z三个方向的剪应力。

S1:1st

Principal stress;S2:2st Principal stress;,S3:3st Principal

stress--第一、二、三主应力。区分:首先把一个微元看成是一个正方体,那么假设三个主应力分别是F1

F2

F3,那么如果三个力中哪个力最大,就是F1,也是最大主应力,也叫第一主应力,第二大的叫第二主应力,最小的叫第三主应力,因此,是根据大小来定的。

SINT:stress intensity--应力强度,是由第三强度理论得到的当量应力,其值为第一主应力减去第三主应力。

SEVQ:Von Mises是一种屈服准则,屈服准则的值我们通常叫等效应力。Ansys后处理中'Von

Mises Stress'我们习惯称Mises等效应力,它遵循材料力学第四强度理论(形状改变比能理论)。

我们分析后查

ansys提取节点应力

? ? ? ?

|

浏览：530

|

更新：2013-12-11 15:02

GUI操作：在General Postproc——Query Results——Subgrid Solu，选择你想显示的节点。 命令流：1. 最简单的办法是使用NSORT，打印出结果，可以通过控制使其输出到文件

2. 使用apdl能复杂一点，下面是以前经常用的一段命令流，参考着修改一下吧 *CREATE,GET_node_inf,mac,

*GET,Nnod,NODE,0,COUNT !获取所选择的节点总数 *DIM,S_Xyz,ARRAY,NNOD,6 !定义1个数组存放数据 *GET,Nd,NODE,0,NUM,MIN !获取最小的节点编号 *DO,I,1,Nnod,1

S_Xyz(I,1)=Nd !将节点列表放数组第1列 S_Xyz(I,2)=NX(Nd) !节点的X坐标放数组第2列 S_Xyz(I,3)=NY(Nd) !节点的Y坐标放数组第3列 S_Xyz(I,4)=NZ(Nd) !节点的Z坐标放数组第4列

*GET,S_Xyz(I,5),NODE,ND,S,EQV !节点的 von mi

188单元作谐响应分析,求解结束后,我想取出模型中的最大应力值作为参数,然后在接下来的优化当中用该最大应力作为状态变量,请问我应该怎么做啊,注意优化时,对应于每组参数值,最大应力点的位置都可能不同.

ANSYS梁单元如何提取应力

这几个月来一直没花时间做论文，以至于现在手忙脚乱，而且发现很多基础的知识都忘光了，抑或是以前压根就没好好学过，我发现自己是个不折不扣的混日子的，只有在最后关头才会花时间做事情，而且常常是一阵忙，忙完就忘光~~

元旦前赶回老家参加妹的婚礼，回来后一直不敢去见导师，因为导师之前要求我们在元旦之前把论文初稿给他，而我和同门都没给，原因很简单，就是没花时间在论文上，大部分时间都是混，除了上网、看电影外就是和同学一起玩~~总值及其堕落~~~

这两天总算花了点时间做论文，当然主要还是为了能回去过年，这次我的口号是：“不做完论文就不回家！”然而很可能我过年前做不完，所以不得不预先面对可能不能回家过年的打算！

好了，也不多说废话了，转入正题！

这几天论文进展缓慢除了没用心外主要还是不知道到底要提取哪些数据，如何提取，模型是否完善，其实都很糊涂，今天拿第3工况大概地进行静力分析，但随后就是提取数据，然而基础不扎实，以至于不知道如何提取数据。

问题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

引 言

岩石同许多金属材