高等桥梁结构论悬索桥课程设计.docVIP

高等桥梁结构理论 学 生 姓 名 学生学号 学院名称 指导老师 悬索桥成桥状态，对于主缆而言，所受荷载为沿弧长均布的主缆自重（包括缠丝及防护）及通过吊索传递的局部荷载，后一部分近似作为集中荷载处理，悬索桥的主缆受力可以简化为承受沿弧长均布的均布荷载加吊索处作用集中力的柔性索。 局部坐标系下，满足边界条件的索曲线方程为： （1） 其中 本桥中吊索将主跨主缆划分为段，每一段的主缆曲线方程由式(1)确定，除此之外还应满足力平衡方程和变形相容方程： ① ，主跨任意点通过给定点，即跨中满足矢跨比为1/10； ，由于主塔塔顶无偏位，边跨水平力等于主跨水平力； ，各局部坐标原点满足力平衡条件： （2） 建立迭代计算过程，通过二分法逼近求解，计算步骤如下： ①，计算主缆索力水平分量的迭代区间，； ②，假定左支座处竖向力的迭代区间，， 由式（2）得,得到和，以及： ， ③，由公式(2)得到下一段的等，依此循环可求得以及……； ④，求，若满足误差限则下步计算否则重复②~④，直至满足条件； ⑤，检验跨中垂度，若不满足则重复①~⑤， ⑥，计算各索段长度，无应力长度，各点的索曲线坐标。 计算得到主缆水平分力以后便可进一步算出悬链线索的形状长度和无应力长度。 (3) 无 （4） 根据以上理论，设未知数为成桥状态下的水平力H0以及鞍座处的竖向力V0，每根吊杆处作用Pi=(60+150)×1000/10=21000KN。建立以下的方程： ① ② 其中，令节点i左边的竖向力为V(i-1,2)，右边的竖向力为V(i,1)， 建立方程V(i,1)=V(i-1,2)-Pi 用matlab程序中的fsolve函数可以迭代解出得到成桥状态下的H0及V0。 缆的最大应力为塔顶支座处的应力，其计算公式为: ， 边跨缆的跨中垂度为: LS--边跨跨径长 1.1.2用matlab编辑子函数【fenduan.m】、【wuyingli2.m】、【xuanlianxian.m】 ①子函数【fenduan.m】 %---------------------------------------------------------------------- function y=fenduan(x) %初始数据输入 p1=53; L1=1000; C1=0; f1=100; Pi=21000; %每根吊杆处作用力 A=0.65; EA=1.2675e8; l=[50 100 100 100 100 100 100 100 100 100 50]; %索的长度向量 %定义三个数组 alpha=zeros(11,1); beta=zeros(11,1); c=zeros(11,1); V=zeros(11,2); H0=x(1); %设水平力为未知数 V0=x(2); %设塔顶竖向力为未知数 alpha(1)=asinh(V0/H0); beta(1)=p1*l(1)/2/H0; c(1)=H0/p1*(cosh(alpha(1))-cosh(2*beta(1)-alpha(1))); V(1,1)=V0; %第二个数1表示节点右边的竖向力 V(1,2)=-H0*sinh(2*beta(1)-alpha(1)); %第二个数2表示节点左边的竖向力 for i=1:10 V(i+1,1)=V(i,2)-Pi; alpha(i+1)=asinh(V(i+1,1)/H0); beta(i+1)=p1*l(i+1)/2/H0; c(i+1)=H0/p1*(cosh(alpha(i+1))-cosh(2*beta(i+1)-alpha(i+1))); V(i+1,2)=-H0*sinh(2*beta(i+1)-alpha(i+1)); %下一个节点左边的竖向力 end y(1)=sum(c(:)); %目标函数1，ci之和为0 y(2)=sum(c(1:5))+H0/p1*(cosh(alpha(6))-1)-100 %目标函数2，经过跨中定点 ②子函数【wuyingli2.m】 %