XFEM的程序

Matlab写的XFEM程序，用于求解裂纹，应力应变位移等变量

% -------------------------------------------------------------------------
clear all
clc
tic;  
format long
% -------------------------------------------------------------------------

disp（‘                         ‘）
disp（‘****************************************************************‘）
disp（‘****                                                        ****‘）
disp（‘****                      2D XFEM                           ****‘）
disp（‘****           2D XFEM for Fracture Analysis                ****‘）
disp（‘****     1. Considering the geometrical enrichment          ****‘）  
disp（‘****     2. Removing the blending elements                  ****‘）
disp（‘****                    双边缺口拉伸试样                     ****‘）
disp（‘****                                                        ****‘）
disp（‘****************************************************************‘）
disp（‘                                                                ‘） 


addpath E:\XFEM_DOUBLE_EDGE\Mesh
addpath E:\XFEM_DOUBLE_EDGE\other
addpath E:\XFEM_DOUBLE_EDGE\quadrature

global node element elemType E nu


% -------------------------------------------------------------------------
%                                 INPUT
% -------------------------------------------------------------------------

% Dimension of the domain （it is simply a rectangular region L x D）
L = 2.0 ;
D = 4.0 ;

% Material properties
E  = 1000.0 ;
nu = 0.3 ;
stressState=‘PLANE_STRESS‘;

% Compliance matrix C
if （ strcmp（stressState‘PLANE_STRESS‘） ）
    C = E/（1-nu^2）*[1nu0;nu10;000.5*（1-nu）];
else
    C = E/（1+nu）/（1-2*nu）*[1-nunu0;nu1-nu0;000.5-nu];
end

% Loading
sigmato = 1.0;
fac_enrich =3.0; % coefficient for geometrical enrichment（几何加强）

% +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%                          MESHING
% +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

% -------------------------------------------------------------------------
% -
disp（[num2str（toc）‘   MESH GENERATION‘]）
% Number of nodes along two directions
nnx = 50 ;
nny = 100 ;

% Four corner points
pt1 = [00];
pt2 = [L0];
pt3 = [LD];
pt4 = [0D];

% Uniform meshing with Q4 elements
% Data structures for nodes and elements
% Node = [x1 y1
%         x2 y2
%         ...
%         xn yn]
% element = [1 3   5 7
%            4 20 35 78
%            ...
%           ]
% 
elemType = ‘Q4‘ ;
[nodeelement] = meshRectangularRegion（pt1 pt2 pt3 pt4 nnxnnyelemType）;

% compute number of nodes of elements
numnode = size（node1）;  % 节点总数
numelem = size（element1）; % 单元总数

% define essential boundaries
uln = nnx*（nny-1）+1;       % upper left node number
urn = nnx*nny;             % upper right node number
lrn = nnx;                 % lower right node number
lln = 1;                   % lower left node number

topEdge  = [ uln:1:（urn-1）; （uln+1）:1:urn ]‘;
botEdge  = [ lln:1:（lrn-

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        906  2014-12-02 21:09  XFEM程序\XFEM_DOUBLE_EDGE\Mesh\make_elem.m

     文件       1350  2015-01-11 22:14  XFEM程序\XFEM_DOUBLE_EDGE\Mesh\meshRectangularRegion.m

     文件        877  2014-12-02 18:07  XFEM程序\XFEM_DOUBLE_EDGE\Mesh\square_node_array.m

     文件        788  2004-12-17 20:25  XFEM程序\XFEM_DOUBLE_EDGE\Mesh\tricheck.m

     文件       1841  2015-03-08 11:40  XFEM程序\XFEM_DOUBLE_EDGE\other\assembly.m

     文件       1928  2015-01-27 11:20  XFEM程序\XFEM_DOUBLE_EDGE\other\branch.m

     文件        651  2015-01-27 16:42  XFEM程序\XFEM_DOUBLE_EDGE\other\branch_node.m

     文件       1383  2015-02-18 08:11  XFEM程序\XFEM_DOUBLE_EDGE\other\discontQ4quad.m

     文件       1523  2015-02-18 08:14  XFEM程序\XFEM_DOUBLE_EDGE\other\disTipQ4quad.m

     文件       1788  2015-03-10 14:37  XFEM程序\XFEM_DOUBLE_EDGE\other\element_disp.m

     文件        486  2015-02-04 22:50  XFEM程序\XFEM_DOUBLE_EDGE\other\ENRdomainf.m

     文件       1257  2015-02-18 08:14  XFEM程序\XFEM_DOUBLE_EDGE\other\gauss_rule.m

     文件       2983  2015-01-25 10:17  XFEM程序\XFEM_DOUBLE_EDGE\other\G_L_coor.m

     文件        942  2015-03-10 17:15  XFEM程序\XFEM_DOUBLE_EDGE\other\Jdomainf.m

     文件        721  2015-02-18 08:08  XFEM程序\XFEM_DOUBLE_EDGE\other\LS.m

     文件      14477  2015-03-07 23:03  XFEM程序\XFEM_DOUBLE_EDGE\other\node_detect.m

     文件       1402  2015-03-07 23:51  XFEM程序\XFEM_DOUBLE_EDGE\other\posi.m

     文件       8935  2015-03-10 16:26  XFEM程序\XFEM_DOUBLE_EDGE\other\SIF.m

     文件        386  2015-02-10 17:39  XFEM程序\XFEM_DOUBLE_EDGE\other\signed_distance.m

     文件       6128  2015-03-09 23:22  XFEM程序\XFEM_DOUBLE_EDGE\other\xfemBmatrix.m

     文件       7370  2015-01-22 21:51  XFEM程序\XFEM_DOUBLE_EDGE\quadrature\lagrange_basis.m

     文件      18049  2015-01-21 20:09  XFEM程序\XFEM_DOUBLE_EDGE\quadrature\quadrature.m

     文件      10387  2015-03-12 11:43  XFEM程序\XFEM_DOUBLE_EDGE\XFEM_DOUBLE_EDGE.m

     文件        906  2014-12-02 21:09  XFEM程序\XFEM_EDGE_CRACK\Mesh\make_elem.m

     文件       1350  2015-01-11 22:14  XFEM程序\XFEM_EDGE_CRACK\Mesh\meshRectangularRegion.m

     文件        877  2014-12-02 18:07  XFEM程序\XFEM_EDGE_CRACK\Mesh\square_node_array.m

     文件        788  2004-12-17 20:25  XFEM程序\XFEM_EDGE_CRACK\Mesh\tricheck.m

     文件       2839  2015-01-26 09:55  XFEM程序\XFEM_EDGE_CRACK\other\assembly.m

     文件       1928  2015-01-27 11:20  XFEM程序\XFEM_EDGE_CRACK\other\branch.m

     文件        651  2015-01-27 16:42  XFEM程序\XFEM_EDGE_CRACK\other\branch_node.m

............此处省略74个文件信息