悬索桥梁工程中数值评价及计算.docVIP

悬索桥梁工程中数值评价及计算摘 要：对悬索桥从安装到最后投入使用过程中的典型受力形式进行了数值分析并为加劲梁的安全施工提供了理论依据。 关键词：数值评价与计算；悬索桥；有限元 有限元分析（FEA）是一种预测结构的偏移与其它应力影响的过程，有限元建模（FEM）将这个结构分割成单元网格以形成实际结构的模型，每个单元具有简单形态（如正方形或三角形）。这样有限元程序就有了可写出在刚度矩阵结构中控制方程方面的信息。每个单元上的未知量就是在节点上的位移，这个点就是单元元的连接点。有限元程序将这些单个单元的刚度矩阵组合起来以形成整个模型的总刚度矩阵，并给予已知力和边界条件来求解该刚度矩阵以得出未知位移，从节点上位移的变化就可以计算出每个单元中的应力。 结构或构件抵抗变形的能力，包括构件刚度和截面刚度，按受力状态不同可分为轴向刚度、弯曲刚度、剪变刚度和扭转刚度等。对于构件刚度，其值为施加于构件上的力（力矩）与它引起的线位移（角位移）之比。对于截面刚度，在弹性阶段，其值为材料弹性模量或剪变模量与截面面积或惯性矩的乘积。一般能满足F＝k△，F为作用力，△为位移，k即为刚度，所以刚度物理意义为单位位移时所产生的力。k可以是某些量的函数，即可为表达式。由F的不同，叫法不同。 另外就是我们要说的刚度叫线刚度，即单位长度上的刚度。 比如，我们在用反弯点法计算多层框架水平荷载作用下内力近似计算时。 计算柱的水平剪力时，剪力与柱层间水平位移△的关系为 V=(12ic/h2)△ 那么d=(12ic/h2)就叫柱的侧移刚度，表示柱上下两端相对有单位侧移时柱中产生的剪力。 其中ic表示柱的线刚度（即ic＝EI/h），h为楼层高，EI是柱的抗弯刚度（M=EI(1/p)，M为弯矩，（1/p）为曲率，也满足F＝k△形式）。 另外还可用D值法，即考虑了梁柱的刚度比变化，因为柱两端梁的刚度不同，即对柱的约束不同，那么它的反弯点，即M＝0的点会随之移动，那端强，反弯点离它越远。而且同层柱剪力分配时也是由柱的线刚度决定，因为同层位移一定，简单讲，由F＝k△，谁的刚度大，谁分得的剪力就大。 1 建模 主缆假定为小应变理想柔性索，其材料满足胡克定律，柏松效应忽略不计，基本公式为： 1=HL0EA0+Hqarsh(VH)-arshV-qL0H(1) h=L0EA0V-qL02+Hq1+(VH)2-1+V-qL0H(2) 吊索两端设有锚头，通过索夹与主缆相连，而锚头以外部分可看做是不可伸长的刚壁。因此，将吊索模拟为两端带刚壁的弹性悬索单元，推到公式如下所示： 1=HL0EA0+Hqarsh(VH)-arshV-qL0H+SSHH2+V2+SXHH2+(V-qL0)2(3) h=L0EA0V-qL02+Hq1+（VH）2-1+V-qL0H2+SSHH2+V2+SXHH2+(V-qL0)(4) 悬索的模拟采用LINK10单元，它具有应力刚化和大变形功能，在每个节点上存在三个平移自由度。LINK10单元带有惯性和阻尼效应，其固有的双线性刚性特性使其成为一个轴向仅受拉或受压的杆单元，广泛应用于悬索结构的静力分析中。加劲梁的模拟采用BEAM3单元，它是单轴承受拉力、压力和力矩作用的单元，每个节点具有X、Y线位移方向及绕Z轴旋转角度3个自由度。 对模型的加劲梁施加简支约束，悬索施加铰约束，和悬索单元相连的所有节点施加Z方向的约束，在悬索桥跨中部位施加Y方向的集中力，具体的加载及约束如下图： 2 找形 在主缆上假定或修正两个切点，则切点和Ip（下图所示）共同形成鞍座面，由此推出鞍座面的法相向量，具体如下所示： 对于施工阶段松弛悬索体系的找形分析.体系的几何不稳定性将引起结构的刚度奇异，使仅适用于稳定结构分析的常规方法失效.针对此问题，基于“力密度法”的基本思想，提出一种给定单元原长的松弛悬索体系的找形方法.从单索特性分析人手，直接推导节点坐标和外荷载之间关系的基本方程，避免了刚度矩阵的建立。 要达到“单元死”的效果，并非将“杀死”的单元从模型中删除，而是将其刚度矩阵乘以一个很小的因子，为了减少求解的方程数和避免病态条件，还要将“死”自由度约束住。此时死单元的载荷、质量、阻尼、比热及其他类似的效果都将被设定为零。“单元活”也并非是在模型中添加单元，而是在当前的载荷步中重新激活已经存在但在前面载荷步里被“杀死”的单元。 通过通用后处理查看悬索桥的跨中挠度、加劲桁架和缆索的最大轴力。比较跨中挠度是否接近零，且加劲桁架的最大轴力是否足够小。如不合适，则需要修改缆索的初应变重新进行计算，直到跨中挠度接近零且加劲桁架的最大轴力最小。由于悬索桥具有明显的非线性，每一次的受力分析，即安装加劲梁以前首先得准确求得索的