请教混凝土损伤塑性的收敛问题【已给出解决方法】

predakings

楼主能方便给个修改后，能运行的INP文件么？想学习下，同样在做这个方面的研究 qq：1215655088，信箱 predakings@163.com，感谢

icestar2005

同求INP    也是在做类似的方面研究    416431350@qq.com     感谢楼主

star2010

据别人经验，初始应力值取0.6的极限应力值.
根据塑性应变=总应变-应力值/弹性模量
我不太明白这个0.6倍怎么使用？塑性应变=总应变-应力值/弹性模量，这个式子中的应力是不是每个应力应变的对应值？这跟厨师应力值有什么关系？我不懂，希望指点

yzuwumengying

4# Analyst

请帮帮忙啊，我出现了类似的错误信息，尽管加载很小，步也很小，可是错误不变
***WARNING: For 134649 beam elements either the average curvature about the
             local 1-direction differs by more than 0.1 degrees per unit length
             as compared to the default curvature or the approximate integrated
             curvature for the entire beam differs by more than 5 degrees as
             compared to the approximate integrated default curvature. This may
             be due to a user-specified normal or due to the nodal averaging
             routine used by Abaqus. This difference may cause unexpected
             behavior of the beam and you may want to verify that the beam
             normals are correct for your problem. The elements have been
             identified in element set WarnBeamCurvature1.
***WARNING: For 115228 beam elements either the average curvature about the
             beam tangent differs by more than 0.1 degrees per unit length as
             compared to the default curvature or the approximate integrated
             curvature for the entire beam differs by more than 5 degrees as
             compared to the approximate integrated default curvature. This may
             be due to a user-specified normal or due to the nodal averaging
             routine used by Abaqus. This difference may cause unexpected
             behavior of the beam and you may want to verify that the beam
             normals are correct for your problem. The elements have been
             identified in element set j_WarnBeamTwist.
怎么办呢？ beam normal 应该没太大问题呢 模型是个网格结构 很大 含有很多单元的

cgyhhu

楼主的C40混凝土参数有问题，fc怎么可能达到39.6呢，你的C60才38左右

foreversun

19# bearandyou
我想问下，你预应力筋用的梁单元，截面是怎么定义的？比如说12.7的钢绞线，截面是圆形，你是把他面积换算成相应的直径建模的，还是直接用的12.7？

青春无虑

学习了，谢谢

zb580504

嗯，粘性系数对收敛影响很大的

lvwei04

Analyst 发表于 2010-6-23 10:24
接着就是材料定义部分。这部分我尝试过四个混凝土本构
因为还是在试算阶段，我将钢筋设置为弹性材料，混凝 ...

求完整C40和C50的本构模型，谢谢，感激。。

OSSCU

lvwei04 发表于 2012-2-15 10:25
求完整C40和C50的本构模型，谢谢，感激。。

同求，楼主你有吗？

zippen

Analyst 发表于 2010-12-28 20:01
13# liyangshark
根据规范得到一组值，其中应变部分是总应变，包含了弹性应变，这个是需要去掉的
据别人 ...

但是弹性模型的取值该怎么办？是取初始弹性模量还是某一点的割线模量，在程序里面输入的弹性模型是不是初始弹性模型呢？对峰值应变和极限压应变有什么要求吗？期待解答。

nashone

正在学习中

wawa_0077

发给inp呗，让我们都学习学习

p2p2l5l5k8k8

输入粘性系数0.0005后比较容易收敛

appiedengqing

很好的学习例子，，，

pbmcts

有意思回帖汇总：

11#
谢谢大家了，后来尝试了五方面的工作，顺利的收敛了——收敛了还要根据力学原理等判断计算结果是否可用。
第一是原先的模型网格画得很细，大约有60万节点的样子，一共有78万多单元，计算耗时很长，而且总是出现不收敛。后来将网格改粗了不少，节点变为20多万——有没再细化网格再算，网格粗细对结果还是有影响的。
第二个是改变了钢筋跟混凝土的接触方式。因为这是一个较为精细的模型，混凝土和钢筋分开建模，在最初的一些处理中，并没有用embed，而是合并节点的方式
。后来检查发现，改用了embed。
第三是混凝土的本构例如
*Concrete Damaged Plasticity
30.,    0.1,   1.16, 0.6667,     0
改为
*Concrete Damaged Plasticity
30.,    0.1,   1.16, 0.6667,     0..0005——粘性参数是依据什么取值的？
第四是把所有的step中的初始时间步改为0.0005
第五是本身是按照混凝土标准取得的本构值，再次按照帮助文档的方式进行了重新计算
收敛效果很明显。
14#
根据规范得到一组值，其中应变部分是总应变，包含了弹性应变，这个是需要去掉的
据别人经验，初始应力值取0.6的极限应力值——预应力筋取值？.
根据塑性应变=总应变-应力值/弹性模量得到塑性应变

16#
钢筋还是要用梁单元，很明显，模型中钢筋还是要承受弯矩的，桁架单元是不考虑弯矩影响，只考虑轴向方向的应力
17#
16# bearandyou truss单元通过力偶对提供弯矩抗力，使用beam单元无非多考虑钢筋截面自身的弯矩抗力，应该很小吧，而且beam单元embed在混凝土中，并不约束转动，意义不大？
19#
这里还是我疏忽了，没有看清楚，普通钢筋用的桁架单元，但是预应力钢筋还是用的梁单元，我看成桁架单元。
   一根钢筋抗弯能力比较小，甚至可以忽略，但是大量的钢筋时，就不能忽略混凝土中钢筋的抗弯能力了，具体问题具体对待。本例是高架桥问题，混凝土中钢筋不仅承受拉力还会提高混凝土抗弯能力，不可忽视其抗弯作用。
20#
19# bearandyou 还望讲得详细些，当用于混凝土结构（主抗弯）中时，二者的优势与劣势。预应力筋与普通筋的区别主要是什么呢？——采用什么单元模拟，是否考虑抗弯？如果使用梁单元的话，转动自由度如何处理呀？学习一下哈


25#
楼主的C40混凝土参数有问题，fc怎么可能达到39.6呢，你的C60才38左右
26#
我想问下，你预应力筋用的梁单元，截面是怎么定义的？比如说12.7的钢绞线，截面是圆形，你是把他面积换算成相应的直径建模的，还是直接用的12.7？
28#
嗯，粘性系数对收敛影响很大的
31#
但是弹性模型的取值该怎么办？是取初始弹性模量还是某一点的割线模量，在程序里面输入的弹性模型是不是初始弹性模型呢？对峰值应变和极限压应变有什么要求吗？期待解答。

xiaokui88

加入粘滞系数能很大提高收敛性

ljxabccba

好东西，楼主！

tangxu1992511

学习了不错不错不错哦

huweizip

楼主方便给个这个算例的INP文件么，谢谢了596949434@qq.com