精选优质文档-倾情为你奉上平面四节点单元的刚度矩阵的fortran编程parameter(E0=2.0e06,u=0.3,h=2) !参数的选定real,dimension(2,2): kii,kij, kji,kjj, kim,kin, kjm,kjn,kmi,kmj,kmm,kmn,kni,knj,knm,knnreal,dimension(3,8):B !单元总应变矩阵real,dimension(2,4):z ！节点的坐标real,dimension(3,2):Bi,Bj,Bm,Bn ！单元应变矩阵real,dimension(8,8):Kee !单元刚度矩阵real,dimension(3,3):D !矩阵 D为 33矩阵的定义 real,x,y data z/xi,yi,xj,xm,ym,xn,yn/ a=z(1,2)-z(1,1)=xj-xi b=z(2,3)-z(2,2)=ym-yj Ni=0.25*(1+xi*x/a)*(1+yi*
