2012结构力学2大作业1.docVIP

PAGE \* MERGEFORMAT 1 2012结构力学2大作业1 编制原理： 本程序的编制原理是基于矩阵位移法的基本思想，首先把斜拉桥整个复杂的结构体系拆分为单个简单的单元，以结点位移为基本未知量，列出各个单元的单元刚度矩阵，再把单元刚度矩阵组装为整体刚度矩阵，同时根据各个结点所受合力算出每个结点的位移，再根据单元刚度方程算出各个单元的弯矩。 程序说明： 27左侧主塔左侧7根拉杆，右侧5跟拉杆，右侧主塔两侧各6根拉杆，把整个体系拆分为51个单元，共30个结点，，采用后处理法，共90个位移编号，对于梁两端结点按照铰支座处理，列出总体刚度矩阵后，用置0置1法处理总刚度矩阵。根据公式算出，同时在根据算出弯矩。整体单元编号和结点编号如图： 27 L L 此处无搭接 ? ? ? ? K=5 eq \o\ac(○,26) 1 1 n=6 2H/3 H/3 L L L 2L L L L eq \o\ac(○,1) eq \o\ac(○,25) eq \o\ac(○,49) 1 26 28 29 30 位移编号为： 结点编号 位移编号 1 （1,2,3） 2 (1,2,3) …… …… 30 (88,89,90) 结点、单元、位移编号结构离散化算法流程： 结点、单元、位移编号 结构离散化 将梁、拉杆、主塔的单元刚度矩阵转换为整体坐标下的单元刚度矩阵单元分析，建立拉杆、梁、主塔的单元刚度矩阵 将梁、拉杆、主塔的单元刚度矩阵转换为整体坐标下的单元刚度矩阵 单元分析，建立拉杆、梁、主塔的单元刚度矩阵 确定结构中各个结点荷载的矩阵对单元刚度矩阵进行组装，形成整体刚度矩阵 确定结构中各个结点荷载的矩阵 对单元刚度矩阵进行组装，形成整体刚度矩阵 利用单元刚度矩阵和所求位移求解内力建立刚度矩阵方程，求解位移 利用单元刚度矩阵和所求位移求解内力 建立刚度矩阵方程，求解位移 源代码： function F=zy(H,L) H=input(塔高：H=); L=input(节点间距：L=); EIc=input(主塔抗弯刚度EIc=); EAc=input(主塔抗压刚度EAc=); EIb=input(主梁抗弯刚度EIb=); EAb=input(主梁抗压刚度EAb=); EA0=input(拉杆抗压刚度EA0=); q=input(均布荷载q=); %输入数据 if H=0||L=0||EIc=0||EAc=0||EIb=0||EAb=0 clear; fprintf(输入错误，请重新输入！\n); hitwh101210232; else k1 = [ EAb/L 0 0 -EAb/L 0 0 0 12*EIb/L^3 6*EIb/L^2 0 -12*EIb/L^3 6*EIb/L^2 0 6*EIb/L^2 4*EIb/L 0 -6*EIb/L^2 2*EIb/L -EAb/L 0 0 EAb/L 0 0 0 -12*EIb/L^3 -6*EIb/L^2 0 12*EIb/L^3 -6*EIb/L^2 0 6*EIb/L^2 2*EIb/L 0 -6*EIb/L^2 4*EIb/L] ;%长L的梁单元刚度矩阵 k2= [ 1/2*EAb/L 0 0 -1/2*EAb/L 0 0 0 3/2*EIb/L^3 3/2*EIb/L^2 0 -3/2*EIb/L^3 3/2*EIb/L^2 0 3/2*EIb/L^2 2*EIb/L 0 -3/2*EIb/L^2 EIb/L -1/2*EAb/L 0 0 1/2*EAb/L 0 0 0 -3/2*EIb/L^3 -3/2*EIb/L^2 0 3/2*EIb/L^3 -3/2*EIb/L^2 0 3/2*EIb/L^2 EIb/L 0 -3/2*EIb/L^2 2*EIb/L];%长2L梁的单元刚度矩阵