function u=u()                    % 本程序采用8结点四边形单元，分析固定位移边界条件的弹性力学平面问题   

fprintf('请输入弹性模量E的值(单位均采用国际单位制，下同) \n ');

E=input('');

fprintf('请输入泊松比NU的值\n ');

NU=input('');

fprintf('请输入结构的厚度hou\n ');

hou=input('');

fprintf('请输入分析类型p，平面应力问题请输入1，平面应变问题请输入2\n ');

p=input('');

while p~=1&p~=2

    fprintf('输入有误。请重新输入分析类型p，平面应力问题请输入1，平面应变问题请输入2\n ');

    p=1;

end

if p==1                           % 根据问题类型确定弹性矩阵D

   D=(E/(1-NU*NU))*[1 NU 0; NU 1 0;0 0 (1-NU)/2];

elseif p==2

   D=(E/(1+NU)/(1-2*NU))*[1-NU NU 0;NU 1-NU 0;0 0 (1-2*NU)/2];

end

fprintf('请输入结点个数nn(nn为不小于8的正整数) \n ');

nn=input('');

while nn<8

    fprintf('输入有误，请重新输入结点个数nn(nn为不小于8的正整数) \n ');

    nn=input('');

end

fprintf('请输入单元个数ne(ne为不小于1的正整数) \n ');

ne=input('');

while ne<1