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

连续梁程序的编制与使用 结构力学中的矩阵位移法是随着电子计算机进入结构力学领域中而产生的一种方法，而Matlab语言正是进行矩阵运算的强大工具，因此，用Matlab语言编写结构力学程序有更大的优越性。本章将详细介绍如何利用Matlab语言编制连续梁结构的计算程序。 矩阵位移法的解题思路是将结构离散为单元（杆件），建立单元杆端力与杆端位移之间的关系－单元刚度方程；再将各单元集成为原结构，在满足变形连续条件和平衡条件时，建立整体刚度方程；在边界条件处理完毕后，由整体刚度方程解出节点位移，进而求出结构内力。 用矩阵位移法计算连续梁的步骤如下： 整理原始数据，如材料性质、荷载条件、约束条件等，并进行编码：单元编码、结点编码、结点位移编码、选取坐标系。 建立局部坐标系下的单元刚度矩阵。 建立整体坐标系下的单元刚度矩阵。 集成总刚。 建立整体结构的等效节点荷载和总荷载矩阵 边界条件处理。 解方程，求出节点位移。 求出各单元的杆端内力。 实际上，上述步骤也是编制Matlab程序的基本步骤，在求出计算结果后，还可以利用Matlab的绘图功能绘制结构图、内力图、变形图等等。 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); －打开输出数据文件，该文件不存在时，通过此命令创建新文件，该文件存在时则将原有内容全部删除。该文件设置为可写格式，可在程序执行过程中向输出文件写入数据。 %----------------------------------------------------------------------------------------------------- 3 读入程序控制信息 NELEM =fscanf(FP1,%d,1); %单元个数（单元编码总数） NPOIN =fscanf(FP1,%d,1); %结点个数（结点编码总数） NVFIX =fscanf(FP1,%d,1); ％约束个数（零位移总数） NFPOIN=fscanf(FP1,%d,1); ％结点荷载个数（作用在结点上集中力偶总数） NFPRES=fscanf(FP1,%d,1); ％非结点荷载数（作用在单元上分布荷载总数） YOUNG=fscanf(FP1,%f,1); ％弹性模量 说明： 从输入文件FP1中读入单元个数，结点个数，约束个数，结点荷载个数，非结点荷载个数，弹性模量； 程序中弹性模量仅输入了一个