双线性四边形等参单元有限元程序.pdfVIP

双线性四边形等参单元有限元程序 本程序采用matlab 编写。程序加载由用户提供的前处理数据，包括网格数据 和载荷数据。采用直接的数值运算和 matlab 符号运算两种方法（可选择）生成 单元刚度矩阵。自动集成结构刚度矩阵，选用直接解法求解线性方程组，解出节 点位移。后处理过程中，程序计算了节点应力值，对于共享节点的应力程序采用 各个单元计算值的平均，程序同时给出了单元最佳应力点的应力值！ 1 使用说明 用户需要在目录中给定的文件中按照既定的格式给出必要的前处理数据（后 面有详细说明）。然后将matlab 的工作目录设为我们提供的“双线性四边形等参 单元程序”目录，然后在命令行中输入main()并按下回车键即可。如图1 所示： 图 1 程序运行说明 1.1 输入数据 程序的输入数据需要在“双线性四边形等参单元程序”目录中给定的文件中 完成。每一个.txt 文件对应一个矩阵，所以数据的输入必要严格按照matlab 矩阵 加载文件的格式完成，即文件的一行对应于矩阵行，数据之间用空格或者逗号隔 开，分号或者换行符表示进入矩阵另一列的输入。如图2：所示输入的矩阵阵为 一4*4 的矩阵： ⎛⎜1 2 4 5⎞⎟ ⎜2 3 5 6 ⎟ ⎜5 6 8 9 ⎟ ⎜ ⎟ ⎝4 5 7 8⎠ 图 2 输入说明 “noteLocation.txt ”为节点坐标数据，a(i,j) 为节点 i 的第 j 个坐 标值，j=1 即为X 坐标，j=2 即为Y 坐标。如图3 所示： 图 3 节点坐标数据示例 “noteNum.txt”为单元节点号矩阵，a(i,j)为第i 个单元的第j 个节点的总体 编号。单元节点按照逆时针编号。对于四节点四边形单元该矩阵为nE*4 的矩阵， nE 为单元数！如图2 所示。 “displacementboundary.txt”为位移边界矩阵，该矩阵为nF*2 的数组，nF 为 总自由度数，dB(i,1)为第i 个节点是否位移边界判据，1 为是，0 为否,dB(i,2)为 第i 个节点常位移值，如果该节点是位移边界则为位移值，否则为0。 “surface load.txt”为表面载荷矩阵，该矩阵为n*6 的矩阵。a(i,1)为第i 个边 界(一个单元的一条边算一个边界条件)所在的单元号，a(i,2)为第i 个边界所在的 边在单元中的编号(编号规则如图4 所示) 。a(i,3)为第i 个边界起始端点（按逆时 ③ ④ ② ① 图 4 边界编号示意图 针方向）上的x 方向的的压力，a(i,4)为第i 个边界起始端点（按逆时针方向）上 的y 方向的的压力，a(i,5)为第i 个边界末端点（按逆时针方向）上的x 方向的的 压力，a(i,6)为第i 个边界末端点（按逆时针方向）上的y 方向的的压力。如图5 所示两个单元的结构，在26 边受均匀横向正压力p ，则其载荷矩阵为： 5 6 ② 3 4 ① 1 2 图 5 载荷矩阵示例