曾攀有限元教材matlab代码

%-----failure----direct method （not deal with reaction force）-----------------------
A1=[0.458330.156250.072917;
    0.156250.45833-0.30208;
    0.041667-0.031250.30208]
B1=[10;10;0]
X1=A1B1
X1 =

   16.5557
   15.7426
   -0.6550
%-----pass--direct method （deal with reaction force）------------------------   
A1=[0.458330.156250.072917;
    0.156250.45833-0.30208;
    0.072917-0.302080.60416]
B1=[10;10;0]
X1=A1B1
X1 =

   11.3448
   25.4280
   11.3448
%----pass-Lagrange multipler-----------------------------------  
A2=[0.458330.156250.0416670.031250;
    0.156250.45833-0.03125-0.270830;
    0.041667-0.031250.45833-0.156251;
    0.03125-0.27083-0.156250.45833-1;
    001-10]
B2=[10;10;0;0;0]
X2=A2B2
X2 =

   11.3448
   25.4280
   11.3448
   11.3448
   -3.1051
%----pass---penalty approach------------------- 
C=1000  
A3=[0.458330.156250.0416670.03125;
    0.156250.45833-0.03125-0.27083;
    0.041667-0.031250.45833+C-0.15625-C;
    0.03125-0.27083-0.15625-C0.45833+C]
B3=[10;10;0;0]
X3=A3B3 
X3 =

   11.3441
   25.4296
   11.3441
   11.3472
%-----pass---penalty approach-------------   
C=100000  
A3=[0.458330.156250.0416670.03125;
    0.156250.45833-0.03125-0.27083;
    0.041667-0.031250.45833+C-0.15625-C;
    0.03125-0.27083-0.15625-C0.45833+C]
B3=[10;10;0;0]
X3=A3B3 
X3 =

   11.3448
   25.4280
   11.3448
   11.3448

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-05-16 20:33  曾攀有限元
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEM
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEMBar1D2Node
     文件         299  2008-10-30 08:20  曾攀有限元MATLABFEMBar1D2NodeBar1D2Node_Assembly.m
     文件         201  2008-09-30 11:42  曾攀有限元MATLABFEMBar1D2NodeBar1D2Node_Forces.m
     文件         202  2008-09-30 11:43  曾攀有限元MATLABFEMBar1D2NodeBar1D2Node_Stiffness.m
     文件         218  2008-09-30 11:43  曾攀有限元MATLABFEMBar1D2NodeBar1D2Node_Stress.m
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEMBar2D2Node
     文件         356  2008-09-30 11:56  曾攀有限元MATLABFEMBar2D2NodeBar2D2Node_Assembly.m
     文件         407  2008-09-30 11:48  曾攀有限元MATLABFEMBar2D2NodeBar2D2Node_Forces.m
     文件         438  2008-09-30 11:49  曾攀有限元MATLABFEMBar2D2NodeBar2D2Node_Stiffness.m
     文件         408  2008-09-30 11:50  曾攀有限元MATLABFEMBar2D2NodeBar2D2Node_Stress.m
     文件        1180  2008-09-30 12:01  曾攀有限元MATLABFEMBar2D2Node[典型例题]3_2_5（1）.m
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEMBeam1D2Node
     文件         324  2008-09-30 18:00  曾攀有限元MATLABFEMBeam1D2NodeBeam1D2Node _Deflection.m
     文件         346  2008-09-30 17:58  曾攀有限元MATLABFEMBeam1D2NodeBeam1D2Node_Assembly.m
     文件         291  2008-09-30 17:57  曾攀有限元MATLABFEMBeam1D2NodeBeam1D2Node_Stiffness.m
     文件         341  2008-09-30 17:58  曾攀有限元MATLABFEMBeam1D2NodeBeam1D2Node_Strain.m
     文件         200  2008-09-30 17:59  曾攀有限元MATLABFEMBeam1D2NodeBeam1D2Node_Stress.m
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEMBeam2D2Node
     文件         376  2008-09-30 17:45  曾攀有限元MATLABFEMBeam2D2NodeBeam2D2Node_Assemble.m
     文件         400  2008-09-30 17:44  曾攀有限元MATLABFEMBeam2D2NodeBeam2D2Node_Stiffness.m
     文件         183  2008-10-30 08:31  曾攀有限元MATLABFEMBeam2D2NodeBeamForces.m
     文件         543  2008-10-03 18:38  曾攀有限元MATLABFEMBeam2D2Node[MATLAB算例]3_3_7（2）.asv
     文件         541  2008-09-30 17:51  曾攀有限元MATLABFEMBeam2D2Node[MATLAB算例]3_3_7（2）.m
     目录           0  2018-04-12 22:06  曾攀有限元MATLABFEMFEM2D
     文件         292  2008-10-19 10:29  曾攀有限元MATLABFEMFEM2Dconstrain.txt
     文件         342  2008-10-19 10:28  曾攀有限元MATLABFEMFEM2Delement.txt
     文件         389  2008-10-22 07:18  曾攀有限元MATLABFEMFEM2Delem_stress.dat
     文件        4722  2008-10-19 22:17  曾攀有限元MATLABFEMFEM2DFEM2D_MAIN.m
     文件         292  2008-10-19 16:35  曾攀有限元MATLABFEMFEM2Dforce.txt
............此处省略61个文件信息