基于MATLAB的四结点等参元的编程.pdfVIP

of InstituteofUrbanConstructionV01．12No．4Dee．2006 JournalTi粕jin 辑F％岛嘞岛怒 铲土木3-程0 叁‘耷4，窜or窜6碡攻 基于MATLAB的四结点等参元的编程 栾小兵，李自林，朱斌 (天津城市建设学院土木工程系，天津300384) 摘要：基于有限元理论，用MATLAB语言编制了四结点等参元计算机程序．计算了悬臂双向 板在均布荷栽与集中力共同作用下的位移、应力和应变，并与大型通用有限元软件ANSYS计 算结果进行了比较，二者吻合较好，验证了本程序编制的正确性．利用MA’眦B编程对结构 进行分析计算比FORTRAN、c语言简便，具有精度高的优点． 关键词：MATLAB；四结点等参元；悬臂双向板；编程 中图分类号：TP314 文献标识码：A 虽然MATLAB¨1具有很强大的功能，但在土木工所示． 程中的应用很少．笔者查阅到几篇有关应用该软件研 究土木工程中问题的文献．东南大学的毛利军、李爱 ③ 群旧刮运用MATLAB中SIMULINK平台进行了基础滑 移隔震结构仿真分析研究，还进行了建筑结构双线性 基础隔震体系地震反应的参数研究．西北工业大学的 雷开贵和李永双Mo应用MATLAB仿真技术进行了建 筑结构的风振控制的研究．重庆大学的李英民哺’应用 MATLAB软件中的SIMULINK平台进行了框架结构非 图1单元尺寸 线性地震反应研究．徐赵东、郭迎庆∞1运用MATLAB 软件进行了被动控制结构的SIMULINK动态仿真研设f，叼是自然坐标．Ⅳ，，Ⅳ2，Ⅳ3，Ⅳ4是自然坐标系 究，并出版了专著“MATLAB语言在建筑抗震工程中 下的4个形函数． 的应用”，该书论述了SIMULINK工具箱、神经网络工因为此单元是等参元，故得 具箱、模糊逻辑工具箱在抗震工程中的应用． 4 本文以编写四节点等参元的实例，说明了在有限 省=∑并i札(亭，叼)1 ”1 I (1) 元编程中，MATLAB比FORTRAN和C语言具有更大 的优越性． Y=∑)，；以(f，叩)j i=1 式中：石和Y分别为迪卡尔坐标系下的横坐标和 1 四节点等参元的理论 纵坐标；搿；和Yi分别为节点i的横坐标和纵坐标． 1．2应变-位移矩阵 1．1 单元简介 在求单元刚度矩阵中，可以根据弹性力学理论和 该单元可用于二维平面单元的模拟，既可解决平 有限元理论推导出单元的应变一位移矩阵曰，具体推导 面应变问题，又可以解决平面应力问题．该单元是4个 过程见文献[7]，则 节点，每个节点有两个自由度．它的几何尺寸如图l 收稿日期：2006-06-28；修订日期：2006-09-06 作者简介：栾小兵(1981一)，男，江苏泰兴人，天津城市建设学院硕士生 万方数据 天津城市建设学院学报 2006年第12卷第