清华大学结构力学下-连续梁内力计算程序.pdfVIP

清华大学结构力学下-连续梁内力计算程序--第1页

清华大学土木工程系连续梁内力计算程序结23李会平2011012208

连续梁内力计算程序

一、连续梁计算图示：

算例1

算例2

算例3

算例4

清华大学结构力学下-连续梁内力计算程序--第1页

清华大学结构力学下-连续梁内力计算程序--第2页

清华大学土木工程系连续梁内力计算程序结23李会平2011012208

二、Matlab程序代码如下所示：

clc

clearall

%输入原始数据

NE=input(输入单元总数NE:);

BL=input(输入单元杆长BL列向量（NE×1）:);

EI=input(输入单元抗弯刚度EI列向量（NE×1）:);

JD=input(输入单元定位向量JD（NE×2）:);

FM=input(输入单元固端弯矩FM（NE×2）:);

%节点转角未知量总数

N=max(max(JD));

%输入直接结点力矩

PJ=input(输入N维直接结点力矩列向量:);

%输出元数据

disp(单元总数),NE

disp(单元杆长),BL

disp(单元抗弯刚度),EI

disp(单元定位向量),JD

disp(单元固端弯矩),FM

disp(直接结点力矩),PJ

%EI改存线刚度

EI=EI./BL;

%形成结点荷载向量

P=zeros(N,1);

清华大学结构力学下-连续梁内力计算程序--第2页

清华大学结构力学下-连续梁内力计算程序--第3页

清华大学土木工程系连续梁内力计算程序结23李会平2011012208

fori=1:NE

forj=1:2

m=JD(i,j);

ifm~=0

P(m)=P(m)-FM(i,j);%结构等效结点荷载向量

end

end

end

F=P+PJ;

%结构整体刚度矩阵的集成

K=zeros(N,N);

fori=1:NE

m=JD(i,1);

n=JD(i,2);

ifm~=0n~=0

K(m,m)=K(m,m)+4*EI(i);

K(n,n)=K(n,n)+4*EI(i);

K(m,n)=K(m,n)+2*EI(i);

K(n,m)=K(n,m)+2*EI(i);

elseifm==0n~=0

K(n,n)=K(n,n)+4*EI(i);