Abaqus 使用中的一些经验 - 图文 下载本文

论坛汇总第一部分

对了,还有一点就是平衡地应力要用*geostatic来做,而不能用*static来做。用static不能让初始位移为0。

十二、有关初始应力的问题

一根悬臂梁,在自重作用下弯曲,求出在自重作用下的单元应力,然后单独把求得的单元应力加于梁上,此时梁的扰曲线应与自重作用下的扰曲线关于水平面对称,如果把自重应力和自重都加上,作为梁的初始状态,则此时梁应处于自重平衡状态,没有变形。

ansys就是利用这种原理施加初始应力场,作为分析的起始步,但在abaqus中,初应力不能从前面的结果文件中读入,在土工问题中,如果用*initial conditions,type=stress,geostatic来定义初应力,按前面所述,如果只加初应力(不加自重),在二维情况下这种与自重作用下的反向变形很明显,此时的地压平衡也很容易实现;在三维情形,用这种方法加上的初始应力就没有这种与自重作用下的反向变形,因此三维的地压平衡就很难实现(重力场作用下产生变形,而初应力作用下土体没有变形)。因课题涉及三维问题,地压平衡这一步总是实现不了,没法往下做,甚为郁闷问题,请指教一二,或是有无这方面的例题。BOW!!!!!!

能否解释解释, 只加初应力(不加自重),产生的反向变形是什么意思?难道说, 只加初应力(不加自重),地表要向上移动?是的,你可以试试,这种情况下地表只产生向上的位移。

三维地压平衡问题已顺利解决(在此感谢版主cafe的关心!!),问题出在对abaqus 坐标轴的理解上,在二维问题中,坐标轴的竖向是2方向(y轴),初始有效应力就是加在2方向上,当然土体的自重也就加在2方向上了: *dload elset_name,grav,10,0,-1,0而在三维问题中,坐标的竖轴就当为3方向(z轴),而不是我们在cae界面中所见的2方向,所以初始有效应力是加在3方向上(而在二维情形,有效应力是加在2方向),因而,在加土体自重时应当加在3方向上: *dload Elset_name,grav,10,0,0,-1

十三、【讨论】再谈初始应力的施加问题

岩土问题的研究都要涉及到一个初始应力的问题,其施加问题一般都要遇到。下面我就此问题谈谈个人的观点:初始应力包括垂直和水平两个方向的力,所以我们在施加的时候要注意两个方向都要有作用的力。一般水平和垂直方向的比是K=1-sinφ .

- 17 -

,各位有无做过这方面的

论坛汇总第一部分 但是我现在有个问题,在abaqus 中施加初始应力的时候我认为是在第一步中施加,可是在create load 对话框中我并没有看到所需要的力,请高手指教。附件中。

176988-a1.rar (0.0k) 如果考虑构造应力,水平应力和垂直应力的比一般为2:1

对于考虑构造应力的话的,水平应力和垂直应力之间的比值并不是2:1。不知你是哪里的结论,请告之,谢谢。因为据我所知,地应力的主要来源是地质构造运动和岩体自重。他们之间的比值是变化的,不是固定的,当然详细的解释有兴趣的朋友可以参阅有关的书籍,比如高等岩石力学,矿山岩体力学,那里有比较详细的解释。希望和楼上的小姐共勉,和我联系,我的邮箱是zxg2018@sina.com.cn我们一起探讨这个问题。

我的说法可能让你误会了,一般情况下,若有构造应力的情况(一般为水平方向),根据各国地应力的测试资料,通常来说,在一定的深度内水平应力大于自重应力,他们比值不是固定的,我们在实际计算时,若考虑构造应力,而且没有实测的地应力资料,可以取2:1

十四、【讨论】紧急求救:geostatic 分析中如何加初始应力和侧压系数? 用AQAQUS CAE时,在step模块选用geostatic 分析时, 如何体现侧压系数? 初始应力如何体现?是不是必须写keywords, 请给一简单例子。 小妹先谢了

cae里不能实现初始应力的施加,可在inp文件里添加,祥见document里abaqus keywords reference manual 里*initial conditions,type=stress,geostatic里的说明。

十五、讨论】再谈ABAQUS中开挖工程的初始应力场问题!

对于涉及开挖问题来说处理好初始应力问题是至关重要的,目前,常常两种方法,其一就是直接把重力场(也可以有其它力存在)始终作用在模型上,进行开挖模拟,这样求得的应力场为实际应力场,而求得的位移场需要减去初始位移场,其二就是想办法读进初始应力场到模型中(没有初始位移),然后进行开挖,这样所求得的位移场和应力场均为实际的位移和应力场。本人认为第一种方法最好,因为我们的初始应力场既可能为弹性也可能局部进行塑性,而第二中方法似乎只能处理弹性初始场,如果真是这样,不知在ABAQUS中能够实现第一种方法,也就是能否方便实现各开挖步中位移场的相减运算,请高手指点!

可以得,在inp里修改或者cae里也行

十六、在inp文件中加地基初始应力(not geostatic),怎么读不进去呢? 采用C3D10MH单元,命令如下: *Initial conditions, type=geostatic 1,1,s11,s22,s33,s12,s13,s23 1,2,s11,s22,s33,s12,s13,s23 1,3,s11,s22,s33,s12,s13,s23 1,4,s11,s22,s33,s12,s13,s23 2,1... 2,2... 2,3,... 2,4,... ... ... ...

该单元四个积分点,不知道是不是应该这么加。abaqus/pre.exe错误提示信息如下:

***ERROR: in keyword *INITIALCONDITIONS, file \

***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO INTERPRET SOME DATA. SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION

***ERROR: AN INITIAL CONDITION HAS BEEN SPECIFIED ON ELEMENT 0 BUT THIS ELEMENT HAS NOT BEEN DEFINED

- 18 -

论坛汇总第一部分

CARD IMAGE: , 1, -5.4809E+04, 27.20, -2.6842E+04, 5.717, 0.9051, 1.552 我试过只按一个积分点加,错误信息同上。 那位帮看看可能是怎么回事?多谢了!

还有一个问题,每个单元是四个积分点,但是abaqus只能取一个,作为初始应力平均加到整个单元上。请问能不能向输出单元应力那样把不同的初始应力分别加到每个积分点上,或是四个节点上也行啊(最起码保证应力在单元内线性)! 你是不是先用abaqus的cae建的模型然后自动剖网格的吧,然后再利用文本编辑器编辑inp文件,是不是,这种情况他都会提示unknown id,他把节点号当成了几何信息里的set,我也不知道更好解决的方法是什么,这个建议你试试吧,你可以把关于inp文件的instance,assembly的信息都删掉,把inp文件里的模型,节点,单元组成形式搞成orphan单元的形势,就像那些example problem manual里的文件一样就可以了

这条经验很宝贵,ABAQUS DOCUMENTATION里面恐怕没有。虽然搞不清原因,但确实我改了ABAQUS自动生成的INSTANCE,ASSEMBLY之后就好用了。这是条值得加分的帖子!真是很感谢inchesnail!

还有一个问题,每个单元是四个积分点,但是abaqus只能取一个,作为初始应力平均加到整个单元上。请问能不能向输出单元应力那样把不同的初始应力分别加到每个积分点上,或是四个节点上也行啊(最起码保证应力在单元内线性)! 十七、【讨论】桩土相互作用中的初始应力场问题处理

在城市地下工程开挖的过程中,常常不能避免近邻的桩基(包括桥梁桩基和建筑物桩基),要解决好近邻桩基的开挖问题,首先需要通过接触面来反映桩土之间的相互作用,其次是要解决好初始应力场问题,接触问题在ABAQUS中相对容易,可是初始应力场问题比较棘手,如果不考虑桩基问题或者考虑桩土完全耦合在一起,初始应力场问题都很好解决。如果通过滑移接触单元来模拟桩土间的相互作用,由于ABAQUS中对用户隐藏了接触单元,而桩单元和土单元中的内应力是重力(假定没有构造应力)和接触力共同引起的,而我们在进行内应力平衡时却只能考虑重力,接触力无法考虑更不能施加,因此,ABAQUS软件处理这种问题似乎有些缺陷,问题要解决需要ABAQUS软件对用户开放接触单元,这样我们变能对桩和土施加接触外力实现真理的平衡,不知各位同行是否同意这种说法! 讲得好!好问题!

这个问题我也注意到啦,其实在土木的这个问题处理上,老板本的界面接触单元,包括ISL家族中的几员,其实可以清晰地做到这点,可惜ABAQUS的趋向是更加的傻瓜和黑箱化(可能这更利于技术支持的赢余),由于对最新版本的疯狂追求,很多可能还没有享用,就扔掉啦。很巧合的是,一个相似的问题在ABAQUS的2004年会上也提出了,对以上问题有同样的置疑。个人认为ABAQUS的5.8-10和6.2-7以及6.4-1都是几个革新性的好版本,其余的很难...,因为我的一个大型蠕变问题的子程序在5.8和6.4都可以顺利运行,在6.3上竟然报出一些错误。可能,功力还差三成啊。 谢谢总管的关注和加分!

目前版本中这个问题的解决需要考虑另外一种方式来模拟有桩土相互作用的开挖问题。我们模拟开挖问题主要有两种方法,一种是在初始位移场为零的真正初始应力场中进行开挖模拟,这种情况下,在随后的载荷步中所得的位移场和应力场均为实际的位移场和应力场;另一种是在初始应力场不为零的真正初始应力场中进行开挖模拟,这种情况下,在随后的载荷步中所得的位移场需要减区初始位移场才为当前步的实际位移场,所得的应力场则是实际的应力场。在ABAQUS中这两种方法均可以采用,故遇到类似问题时,我认为可以采用后一种办法。 有没有监测资料,可以反演出来初始应力场

十八、【原创】如何彻底解决地表不水平下的初始应力场问题?

我以前一直用ANSYS,我是做地下工程的,最近觉得ABAQUS是一个好东西,但是对于地面不水平的情况如何通过文件读取初始应力的问题,是岩土开挖工程中至关重要的问题,一般是通过先加重力计算一次,然后通过后处理获取单元中心的应力分量,然后把每个单元的应力分量放入inp文件中,我的表达方式为: *initial conditionds,type=stress,input=xx.inp 其中文件xx.inp每行的格式比如为: 1(单元号),s1,s2,s3,s4,s5,s6

在运行中程序似乎只认单元集名(ESET),单元号难道不行还是有别的原因,请多多指教,谢谢回复. 单元号也是可以的。看看下面的例子:

*INITIAL CONDITIONS,TYPE=STRESS,GEOSTATIC 1,-2.12E+07,0.0 ,-2.12E+07,-200000,2.094E-01,6.071E-01

- 19 -

论坛汇总第一部分

2,-6.35E+07,0.0 ,-6.35E+07,-200000,2.629E-01,4.779E-01 3,-1.06E+08,0.0 ,-1.06E+08,-200000,2.931E-01,4.330E-01 4,-1.48E+08,0.0 ,-1.48E+08,-200000,3.124E-01,4.116E-01 5,-1.91E+08,0.0 ,-1.91E+08,-200000,3.274E-01,4.011E-01 6,-2.33E+08,0.0 ,-2.33E+08,-200000,3.389E-01,3.954E-01 ...

但是你这种处理方法可能效果不是很好。尤其是临空面存在的情况下。关键是通过GRAV建立起来的应力场处理方法不是很好。ABAQUS在Geotech上明显存在不足。

10节点四面体单元,加初始应力,能否4个高斯点上加不同应力值?

我用*initial conditions, type=stress,采用C3D10MH四面体单元,10节点,4积分点 帮助上这么说: Stress values given on data lines are applied uniformly over the element.

abaqus输出单元应力的时候是输出每个单元的每个高斯点上的六个应力分量,但加初始应力的时候为什么每个单元只用一组应力分量呢?这样由原来的应力线性变成应力常量了,降低了精确度。能否克服这个缺点? 积分点是计算中做数值积分时才用到,在初始条件中怎么给积分点上加东西呀 十九再谈ABAQUS中开挖工程的初始应力场问题!

对于涉及开挖问题来说处理好初始应力问题是至关重要的,目前,常常两种方法,其一就是直接把重力场(也可以有其它力存在)始终作用在模型上,进行开挖模拟,这样求得的应力场为实际应力场,而求得的位移场需要减去初始位移场,其二就是想办法读进初始应力场到模型中(没有初始位移),然后进行开挖,这样所求得的位移场和应力场均为实际的位移和应力场。本人认为第一种方法最好,因为我们的初始应力场既可能为弹性也可能局部进行塑性,而第二中方法似乎只能处理弹性初始场,如果真是这样,不知在ABAQUS中能够实现第一种方法,也就是能否方便实现各开挖步中位移场的相减运算,请高手指点!

可以得,在inp里修改或者cae里也行

二十、【讨论】ABAQUS中能否做岩土开挖问题的时空效应分析? gremms wrote:

对于开挖问题来说,要真正实现时空效应分析,只有采用三维粘弹-粘塑性本构关系或粘弹-粘塑性本构,不知ABAQUS中本身是否有这些本构关系,如果没有是否可以通过用户程序来实现,新手上路,请各位高手发表一点高见! ABAQUS可以实现三维粘弹-粘塑性本构关系或粘弹-粘塑性本构

二十一、【讨论】在ABAQUS中如何直接建立线或面单元与实体单元间的耦合关系?

由于本人做隧道开挖问题比较多,用了ABAQUS后,最大的感受就是建隧道开挖问题的计算模型非常麻烦,比如我们常常用壳单元(三维问题)或梁单元(二维问题)来模拟初期支护结构,而梁和壳单元与他们相连的实体单元可以处理为完全耦合关系,这个问题在ansys、marc、adina等软件中非常好处理,可以直接在实体单元的所在的几何面或边上划分壳或梁单元,而在abaqus中似乎需要对梁和壳、以及实体单元分别建立Part,再通过耦合关系把他们在节点处耦合在一起,这样处理起来非常麻烦,特别是涉及到大量的开挖步时,更是如此,难道我得放弃abaqus吗,不知各位高手似乎有更好的办法(本人接触 ABAQUS的时间较短)

其实我觉得一样的,你可以建立不同的part再assem一起也可以,至于梁和体连接可以用多点约束,壳和体有shell to solid couple的选项可以用

如何利用abaqus建立组合有限元模型,请大家发表意见。

如果是兩個物件相黏接的話我都是用tie ,並不需要去切成另一個part除非你要接合的位置是比較特殊,才會特別去切 是2个不同的部分之间设置接触面

在ABAQUS中建立接触面的方式只能通过不同的PART来进行,不过处理的方法有两种,一种是通过网格PART来进行,这种办法用来建立一致网格的接触关系时非常好,另一种就是通过几何PART。此外,在ABAQUS中建立壳、体耦合关系时,通过其蒙皮技术能够很好解决。

“shell to solid couple”这种办法只能用于壳的边和实体的耦合,不能用于壳面和实体面的藕荷”其蒙皮技术,是否能详细说明?谢谢!

- 20 -