程序说明.txt

四节点矩形单元matlab有限元程序

                                           程序说明


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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


结算结果有应力值，节点变形值，网格变形图