matlab 矩阵位移法编程 结构力学

矩阵位移法编程大作业

（091210211）

一、编制原理

本程序的原理是基于结构力学矩阵位移法原理，以结构结点位移作基本未知量，将要分析的结构拆成已知节点力—结点力位移关系的单跨梁集合，通过强令结构发生待定的基本未知位移，在各个单跨梁受力分析结果的基础上通过保证结构平衡建立位移法的线性方程组，从而求得基本未知量。

二、程序说明

本程序是计算3层11跨框架右侧结点的位移和弯矩的程序，编译过程是按照矩阵位移法的先处理法进行的。首先将结构杆件的交汇点作为结点，共有36个结点和108个位移编号，然后根据梁、柱、斜杆的不同分别建立单元刚度矩阵，然后转换为整体坐标系下的刚度矩阵，然后将所有杆件的单元刚度矩阵整合成为总体刚度矩阵，在进行整合时连续运用for函数，最终形成108阶的总体刚度矩阵。然后通过对荷载的分析自己确定出荷载矩阵，直接写进程序。这样就可以把36个结点的108个位移求得，然后再利用各个单元的单元刚度矩阵和所得的位移求得单元杆件的内力。

离散化编号如下图：

三、算法流程

四、源代码

%结构力学大作业 3层11跨框架矩阵位移法编程王贝 091210211

h=input('输入单层高h：');

L=input('输入单跨度L：');

EIc=input('输入柱子的抗弯刚度EIc：');

EAc=input('输入柱子的抗压刚度EAc：');

EIb=input('输入梁的抗弯刚度EIb：');

EAb=input('输入梁的抗压刚度EAb：');

EIo=input('输入斜杆的抗弯刚度EIo：');

EAo=input('输入斜杆的抗压刚度EAo：');

q=input('输入侧向均布荷载集度q:');

T1=[1,0,0,0,0,0;

0,1,0,0,0,0;

0,0,1,0,0,0;

0,0,0,1,0,0;

0,0,0,0,1,0;

0,0,0,0,0,1];%角度为0°的转换矩阵

T2=[0,1,0,0,0,0;

-1,0,0,0,0,0;

0,0,1,0,0,0;

0,0,0,0,1,0;

0,0,0,-1,0,0;

0,0,0,0,0,1];%角度为90°的转换矩阵

x=atan(h/L);

T=[cos(x),sin(x),0,0,0,0;

-sin(x),cos(x),0,0,0,0;

0,0,1,0,0,0;

0,0,0,cos(x),sin(x),0;

0,0,0,-sin(x),cos(x),0;

0,0,0,0,0,1];%斜杆的转换矩阵

T3=T;

%梁的单元刚度矩阵

kb0=[EAb/L 0 0 -EAb/L 0 0;

0 12*EIb/(L*L*L) 6*EIb/(L*L) 0 -12*EIb/(L*L*L) 6*EIb/(L*L);

0 6*EIb/(L*L) 4*EIb/L 0 -6*EIb/(L*L) 2*EIb/L;

-EAb/L 0 0 EAb/L 0 0;

0 -12*EIb/(L*L*L) -6*EIb/(L*L) 0 12*EIb/(L*L*L) -6*EIb/(L*L);

0 6*EIb/(L*L) 2*EIb/L 0 -6*EIb/(L*L) 4*EIb/L];

%柱子的单元刚度矩阵

kc0=[EAc/h 0 0 -EAc/h 0 0;

0 12*EIc/(h*h*h) 6*EIc/(h*h) 0 -12*EIc/(h*h*h) 6*EIc/(h*h);

0 6*EIc/(h*h) 4*EIc/h 0 -6*EIc/(h*h) 2*EIc/h;

-EAc/h 0 0 EAc/h 0 0;

0 -12*EIc/(h*h*h) -6*EIc/(h*h) 0 12*EIc/(h*h*h) -6*EIc/(h*h);

0 6*EIc/(h*h) 2*EIc/h 0 -6*EIc/(h*h) 4*EIc/h;];

%斜杆的单元刚度矩阵

H=sqrt(h*h+L*L);

ko0=[EAo/H 0 0 -EAo/H 0 0;

0 12*EIo/(H*H*H) 6*EIo/(H*H) 0 -12*EIo/(H*H*H) 6*EIo/(H*H);

0 6*EIo/(H*H) 4*EIo/H 0 -6*EIo/(H*H) 2*EIo/H;

-EAo/H 0 0 EAo/H 0 0;

0 -12*EIo/(H*H*H) -6*EIo/(H*H) 0 12*EIo/(H*H*H) -6*EIo/(H*H);

0 6*EIo/(H*H) 2*EIo/H 0 -6*EIo/(H*H) 4*EIo/H];

kb=T1'*kb0*T1;%总体坐标下梁的单元刚度矩阵

kc=T2'*kc0*T2;%总体坐标下柱子的单元刚度矩阵

ko=T3'*ko0*T3;%总体坐标斜杆的单元刚度矩阵

X=zeros(108,108);Y=zeros(108,108);Z=zeros(108,108);%定义108阶0矩阵K1=zeros(108,108);K2=zeros(108,108);K3=zeros(108,108);

K4=zeros(108,108);K5=zeros(108,108);K6=zeros(108,108);

K7=zeros(108,108);K8=zeros(108,108);K9=zeros(108,108);

%把梁杆单元矩阵整合到总体刚度矩阵的循环语句

for ii=1:11

X(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb;

K1=K1+X;X=zeros(108,108);

end

for ii=13:23

Y(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb;

K1=K1+Y;Y=zeros(108,108);

end

for ii=25:35

Z(3*ii-2:3*ii+3,3*ii-2:3*ii+3)=kb;

K1=K1+Z;Z=zeros(108,108);

end

%把柱杆单元矩阵整合到总体刚度矩阵的循环语句

for jj=1:36

K2(3*jj-2:3*jj,3*jj-2:3*jj)=kc(4:6,4:6);

end

for jj=1:24

K3(3*jj-2:3*jj,3*jj-2:3*jj)=kc(1:3,1:3);

end

for jj=1:24

K4(3*jj-2:3*jj,3*jj+34:3*jj+36)=kc(1:3,4:6);

end

for jj=1:24

K5(3*jj+34:3*jj+36,3*jj-2:3*jj)=kc(4:6,1:3);

end

%把斜杆杆单元矩阵整合到总体刚度矩阵的循环语句

for gg=3:12:27

K6(3*gg-2:3*gg,3*gg-2:3*gg)=ko(4:6,4:6);

end

for gg=2:12:14

K7(3*gg-2:3*gg,3*gg-2:3*gg)=ko(1:3,1:3);

end

for gg=2:12:14

K8(3*gg-2:3*gg,3*gg+37:3*gg+39)=ko(1:3,4:6);

end

for gg=2:12:14

K9(3*gg+37:3*gg+39,3*gg-2:3*gg)=ko(4:6,1:3);

end

K=K1+K2+K3+K4+K5+K6+K7+K8+K9;%总体刚度矩阵

P=zeros(108,1);P(1,1)=h*q;P(37,1)=h*q;P(73,1)=h*q/2;P(75,1)=q*h*h/12; A=K\P;%结构位移

B1=kb*A(103:108,1);

B2=kb*A(67:72,1);

B3=kb*A(31:36,1);

D1=zeros(6,1);D1(1:3,1)=A(70:72,1);D1(4:6,1)=A(106:108,1);

D2=zeros(6,1);D2(1:3,1)=A(34:36,1);D2(4:6,1)=A(70:72,1);

C1=kc*D1;C2=kc*D2;C3=kc(4:6,4:6)*A(34:36,1);

M1(1,1)=B1(6,1);M2(1,1)=C1(6,1);M3(1,1)=B2(6,1);M4(1,1)=C2(6,1);

M5(1,1)=C1(3,1);M6(1,1)=B3(6,1);M7(1,1)=C3(3,1);M8(1,1)=C2(3,1);

for i=1:3

m=36*i-2;

fprintf('第%d层最右侧节点的位移是%d\n',i,A(m,1))

end

fprintf('第1层最右侧节点左侧杆的弯矩是%f\n',M6)

fprintf('第1层最右侧节点下侧杆的弯矩是%f\n',M7)

fprintf('第1层最右侧节点上侧杆的弯矩是%f\n',M8)

fprintf('第2层最右侧节点左侧杆的弯矩是%f\n',M3)

fprintf('第2层最右侧节点下侧杆的弯矩是%f\n',M4)

fprintf('第2层最右侧节点上侧杆的弯矩是%f\n',M5)

fprintf('第3层最右侧节点左侧杆的弯矩是%f\n',M1)

fprintf('第3层最右侧节点下侧杆的弯矩是%f\n',M2)

五、试算算例

输入数据：输入单层高h：1

输入单跨度L：1

输入柱子的抗弯刚度EIc：1

输入柱子的抗压刚度EAc：1

输入梁的抗弯刚度EIb：1

输入梁的抗压刚度EAb：1

输入斜杆的抗弯刚度EIo：1

输入斜杆的抗压刚度EAo：1

输入侧向均布荷载集度q:1

计算结果：第1层最右侧节点的位移是-6.219850e-003 第2层最右侧节点的位移是-2.152659e-002

第3层最右侧节点的位移是-4.131873e-002

第1层最右侧节点左侧杆的弯矩是0.000729

第1层最右侧节点下侧杆的弯矩是0.008642

第1层最右侧节点上侧杆的弯矩是-0.009371

第2层最右侧节点左侧杆的弯矩是0.000074

第2层最右侧节点下侧杆的弯矩是0.004158

第2层最右侧节点上侧杆的弯矩是-0.004232

第3层最右侧节点左侧杆的弯矩是-0.000761

第3层最右侧节点下侧杆的弯矩是0.000761

矩阵位移法编程大作业王贝091210211

MATLAB在结构力学分析中的应用

MATLAB在结构力学分析中的应用 摘要:传统的手算方法解超静定结构工作量繁重，有时甚至是不可能，运用结构有限元编程的一般方法，通过两个实例的对照，展示MATLAB在结构力学分析中的应用，MATLAB具有高性能，方法具有普遍的适用性,实现弯矩图自动绘制。 关键词: MATLAB结构有限元弯矩图 Abstract：While using the traditional manual method to resolve complex statically indeterminate structures, it is heavy workloads, sometimes even impossible，using finite element programming of the general method, Based on two examples, This paper introduces a method of application of MATLAB in structure mechanics, MATLAB has the advantages of high performance, it can be applied to many kinds of structures, realization of automatic drawing bending moment diagram. Key words: MATLAB; Finite element; Bend moment diagram 引言 结构力学[3]中，常利用传统的力法与位移法求解超静定结构，力法是几何问题，位移法把复杂的几何图乘转化为代数运算，但它们基本未知量很多时，系数构成的矩阵计算巨大，两者都不能满足科研工作者的需要。应用MATLAB 软件丰富可靠的矩阵运算、数据处理、图形绘制等便利工具，可使得计算和图象一体化。对于结构力学计算是十分有利的工具。 1基本方法 MATLAB结构有限元编程的基本思路是先分后合，即将结构分成各个单元和节点，桁架与刚架已经离散化，对于连续系统这一步极其重要，然后进行单元分析，集成整体刚度矩阵，引入边界条件，最后解方程。在求解平面桁架结构，虽然结构简单，用手算可得各杆件的轴力，但重复的过程太多，现在使用MATLAB语言来编制有限元位移法的程序时，则编程的难度明显降低，对有限元位移法的概念的理解更加深入，编程所需时间也大大减少。 图1为一平面桁架，各杆E=70GPaA=0.004，试用矩阵位移法求解各杆轴力 图1 解：平面桁架元是既有局部坐标又有总体坐标的二维有限元；对各结点

有限元MATLABWord版

MATLAB报告 Matlab程序求解简要过程如下： （1）求取单元节点位移提取矩阵T 单元节点位移提取矩阵T本质上是置换矩阵群中的一个，结果可将任意杂乱的节点顺序置换成统一的顺序。另一方面其作用是对单元刚度矩阵进行“升维操作”，将单元刚度矩阵统筹到整体刚度矩阵上来，便于对总体节点位移矩阵和支座反力进行求取。 本程序分析过程中对单元1的节点提取是按顺序编号1-2-3，对单元2的节点提取是按顺序编号2-3-4。单元1的节点位移提取矩阵如下：

单元2的节点位移提取矩阵如下： （2）求取单元几何矩阵B 单元1的节点按编号顺序1-2-3分别进行对几何函数矩阵或算子矩阵的bi逆时针操作，对ci顺时针操作；单元2的节点按编号顺序2-3-4分别进行对几何函数矩阵的bi顺时针操作，对ci逆时针操作.在MATLAB程序中通过mod()取模函数来达到对节点的顺时针或逆时针循环操作。 单元1的几何矩阵如下：

单元2的几何矩阵如下： （3）求取应力矩阵S 单元应力矩阵满足S=D*B,其中D为弹性矩阵，B为单元几何矩阵 各单元的弹性矩阵如下： 单元1的应力矩阵如下： 单元2的应力矩阵如下： （4）求取单元刚度矩阵K

单元刚度矩阵K满足公式K=B’*D*B*t*A,其中t为平面板的厚度，A为单元面积，且单元刚度矩阵为对称矩阵。 单元1的刚度矩阵如下： 单元2的刚度矩阵如下： （5）求取总体刚度矩阵sumKK 由上述步骤求得的单元刚度矩阵K利用单元虚功原理和刚度方程可导出K’*δ=f，其中δ为单元节点位移列阵，f为单元等效节点载荷列阵，为了能将各个单元刚度方程统一到一个整体，便需要步骤（1）的单元节点提取矩阵对单元刚度方程进行变换，将两个变换结果联立便得到总体刚度方程，其中也可得到总体刚度矩阵sumKK，且总体刚度矩阵可由sumKK=ΣT’*K*T求得。 总体刚度矩阵如下：

《结构力学习题集》下矩阵位移法习题及答案 2

第七章 矩阵位移法 一、就是非题 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性与奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 就是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。 6、结 构 刚 度 矩 阵 就是 对 称 矩 阵 ,即 有K i j = K j i ,这 可 由 位 移 互 等 定 理 得 到 证 明 。 7、结构刚度方程矩阵形式为:[]{}{}K P ?=,它就是整个结构所应满足的变形条件。 8、在直接刚度法的先处理法中,定位向量的物理意义就是变形连续条件与位移边界条件。 9、等效结点荷载数值等于汇交于该结点所有固端力的代数与。 10、矩阵位移法中,等效结点荷载的“等效原则”就是指与非结点荷载的结点位移相等。 11、矩阵位移法既能计算超静定结构,也能计算静定结构。 二、选择题 1、已知图示刚架各杆EI = 常数,当只考虑弯曲变形,且各杆单元类型相同时,采用先处理法进行结点位移编号,其正确编号就是: (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 2、平面杆件结构一般情况下的单元刚度矩阵[]k 66?,就其性质而言,就是: A.非对称、奇异矩阵; B.对称、奇异矩阵; C.对称、非奇异矩阵; D.非对称、非奇异矩阵。 3、单元i j 在图示两种坐标系中的刚度矩阵相比: A.完全相同; B.第2、3、5、6行(列)等值异号;

MATLAB在结构力学分析中应用

MATLAB在结构力学分析中的应用摘要:传统的手算方法解超静定结构工作量繁重，有时甚至是不可能，运用结构有限元编程的一般方法，通过两个实例的对照，展示matlab在结构力学分析中的应用，matlab具有高性能，方法具有普遍的适用性,实现弯矩图自动绘制。 关键词: matlab结构有限元弯矩图 abstract：while using the traditional manual method to resolve complex statically indeterminate structures, it is heavy workloads, sometimes even impossible，using finite element programming of the general method, based on two examples, this paper introduces a method of application of matlab in structure mechanics, matlab has the advantages of high performance, it can be applied to many kinds of structures, realization of automatic drawing bending moment diagram. key words: matlab; finite element; bend moment diagram 引言 结构力学[3]中，常利用传统的力法与位移法求解超静定结构，力法是几何问题，位移法把复杂的几何图乘转化为代数运算，但它们基本未知量很多时，系数构成的矩阵计算巨大，两者都不能满足科研工作者的需要。应用matlab软件丰富可靠的矩阵运算、数据处理、图形绘制等便利工具，可使得计算和图象一体化。对于结构

结构力学位移法题与答案解析

超静定结构计算一S移法 —.判断题： Is判断下列结构用位移法计算时基本未知呈的数目。 2、位移法求解结构力时如果Mp图为零，则自由项血一走为零。 3、位移法未知呈的数目与结构的超静定次数有关。 4、位移法的基本结构可以是静定的，也可以是超静走的。 5、位移法典型方程的物理意义反映了原结构的位移协调条件。 二计算题： (2) (3) (1) (6) ￡/=■ El El EA 2EI 、b EA E/=oc d 4EI一— J E/=oo 2E1 4A7 2EI 4 El

12.用位移法计算图示结构并作〃图，横梁刚度EA -8 ,两柱线刚度/相同。 13、用位移法计算图示结构并作〃图。F/二常数。 14、求对应的荷载集度g。图示结构横梁刚度无限大。已知柱顶的水平位移为512/（3 曰）（T）。 15、用位移法计算图示结构州乍M图。曰=常数。

16、用位移法计算图示结构r求出未知呈，各杆曰相同。 4m 4m 19、用位移法计算图示结构并作〃图。 -2/ 2f q 二i i 20、用位移法计算图示结构并作〃图。各杆日=営数r q = 20kN/m o 6m 4 ------- B 6m 6m R --- k ----- 1 23、用位移法计算图示结构州乍M图。曰=常数。 7T7F 24、用位移法计算图示结构州乍M图。曰=常数。

°^=ZJ 週AV 酔辭圍闕￥觀⑨由、充 。回申Z7阴甘县欲 遍如士星與莎竺园蔑44辛觀⑨由、6 乙 Ic n n M M I Z M f c/i in

38、用位移法计算图示结构并作〃图。曰=常数。 42、用位移法计算图示结构州乍〃图。 43、用位移法计算图示结构州乍〃图。曰=常数。 48、已知0点的位移0，求几

结构力学习题集(下)-矩阵位移法习题及答案

第七章 矩阵位移法 一、是非题 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。 6、结 构 刚 度 矩 阵 是 对 称 矩 阵 ，即 有K i j = K j i ，这 可 由 位 移 互 等 定 理 得 到 证 明 。 7、结构刚度程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 8、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 9、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 10、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 11、矩阵位移法既能计算超静定结构，也能计算静定结构。 二、选择题 1、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 2、平面杆件结构一般情况下的单元刚度矩阵[]k 66?，就其性质而言，是： A ．非对称、奇异矩阵； B ．对称、奇异矩阵； C ．对称、非奇异矩阵； D ．非对称、非奇异矩阵。 3、单元i j 在图示两种坐标系中的刚度矩阵相比： A ．完全相同；

MATLAB程序：已知三个位置设计平面四杆机构求解程序(位移矩阵法)

%MATLAB程序：已知三个位置设计平面四杆机构求解程序（位移矩阵法） clear;clc; %凡是变量名前带v的为数值变量，不带的是符号变量 vxp1=0; vyp1=0; vsita1=0*pi/180; vxp2=-2; vyp2=6; vsita2=40*pi/180; vxp3=-10; vyp3=8; vsita3=90*pi/180; %精确位置P1，P2，P3及各角度 vsita12=vsita2-vsita1; vsita13=vsita3-vsita1; vxa=-10; vya=-2; vxd=-5; vyd=-2; %选定A,D点 %所有数值均在此确定，更改此处即可解出不同数值的四杆机构位移矩阵方程 syms xp1 yp1 xp2 yp2 xp3 yp3 sita12 sita13; syms xa ya xb1 yb1 xb2 yb2 xb3 yb3; f1='(xb2-xa)^2+(yb2-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; f2='(xb3-xa)^2+(yb3-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; %前两个机构方程 f3='xb2=cos(sita12)*xb1-sin(sita12)*yb1+xp2-xp1*cos(sita12)+yp1*sin(sita12)'; f4='yb2=sin(sita12)*xb1+cos(sita12)*yb1+yp2-xp1*sin(sita12)-yp1*cos(sita12)'; %由第一个位移矩阵方程得出 f5='xb3=cos(sita13)*xb1-sin(sita13)*yb1+xp3-xp1*cos(sita13)+yp1*sin(sita13)'; f6='yb3=sin(sita13)*xb1+cos(sita13)*yb1+yp3-xp1*sin(sita13)-yp1*cos(sita13)'; %由第二个位移矩阵方程得出 f1=subs(f1,{xa,ya},{vxa,vya}); f2=subs(f2,{xa,ya},{vxa,vya}); f3=subs(f3,{xp1,xp2,yp1,sita12},{vxp1,vxp2,vyp1,vsita12}); f4=subs(f4,{xp1,yp1,yp2,sita12},{vxp1,vyp1,vyp2,vsita12}); f5=subs(f5,{xp1,xp3,yp1,sita13},{vxp1,vxp3,vyp1,vsita13}); f6=subs(f6,{xp1,yp1,yp3,sita13},{vxp1,vyp1,vyp3,vsita13}); %代入具体数值 [xb1,xb2,xb3,yb1,yb2,yb3]=solve(f1,f2,f3,f4,f5,f6); %解方程 vxb1=vpa(xb1); vyb1=vpa(yb1); vxb2=vpa(xb2); vyb2=vpa(yb2); vxb3=vpa(xb3); vyb3=vpa(yb3); (vxb1-vxa)^2+(vyb1-vya)^2; (vxb2-vxa)^2+(vyb2-vya)^2; (vxb3-vxa)^2+(vyb3-vya)^2; %去掉这三行分号可验证B点三个位置是否距离A点相等 syms xd yd xc1 yc1 xc2 yc2 xc3 yc3;

《结构力学》典型习题与解答

《结构力学》经典习题及详解 一、判断题（将判断结果填入括弧内，以 √表示正确 ，以 × 表示错误。） 1．图示桁架结构中有3个杆件轴力为0 。（×） 2．图示悬臂梁截面A 的弯矩值是ql 2。 (×) l l 3．静定多跨梁中基本部分、附属部分的划分与所承受的荷载无关。（√ ） 4．一般来说静定多跨梁的计算是先计算基本部分后计算附属部分。（× ） 5．用平衡条件能求出全部内力的结构是静定结构。（ √ ） 6．求桁架内力时截面法所截取的隔离体包含两个或两个以上的结点。（√ ） 7．超静定结构的力法基本结构不是唯一的。（√） 8．在桁架结构中，杆件内力不是只有轴力。（×） 9．超静定结构由于支座位移可以产生内力。 （√ ） 10．超静定结构的内力与材料的性质无关。（× ） 11．力法典型方程的等号右端项不一定为0。 （√ ） 12．计算超静定结构的位移时，虚设力状态可以在力法的基本结构上设。（√） 13．用力矩分配法计算结构时，汇交于每一结点各杆端分配系数总和为1，则表明分配系 数的计算无错误。 （× ） 14．力矩分配法适用于所有超静定结构的计算。（×） 15．当AB 杆件刚度系数i S AB 3 时，杆件的B 端为定向支座。 (×)

二、单项选择题（在每小题的四个备选答案中选出一个正确答案，并将其代号填在题干后面的括号内。不选、错选或多选者，该题无分。） 1．图示简支梁中间截面的弯矩为（ A ） q l A ． 82ql B ． 42ql C ． 22 ql D ． 2ql 2．超静定结构在荷载作用下产生的内力与刚度（B ） A ． 无关 B ． 相对值有关 C ． 绝对值有关 D ． 相对值绝对值都有关 3．超静定结构的超静定次数等于结构中（B ） A ．约束的数目 B ．多余约束的数目 C ．结点数 D ．杆件数 4．力法典型方程是根据以下哪个条件得到的（C ）。 A ．结构的平衡条件 B ．结构的物理条件 C ．多余约束处的位移协调条件 D ．同时满足A 、B 两个条件 5． 图示对称结构作用反对称荷载，杆件EI 为常量，利用对称性简化后的一半结构为（A ）。 6．超静定结构产生内力的原因有（D ） A ．荷载作用与温度变化 B ．支座位移 C ．制造误差 D ．以上四种原因

matlab连续梁程序的编制与使用

第三章连续梁程序的编制与使用 入结构力学领域中而产生的一种方法，而Matlab 语言正是进行矩阵运算的强大工具，因此，用 Matlab语言编写结构力学程序有更大的优越性。本 章将详细介绍如何利用Matlab语言编制连续梁结 构的计算程序。 矩阵位移法的解题思路是将结构离散为单元 （杆件），建立单元杆端力与杆端位移之间的关系 －单元刚度方程；再将各单元集成为原结构，在满 足变形连续条件和平衡条件时，建立整体刚度方 程；在边界条件处理完毕后，由整体刚度方程解出 节点位移，进而求出结构内力。 用矩阵位移法计算连续梁的步骤如下： 1）整理原始数据，如材料性质、荷载条件、约 束条件等，并进行编码：单元编码、结点编 码、结点位移编码、选取坐标系。 2）建立局部坐标系下的单元刚度矩阵。 3）建立整体坐标系下的单元刚度矩阵。 4）集成总刚。 5）建立整体结构的等效节点荷载和总荷载矩 阵 6）边界条件处理。 7）解方程，求出节点位移。 8）求出各单元的杆端内力。 实际上，上述步骤也是编制Matlab程序的基本 步骤，在求出计算结果后，还可以利用Matlab的绘图功能绘制结构图、内力图、 变形图等等。 图3-1程序流程图

3.1 程序说明 %******************************************************************* % 矩阵位移法解连续梁主程序 %******************************************************************* ●功能：运用矩阵位移法解连续梁的基本原理编制的计算主程序。 ●基本思想：结点（结点位移）编码默认为从左至右，从1开始顺序进行；杆 端弯矩的方向默认为逆时针。 ●荷载类型：可计算结点荷载，每单元作用的跨中集中力和均布荷载。 ●说明：主程序的作用是通过赋值语句、读取和写入文件、函数调用等完成算 法的全过程，即实现程序流程图的程序表达。 %----------------------------------------------------------------------------------------------------- 1 程序准备 format short e %设定输出类型 clear all %清除所有已定义变量 clc %清屏 ●说明： format short e －设定计算过程中显示在屏幕上的数字类型为短格式、科学计数法； clear all －清除所有已定义变量，目的是在本程序的运行过程中，不会发生变量名相同等可能使计算出错的情况； clc －清屏，使屏幕在本程序运行开始时 %----------------------------------------------------------------------------------------------------- 2 打开文件 FP1=fopen('input.txt','rt'); %打开输入数据文件存放初始数据FP2=fopen('output.txt','wt'); %打开输出数据文件存放计算结果 ●说明： FP1=fopen('input.txt','rt'); －打开已存在的输入数据文件input.txt，且设置其为只读格式，使程序在执行过程中不能改变输入文件中的数值，并用文件句柄FP1来 FP2=fopen('output.txt','wt'); －打开输出数据文件，该文件不存在时，通过此

《结构力学习题集》-矩阵位移法习题及答案

第八章 矩阵位移法 – 老八校 一、判断题： 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、结构刚度方程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。 7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 10、矩阵位移法既能计算超静定结构，也能计算静定结构。 11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 ( )

二、计算题： 12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。 12 3l l 4 l 5EI 2EI EA (0,0,0) (0,0,1) (0,2,3) (0,0,0) (0,2,4)(0,0,0) EI 13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。EI ，EA 均为常数。 l 14、计算图示结构整体刚度矩阵的元素665544,,K K K 。E 为常数。 l l 1 3 4 2 A , I A A /222A I , 2A 15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵 [][]K K 22 24 ,。 [][]k k 1112 [][] k k 2122 [] k = i i i i i 单刚分块形式为 ：

基于MATLAB的平面刚架静力分析

基于MATLAB 的平面刚架静力分析 为了进一步理解有限元方法计算的过程，本文根据矩阵位移法的基本原理应用MATLAB 编制计算程序对以平面刚架结构进行了静力分析。本文还利用ANSYS 大型商用有限元分析软件对矩阵位移法的计算结果进行校核，发现两者计算结果相当吻合，验证了计算结果的可靠性。 一、 问题描述 如图1所示的平面刚架，各杆件的材料及截面均相同，E=210GPa ，截面为0.12×0.2m 的实心矩形，现要求解荷载作用下刚架的位移和内力。 5m 4m 3m 图1 二、矩阵位移法计算程序编制 为编制程序方便考虑，本文计算中采用“先处理法”。具体的计算步骤如下。

（1） 对结构进行离散化，对结点和单元进行编号，建立结构（整体）坐标系 和单元（局部）坐标系，并对结点位移进行编号； （2） 对结点位移分量进行编码，形成单元定位向量e λ； （3） 建立按结构整体编码顺序排列的结点位移列向量δ，计算固端力e F P 、等 效结点荷载E P 及综合结点荷载列向量P ； （4） 计算个单元局部坐标系的刚度矩阵，通过坐标变换矩阵T 形成整体坐标 系下的单元刚度矩阵e T e K T K T = ； （5） 利用单元定位向量形成结构刚度矩阵K ； （6） 按式1=K P δ- 求解未知结点位移； （7） 计算各单元的杆端力e F 。 根据上述步骤编制了平面刚架的分析程序。程序中单元刚度矩阵按下式计算。 32322 23 2 32 22 0000 1261260 064620 00001261260062640 EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l K EA EA l l EI EI EI EI l l l l EI EI EI EI l l l l ??- ??? ???- ?? ? ???- ??? ?=??-?? ? ???---??? ???-??? ?

《结构力学习题集》-矩阵位移法习题及标准答案

第八章 矩阵位移法 一、判断题： 1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。 2、单元刚度矩阵均具有对称性和奇异性。 3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。 4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。 5、结构刚度方程矩阵形式为：[]{}{}K P ?=，它是整个结构所应满足的变形条件。 6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。 7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。 8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 10、矩阵位移法既能计算超静定结构，也能计算静定结构。 11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： (0,1,2) (0,0,0) (0,0,0) (0,1,3) (0,0,0)(1,2,0) (0,0,0)(0,0,3) (1,0,2) (0,0,0) (0,0,0)(1,0,3) (0,0,0) (0,1,2) (0,0,0)(0,3,4) A. B. C. D. 2134123412341234 ( )

二、计算题： 12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。 12 3l l 4 l 5EI 2EI EA (0,0,0) (0,0,1) (0,2,3) (0,0,0) (0,2,4)(0,0,0) EI 13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。EI ，EA 均为常数。 l 14、计算图示结构整体刚度矩阵的元素665544,,K K K 。E 为常数。 l l 1 3 4 2 A , I A A /222A I , 2A 15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵 [][]K K 22 24 ,。 [][]k k 1112 [][] k k 2122 [] k = i i i i i 单刚分块形式为 ：

矩阵位移法大作业

矩阵位移法大作业 学号：151210122 姓名：谭逸天 班级：土木一班

编制原理： 使用Math Work公司开发的科学与工程计算机软件——MATLAB， 利用其矩阵运算的便利性，将题目要求结构的基本信息编入脚本命令文件中，并编入求解步骤。加上刚度信息的输入指令，以及提取解答要求信息并输出的指令。令使用者只需输入结构材料相关信息便可计算题目对应悬索—拱组合体系的信息，并直接在命令窗口输出。 利用计算套路的重复性，程序开发时进行模块化设计。再由重复单元完成多次、重复的运算。 从整体性考虑，数据储存采用“算后集装，装后回收”对变量及数组重复使用，由配音进行简单命名，提高可辨识度。由于计算套路及程序本身高度模块化，并且题目所需个体信息相对于整体极少，提取个体化的信息只需简单改造命令模块，从整体信息中提取处理得出。编程所需的“数据化”“编码”等预处理由人工在编程开始前完成，由左下斜索基座作原点，正右向为X轴正向，正上为Y轴正向，建立右手系。编码顺序从左倒右由上及下，并用先处理法处理基座。（如下图所示）

6 7 共45个单元，32个结点编号，71个位移编号。 本人学号对应节间数m=14;f1=7L/4;f2=7L/10;h=7L/2;以上数据 为编程中人工设定值，结构的其余信息根据用户的输入进行计算得出。

程序说明： 初始计算结构在坐标系中的坐标信息，手动编入悬索与拱的曲线关键点信息，代入方程求解。随后由循环语句模块计算并存储结构中各类杆件的角度、长度信息，采用以直代曲的方法处理曲线。 由于先处理法，两端各四个单元不与其余单元通用编码递进规律，采用单独的语句进行计算并集装入总体信息储存矩阵中，其余规律性单元信息由循环的语句模块进行集装，便于之后的计算。定位向量统一装至71行6列的矩阵“dingwei”中，单元的长度与夹角信息统一装至71行2列的矩阵“danyuan”中，第一列为长度，第二列为角度。使两个信息矩阵的行序号对应单元序号，便于之后使用。 之后进入单元分析部分。先是对上部悬索进行单元分析，此部分为桁架单元，从“danyuan”矩阵中提取长度信息与角度信息，结合 开始时输入的刚度信息组装单刚矩阵与坐标变换矩阵，进行坐标变换后直接提取定位向量进行集装部分总刚矩阵的步骤。集装命令通过循环嵌套配合判断语句，对单刚矩阵进行二维遍历，并提取合格的元素填充至对应位置。随后，通过少量改动实现对斜索、吊杆、拱、主塔的处理。 之后保留基本结构，进行单元结点荷载的分析，并集装出结构结点荷载矩阵。 之后通过简单矩阵运算即得结构结点位移列阵。 进入单元后处理。将集装循环语句进行改造，达成逆向提取单元结点位移的功能。提取之前存储的单元信息进行坐标变换。最后算出

结构力学位移法题及答案

超静定结构计算——位移法 一、判断题： 1、判断下列结构用位移法计算时基本未知量的数目。 （1） （2） （3） （4） （5） （6） EI EI EI EI 2EI EI EI EI EA EA a b EI= EI=EI= 24442 2、位移法求解结构内力时如果P M 图为零，则自由项1P R 一定为零。 3、位移法未知量的数目与结构的超静定次数有关。 4、位移法的基本结构可以是静定的，也可以是超静定的。 5、位移法典型方程的物理意义反映了原结构的位移协调条件。 二、计算题： 12、用位移法计算图示结构并作M 图，横梁刚度EA →∞，两柱线刚度 i 相同。 2 13、用位移法计算图示结构并作M 图。E I =常数。

l l l/2l/2 14、求对应的荷载集度q。图示结构横梁刚度无限大。已知柱顶的水平位移为 () 5123 /() EI→。 12m12m 8m q 15、用位移法计算图示结构并作M图。EI =常数。 l l l 16、用位移法计算图示结构，求出未知量，各杆EI相同。 4m 19、用位移法计算图示结构并作M图。 q l l

20、用位移法计算图示结构并作M 图。各杆EI =常数，q = 20kN/m 。 6m 6m 23、用位移法计算图示结构并作M 图。EI =常数。 l l 2 24、用位移法计算图示结构并作M 图。EI =常数。 q 29、用位移法计算图示结构并作M 图。设各杆的EI 相同。 q q l l /2/2 32、用位移法作图示结构M 图。 E I =常数。

q l l /2 l /2l 36、用位移法计算图示对称刚架并作M 图。各杆EI =常数。 l l 38、用位移法计算图示结构并作M 图。EI =常数。 q l l l l 42、用位移法计算图示结构并作M 图。 2m 2m 43、用位移法计算图示结构并作M 图。EI =常数。

结构静力计算和动力计算的对比分析

结构静力计算与动力计算的对比分析 结构精力计算和结构动力计算是一个比较理论化和深度比较广的论述题目，在此，我仅凭本人有限的学识来展开对两者内容及关系的介绍和论述。也藉此契机，对结构力学上下册作一个比较系统的梳理和总结，为以后的学习以及工作打下坚实的基础。 首先，我想先介绍一下有关结构力学的基本概念，让读者可以带着一个整体、宏观的概念去深入理解具体的内部结构内容。 那么，我想从静力荷载和动力荷载的含义入题。静力荷载是指其大小、方向和位置不随时间变化或变化很缓慢的荷载，它不致使结构产生显著的加速度，因而可以略去惯性力的影响。结构的自重及其他恒荷载即属于静力荷载。动力荷载是指随时间迅速变化的荷载，它将引起结构振动，使结构产生不容忽视的及速度，因而必须考虑惯性力的影响。除荷载外，还有其他一些非荷载因素作用也可使结构产生内力和位移，例如温度变化、制造误差、材料收缩以及松弛、徐变等。 在结构静力计算中，最核心的内容就是计算结构的位移，而一切都要从虚功原理说起。虚功原理的两种表述：1、对于刚体体系，刚体体系处于平衡的必要和充分条件是，对于任何虚位移，所有外力所作虚功总和为零；2、对于变形体系，其处于平衡的必要和充分条件是，对于任何虚位移，外力所作虚功总和等于各微段上的内力在其变形上所作的虚功总和，简单说，外力虚功等于内力虚功。 虚功方程： 由于力状态与位移状态是彼此独立无关的，因此运用单位荷载法： 由： 得位移计算一般公式： 同过几何关系可得弯矩图乘法便捷计算公式（为计算带来极大的方便）: 力法： 力法典型方程： （系数δ?、的求解方法如同上述虚功原理的原理。） 该方程的物理意义为：基本结构在全部多余未知力和荷载共同作用下，在去掉各多余联系处沿各多余未知力方向的位移，应与原结构相应的位移相等。可见，力法可以求解出超静N u s s W F d Md F d ?γ=++∑∑∑???1k R N u s s F c F d Md F d ?γ?+=++∑∑∑∑??? N S S s k N s R F ds Md F d F M F F c EA EI GA γ?=++-∑∑∑∑???S w c Md A y M EI EI =∑?1111221211222200P P X X X X δδδδ++?=??++?=?基本体系 1 X 结 构

结构力学位移法解析

第十章位移法 §10－1 概述 位移法——以结点位移（线位移，转角）为基本未知量的方法。 基本概念：以刚架为例（图10-1） 基本思路：以角位移Z1为基本未知量 平衡条件——结点1的力矩平衡 位移法要点：一分一合 ①确定基本未知量（变形协调）基本体系－独立受力变形的杆件 ②将结构拆成杆件－杆件分析（刚度方程－位移产生内力、荷载产生内力） ③将结构杆件合成结构：整体分析——平衡条件——建立方程 §10－2 等截面直杆的转角位移方程 单跨超静定梁——由杆端位移求杆端力——转角位移方程 矩阵形式 一、端（B端）有不同支座时的刚度方程 （1）B端固定支座 （2）B端饺支座 （3）B端滑动支座 二、由荷载求固端力（3*，4，11*，12，19，20） （1）两端固定 （2）一端固定，一端简支 （3）一端固定，一端滑动（可由两端固定导出） 三、一般公式 叠加原理杆端位移与荷载共同作用 杆端弯矩：（10-1） 位移法意义（对于静定、超静定解法相同） 基本未知量－被动（由荷载等因素引起） →按主动计算——位移引起杆端力＋荷载的固端力 →结点满足平衡 正负号规则——结点转角（杆端转角） 弦转角——顺时针为正 杆端弯矩 位移法三要素： 1.基本未知量－独立的结点位移 2.基本体系－原结构附加约束，分隔成独立变力变形的杆件体系。 3.基本方程－基本体系在附加约束上的约束力（矩）与原结构一致 （平衡条件）

§10－3基本未知量的确定 角位移数＝刚结点数（不计固定端） 线位移数＝独立的结点线位移 观察 几何构造分析方法——结点包括固定支座）变铰结点 铰结体系的自由度数＝线位移数 ――即使其成为几何不变所需添加的链杆数。 §10－4典型方程及计算步骤 典型方程（10－5、6） 无侧移刚架的计算 无侧移刚架－只有未知结点角位移的刚架（包括连续梁）（△＝0） 有侧移刚架计算 有侧移刚架――除结点有位移外还有结点线位移 求解步骤： （1）确定基本未知量：Z i （按正方向设基本未知量）——基本体系， （2）作荷载、Z i = 1 —— ()()01i P i i M M ??==、图 （3）求结点约束力矩：荷载 —— 自由项R Ip ，及ΔJ = 1 —— 刚度系数 k IJ （4）建立基本方程：[k IJ ]{ Z i } + { R Ip } = {0} —— 附加约束的平衡条件 求解Z i （Δi ） (5) 叠加法作i i P Z M M M ∑+= §10－5 直接建立位移法方程 求解步骤： （1）确定基本未知量：Z i （按正方向设基本未知量）——基本体系， （2）写杆端弯矩（转角位移方程） （3）建立位移法方程—— 附加约束的平衡，求解Z i (4) 叠加法作i i P Z M M M ∑+= §10－6 对称性利用 对称结构 对称荷载作用 —— 变形对称，内力对称 （M 、N 图对称，Q 图反对称——Q 对称） 反对称荷载作用 —— 变形反对称，内力反对称 （M 、N 图反对称，Q 图对称——Q 反对称） —— 取半跨 对称结构上的任意荷载 ——对称荷载＋反对称荷载

结构力学-第7章 位移法解析

第7章位移法 一. 教学目的 掌握位移法的基本概念； 正确的判断位移法基本未知量的个数； 熟悉等截面杆件的转角位移方程； 熟练掌握用位移法计算荷载作用下的刚架的方法 了解位移法基本体系与典型方程的物理概念和解法。 二. 主要章节 §7-1 位移法的基本概念 §7-2 杆件单元的形常数和载常数—位移法的前期工作 §7-3 位移法解无侧移刚架 §7-4 位移法解有侧移刚架 §7-5 位移法的基本体系 §7-6 对称结构的计算 ＊§7-7支座位移和温度改变时的位移法分析（选学内容） §7-8小结 §7-9思考与讨论 三. 学习指导 位移法解超静定结构的基础是确定结构的基本未知量以及各个杆件的转角位移方程，它不仅可以解超静定结构，同时还可以求解静定结构，另外，要注意杆端弯矩的正负号有新规定。 四. 参考资料 《结构力学(Ⅰ)-基本教程第3版》P224～P257 第六章我们学习了力法，力法和位移法是计算超静定结构的两个基本方法，力法发展较早，位移法稍晚一些。力法把结构的多余力作为基本未知量，将超静定结构转变为将定结构，按照位移条件建立力法方程求解的；而我们今天开始学的这一章位移法则是以结构的某些位移作为未知量，先设法求出他们，在据以求出结构的内力和其他位移。由位移法的基本原理可以衍生出其他几种在工程实际中应用十分普遍的计算方法，例如力矩分配法和迭代法等。因此学习本章内容，不仅为了掌握位移法的基本原理，还未以后学习其他的计算方法打下良好的基础。此外，应用微机计算所用的直接刚度法也是由位移法而来的，所以本章的内容也是学习电算应用的一个基础。

本章讨论位移法的原理和应用位移法计算刚架，取刚架的结点位移做为基本未知量，由结点的平衡条件建立位移法方程。位移法方程有两种表现形式：①直接写平衡返程的形式（便于了解和计算）② 基本体系典型方程的形式（利于与力法及后面的计算机计算为基础的矩阵位移法相对比，加深理解） §7-1 位移法的基本概念 1.关于位移法的简例 为了具体的了解位移法的基本思路，我们先看一个简单的桁架的例子：课本P225。图7-1和图7-2所示。 (a) (a) (b) (b) 图7-1 图7-2 第一步：从结构中取出一个杆件进行分析。（杆件分析） 图7-2中杆件AB 如已知杆端B 沿杆轴向的位移为i u （即杆件的伸长）则杆端力Ni F 为： i i i Ni u l EA F （7-1） E-为弹性模量，A-为杆件截面面积，i l -为杆件长度

结构力学之矩阵位移法

第十二章 矩阵位移法 【例12-1】 图 a 所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。分别用位移法和矩阵位移法计算。 图12-1 解：（1）位移法解 ?基本未知量和基本结构的确定 用位移法解的基本结构如图c 所示。这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。 ?位移法基本方程的建立 ?? ? ?? =+θ+θ+θ=+θ+θ+θ=+θ+θ+θ000333323213123232221211313212111P P P R K K K R K K K R K K K 将上式写成矩阵形式

?? ??? ?????=??????????+??????????θθθ?? ????????0003213213332 31 232221131211P P P R R R K K K K K K K K K ?系数项和自由项 计算（须绘出单位弯矩图和荷载弯矩图） 由图d ，结点力矩平衡条件 ∑=0M ，得 EI K 411=，l EI K 221=，031=K 由图e ，结点力矩平衡条件 ∑=0M ，得 l EI K 212=，l EI l EI l EI K 84422=+=，l EI K 232= 由图f ，结点力矩平衡条件 ∑=0M ，得 013=K ，l EI K 223=，l EI EI EI K 84433=+= 由图g ，结点力矩平衡条件 ∑=0M ，得 81Pl R p -=，2Pl R P -=，03=P R 将系数项和自由项代入位移法基本方程，得 ??? ???????=??????????--+?? ??? ?????θθθ??????????0000118820282024321Pl l EI ?解方程，得?? ????????-= ?? ? ?? ?????θθθ14114162321EI Pl ?由叠加法绘弯矩图，如图h 所示。 （2）矩阵位移法解 ?对单元和结点编号（图a ） 本题只考虑弯曲变形的影响，故连续梁每个结点只有一个角位移未知数。若用后处理法原始结构刚度阵为44?阶；用先处理法结构刚度阵为33?阶（已知角位移04=θ）。下面采用先处理法来说明矩阵位移法计算过程。 单元标准形式为（图b ） )(e k ?? ????=?? ?? ??????=)()()()() (4224e jj e ji e ij e ii e k k k k l EI l EI l EI l EI