%?-------------------------------------------------------------------------
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（幾何加強(qiáng)）

%?+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%??????????????????????????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）;??%?節(jié)點(diǎn)總數(shù)
numelem?=?size（element1）;?%?單元總數(shù)

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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????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個(gè)文件信息