四节点矩形单元的matlab原程序.rar_LX9M_knowledgehnd_neare77_四节点单元_四节点有限元程序

程序说明 plane4 主程序 Coordinate 节点坐标 ME 单元编号 Showmesh 网格图 IS 求单刚在总刚中的位置 Stiffness 求单元刚度矩阵 Assemble 单刚组结成总刚 RestictU 输入位移边界条件 RestrictP 输入力的边界条件 Caculate 边界条件处理，求出力和位移 ShowDistortion 画变形图 Stresses 求单元应力 PStresses 求主应力 Mesh 网格图 Distortion 变形图 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 程序可求解任意边界条件任意力（需手算或增加m文件将力等效到单元节点）的矩形板的单元受力及应力。 运行主程序后，首先弹出网格图，根据网格情况输入位移边界条件及力的边界条件。如提示有[]表示输入一个矩阵。 以一个3*4网格的板为例，左边固支，右边受力。先输入相同位移约束个数5，再输入编号矩阵[1 5 9 13 17]， 接着按提示输入约束矩阵[0 0]（未知的用888表示），全部输入完成后输入0表示完成。力的边界条件输入方式一样， 可多次输入。输入0表示完成。 此次配置相同位移的结点个数　NR=5 配置相同结点位移的结点编号　NU[NR]=[1 5 9 13 17] 结点位移大小　U[NF]=[0 0] NR=0 此次配置结点力的结点个数　NL=2 加结点力的结点编号　NFP[NL]=[4 20] 结点力大小　p[NF]=[12500 0] NL=3 加结点力的结点编号　NFP[NL]=[8 12 16] 结点力大小　p[NF]=[3125 0] NL=0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 结算结果有应力值，节点变形值，网格变形图