Abaqus速成技巧 - 图文(9)

这样看来，无应力状态法似乎是完美的。精确计算出各单元无应力尺寸后，剩余的工作只是“搭积木”。但从以下的分析表明，无应力状态法对误差很敏感，实际适用时要慎重。

（Marchetti M，1999）经验表明，拉索制作精度大约为±10 mm/100 m，相对误差约10-4。这个精度看起来是非常高的。然而，斜拉索的永久应变范围约为2.5×103~3.5×103。因此，100 m的拉索的弹性伸长约为0.3 m。所以，用无应力索长控制索力的精度≤0.01/0.3＝3.3％。以上没有考虑锚固点坐标误差。对现浇的混凝土斜拉桥，这个误差也将达到厘米级。所以，可以得出以下结论：对短索用无应力索长控制张拉要慎重；对长索，由于索长改变绝大部分由垂度变化引起，只要计算精确，用无应力索长控制张拉的精度很高。

（对南京长江二桥，最长索下料长度335.911 m，其中，垂度2.518 m，弹性伸长0.742 m，垂度修正0.046 m。）

（Paik，2001）韩国SeoHae桥是主跨470 m的双塔双索面叠合梁斜拉桥，采用索长控制为主，索力控制为辅。

由于最终构形对索长非常敏感，为了理论考虑拉索锚固点与现场锚固点之间的差别（difference between the locations of the as-built and theoretical cable work points），计算得到的无应力长度必须调整。为了建立as-built work point locations，在工厂仔细检查预制主梁和拉索之间的连接以及主塔内拉索锚固点。根据这些实测数据，拉索安装前作出索长修正。这种方法应用于美国William H. Harsha桥和William Natcher桥（Felber A，2001）。

- 27 -

论坛汇总第一部分

另外还需要考虑的一点是，无应力状态法对计算参数变化是否敏感，即系统的鲁棒性如何。这一点至关重要，至今未见有文献论述。

看了上面的贴子,大开眼界.问题是这几种算法如何用程序实现.用ansys的生死单元可实现否? ADINA应该可以,你们也可以去ASCE的网站去查查,有需要的文章告诉我,我帮大家弄来. 二十七、〖讨论〗abaqus中断裂参量的计算 1、畸变网格怎么划分？

我目前采用的是icem划分网格的方法，在abaqus cae中如何实现不知道，直接修改inp文件肯定可以，但是inp文件修改真的是觉得比较麻烦。2、断裂参量如果计算？ KI和J怎么计算？cae中有办法实现？ Blast：你怎么弄的？

首先你必须有abaqus65

在65版前，abaqus的CAE不支持断裂特性计算装好后在的定义接触、连接的那部分模块可以定义裂纹前缘、或者体内部的隐含裂纹进入后，直接在菜单里面可以找到crack子菜单，定义奇异单元、隐含裂纹都在那里 这样mesh完后，会得到奇异单元，但是要注意，裂尖要用三角形或三棱柱的单元，这样系统才会智能的给你产生所需奇异单元定义好裂纹轮廓后进入step设定那里直接定义历史的输出，可以找到定义J\\K\\T\\G的选项。这样定义后就行了注意塑性材料不能定义K什么的，弹性材料不能定义J 二十八、〖讨论〗关于空隙介质的地应力平衡问题

在STANDARD手册中，提到考虑水时的地应力平衡时，土的密度为DRY density,干密度可是不同于浮容重换算得到的“浮密度”呀！不知道大家怎么理解？另外，定义渗透性时，需要定义水的specific weight，是比重还是水的容重/不太明白，敬请指点！

总应力当然是由干密度和水的重量来平衡（主要用来考虑非饱和区的平衡），这样描述比较方便．specific weight是水的容重，密度乘以重力加速度．

specific weight是水的比重，当用kg时，值为10。无论计算地下水位上、下时，采用gravity加重力荷载时，密度为干密度，由于考虑initial condition=pore pressure,程序算得gestatic步中应力皆为有效应力。具体可参见手册

迷糊了，究竟哪一位大哥说得正确？关于specofic weight!比重是指物质的重量与同体积的纯水在4℃时的重量之比。液体比重是指相同体积的液体重量与水的重量之 比，那究竟是1还是10？另外，我们一般在土工分析中，计算有效应力，可都是采用浮容重呀，没有听说用干容重

建议你看一下帮助文档，给你提个醒，abaqus中采用拉应力为正的规定在多孔介质分析中仍然适用．如果再不明白随便找个土力学．肯定能搞明白其中一些基本概念．明白了再做，不误事． 二十九〖讨论〗接触面应力异常的原因是什么

我在计算接触的时候发现,接触面上总是出现应力异常的情况,比如某个单元的应力很大,和周围单元的应力值不连续,出现断层,而实际上是不会这样的.如图所示.请问大侠这是什么原因造成的?如何避免这种情况? 问题已解决.该面(从面)网格加密后,求解的结果就比较满意了.

此主题相关图片如下：

- 28 -

论坛汇总第一部分

三十、〖讨论〗如果处理构造地应力场

如果地应力场是自重应力场＋构造应力场，应如何处理？

是否可以将材料的强度放大（目的是在施加重力和边界应力条件时不发生塑性变形），使用geostatic求解，然后将节点应力读出，作为ini应力条件再导回原模型。这种方法的关键是确定好K0和构造应力边界条件。

三十一、〖讨论〗如何实现应力驱动的本构关系？

由于现有的有限元程序一般均基于位移法（或混合法），导致整体非线性方程组的一般均针对位移求解（或包括位移变量），于是单元层次的应力－应变关系一般都是应变驱动的，这也是UMAT@ABAQUS或其他自定义本构关系接口一般均是给出应变和应变增量，而要求在接口中更新应力的原因所在。实际上，材料的本构行为往往与荷载相关，如材料的各向异性往往由荷载历史决定，即本构关系ε＝C(σ):σ中的柔度矩阵C是当前应力σ的函数，这就导致一个问题：如何实现此类应力驱动的本构关系，尤其是涉及到应变软化段？迭代法可能是一个办法，但总应变ε中除了包括弹性应变、塑性应变外，往往还包括微裂缝引起的损伤应变，如何区分开来似乎还存在一定问题。 自编程序虽然可以在单元层次实现一些简单的问题，但应用到结构的非线性分析似乎还 是会存在问题，因为基于位移法的有限元不是提供应力增量来求解应变的。望各位就此问题提出自己的意见，谢谢！！

我觉得可以结合uel呀， 就是杂交元或者混合元等等，变量是应力，是不是可以实现UMAT◎ABAQUS接口传递应力增量，而不是应变增量？

就像经典塑性力学的屈服函数是应力的函数，但可以通过returning-mapping算法来转化为应变驱动算法一样，看来可能还得从算法上来找解决办法。在清华BBS上看到有人讨论－－全局积分与局部积分/隐式积分与显式积分，还有人说陈惠发的m叠代就是显式积分（在朱伯芳的书上也讲了），但那个好像不是本构积分把，而是过渡单元的处理，请斑竹能不吝赐教！另，三维的本构可以用在平面应变以及轴对称问题吗？（平面应力肯定不行）就是我的NDI,NSHR尽管可能由主程序传过来，但是我的计算应力更新时全部照着6×6来写？ SMTH上不了,不知道那个m叠代是个什么咚咚？

是Chen 1994年那本书第二卷上的 “塑性状态的确定”一

- 29 -

论坛汇总第一部分

节么？如果是，感觉那个什么比例因子r的求解好像还是要用迭代法来搞定阿？ 我的理解是，首先你要知道分析的问题NDI,NSHR和NTENS的值，然后根据相应的值来编写UMAT，所以很多例子中都有一句判断： IF (NDI.NE.3) THEN WRITE(6,1)

1 FORMAT(//,30X,'***ERROR - THIS UMAT MAY ONLY BE USED FOR ', 1 'ELEMENTS WITH THREE DIRECT S …… 此处隐藏：3513字，全部文档内容请下载后查看。喜欢就下载吧 ……