matab 在结构力学中的应用.doc

PAGE Matlab在结构力学中的应用 马晓光，于国清 （白城师范学院　机电系，白城师范学院　土木系） 摘　要　　本文结合实例，介绍了Matlab在结构力学计算中的应用方法，展示了Matlab在结构有限元编程中的简捷性与高效性。本文介绍的方法，在结构有限元计算中具有普遍的适用性。 关键词　　Matlab　 有限元 　编程 APPLICATION OF MATLAB IN STRUCTURAL MECHANICS Ma Xiaoguang，Yu Guoqing ( Baicheng Normal College, Baicheng Jilin 137000, China ) Abstract This paper introduces a method of application of Matlab in structural FEM, taking a continuous beam as an example. The method indicates the neatness and high efficiency of Matlab in programming a structural FEM and can be applied to many kinds of structures. Keywords FEM; Matlab; programming; １　引　言 　　Matlab既是强大的数学软件，又是一种简洁、高效的科技编程语言，这一点在结构力学有限元计算中得到了充分的体现。Matlab最突出的特点就是处理矩阵的各种运算，而大量的矩阵运算正是结构有限元计算非常需要的。此外，利用Matlab的绘图功能，还可直接绘出结构的内力图。因此，Matlab对于结构力学计算是十分有利的工具。 ２　MATLAB编程实例 单元编号 节点i 节点j ① 1 2 ② 2 3 ③ 3 4 10kN/m30kN3m3m2m2m图1１　　　　　　　　　　２　　　　　　　　 ３　　　　　　　　４ 10kN/m 30kN 3m 3m 2m 2m 图1 １　　　　　　　　　　２　　　　　　　　 ３　　　　　　　　４ ① ② ③ 　　该结构可以用三弯矩方程求解，也可以用矩阵位移法求解。前者仅适用于连续梁结构，而后者除了适用于连续梁之外，还能适用于其他类型的结构，应用范围很广，更具有普遍性。因此，以下根据矩阵位移法编写计算程序，以展示Matlab在结构力学中应用的一般方法。 　　第一步，给整体各节点编号，如图1所示。各单元的编号及其连通性见表格。 　　第二步，利用Matlab函数建立各单元的刚度阵。该函数为 function y=BeamES(x) % x为梁的原始参数矩阵 y = x(1)*x(2)/(x(3)*x(3)*x(3))*[12 6*x(3) -12 6*x(3) ; 6*x(3) 4*x(3)*x(3) -6*x(3) 2*x(3)*x(3) ; -12 -6*x(3) 12 -6*x(3) ; 6*x(3) 2*x(3)*x(3) -6*x(3) 4*x(3)*x(3)]; 　　第三步，建立整体刚度阵。因为该结构有4个节点，每个节点有两个未知数（可考虑支座沉降），所以，整体刚度阵为8×8的方阵。首先利用下面的语句生成一个8×8的0矩阵： K=zeros(8), 然后，反复调用下面的函数，便可由单元刚度阵生成整体刚度阵： function y = BeamAssemble(K,k,i,j) % 其中K为整体刚度阵，k为单元刚度阵， ％ i，j为单元两端在整体节点上的编号。 K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1); K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2); K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3); K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4); K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1); K(2*i,2*i) = K(2*i,2*i) + k(2,2); K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3); K(2*i,2*j) = K(2*i,2*j) + k(2,4); K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1); K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2); K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3); K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4); K(2