曾攀有限元教材matlab代码

清华大学曾攀版本有限元教材以及书中例题对应的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=A1\B1
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=A1\B1
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=A2\B2
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=A3\B3 
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=A3\B3 
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  曾攀有限元\MATLABFEM\Bar1D2Node\
     文件         299  2008-10-30 08:20  曾攀有限元\MATLABFEM\Bar1D2Node\Bar1D2Node_Assembly.m
     文件         201  2008-09-30 11:42  曾攀有限元\MATLABFEM\Bar1D2Node\Bar1D2Node_Forces.m
     文件         202  2008-09-30 11:43  曾攀有限元\MATLABFEM\Bar1D2Node\Bar1D2Node_Stiffness.m
     文件         218  2008-09-30 11:43  曾攀有限元\MATLABFEM\Bar1D2Node\Bar1D2Node_Stress.m
     目录           0  2018-04-12 22:06  曾攀有限元\MATLABFEM\Bar2D2Node\
     文件         356  2008-09-30 11:56  曾攀有限元\MATLABFEM\Bar2D2Node\Bar2D2Node_Assembly.m
     文件         407  2008-09-30 11:48  曾攀有限元\MATLABFEM\Bar2D2Node\Bar2D2Node_Forces.m
     文件         438  2008-09-30 11:49  曾攀有限元\MATLABFEM\Bar2D2Node\Bar2D2Node_Stiffness.m
     文件         408  2008-09-30 11:50  曾攀有限元\MATLABFEM\Bar2D2Node\Bar2D2Node_Stress.m
     文件        1180  2008-09-30 12:01  曾攀有限元\MATLABFEM\Bar2D2Node\[典型例题]3_2_5（1）.m
     目录           0  2018-04-12 22:06  曾攀有限元\MATLABFEM\Beam1D2Node\
     文件         324  2008-09-30 18:00  曾攀有限元\MATLABFEM\Beam1D2Node\Beam1D2Node _Deflection.m
     文件         346  2008-09-30 17:58  曾攀有限元\MATLABFEM\Beam1D2Node\Beam1D2Node_Assembly.m
     文件         291  2008-09-30 17:57  曾攀有限元\MATLABFEM\Beam1D2Node\Beam1D2Node_Stiffness.m
     文件         341  2008-09-30 17:58  曾攀有限元\MATLABFEM\Beam1D2Node\Beam1D2Node_Strain.m
     文件         200  2008-09-30 17:59  曾攀有限元\MATLABFEM\Beam1D2Node\Beam1D2Node_Stress.m
     目录           0  2018-04-12 22:06  曾攀有限元\MATLABFEM\Beam2D2Node\
     文件         376  2008-09-30 17:45  曾攀有限元\MATLABFEM\Beam2D2Node\Beam2D2Node_Assemble.m
     文件         400  2008-09-30 17:44  曾攀有限元\MATLABFEM\Beam2D2Node\Beam2D2Node_Stiffness.m
     文件         183  2008-10-30 08:31  曾攀有限元\MATLABFEM\Beam2D2Node\BeamForces.m
     文件         543  2008-10-03 18:38  曾攀有限元\MATLABFEM\Beam2D2Node\[MATLAB算例]3_3_7（2）.asv
     文件         541  2008-09-30 17:51  曾攀有限元\MATLABFEM\Beam2D2Node\[MATLAB算例]3_3_7（2）.m
     目录           0  2018-04-12 22:06  曾攀有限元\MATLABFEM\FEM2D\
     文件         292  2008-10-19 10:29  曾攀有限元\MATLABFEM\FEM2D\constrain.txt
     文件         342  2008-10-19 10:28  曾攀有限元\MATLABFEM\FEM2D\element.txt
     文件         389  2008-10-22 07:18  曾攀有限元\MATLABFEM\FEM2D\elem_stress.dat
     文件        4722  2008-10-19 22:17  曾攀有限元\MATLABFEM\FEM2D\FEM2D_MAIN.m
     文件         292  2008-10-19 16:35  曾攀有限元\MATLABFEM\FEM2D\force.txt
............此处省略61个文件信息