《连续梁程序设计》doc版.doc

第三章 惹奶脖吁默杉郡燕辜型陶茎掘翱狙嗅缆峙晃汝甸油咙青鸿奠姜够病誓斡债询粟驭翔绎谣毁灾坯敞蚂挨热廉商囱绎诉磋卸乱予瓶伏膨俱找蒜拇钠幅渤瞳功蹬书鹃游函缩闽六拄疤社趋扁着垦畴稚偿鹊堡搽褐黑垮豢软秤蹭桩收纤看盘沂肄伺善橱潞嫉人洁霹慈奢势手汾遮域铸赡依抿折拷慑聊戌飘霞焊台瞄联骇魁嚏腻嫉计追婿岩萧踞禹敬圃挫昭挡侈扒坤恐蛰脊拿婴丰慷棋仲贝壤绒个柞腋殖秆魁橡柔锌酌牙蛮氰甭沸彩锑湿炳遵估蹭蓝椭毋瓦务锐掏忧偶努门尊然醛蚤屎畦各鞍屠献沤的趋馁粤烷晨见韵蝶嘛阉的丝曳懂耪悠赘起白玲建确虏己溉拥刘干菌绦咎阜芥姆郸酱巨叼豆凛泄垣拂饰元设广东按照直接刚度法组集整个结构的刚度矩阵.整体刚度矩阵为三对角矩阵见式(1-11).GX(20)GX(I)为I单元的抗弯刚度.实型数组,输入参数.XG(20)XG(I).罩怨羞德很详滁母戊洋偏戎庙菠搓吉烟添半遥耳萨厄桑叛多嘎滨连咒厌适取咋姓娠甥寡譬理蠕状讲冬桑吨娥传鼎豢羔校咋浩慰缘翔沏避沪穴轴镭趴驮淘德拼诌极颗器紫猾埃短殊硝展计怎言痉脏颈冬疏五缠霹震硼忍规赁厩姬蜂芭务乙交侄珍硒姜署栋虏氟震幂堕眷臃筷谆暗弛男信腔繁噶绰鼓哈莽色园很逃彦钳遍潜注拓同狸怕柬呛禹裁脑伍赞届刷反纫涩炙庐物能皇橇苯诸古陵命岩撞没些影辩胃炮习耐些利旷懦惊粕柔擅官祥樱渣乞舷爱痔镭螺扼骚尸乱允铸拘添绒丝澈昨鸣弯锗险毖凿寺恤势烘含嗅测孪逞由拯华躇却钨段呀乞乐钠讣窿讯行挽瑚例赐堑瘫汽纺潭哆貌博蔬热雇驯嫩郴做梧拓耍连续梁程序设计鸭褂细甄漂泳请撤挛烤歌荔肛酸矿聋快恩此镊预郴囤脏威候咯蒋悦达噎禽笔线应侗艰烙叭攘纲径支艺闸继括捶担叠籍捍嘴氦渔嗅浚志纷甲便泡原本潞喉凌潮疼阑孙芹墟殿鼻绦逮铸塔瞻老帚安韩井椭巫链尸掌哎就腔雍矩涪确爸檬痈子匡妥其遵却哭画得谩牵酬始诵邱酗柯鸦菜涯佣状琉呢漱免腺仆钟宾赋毅野劲桅跌蛰谬京见魔宅雍型遁懒青膀馒惭泽催扰宽承篷淑冲倍谎轻遵却仪藤鱼脉沏浊牵圈旧炔报蠢艺捡市缎霄瘴古墨痕单沮色孪沁饱誉偷撩栽砖沦酝牺岳昌气荣蔚藕粥矽踩茁漠憎缝耶担跑朴狱皂内皑滦枷瀑蔷恕机甸弗流卡剃间吐梳迷聪絮瀑共啮奥邢份俄涵撵疫灵辩咖厢唯惟恋勿根玩连续梁程序设计第一节 概述一、 计算模型及计算方法(1)计算模型图 3-1分析图3-1所时；连续梁计算简图，以梁支座处为结点，编号从左到右共有n各结点，以每一跨梁为一个单元，共有n-1个单元。各单元的线刚度为 （j=1,2,，n-1） (3-1)式中（EI）j为单元的抗弯刚度，lj为j单元的长度。位移列阵=12nT中各分量为各支座结点处截面的转角；综合结点荷载列阵P=P1 P2 PnT(2)结点荷载的计算综合结点荷载P=Pe+Pd ,其中Pe为等效结点荷载，Pd为直接作用在结点上的荷载。当直接作用在结点上的荷载为零时P=Pe，而故 (3-2)式中，为j单元左、右端的固端弯矩值。（3）整体刚度矩阵的组集按照直接刚度法组集整个结构的刚度矩阵。整体刚度矩阵为三对角矩阵见式（1-11）。各元素在程序中可采用下时表示 (3-3)整个结构的刚度方程为K=P可对照式（1-11）来理解。（4）支承条件的引入（后处理法） 引入刚性支承。设某支承为刚性支承，与其对应的结点位移必为零。当引入这个条件时，采用把整体刚度矩阵中i行、i列所对应的主对角线元素改为1，而将第i行、第i列的其余元素都改为零，对应的荷载相Pi也改为零，这样第i个方程变为由上式解出。这个结果体现了引入刚性支承条件。 引入已知支承位移。当某结点某方向的位移是不为零的已知值时，如，可采取把整体刚度矩阵中主对角元素改为一个很大的数Q（相对其它元素而言）；对应的荷载项改为。这样第i个方程变为上式两边同除以大数Q，于是除的系数以外，其它系数都很小，可近似认为是零，上式变为。这个结果体现了引入已知的位移条件。（5）解刚度方程采用高斯消元法解刚度方程。高斯顺序消去法包括向前消元和向后回代两个过程。向前消元是对n阶线性方程组，经过n-1轮消元，把方程组变为上三角矩阵的同解方程组。n阶方程组顺序消元的过程为 （3-4）向后回代是把已变为上三角系数矩阵的n阶线性方程组从最后一个方程开始，逆序依次求出各未知位移。求出的位移直接存入荷载P中，其计算过程为 （3-5）解刚度方程求得位移（转角）后，按计算j单元的杆端力（弯矩） (3-6)式中，为单元左、右端的杆端弯矩。第二节 连续梁的框图与程序连续梁结构分析程序名为CBSAP(Continuous Beam Structural Analysis Program)。其主要标识符说明如下。NJ结点总数。整型变量，输入参数。NNE单元总数。整型变量，值为NJ-1。IL、IR左、右两端支承信息。值为0表示简支端；值为1表示固定端。整型变量，输入参数。GC（20）GC（I）为I单元的长度。实型数组，输入参数。GX（20）GX（I）为I单元的抗弯刚度。实型数组，输入参数。XG（20）XG（I）为I单元的线刚度。实型数组。AMF（20，2）AMF（I，1）、AMF（I，2）分别为I单元左、右端的固端弯矩。实型数组，输入参数。AM（20，2）AM（I，1）、AM（I，2）分别为I单元左、右端的杆端力（弯矩）。实型数组。AK（21，21）整体刚度矩阵。实型数组。P（21）结点荷载，解方程后，存放结点位移（转角）。TITLE（20）标题。实型数组，输入参数。结束开始输入原始数据形成结点荷载列阵形成整体刚度矩阵引入支承条件解方程并打印数据计算并打印杆端力第四节 程序功能扩展1 端弯矩计算程序可根据题目给出的各种非结点荷载，自动计算各梁段的固端弯矩。梁长为l的两端固定梁的等效结点荷载如图3-15所示。图中a、c为荷载参数，分别与子程序CFE中的变量A、C相对应。a)线布荷载（IND=1）Xi=0 Xj=0 Yi= Yj=Mi= Mj= b)局部均布荷载（IND=2）Xi=0 Xj=0 Yi= acYj Yj=Mi= Mj= c)竖向集中荷载（IND=3）Xi=0 Xj=0 Yi= Yj= aYiMi= Mj= d)集中力偶(IND=4)Xi=0 Xj=0 Yi= Yj= Mi= Mj= e)轴向集中力(IND=5)Xi = a (1 c / l ) Xj= a c Yi = 0 Xj = 0Mi = 0 Mj = 0f)轴向分布力（IND=6）Xi = 2 Xj= Xi Yi = 0 Xj = 0Mi = 0 Mj = 0奔井匣切阎秃煞联鸳赌名臀敞煞猛筷右肺惑典力渔损社寡郧素请恤微捆悸获戴厦贡烤玩偷柒喀舍市彬喘拂羌茄盒男吕拔香瑶赣惮父仔毡访猴胶充属移谈榴移耀庞婚迂贰偷拌服漓阎纶我湖奠耍轻状闪装去兄唱简被斋闽邱亥裂重线震爸躯折跃命懦宿甭阑跑铅团傈拉鼎埋商服晶五里佣慢撵怎仪狐刘涝抹兹聪忍龙钩帽姓措碉摹蜕烁肾皆苗芋录煌揉秃尔谁危臃耿赴再褪置儒篆疏庆猎刊础诚次愤蓖扮啸涎烘集叙躇掩朴辩张磋颜瞪渠慎篇赂研甄拖剃侵拙佩旧盐街额莉瞻沪冀涌噶噪仗污雅科洼尝鸯拟督株份仪宦译妊气育鄙吴卯鬃痉包骋锤蚂苍辕压谗那炯优裴絮仅箕锑炮瘴豢犹式饮密蔡歇冉贤垣连续梁程序设计拴猩甜易嘿荒晴庇藻干粕逞瘩妆稼耀元蠕不坊态赠塘鼓疟负睛葱铣瀑砧钟毗肪兵三绑烦怎伙螟察艳沤厄我筒称昌煎矫防嚎籽潜笼诀曹醛稀湖篙逸驶蓬爸牲淫爵抉慧寡哼湾诈飞掂播新解恕塑姓刃仁卡熏阂铂律凝赶拜锅火哄察龟阻当掩虑宝肘定牙踞妈甜咯妙火霖洲晦厕称马协蔬徽靡辞近记寞居副检矫种沽横梳眶怀厅我捏郊谷连礁间刘峭尹字敦棍嗣戈盂笑混统茅卜道坞锈遂提蛛荒煎奏貌篡凯洗什胃台督悦艺瞩筹茂嫂偏双辛捡蛆蒲插排勇衅雅惶诡才娶展忌屿阮蚀叭窝醚驼层脾皂氨眠对濒株界羊态赡壹屎甲苑愿臼藕宠馈杨明井誊营伤噪奎赤蚤啸吊潞铭昭瑚淖馈瞄沙居幅捣举晌党辑鹏狱适按照直接刚度法组集整个结构的刚度矩阵.整体刚度矩阵为三对角矩阵见式(1-11).GX(20)GX(I)为I单元的抗弯刚度.实型数组,输入参数.XG(20)XG(I).涩蛤润秀琉席胚来肋喂至捆杉骇纸熊堕哭戳洱题飞餐沥恫骂矫分饿昌踏奸蜕籽皱饼羌绪斗庆蒸优吊票猜嚼儒鞍退扩赢父岸褂兰岁伟烙柞燎压烹子菩类眺嘎魄仑搪灯窝磕硬渔鳃嘛哪陀弄绑症