ANSYS桥梁建模经验1

1.从连续刚构桥纵向来看，保证各个截面具有几何拓扑一致性，对于箱形截面来说，只有空间倒角部位可能出现截面拓扑不一致，在建模时人为将其定义为几何拓扑一致。对于横隔板建模采用后补法来实现。
2.从横截面来看，一般单箱室连续刚构桥具有外轮廓的几何拓扑，因此将其作为截面的拓扑。对于和其拓扑一致或者可以调整为该拓扑结构的连续刚构都能分析。对于箱内倒角为双折线的连续刚构只能近似模拟。具体建模时，我们只需输入外轮廓尺寸，截面纵向位置，APDL命令流即可自动计算其他关键点位置，用循环语句生成实体模型。
442钢筋与混凝土的连接如果采用线单元初应力法模拟预应力有限元模型中钢筋和混凝土的连接主要分为三种即分离式整体式和组合式4252分离式模型把混凝土和钢筋作为不同的单元来处理即混凝土和钢筋各自被划分为足够小的单元两者的刚度矩阵是分开来求解的考虑到钢筋是一种细长的材料通常可以忽略其横向抗剪强度因此可以将钢筋作为线单元处理
3.对于每一根预应力束，定义张拉端为局部坐标原点，通过指定导线点局部坐标，和局部坐标在整体坐标中的位置进行定位。对于双向张拉的预应力束将其分为两根来建模，分割点位于0.5L处（L预应力束总长）。这样做看上去繁琐，但却解决了所有预应力束都可以通过一个宏命令完成预应力束损失的计算，单元的生成，初应变的赋值功能，没有任何局限性，为参数化建模提供了保证。参数的输入只有各导线点局部坐标，局部坐标在整体坐标中坐标，张拉控制力，波纹管类型等计算损失的参数。
1）输入参量： ； ； ； ；
2）计算建模控制点参量 ： ；
3）计算上层配筋层纵向配筋率 ：
4）返回计算数值 和 进行建模和单元属性赋值。
图4-1配筋率计算示意图
Fig.4-1Schematic plan of ratio of reinforcement calculation
4
从上一节我们可以看出，本文采用的连续刚构实体有限元模型，从真实反应结构的力学特性方面做的比较细致，采用APDL建模技术后，可操作性也大大提高，但尽管如此，模型的规模还是比较庞大的，因此必须考虑计算效率的问题。
4
APDL也就是ANSYS参数化设计语言，是一种类似FORTRAN的解释性语言。它具有一般程序语言所具有的功能，如参数、宏、标量、向量及矩阵运算、分支、循环、重复以及访问ANSYS有限元数据库等功能。是优化设计的基础，也是参数化设计的最高技术[52]。APDL命令流通常具有以下优点：
1．模型文件小，不同版本间通用性强。
1.为了大量降低单元数，纵向仅对一个T进行模拟（即：只模拟一个边跨和半个中跨）；在跨中处施加对称约束；横向上仅模拟半个箱梁，在对称截面处施加对称约束。也就是采用四分之一模型。单元总数控制在10万以内，单次求解时间小于15分钟。这样做只能模拟对称荷载，近似模拟部分非对称荷载效应，当然扭转效应无法模拟，这样做对于本文的分析研究基本满足，如果有必要，也可以采用单元镜像功能对APDL命令流稍作土和钢筋作为不同的单元来处理，即混凝土和钢筋各自被划分为足够小的单元，两者的刚度矩阵是分开来求解的，考虑到钢筋是一种细长的材料，通常可以忽略其横向抗剪强度，因此可以将钢筋作为线单元处理。钢筋和混凝土之间可以插入粘结单元来模拟钢筋和混凝土之间的粘结和滑移。一般钢筋混凝土是存在裂缝的，而开裂必然导致钢筋和混凝土变形的不协调，也就是说要发生粘结的失效与滑移，所以此种模型的应用最为广泛。
1.施加预加力
如果是单纯施加一定的预加应力,通常有等效荷载法、线单元初应力法。前者的弊端是不能模拟后期预应力损失，以及被动受力特性。
2.预应力张拉损失
预应力张拉损失包括管道摩擦损失；锚具变形、钢筋回缩损失；混凝土加热养护时，预应力钢筋与台座之间的温差损失；混凝土弹性压缩损失。这些损失有的可以通过事先计算然后采用等效荷载法、线单元初应力法对其进行模拟，然而弹性压缩损失按规范计算难以对其真实模拟，必须通过施工顺序的模拟来实现其真实的模型预应力赋予值。
值得一提的是补偿弹性压缩损失。关于弹性压缩损失，有两点需要澄清一下：a.弹性压缩损失有时指实际工程的弹性压缩损失，有时指有限元模型的弹性压缩损失。前者是指先张法预应力混凝土放张后产生的弹性压缩损失和后张法预应力混凝土分批张拉引起先张拉的预应力束，在其他预应力束张拉时可能引起的弹性压缩损失；后者是由于有限元采用线单元初应力法模拟预应力所带来的；本文所说的补偿弹性压缩损失是针对于后者而言，也是有限元正确模拟预应力分批张拉的必经之路。b.对于先张法和后张法预应力的工程弹性压缩损失，有限元模拟和手算是截然相反的。先张法工程弹性压缩损失手算是在控制力乘以一个折减系数，而有限元模拟则是输入与控制力等价的初始应变，通过预应力和实体元变形协调真实得模拟混凝土的受力状态，也就是说不需要、也应该再调整初始应变。后张法预应力，对于张拉过程中产生的弹性压缩，手算不折减，对于分批张拉引起的弹性损失进行折减；而有限元模拟恰恰相反，张拉过程的弹性压缩要通过迭代计算调整初应变来真实模拟张拉，而对于分批张拉的弹性压缩损失，通过预应力和实体元变形协调自行损失，不得再次调整。图4-2是后张法悬臂施工预应力初应变调整流程图：
4.3.3
整个模型共采用三种单元、可定义多种材料属性。其中梁体钢筋混凝土采用SOLID65实体元，支座以及双薄壁墩采用SOLID45实体元，预应力钢筋采用LINK8杆单元。其中实体单元的材料特性，普通钢筋和预应力钢筋的材料特性，均作为参数化变量可以调整；实体元的实常数（配筋率、配筋角度）、以及预应力单元的实常数（初应变、面积）由APDL命令流根据所输入的控制参数（普通钢筋直径、间距；预应力钢筋直径、锚下控制力、材料特性）自行计算。
3.预应力在整个预应力钢筋中的分布
实际的预应力值由于张拉预应力钢筋与混凝土的摩擦导致预应力值的分布不均一。如果采用线单元初应力法，则必须采用多单元赋予多个初应变来模拟，当然，带来了模型的复杂化。
4.混凝土收缩徐变带来后期预应力损失
混凝土的收缩徐变可以采用调整弹性模量法、和定义蠕变材料非线性来实现。其带来的预应力损失效应最好用线单元初应力法模拟，程序进行自动应变协调计算，不需要人工干预，或者手算。
具体建模采用几何建模法和单元建模法的组合建模，钢筋混凝土采用建立几何模型后扫略成SOLID65实体单元，预应力束采用直接连接钢筋混凝土单元上的节点形成LINK8单元。这样建模的前提是实体单元的建立必须在预应力束所在位置生成节点。通过ANSYS的体切割命令,可以在实体几何模型上切出预应力筋的节点位置,这对于几何形状不太复杂的模型是比较方便的，如果模型较为复杂时，（如实体不太规律，含有曲线预应力钢筋等），切割法会产生形状不符合求解要求的实体单元，即坏单元。同时切割的APDL编程不太容易控制。本文采用预应力筋处预留关键点，以及适当加密腹板单元的竖向划分数来保证模型的预应力束位置与实际工程吻合（见图4-2，4-4）。这样做一方面增加了分析命令流程序的通用性和APDL编程的可操作性；另一方面有利于模型的生成，即保证了单元的正确性，又大大降低了单元的数目（见图4-4，4-5），这对于求解是非常有利的，但带来了大量的模型建立前期工作。要使模型具有一定的适用性和可调性，必须进行容错判断，需要进行合理的几何数学计算和严密的逻辑判断。
4.4
预应力混凝土的建模是整个模型建立的关键所在，下面对其进行单独论述。预应力混凝土模型的建立有两个重点：1.预应力效果的模拟；2.钢筋（包括普通钢筋和预应力钢筋）与混凝土的连接。
4.4.1
所谓预应力效果，是指预应力的力学行为，包括对混凝土施加预加力、预应力张拉损失、预应力在整个预应力钢筋中的分布、混凝土收缩徐变带来后期预应力损失、预应力的被动受力效应、预应力钢筋与混凝土之间的滑移。当然如果将这些预应力效果在一个模型里同时考虑是很不现实的。以下这几种预应力力学行为的模拟做一简单分析。
当不考虑混凝土和钢筋二者之间的滑移时，三种模型都可以。
本文建立的模型既考虑普通钢筋、又考虑预应力筋，因此采用预应力钢筋与混凝土的连接采用整体式模型，不考虑预应力钢筋与混凝土之间的滑移；对于普通钢筋与混凝土的连接采用分层式带筋预应力混凝土单元SOLID65。这样建立的模型既可以分析钢筋混凝土梁桥，也可以分析预应力混凝土梁桥。
5.预应力的被动受力效应以及预应力钢筋与混凝土之间的滑移的模拟
这种效应的模拟必须采用线单元初应力法模拟，钢筋与混凝土之间的滑移还必须引进接触单元。
本文所建立的模型采用线单元初应力法，分段模拟预应力值，不考虑钢筋与混凝土之间的滑移和混凝土收缩徐变带来后期预应力损失。管道摩擦损失、锚具变形、钢筋回缩损失按照现行规范编写预应力损失计算程序，在预应力钢筋单元生成的同时，自动计算相应的初应变赋予单元，并按顺序将实常数编号、相应初应变返回预定的数组，便于下一步根据施工阶段的模拟补偿弹性压缩损失。
2.整体式模型
将钢筋分布于整个单元中，假定混凝土和钢筋粘结很好，并把单元视为连续均匀材料，与分离式不同的是，它求出的是综合了混凝土和钢筋单元的刚度矩阵；与组合式不同之点在于它不是分别求出混凝土与钢筋对单元刚度的贡献然后再组合，而是一次求得综合的刚度矩阵。
3.组合式模型
组合式模型又分为两种：一种是分层组合式，在横截面上分成许多混凝土层和若干钢筋层，并对截面的应变作出某些假设，这种组合方式在钢筋混凝土板、壳结构中的应用较广；另一种组合方法是采用带钢筋膜的等参单元。
—i批预应力单元的控制应力向量；
—允许误差。
通过以上预应力单元的控制模拟，前三项预应力效果都得到了很好的模拟，通过APDL的功能，使得预应力的模拟得到了很好的实现，具有可推广性。
4.4.2
如果采用线单元初应力法模拟预应力，有限元模型中钢筋和混凝土的连接主要分为三种,即分离式、整体式和组合式[42~52]。
图4-2初应变调整流程图
Fig. 4-2 Flow diagram of initial strain adjustment
其中： —施工阶段总数；
—初应变数组；
—1到i施工阶段实体单元；
—1到i施工阶段预应力单元；
—1到i预应力单元的初应变向量；