对MIDAS等软件的学习、理解
Y_n3O@, -------继续探讨笔者“曲线梁桥”程序与
u' Q82l&Y MIDAS、ANSYS、“桥梁博士”对比计算中
s2wDJ| 出现显著差异的原因------
$eG_LY 1v 2006-02-19
NsSZ?ky 珠海东部久远科技有限公司 孙广华博士、副研究员
#;W4$ q 最近笔者进入MIDAS官方网站
http://cn.midasuser.com，阅读并下载了一些文件，仔细地进行了研究。有以下了解：
U s!ZQ#pP 1． 证实了笔者以前的推测：MIDAS计算梁的时候，是用形心代替剪力中心的。
@ =~k[o 进入网站的“技术资料”——“变截面梁”，出现关于变截面梁单元基本假定、单元刚度矩阵生成方法、公式等资料（英文）。其中基本假定的第3条是：变截面梁单元的弯曲中心和剪力中心是重合的（The element cross_section bending and shear centers are coincited）； “梁单元”（即等截面梁）也有同样的规定。
.P/0 `A{& 2．“技术资料”——“变截面梁”，基本假定的第4条是：梁单元可以承受轴向力、绕横截面两主轴的弯矩和绕其重心轴的扭矩（The element is capable of resisting axial forces,bending moments in the plane of its cross_section, and twisting moments about its centroidal axis）,注意：它不考虑与梁扭转有关的另一项内力——双力矩。双力矩是在梁的横截面上弯矩的反对称分配。如果把箱梁简化为单根梁单元，由于不考虑双力矩，计算出来的受扭箱梁左右主梁的弯矩必然是相等的，仅剪力有差异。这显然不对。实际情况应当是：箱梁左右主梁的剪力、弯矩都会有差异。
:KA)4[#;W 3．剪力中心是弹性薄壁杆理论里的一个概念。形心是所有杆（梁）理论里都有的概念。弹性薄壁杆受到的垂直于杆轴线的外力，要等效平移到剪力中心轴线上，成为垂直于剪力中心轴线的力和绕剪力中心轴线的力矩，然后进行包括扭转在内的变形计算。如果杆单元不受垂直于杆轴线的外力，不考虑剪力中心也不妨碍进行正确的扭转内力计算。当采用更精确的力学模型时，形心、剪力中心的概念就不存在了。比如，用板壳单元模拟薄壁杆结构，就不存在整体截面的形心、剪力中心问题。 箱形梁桥这样的结构简化为弹性薄壁杆进行分析，就必须找出它横截面的形心、剪力中心，才能正确计算它在横向力作用下的弯曲、扭转。如果用梁格体系模拟箱形梁桥，形心、剪力中心的横桥向位置都不需要，但其竖向位置是必不可少要知道的。
B{D4.!a 4．梁桥的形心、剪力中心绝大多数情况下都不会重合。单箱多室箱形梁桥的剪力中心一般在形心下方，多箱多室箱形梁桥、开口截面梁桥的剪力中心一般在形心上方。只有实心板桥，如果它两侧没有悬臂板，或者悬臂板影响可以忽略，其剪力中心与形心是重合的。（笔者“曲线梁桥”程序可以用图形显示梁截面上的主梁自动划分、翼板有效宽度、形心和剪力中心位置）
WZPj?ou`G 5．用MIDAS的板壳单元模拟箱形梁桥，可以避开它不计算剪力中心的缺点，但在应用上有相当的困难。最主要的是各国土木结构设计规范有关混凝土梁强度计算的条文，都是以弯矩、剪力、扭矩和梁横截面尺寸为基础制定的。板壳单元输出的是梁横截面上若干点的应力，不能直接用于强度计算。对于位置固定的静力荷载，当然可以把若干点的应力换算成横截面上的内力。对于位置不固定的车辆荷载，理论上必须采用影响面方法求最大、最小效应。板壳单元输出的只能是各点的应力影响面。把各点的应力影响面重新合成为横截面的内力影响面，要另外附加大量工作，不知道MIDAS有无这项功能。因此，在计算梁式桥时，最实用的方法，应当是梁格法（MIDAS的文档资料里，也有这句话）。
W=w]`' 6．MIDAS是不能按照输入的桥梁结构形状自动建立梁格模型。网站/FAQ里有一资料“（54） 梁格法建模注意事项”，见附件2。 “在梁桥中会经常使用梁格法建立模型，因为不同的设计人员对横向联系的模拟(虚梁的设置)不尽相同，所以分析结果会略有差异。”不能自动建立梁格模型，不仅大大增加了用户的工作量，也使计算者主观因素对计算结果发生影响，不是“略有差异”，而是差别很大，不能说不是一个缺憾。注意事项的第4项，“虚拟的横向联系梁的刚度可按一字或二字形矩形截面计算”，横向联系梁的抗弯、抗扭刚度可以这么算，但箱形梁横向抗剪刚度计算就不是这样了，正确算法见E.C.汉勃利《桥梁上部构造性能》（郭文辉译, 人民交通出版社, 1982 年），孙广华的《曲线梁桥计算》（人民交通出版社, 1995 年）也引述了汉勃利的算法。
kWV k^ , 7．用MIDAS计算梁式桥，在什么情况下结果可能是可靠的？
dE_"|,: 情况一、预应力、非预应力直梁桥的静力计算
!<HMMf,-D 情况二、非预应力曲线梁桥的静力计算
w Wfj#IB;R 情况三、预应力、非预应力的直梁、曲线实心板梁桥的静力计算、动力计算
J -Lynvqm 还有一个前提：要按照前面6里的正确方法建立梁格模型。
2[uFAgf@ 说“可能是可靠的”，因为要得到正确结果，相关的因素太多太多，有程序编制方面的因素，也有使用者的人为因素。
9d- '%Q>+ 8．MIDAS和ANSYS、SUP2000及国内桥梁分析软件
+ xO3<u MIDAS不考虑剪力中心的问题，在ANSYS、SUP2000、“桥梁博士”也同样存在。在用它们建立梁格法时候，同样有工作量大和在什么情况下结果才可能是可靠的的问题。但这毕竟是小缺点，绝对无损于它们优秀软件的名誉。笔者查阅了ANSYS、SUP2000的杆单元资料，其中有不少横截面非中心对称的杆单元，这样的截面，形心和剪力中心是不重合的。为什么它们都不提剪力中心（MIDAS倒是声明了“形心和剪力中心重合”）？应当不是由于疏忽，可能有他们的思考。笔者还想钻这个牛角尖，争取找到答案。（后注：笔者已经查到，ANSYS有一种薄壁杆单元，可以计算单室截面的剪力中心，SUP2000、“桥梁博士”肯定是不考虑剪力中心的）
j/=Tj'S?D 9．相比之下，笔者的“曲线梁桥”程序显示出了突出的优点。（1）它计算梁截面的形心和剪力中心，形心位置考虑了翼板有效宽度，钢索水平分力对剪力中心取矩；（2）程序在输入了结构外部尺寸、材料标号、顺桥向划分的单元数以后，自动建立梁格力学模型，计算者不能干预；（3）程序在输入了最大车道数以后，即可以完成动态布载，输出各主梁各截面的组合内力、应力、钢筋数量、预应力钢筋应力；（4）除了顺桥向划分的单元数的多少是由计算者决定、对精度有一定影响外，整个计算基本不受计算者主观因素影响。如果和MIDAS“梁格法建模注意事项”比较一下，“曲线梁桥”程序的应用实在是够简单了。
v[{7\Hha Xq03o#-p+ 附件： MIDAS梁格法建模注意事项
oy5K*
} 在梁桥中会经常使用梁格法建立模型，因为不同的设计人员对横向联系的模拟(虚梁的设置)不尽相同，所以分析结果会略有差异。下面就一些注意事项供设计人员参考。
GD]epr%V 1. 将多室箱梁分割为梁格时，注意纵梁的中和轴位置应尽量一致。
,d 38TN 2. 每跨内的虚拟的横向联系梁数量不应过少(划分为1.5m左右一个在精度上应能满足要求)。
:{iH(ae; 3. 虚拟的横向联系梁之间尽量要设为铰接(可将纵梁之间的虚拟横梁分割为两个单元，将其中一个释放梁端约束)。
R6o D 4. 虚拟的横向联系梁的刚度可按一字或二字形矩形截面计算。
f_[d FKoX 5. 虚拟的横向联系梁的重量应设为零(可在截面刚度调整系数中调整)。
buN@O 7\ 6. 当虚拟的横向联系梁悬挑出边梁外时，应设置虚拟的边纵梁(为了准确地计算自振周期和分配荷载)，此时可将虚拟的边纵梁作为一个梁格进行划分。
N {Z+ 7. 定义移动荷载的车道时，应尽量选择按“横向联系梁”方法分布移动荷载，此时应将所有的横向联系梁定义为一个结构组，并在定义车道时选择该结构组。
kk6
!krZ 8. 定义车道时最好定义两次车道，一次按横向偏载定义，一次按横向中间向两边定义。定义移动荷载工况时可定义偏载和居中两个工况(荷载组合中会自动找到包络结果)。
~mXzQ be
p 9. 定义支座时尽量遵循一排支座中只约束其中一个支座在X, Y方向的自由度的原则(否则温度荷载结果会偏大)。另外，多支座时一般可不约束旋转自由度。
Nc{&AV8Y_v 10. 注意输入梁截面温度荷载时宽度B的取值为实际翼缘宽度(或腹板宽度之和)。
y&5
O) 11. 弯桥时应注意支座的约束方向(设置节点局部坐标系)。
M'< % d[ 以上注意事项仅供参
