一 问题描述如图所示的平面矩形结构，设E=1，NU=0.25，h=1，考虑以下约束和外载：位移边界条件BC（u）：UA=0，VA=0，UD=0，力边界条件BC（p）：在CD边上有均布载荷q=1，建模情形：使用四个四节点矩形单元，试在该建模情形下，求各节点的位移以及各个单元的应力分布。二 Matlab程序(1).函数定义：function k= Quad2D4Node_Stiffness(E,NU,h,xi,yi,xj,yj,xm,ym,xp,yp,ID) syms s t; a = (yi*(s-1)+yj*(-1-s)+ym*(1+s)+yp*(1-s)/4; b = (yi*(t-1)+yj*(1-t)+ym*(1+t)+yp*(-1-t)/4; c = (xi*(t-1)+xj*(1-t)+xm*(1+t)+xp*(-1-t)/4; d = (xi*(s-1)+xj*(-1-s)+xm*(1+s)+xp*(1-s)/4; B1 = a*(t-1)/4-b*(s-1)/4 0 ; 0 c*(s-1)/4-
