【讨论】有人做过混凝土内部的预应力钢绞线吗?

jinjinhao

天啊,我做的 U 形渡槽,在沿弧形应力损失是变化的, rebar单元有上万个,难道只能用修改 dat 文件的方法考虑损失吗? 那岂不是要累死我!

jinjinhao

还不明白的是,在 Model->Geometry Atrributes->Simple Geometry->Truss 和 Model->Initial Conditions 里设置初始应变有何不同?后者可以按空间函数加载,是否可以解决应力损失的问题?

ncwuyangguang

我今天也做了个truss和rebar，但结果好像不太一样。？？
请大家指教？？
我把两个in文件都贴上来了！！！

iamperisher

to ncwuyangguang:带这么多垃圾的in文件别人怎么看啊，应该稍加整理。truss文件最后一个线号错误，reabr文件也有点问题。
与同样的模型，整理后的reabr和truss命令文件如下，自己对照看看是哪里有问题。
COORDINATES POINT
entries name y z
1  0   0  
2  0   1
3  1   1
4  1   0
5  0.5 0
6  0.5 1
surface vertex name=1 1 2 3 4
line stra name=5 5 6
material elastic name=1 e=52452 nu=0.2
material elastic name=2 e=1.235541e8 nu=0.4
master idof=100111
load pressure name=1 magn=112
apply-load
1 pressure 1 line 2  
fixity name=hg
'y-tran'
fixboundary lines
1 'hg'
3 'hg'
fixboundary lines
4 'all'
rebar-line name=1
5
egroup twodsolid name=1 subtype=strain
egroup truss name=2 option=rebar rebar-line=1
*对平面问题面积0.005没有意义
line-elemdata truss
5 2 0.005
surface-elemdata twodsolid
1 1
subdivide line name=1 mode=division ndiv=10
2
3
4
gsurface name=1 nodes=8 group=1 meshing=mapped
  
*nodes=4/2 truss 1.29547e-3  
*nodes=4 rebar   1.29547e-3 误差0%
*nodes=8/3 truss 1.30169e-3
*nodes=8 rebar   1.31137e-3 误差0.7%  4/8误差1.2%
*--------------------truss-----------------------------------------
COORDINATES POINT
entries name y z
1  0   0  
2  0.5 0
3  1   0
4  0   1
5  0.5 1
6  1   1
surface vertex name=1 1 2 5 4
surface vertex name=2 2 3 6 5
material elastic name=1 e=52452 nu=0.2
material elastic name=2 e=1.235541e8 nu=0.4
master idof=100111
load pressure name=1 magn=-112
apply-load
1 pressure 1 line 3
2 pressure 1 line 7
fixity name=hg
'y-tran'
fixboundary lines
4 'hg'
6 'hg'
fixboundary lines
1 'all'
5 'all'
egroup twodsolid name=1 subtype=strain
egroup truss name=2
*对平面问题面积0.005没有意义
line-elemdata truss
2 2 0.005
surface-elemdata twodsolid
1 1
2 1
subdivide line name=2 mode=division ndiv=10
4
6
subdivide line name=1 mode=division ndiv=5
5
3
7
gsurface name=1 nodes=8 group=1 meshing=mapped
gsurface name=2 nodes=8 group=1 meshing=mapped
gline name=2 nodes=3 group=2 ncoincid=all

dali

jinjinhao wrote:
在 Model->Geometry Atrributes->Simple Geometry->Truss 和 Model->Initial Conditions 里设置初始应变有何不同?后者可以按空间函数加载,是否可以解决应力损失的问题?

  
jinjinhao 兄有没有试过这个功能？如果能作，简直太好了，我也有这样的问题要做。

jinjinhao

dali wrote:
jinjinhao 兄有没有试过这个功能？如果能作，简直太好了，我也有这样的问题要做。

  
正在摸索,暂时还没弄明白

iamperisher

再以平面问题的那个离题考察发现，当使用节点初应变(initial-conditiuon) 和单元初应变(line-elemdata)计算时，输入同样的初应变数值，但结果确相差2倍！是否都是这样？算法的理论根据在哪里呢？
jinjinhao兄可以试试看是否是这样的？

iamperisher

考虑并试验了一下，个人认为是这样的：
1.adina是将节点初应变转换成单元初应变进行计算的。并且在理论手册的10.6中也谈到了节点初应变和单元初应变可能不同，因此二者存在差异是可能的。
2.对truss可以采用两种方法施加，即单元初应变和节点初应变，而且节点初应变可利用空间函数施加。
3.对于rebar则只能施加单元初应变，当为不均匀分布初应变时，只能通过修改dat文件进行。

catkin

根据各位的讲解我做了一个简单的预应力混凝土梁的例子：
材料：C60混凝土；Φ7的刻痕钢丝,总张力为338kN
压力为0.35N/mm*mm
图如下

catkin

这是*.in文件，请各位大侠看看对不对？请指教，谢谢！！支撑应该是个面，现在简化成线了，局部应力可能大了！
混凝土材料是用的弹性材料代替的，如果用混凝土材料应该怎么做，参数怎样输入？

weikaihere

虽然很久前的东西了，不是需要学习一下