基于MATLAB语言的结构重分析技术的应用.pdfVIP

工程地质计算机应用 2012年第4期 总68期 ※ 编程技术 应用实践※ 基于MATLAB语言的结构重分析技术的应用 詹顺 (深圳华森建筑与工程设计顾问有限公司) 【摘 要】在结构设计过程中大量的重分析问题是导致计算效率低下的重要原因，构造高效实用的重分析算法是解决这 一 问题的有效途径。本文主要围绕结构分析、重分析方法展开研究，并探讨了相应程序的实施。经算例演示表明，采用本文 提供的结构分析和重分析方法及编制的程序进行桁架结构的计算分析是可行有效的，可以用于结构优化设计等后继的系列研 究中。 【关键词】桁架 结构重分析 结构修改 基向量 精确度 1．静力重分析 不管是传统的结构设计，还是现代的结构优化设计，都需要对结构进行计算分析。结构优化问题变量 多，约束也多，并且大都是复杂的隐式函数，每做一次重分析的工作量浩大。为结构优化设计增加了难度。 对此，一个高效准确的重分析程序对于结构优化是相当重要的，故而，对结构重分析方法的研究是很有价 值的。 本文主要研究的是结构静力重分析问题的迭代算法，这里只给出结构静力重分析问题的简单描述： 设初始设计变量为X ，相应的初始整体劲度矩阵K ，相应的初始位移r牛，结构所受荷载为P。可有 公式： K r}=P (1．1) 假设设计变量的变化为△x，则修改后的设计x： X=X+AX (1．2) 而对应的整体劲度矩阵变化△K，相应的刚度矩阵K： K=K+AK (1_3) 显然，整体劲度矩阵K是随着x的改变而被决定的。 若设r为结构有x改变后的位移，因为有： Kr= (K +△K)r=p (1．4) 显然，结构修改后，因为K的改变，位移r也随着改变。根本来说是由于x的改变导致了位移r的改 变 。 一 旦位移r确定后，结构个单元的内力f便被确定了，S为结构应力变换矩阵，即： f=sr (1．5) 利用降基法进行结构的重分析，不需直接求解式 (1．4)，而是通过构造 q个基向量的线性组合来近似 逼近结构改变后的位移向量，具体思路如下： 基于 q个位移数组的计算评估下，存在的 ，r2，r3，··．， ，q个位移基 向量 。 r=Y1-I-Y2r2+… 一-I- r =rby (1．6) 式子右边矩阵rb的每一列是由我们选择的唯一基向量：rl，r2，·“， 组成的。 rb==={rl,r2…，rq} (1．7) Y为位移基向量的系数矩阵，表示如下： = {1，2，…．．．，g} (1．8) 将式 (1．1)代入式 (1．4)，等号两边都左乘 ，有： 15 工程地质计算机应用 2012年第4期 总68期 K巾y： P (1．9) 可知，系数矩阵y