连续刚构的空间仿真分析

斑竹说的很好，鱼与熊掌不可兼的啊！38万节点的模型，内存没有2G，恐怕算不动吧。

     其实ansys的梁单元没法加预应力（主要是预应力损失），算蠕变又要自定义材料性质，比起其它的桥梁专用程序显得很麻烦，全桥分析很多设计院都开始用midas了。但是在局部分析上，ansys还是有它的优势的。

关于查看主拉应力的经验：
S1,S2,S3的等值线图，根据材料力学知识按大小排序s1>s2>s3
就是结构中任一结点在剪应力为零的三个主平面上的正应力
我们根据S1的等值线图来判断拉应力
根据S3来判断压应力
不知道各位还有什么高见

我们最近也在做连续刚构的局部分析,跟楼主的问题差不多,只是我们把三向预应力分开分析,现在横向和竖向都做好了,就差齿板分析,请问楼主是怎么处理既有竖弯又有平弯的预应力筋的?如果用link单元,直接降温或给一个初应变,那样算出来是不对的,那样只适用于直线配筋,请问楼主的处理方式,我们现在在尝试把预应力筋作成实体模型,然后用接触分析.谢谢.

下面是我们有关齿板的模型

初学者，问题见笑了，问下楼上
“既有竖弯又有平弯的预应力筋，如果用link单元,直接降温或给一个初应变,那样算出来是不对的,那样只适用于直线配筋”
为什么？能给小弟解释下吗？

如果你考虑孔道的影响（局部分析时的最不利状态，未灌浆），那么直接降温或给一个初应变，由于预应力筋与混凝土之间没有任何联系，那样做的结果是预应力筋变形后直接进入到混凝土中，这从后处理中可以看出来，这与实际情况不符。本人拙见。

我没有考虑过孔洞的效应，但是我想
可以考虑采用先建立模型，然后分网，建立单元和节点，通过节点连接钢筋单元，建立出曲线钢筋，然后采用接触分析，应该可以解决，过段时间我研究下，有结果的话我贴出来，不知道大家还有什么高见

空间的局部分析是必不可少的，如果能够进行全桥分析就更好了，水兄能否再谈谈具体的建立模型的思路

[ 本帖最后由 沧海学生小 于 2006-10-24 21:52 编辑 ]

这是我发在okok上的帖子，转了过来
可以采用solid单元或者shell单元，其各有优缺点。
solid单元实体建模型可以更加真实的考虑桥梁的钢筋张拉，挂篮的移动，模板的架设，混凝土的浇注等过程，但缺点是单元数量过多，节点数量过大，计算时间过长。
shell单元克服了这些缺点，但是在考虑空间问题上又没有solid单元来的精确。
前段时间我们刚刚完成了一个连续刚构的分析，把思路贴出来供你参考：
1、材料：混凝土采用solid45单元，钢筋采用link8单元
2、关键技术：采用单元生死和迭加原理模拟施工过程，采用节点力的转换模拟挂篮施工过程
3、建模型过程
1）确定各断面的坐标
2）通过关键点沿抛物线复制的方法组成小体积，再将小体积粘合。
3）用工作平面切分出钢筋的具体位置
4）将线和体赋予不同的特性
5）划分钢筋和混凝土的单元
6）定义不同的组，比如0#块，1#块等
7）按施工阶段，分别杀死不需要的单元，等到需要的时候再激活它

楼主能否提供一下局部分析的命令流？或者告知局部分析的一些情况，如采用的单元类型，约束的施加，荷载的施加，建立局部模型目的是什么，查看应力分布吗？
另外，请楼主介绍一下活载是如何考虑加载的及全桥模型的考察目标，如楼主提到的腹板主拉应力是其考察目标之一，另外还考虑那些目标呢？盼复，谢谢。

[ 本帖最后由 haporth 于 2006-10-26 12:51 编辑 ]

1、局部分析混凝土应该采用实体solid单元（45、65、还是95酌情处理）,钢筋采用link8单元。
比如分析0#块的局部应力我的思路是这样的
我们至少要做到2或者3号块，在其上施加相应边界条件。如果只建模到0#块的话，在其上施加边界条件，那么边界条件造成局部应力的影响会很大，造成分析的不准确。这从《弹性力学》的圣维南原理可以得出结论！（圣维南原理是面对杆的，不知道能否适应实体单元，高手指正）
至于边界条件的施加，我们先要进行平面计算，算出相应位置的力、弯矩等，施加到实体模型上（实体单元不能施加弯矩，可以考虑采用对称力）
2、活载加载当时我们并没有考虑，但是我想可以参照《algor、ansys在桥梁工程种应用》这本书
全桥模型的考察目标主要有：箱梁的纵向应力分布规律、剪力滞效应、腹板主拉应力的分布规律、扭转效应等
这是本人的浅见，如有问题请指教

我用beam１８８建的刚构桥模型，有很多问题要请教下

用beam188单元建连续刚构，其截面处理有两种方法
1、在cad内画出截面，然后定义成面域，再倒如到ansys里面
2、直接在ansys里面用板单元定义截面
不知道你用的是哪种，需要什么方面的帮助？

这是Algor、Ansys在桥梁工程中的应用---------------连续刚构桥的建模及分析的例子，详见附件
/prep7
/title,the analsys of ganggou bridge
et,1,solid95
et,2,link8
mp,prxy,1,0.1667
mp,dens,1,2600
mp,ex,1,3.5e10
mp,ex,2,1.95e11
mp,dens,2,7800
mp,prxy,2,0.3
areagjx=1.4e-4
*dim,gjarea,array,4
gjarea(1)=15*areagjx
gjarea(2)=18*areagjx
gjarea(3)=20*areagjx
gjarea(4)=22*areagjx
F=195.5e3
E=1.95e11
L=1
detl=F*L/(E*areagjx)
istrn=detL/(detL+L)
r,2,gjarea(1),istrn
r,3,gjarea(2),istrn
r,4,gjarea(3),istrn
r,5,gjarea(4),istrn
acel,,9.8
k,2,-2.45,0
k,3,-3.25,0
k,4,-6,0
k,5,-6,-0.18
k,6,-3.25,-0.45
k,7,-2.45,-0.45
k,11,-3.25,-3.9
k,12,-2.45,-3.9
k,13,-2.45,-4.5
k,14,-3.25,-4.5
ksym,x,2,14,1,100
kgen,2,9,14,1,0,-0.00125*(40-10)*(40-10)+2,10,2000
kgen,2,109,114,1,0,-0.00125*(40-10)*(40-10)+2,10,2000
kgen,2,9,14,1,0,-0.00125*(40-20)*(40-20)+2,20,3000
kgen,2,109,114,1,0,-0.00125*(40-20)*(40-20)+2,20,3000
kgen,2,9,14,1,0,-0.00125*(40-30)*(40-30)+2,30,4000
kgen,2,109,114,1,0,-0.00125*(40-30)*(40-30)+2,30,4000
kgen,2,9,14,1,0,-0.00125*(40-40)*(40-40)+2,40,5000
kgen,2,109,114,1,0,-0.00125*(40-40)*(40-40)+2,40,5000
kgen,2,2,14,1,0,0,1,1000
kgen,2,102,114,1,0,0,1,1000
kgen,2,2,7,1,0,0,10,2000
kgen,2,102,107,1,0,0,10,2000
kgen,2,2,7,1,0,0,20,3000
kgen,2,102,107,1,0,0,20,3000
kgen,2,2,7,1,0,0,30,4000
kgen,2,102,107,1,0,0,30,4000
kgen,2,2,7,1,0,0,40,5000
kgen,2,102,107,1,0,0,40,5000
k,100,0,0,0
k,200,1,0,0
k,300,0,1,0
k,400,0,-4.5,0
k,500,0,-14.5,0
l,400,500
k,21,-3.75,-0.25,-40
k,22,3.75,-0.25,-40
k,23,3.75,-0.25,40
k,24,-3.75,-0.25,40
k,25,-6,-0.1,-40
k,26,-3.75,-0.25,-40
k,27,-3.75,-0.25,40
k,28,-6,-0.1,40
k,29,6,-0.1,-40
k,30,3.75,-0.25,-40
k,31,3.75,-0.25,40
k,32,6,-0.1,40
a,21,22,23,24
a,25,26,27,28
a,29,30,31,32
*do,i,1000,5000,1000
v,2+i-1000,i-1000+7,i-1000+107,i-1000+102,2+i,7+i,107+i,102+i
v,2+i-1000,i-1000+3,i-1000+6,i-1000+7,2+i,3+i,6+i,7+i
v,3+i-1000,i-1000+4,i-1000+5,i-1000+6,3+i,4+i,5+i,6+i
v,7+i-1000,i-1000+6,i-1000+11,i-1000+12,7+i,6+i,11+i,12+i
v,11+i-1000,i-1000+14,i-1000+13,i-1000+12,11+i,i+14,i+13,i+12
v,12+i-1000,i-1000+13,i-1000+113,i-1000+112,12+i,i+13,i+113,i+112
v,100+2+i-1000,100+i-1000+3,100+i-1000+6,100+i-1000+7,100+2+i,100+3+i,100+6+i,100+7+i
v,100+3+i-1000,100+i-1000+4,100+i-1000+5,100+i-1000+6,100+3+i,100+4+i,100+5+i,100+6+i
v,100+7+i-1000,100+i-1000+6,100+i-1000+11,100+i-1000+12,100+7+i,100+6+i,100+11+i,100+12+i
v,100+11+i-1000,100+i-1000+14,100+i-1000+13,100+i-1000+12,100+11+i,100+i+14,100+i+13,100+i+12
*enddo
v,7,12,112,107,1007,1012,1112,1107
cskp,100,0,100,200,300
vsel,all
csys,100
vsymm,z,1,100,1,10000,0,0
csys,0
asel,s,loc,y,-4.5
vdrag,all,,,,,,1
allsel
vsel,all
wprota,0,0,90
wpave,2.7
vsbw,all
wpave,3
vsbw,all
wpave,-3
vsbw,all
wpave,-2.7
vsbw,all
wpave,2.45
vsbw,all
wpave,3.25
vsbw,all
wpave,-3.25
vsbw,all
wpave,-2.45
vsbw,all
wpave,2.2
vsbw,all
wpave,3.5
vsbw,all
wpave,-3.5
vsbw,all
wpave,-2.2
vsbw,all
wpave,1.95
vsbw,all
wpave,3.75
vsbw,all
wpave,-3.75
vsbw,all
wpave,-1.95
vsbw,all
vsel,s,loc,y,-0.45,0
vsba,all,1
vsba,all,2
vsba,all,3
allsel
vglue,all

对的，判断应力就是这样判断的：
S1,S2,S3的等值线图，根据材料力学知识按大小排序s1>s2>s3
就是结构中任一结点在剪应力为零的三个主平面上的正应力
我们根据S1的等值线图来判断拉应力
根据S3来判断压应力
38万的节点，如果是静力分析的话，1G内存足够了。
我建立了一个局部模型，6万多个节点，静力分析，1G内存，双核cpu，5分钟就算完了。
但如果考虑了单元生死的话，就相当于非线性分析了，算了2个小时的。