令人头疼的纤维梁模拟结果

wei1012

一、模拟纤维梁的目的
了解纤维梁中混凝土的应力应变关系；
二、算例设定
对长度为0.2m的梁单元，设置纤维截面，采用*integration beam设定，单元截面为1m*1m，混凝土材料假定为C45，受压强度fc=30e6 Pa，受拉强度ft=3e6 Pa，梁单元内不含钢筋，仅考虑混凝土（即fracr=0）。梁一端固定约束，一端施加受压位移｛（0，0），（1，-0.001）｝，对应最大应变为-0.005。
混凝土材料定义：
*MAT_RC_BEAM_TITLE
concrete
$#     mid        ro      eunl        pr        fc       ec1      ec50     resid
         1 2700.00003.0000E+10  0.300000 3.0000E+7  0.002200  0.005000  0.200000
$#      ft     unitc                                   esoft     lchar    output
3.0000E+6  1.000000     0.000     0.000     0.000 3.0000E+9     0.000  1.000000
$#   fracr   ymreinf   prreinf   syreinf   sureinf      eshr       eur    rreinf
     0.0002.0000E+11  0.300000 4.0000E+8 4.7000E+8  0.008000  0.016000  4.000000
纤维梁定义：
*INTEGRATION_BEAM
$#    irid       nip        ra      icst         k
         1         0     1.000        11         2
$#      d1        d2        d3        d4      sref      tref        d5        d6
  1.000000  1.000000     0.000     0.000  1.000000     0.000     0.000     0.000
施加的位移曲线：
*DEFINE_CURVE_TITLE
diaplacement -0.005
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         2         0  1.000000  1.000000     0.000     0.000         0
$#                a1                  o1
               0.000               0.000
            1.000000           -0.001000
三、令人头疼的计算结果

位移结果是对的！

应变曲线图

应力曲线图
结果有正负，是震荡吗???
为什么？

四、不好掌握的LSPP操作
看帮助里面提到查看History结果，需要打开D3THT，结果我打开这个文件反而没有结果，没打开的话却有。

更为关键的是，该如何调整，让dyna的结果合理？


http://blog.sina.com.cn/s/blog_642b4b710101dnvb.html

yuminhust2005

将应力和应变时程曲线合并成应力-应变曲线。看是否和你定义的一样。
因为ls-dyna默认计算的是动力问题，有振荡也不奇怪。
如果只是为了测试应力-应变关系，建议梁只划分一个单元。
或是加大计算时间消去动力效应，或是采用隐式计算。

wei1012

k文件中就只有一个单元。位移结果是对的，应变确是有数量级的差！

wei1012

我把纤维梁的截面尺寸改小为0.001*0.001，结果显示所有积分点的应变均相同，但还是存在正负问题，已经有点进步了！

wei1012

继续给模型增加阻尼，阻尼值为0.99，准静力加载，得到如下图的应变结果：

又有进步！！！

yimengweima

学习一下纤维模型

orange_yz

请问，哪里有教模拟纤维梁的书呀？或者例子？我刚刚学

orange_yz

请问怎样定义包含钢筋的纤维梁呢。我正在做这个模拟。卡住了。打算用type11,solid box截面形式。可是怎么也没有找到把其中几个纤维定义成钢筋的地方。