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 (责任编辑:泉水) |