基于MATLAB中SSI函数的杆件受力分析.docx

2015-6-1 2015-6-1 姓 名： 王 博 姓 名： 王 博 班 级：工程造价11043 学 号： 1416220327 指导老师: 詹志兰 zz 基于 基于matlab中SSI函数的杆件受力分析 摘要:我们对杆件组合的受力分析是通过对关键点进行分析得到的，是对各个点组合的分析，本文以图论的思维来来分析杆件受力，并通过matlab中的ssi程序实现各个受力计算，并进行绘图。 关键字:杆件受力分析；matlab；ssi函数 引言 MATLAB是美国MathWorks公司出品的商业数学软件，主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言的编辑模式，代表了当今国际科学计算软件的先进水平。 本文利用matlab中的ssi和绘图程序进行杆件组合的静力学分析。 SSI函数 原理 由静力学的静力平衡公式可以知道，我们对一个杆件组合的受力分析是通过对关键面和关键点得到的，各个点的数值在X方向，Y方向和Z方向的加总之和是为0的。由此我们只需找准计算起止点，并对过程中相互作用的各个点进行分析就可以了。 2）程序 function [ Elem ] = SSI( Joint,Elem,SuptUnknow,varargin) %SSI 使用联立法计算单元坐标系中单元杆端内力 % %调用：[ Elem ] = SSI( Joint,Elem,SuptUnknow) % 计算刚架单元杆端内力，无内力约束条件 % [ Elem ] = SSI( Joint,Elem,SuptUnknow,Rest) % 计算刚架单元杆端内力，有内力约束条件 % [ Elem ] = SSI( Joint,Elem,Rest,SuptUnknow,varargin) % 计算结构单元杆端内力 % %输入： Joint = 结点信息 % Elem = 单元信息 % SuptUnknow = 支座反力向量 % varargin{1} = Rest,单元杆端内力约束条件矩阵 % 每行为：[单元编号，单元杆端内力分量编号(1～6)] % varargin{2} = JointVec：按刚架结点和桁架结点分类存放结点编号的元胞数组 % 默认：JointVec = 1:NJointNo，NJointNo为最大结点号 % varargin{3} = isCom：存放刚架结点与桁架结点逻辑变量的向量 % isCom = [isCom1，isCom2] % isCom1或isCom2 = true：刚架结点 % isCom1或isCom2 = false：桁架结点 % 默认：isTruss = false % varargin{4} = ElemVec：按刚架单元和桁架单元分类存放单元编号的元胞数组 % 默认：ElemVec = 1:NElemNo，NElemNo为最大单元号 % varargin{5} = isTruss：存放刚架或桁架单元的逻辑变量的向量 % isTruss = [isTruss1，isTruss2] % isTruss1或isTruss2 = true：桁架单元 % isTruss1或isTruss2 = false：刚架单元 % 默认：isTruss = false % %输出：Elem.InFL = 单元坐标系中单元杆端内力 案例分析 图1 杆件组合受力及其分析 对图中的案例分析， 单元及力的表示 图2 单元力分析图 如图2所示，我们将我们将杆件单元中的杆段内力写成6v1的矩阵 Fe=[FX1 0 0 FX2 0 0]T 对结构坐标中单元杆端内力写为 FE=[FX1 FY1 0 FX2 FY2 0 0]T 表示杆端弯矩为0，杆件单元