平面四节点等参单元matlab实现(3)

图1

三、有限元程序及求解

程序求解使用了matlab语言。具体如下： 程序： clc clear

E=2e11; %弹性模量 NU=0.3; %泊松比 t=0.1; %厚度

X=xlsread('D:\data','nodes'); %读取节点坐标 elem=xlsread('D:\data','elements'); %读取单元编号

w=[1,2,3,4,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]; %有位移