the mechanical response of materials with discrete micro- or meso - structures such as 3d - printed structures , woven textiles , paper , or foams can be modelled with dissipative lattice networks , cf .
e.g.  @xcite .
the main advantage of these models consists in their conceptual simplicity , because lattice springs or beams can be identified as individual fibres or yarns of the underlying structure .
therefore , material parameters such as young s or hardening moduli , and constitutive damage or plasticity laws can be determined in a relatively straightforward manner .
furthermore , lattice networks incorporate large deformations and yarn reorientations rather easily compared to phenomenological continuum models , see e.g.  @xcite .
for heterogeneous cohesive - frictional materials such as concrete , lattice networks are capable of capturing distributed microcracking , the heterogeneity at the microscale , and size effects ; applications and further discussions can be found , e.g. , in  @xcite .
the main drawback of lattice structures is their considerable computational cost  which may well be prohibitive for engineering applications , in which the lattice spacing is generally several orders of magnitude smaller than the problem size .
ergo , multiscale , or reduced - order modelling methods are required to tackle realistic applications .
a reduced - order modelling method with a concurrent multiscale character is the quasicontinuum  ( qc ) method .
it was originally designed for conservative atomistic systems at the nano - scale level , see  @xcite for its initial formulation and  @xcite for various extensions .
this approach was further generalized to tackle materials at the meso - scale level by introducing dissipation along with internal variables in a virtual - power - based format by  @xcite and in a variational format by  @xcite .
in its essence , the qc methodology resolves the underlying lattice only in small regions of interest , whereas it coarsens it elsewhere .
this results in a considerable reduction of the number of degrees of freedom  ( dofs ) , internal variables , and effort to construct the governing equations .
the problem of interest in this contribution is the combination of the qc method and crack propagation in damageable lattice structures , as depicted in fig .
[ fig : motivation ] . because most of the dissipation and fibre reorientation occurs near the crack tip
, this region must be fully resolved . to this end , a mesh indicator that follows the evolution of the crack tip needs to be provided .
such an approach typically leaves a trail of the fully resolved region in the wake of the crack if no coarsening is included , cf .
it may be clear that , apart from allowing the crack to open up , these fully resolved regions do not contribute to the overall physics and accuracy .
it is therefore desirable to coarsen them .
however , this coarsening should take into account the effect which the crack has on the local , coarse - scale kinematics of the problem
i.e . it should allow an arbitrary amount of crack opening without any mechanical resistance .
the aim of this contribution is to provide a refinement _ and _ coarsening strategy that allows for the removal of nodes in the wake of the crack tip to increase the efficiency of the variational qc method .
because the damage in the crack wake is highly localized , it can be treated as a `` discontinuity '' , although the system is obviously discrete .
it is reasonable , therefore , to introduce techniques known from the extended finite element method  ( xfem ) , @xcite , or generalized finite element method  ( gfem ) , @xcite , that rely on the partition of unity  ( pu ) concept , @xcite .
a certain similarity to the particle partition of unity method  ( ppum ) , @xcite , can be also observed .
namely , the space of interpolation functions is enriched to also include `` jumps '' in kinematic as well as in internal variables . as will be explained below , and demonstrated in the examples section , such an enrichment provides an accurate and efficient framework to model cracks in damageable lattices . in
what follows , this approach will be referred to as the _ extended qc _ ( in analogy to continuum systems ) , or _
x - qc _ for short .    several previous works have focused on adaptivity and coarsening in qc methodologies as well , but these mainly dealt with conservative , atomistic lattice models , cf .
an xfem - type qc framework was reported in  @xcite , which treated atomistic models within the scope of molecular dynamics and the bridging domain method ( concurrent coupling between the atomistic and continuum regions ) .
in contrast to these previous works , this contribution focuses on structural lattice networks at the meso - scale level .
moreover , we employ the _ dissipative _ variational qc with damage , including a _ fully non - local _ formulation .
the paper is organized as follows . in section  [
sect : varform ] , we briefly recall the theory of rate - independent systems on which the variational qc is built . for simplicity ,
only time - discrete versions of all equations are provided .
having specified the theoretical background , we can introduce the model at hand , its geometry , state variables , and energies driving the evolution .
the main developments of this contribution , i.e. the interpolation and summation qc steps revisited from the point of view of lattice structures with localized damage , are discussed in section  [ sect : qc ] .
namely , the enrichment functions together with the required changes in the summation rules are introduced in detail .
the numerical solution strategy is addressed in section  [ sect : sol ] .
two examples demonstrate the additional efficiency of the proposed method .
the results show that the coarsening strategy , in particular , allows one to reduce the number of dofs to approximately 25   75% of that of an adaptive qc without coarsening , with no significant increase in error .
the number of dofs reduces to 1
@xmath0 of that of the direct numerical simulations ( dns ) .
the corresponding computing times are reduced by a factor of 4   20 , depending on the particular example .
the paper closes with a summary and conclusions in section  [ sect : conclusion ] .
in this section , the basic principles of the variational formulation of rate - independent systems in the context of lattice networks are recalled . for the sake of brevity ,
we limit ourselves to the time - discrete setting ; the general theory is discussed in  @xcite , whereas applications to continuous systems can be found , e.g. , in  @xcite .
applications to lattice networks and qc are presented in  @xcite for isotropic hardening plasticity and localized damage , respectively .
an admissible configuration of a rate - independent system of interest is specified by a state variable  @xmath1 , where  @xmath2 stores the kinematic variables and  @xmath3 the internal variables .
all admissible configurations are specified by the state space  @xmath4 , @xmath5 , @xmath6 , which incorporates , e.g. , prescribed dirichlet boundary conditions or constraints on damage variables .
the energetic solution  @xmath7 is determined using the following incremental minimization problem at time  @xmath8 @xmath9 with an initial condition  @xmath10 , where  @xmath11 describes a discrete time horizon  @xmath12 $ ] , @xmath13 denotes an incremental energy of the form @xmath14 and the inclusion sign  @xmath15 indicates that the potential  @xmath13 is in general nonsmooth or may have multiple minima .
the incremental energy consists of the total ( helmholtz type ) energy  @xmath16 \times \mathscr{q } \rightarrow \mathbb{r}$ ] , and the dissipation distance  @xmath17 .
the dissipation distance  @xmath18 measures the minimum dissipation by a continuous transition between two consecutive states  @xmath19 and  @xmath20 ; for further details see  @xcite , section  3.2 .
the potential energy reads @xmath21 where  @xmath22 \times \mathscr{q } \rightarrow \mathbb{r}$ ] is the internally stored energy , @xmath23 \rightarrow \mathscr{r}^*$ ] represents the column matrix with external loads , and  @xmath24 is the space dual to  @xmath25 .
each time the minimization problem   is solved in this contribution , a local minimum that satisfies the energy balance @xmath26 is searched .
the energy balance   equates the internally stored energy plus the _ dissipated energy _
@xmath27 with the work performed by the external forces @xmath28^\mathsf{t}[{{\boldsymbol{r}}}(t_l)-{{\boldsymbol{r}}}(t_{l-1 } ) ] .
\label{sect : varform : eq:2}\ ] ] in eqs .  , , and
, the symbol  @xmath29 indicates the dependence on  @xmath30 for  @xmath31 .      in this section
, we briefly specify the geometry , underlying lattice , and the kinematic as well as the internal variables .
our considerations below are limited to 2d lattices for notational simplicity , and lattice nodes are called `` atoms '' .    first , the reader is referred to fig .
[ subsect : geometry : fig:1 ] , where a sketch of the kinematic as well as internal variables is presented .
as indicated , we assume an x - braced lattice with the nearest - neighbour interactions .
the original location vector of atom  @xmath32 is denoted as  @xmath33 .
for all  @xmath34 atoms ( collected in an index set  @xmath35 ) , they are stored in the column  @xmath36^\mathsf{t}$ ] , @xmath37 . by analogy ,
the current location vectors  @xmath38 are stored in  @xmath39^\mathsf{t}$ ] .
the entire system contains  @xmath40 interactions stored in an index set  @xmath41 .
each interaction  @xmath42 that connects atoms  @xmath32 and  @xmath43 is of the length @xmath44 where  @xmath45 denotes the euclidean norm .    for _ damaging _ interactions , cf .
[ subsect : geometry : fig:1b ] , the only internal variable is the damage variable  @xmath46 $ ] .
all of these damage variables are stored in the column  @xmath47 , which is hence of length  @xmath40 .
the internally stored energy associated with damaging lattices reads @xmath48 , \label{subsect : lndamage : eq:1}\ ] ] where  @xmath49 denotes the set of the nearest neighbours associated with atom  @xmath32 .
the elastic energy of a bond  @xmath50 is reduced by the damage variable  @xmath51 only if the bond is loaded in tension to better reflect physical behaviour of damageable materials . in eq .
, we have introduced  @xmath52 , @xmath53 , and assumed  @xmath54 .    the dissipation distance for bond  @xmath55 for two consecutive states  @xmath56 and  @xmath57
is defined as @xmath58 and the total dissipation distance is specified as @xmath59 in definition  , @xmath60 is the energy dissipation of a single bond during a unidirectional damage process up to a damage level  @xmath61 .
this function therefore increases from  @xmath62 to  @xmath63 , where  @xmath64 is the energy dissipated at complete failure .    the incremental interaction energy ,
@xmath65 , and incremental site energy , @xmath66 , then read :    [ subsect : lndamage : eq:6 ] @xmath67
in this section we specify the two steps of the standard qc methodology , interpolation and summation , as well as additional procedures that must be adopted to extend it to an efficient description for localized damage .
we start with a geometric crack description for lattice systems in section  [ subsect : crack ] , proceeding to the interpolation step in section  [ subsect : interpolation ] , where the standard framework as well as the mesh refinement and the mesh coarsening are treated .
the first part closes with the description of the enrichment functions .
the second part is devoted to summation , starting with the standard framework in section  [ subsect : summation ] , which is later extended to incorporate also the enrichment functions .
finally , the complete x - qc procedure is summarized in section  [ subsect : xqc ] in terms of a general algorithm , and energy implications are discussed in section  [ subsect : enimplications ] in order to allow for proper verification of the energy equality  .      in
what follows , a crack is defined as an ordered set of points collected in a set  @xmath68 , cf .
also  @xcite : @xmath69 the set  @xmath70 collects all crack wake atoms , whereas  @xmath68 stores the midpoints of all fully damaged bonds . upon ordering , the set  @xmath68 defines an oriented polygon  @xmath71 ( a curve parametrized by a curvilinear coordinate  @xmath72 ) geometrically representing the crack .
the threshold  @xmath73 should be close enough to  @xmath74 , e.g.  @xmath75 .
naturally , any duplicates in   are eliminated , and the ordering may be implemented in various ways , e.g. by sorting the points in  @xmath68 with respect to their @xmath76-coordinates .    for future use , it is convenient to introduce the signed distance function  @xmath77 , cf . e.g.  @xcite , section  3 .
its definition reads @xmath78 where  @xmath79 denotes the column storing the unit normal to the crack  @xmath80 at the closest point  @xmath81 .
for instance , @xmath79 can be defined such that  @xmath82 , where  @xmath83 is the tangent vector to  @xmath80 at  @xmath81 , and  @xmath84 is the unit basis vector pointing in the @xmath85-direction of the adopted coordinate system . if  @xmath81 is situated at a kink of  @xmath80 , entire cone of normals must be considered ; then , the sign of  @xmath86 is positive if the vector  @xmath87 belongs to the cone of normals and negative otherwise .
the function  @xmath77 is therefore positive on one side of the crack and negative on the other . because the crack tip is always located inside the fully refined region , no second level set function is required to describe the crack tip , cf .
@xcite , section  3.2 . for a pictorial representation of the introduced quantities
see fig .
[ subsect : crack : fig:1 ] .    ,
definition  ) are depicted by red lines ( other interactions are omitted for clarity ) and set  @xmath68 is depicted using green dots .
the signed distance function  @xmath77 is represented by contour lines , and the crack polygon  @xmath80 is presented as a green line .
atom sites are shown as black dots . ]      in the four following subsections , the interpolation step of the qc method and its effect on the incremental energy of the full system ( @xmath13 in  ) are detailed .
interpolation introduces to the minimization problem   the following equality constraints : @xmath88 where  @xmath89 denotes a column of generalized dofs located in a kinematically constrained subspace  @xmath90 of the fully dimensional space  @xmath91 .
after substitution of eq .  , the incremental energy   becomes a function of  @xmath89 , i.e. @xmath92 which can be minimized over the subspace  @xmath90 .
this reduces the computational effort if  @xmath93 .    in qc approaches ,
the generalised dofs  @xmath89 are chosen as the dofs of a small set of atoms , the so - called _
repatoms_. usually , repatoms are chosen as nodes of a mesh  @xmath94 that triangulates the domain  @xmath95 , and the interpolation is implemented via piecewise affine ( @xmath96 ) fe shape functions constructed between the repatoms .
consequently , @xmath97 contains standard fem shape function evaluations at all atom positions , cf .
e.g.  @xcite .
note that higher - order polynomial approximations inside elements can be adopted as well , see e.g.  @xcite , @xcite , or  @xcite .    in order to distinguish classical qc approach from its extended version presented below
, we use the following notation for the classical qc : all  @xmath98 repatoms are collected in an index set  @xmath99 , whereas their admissible positions  @xmath100 , @xmath101 , are stored in the column  @xmath102 by analogy to  @xmath2 .
associated interpolation matrix is denoted  @xmath103 .      in this contribution , the mesh refinement indicator of  @xcite , is adopted with minor changes .
it proceeds as follows : a triangle  @xmath104 of the current triangulation  @xmath105 is endowed with a set of sampling interactions  @xmath106 specified as @xmath107 then , the following energy condition is introduced @xmath108 which specifies interactions that are likely to be damaged . in eq .
, @xmath50 denotes the pair potential , @xmath109 a refinement safety parameter , and  @xmath110 the stored threshold energy at which the internal variable starts to evolve .
the mesh refinement indicator then reads @xmath111 having specified a set of triangles marked for refinement , @xmath112 , the current triangulation  @xmath105 must be refined . to this end
, the backward - logest - edge - bisection algorithm by  @xcite is used .      before coarsening the current triangulation  @xmath105 ,
first an indicator is presented that marks which triangles need to be coarsened . in analogy to eq .
we introduce the following energy condition @xmath113 where the compressive part is now included as well .
the rationale behind condition   is that in the crack wake the stress is low , and hence also the energy .
the mesh coarsening indicator then reads @xmath114 verifying condition   for all elements yields a set of triangles marked for coarsening , @xmath115 .
note that because the damage remains localized , it yields steep spatial ( interaction ) energy variations in the vicinity of the crack tip .
for instance , in fig .
[ subsect : nodeprotection : fig:1 ] , the normalized interaction energy can be seen to drop from  @xmath116 to  @xmath117 over just two to three lattice spacings . from this
it may be clear that if the values of  @xmath118 and  @xmath109 are too close , or if  @xmath118 is too high , mesh oscillations may occur in subsequent time steps . in order to achieve good performance , i.e. to avoid mesh oscillations , @xmath118 should be relatively small compared to  @xmath109 .
in particular , we choose  @xmath119 or  @xmath120 and  @xmath121 in both examples of section  [ sect : examples ] .     around a crack tip ( based on linearly interpolated data ) of a qc system .
regular atom sites are shown as black dots and the existing crack as a green line . ]    in order to coarsen the current triangulation  @xmath105 , a similar algorithm to those proposed by  @xcite and  @xcite is adopted . because we use the two - triangle refinement scheme , the entire binary - tree structure of the triangulation history must be stored .
coarsening of a patch of triangles is decided based on their nodes , cf .
[ subsect : coarsening : fig:1 ] .
all nodes ( repatoms ) of the current triangulation are stored in a set  @xmath122 , which is divided into three disjoint sets : @xmath123 , @xmath124 , and  @xmath125 .
the set  @xmath123 stores nodes associated with triangles in  @xmath126 ( and also other protected nodes , see alg .
[ subsect : xqc : alg:1 ] ) .
the set  @xmath124 collects nodes associated with triangles in  @xmath127 , and  @xmath125 complements  @xmath122 . as a consequence of the mesh hierarchy , only a subset of  @xmath124 can be removed .
these are identified through their valences  @xmath128 ( the number of triangles connected to a node ) : @xmath129 only nodes with valence  @xmath130 ( boundary nodes ) or  @xmath131 ( internal nodes ) can be removed ( the blue points in fig .
[ subsect : coarsening : fig:1 ] ) .
the corresponding two or four child elements are subsequently replaced with one or two parent elements , shifting locally the binary - tree hierarchy of the mesh one level up .
as valences change during the coarsening process , this procedure is repeated until convergence .    .
they are found at , and ahead of , the crack tip , and are protected from coarsening .
the yellow along with the blue points define the set  @xmath124 , and can be coarsened .
the green points remain unchanged and correspond to the set  @xmath125 .
only the blue points satisfy the valence condition and can be removed .
regular atoms are shown as black dots and the current triangulation  @xmath105 is depicted in black thin lines .
the existing crack is indicated as a thick black line .
all interactions are omitted for clarity . ]    in what follows , we use a right - angled initial triangulation  @xmath132 , as already indicated in fig .
[ subsect : coarsening : fig:1 ] . because the above described refinement and coarsening algorithms yield self - similar triangulations , all meshes  @xmath105 are right - angled as well .
this is desirable because the summation error is minimized , cf .
@xcite , and poorly shaped triangles are avoided .    an alternative approach to the mesh refinement and coarsening indicators , proposed in eqs .   and  , would be to use _ error estimators _ such as the _ goal - oriented _ error estimator reported e.g. in  @xcite .
generalizations of this type will not be developed in this contribution and are left as a possible future challenge .      in sections  [ subsect :
genint ]   [ subsect : coarsening ] we discussed the qc interpolation and mesh coarsening independently of an existing crack  the crack served merely to identify the fully resolved region at its tip and a coarse region elsewhere . in this subsection , we describe how to extend the qc interpolation when the current triangulation  @xmath105 is such that a crack cuts through a coarse triangle  @xmath104 , and how it influences  @xmath97 , @xmath133 , and  @xmath89 . similarly to the previous sections , the discussion considers the system at time step  @xmath8 .
we first specify the cut elements ( also called reproducing , @xcite ) , cf . fig .
[ subsect : enrichment : fig:1a ] .
a coarse element  @xmath134 such that @xmath135 is considered to be cut by a crack .
conditions   require that at least one point in  @xmath68 belongs to  @xmath104 and that the signed distance function  @xmath77 changes its sign inside  @xmath104 .
let us assume in what follows that a crack can not cut directly through any of the nodes , atoms , or along domain boundaries .
assume then that our system is approximated with a standard qc framework in which the positions of the repatoms are collected in a column  @xmath133 . then , the reconstruction in eq .   can be rewritten as @xmath136 where  @xmath137 represents the @xmath96 fe shape function associated with the repatom  @xmath43 ( this function is evaluated at the undeformed configuration of atom  @xmath32 ) , and  @xmath138 denotes an admissible position of a repatom  @xmath43 .
this definition does not include an existing crack yet .    in order to incorporate a crack ,
an enrichment using the local pu concept ( standard local xfem approximation ) is introduced .
enrichment terms are included as follows @xmath139}_{\varphi_j^\star({{\boldsymbol{r}}}_0^\alpha ) } \ , \widehat{{{\boldsymbol{g}}}}^\star_j , \quad \alpha \in n_\mathrm{ato } , \quad \beta_j \in n^\star , \label{subsect : enrichment : eq:3}\ ] ] where the same @xmath96 fe shape functions for the extrinsic enrichment terms are used as well .    in eq .
, @xmath140 is a set of  @xmath141 enriched repatoms associated with the cut elements , @xmath142 denotes the enrichment function of the  @xmath143-th repatom , and  @xmath144 denotes an admissible vector of generalized dofs associated with @xmath145-th enriched repatom .
let us emphasize that the only enrichment function is  @xmath146 , depicted in fig .
[ subsect : enrichment : fig:1b ] , and that the approximation in eq .
yields continuous , piecewise linear interpolation of the displacement jump .
note also that no crack tip enrichments are required , since the crack tip is always located inside the fully resolved region , and that the shifting adopted in eq .
guarantees the kronecker-@xmath147 property of the resulting approximation , cf .
@xcite or @xcite .
furthermore , the shifting implies that there are actually no blending elements , as shifted enrichment functions vanish outside the cut elements .
even when shifting is not adopted , no difficulties with blending elements occur , because the interpolation is based on @xmath96 fe shape functions and the enrichment function is piecewise constant .
for further details see  @xcite , sections  4.3 and  5.1.2 .    in continuum xfem ,
poorly conditioned stiffness matrices occur if the crack path is close to triangle edges and nodes , cf .
e.g.  @xcite , or  @xcite . for lattice networks ,
the measure for closeness is the lattice spacing  @xmath148 .
hence , if a crack path is closer than  @xmath148 to an element edge or node , the resulting system matrix is rendered singular , cf .
[ subsect : enrichment : fig:2 ] .
two approaches are used for continua to remedy this situation .
the first approach is the use of clustering ( or node gathering ) implemented through flat - top functions , cf .
e.g.  @xcite .
these techniques usually increase the polynomial order of the resulting approximations if the kronecker-@xmath147 property is required .
hence , this is not desirable within the qc framework , as more extensive changes in the summation rule would be needed .
the second approach is to neglect those  @xmath142 that would render the stiffness matrix singular altogether  ( * ? ? ?
* section  11 and references therein ) , which is adopted in this contribution .
the reason is that such shifted enrichments are identical to zero vector for lattice networks and hence , they provide no additional information to the system .
finally , upon introducing @xmath149 along with @xmath150 \quad \mbox{and } \quad \widehat{{{\boldsymbol{g } } } } = [ \widehat{{{\boldsymbol{r}}}}_\mathrm{rep},\widehat{{{\boldsymbol{g}}}}_1^\star,\dots,\widehat{{{\boldsymbol{g}}}}_{n^\star}^\star]^\mathsf{t } , \label{subsect : enrichment : eq:8}\ ] ] the approximation in eq
.   takes the form of eq .  .
the second step of qc frameworks , summation , is the qc counterpart to numerical integration in the finite element method .
it serves two purposes .
the first one is to estimate the exact system s incremental energy  @xmath13 in an efficient , yet accurate way .
this reduces the computational effort related to determining the energies , gradients , and hessians .
the second purpose is to reduce the dimensionality of the internal variable  @xmath151 , which in combination with the interpolation step realises the reduction of the entire state variable  @xmath152 .
the exact incremental energy  @xmath13 is approximated using a limited number of site- or interaction - energies and their weight factors .
these are typically chosen such that the error of the induced approximation is limited , but a low number of sampling atoms / interactions is required .
the interested reader is referred for further details to  @xcite for atom - based summation rules , to  @xcite for an interaction - based summation rule , or to  @xcite for an overview of summation rules in atomistic systems .    here
, we employ the central site - energy - based summation rule  @xcite , which considers only the atoms at the element corners plus one near the center .
the central one is taken to be representative also of atoms which interact across the element boundaries .
if no internal atoms exist , all boundary atoms are considered .
the central rule involves a certain degree of approximation , but is significantly more efficient than the exact summation rule not discussed here ; for details see  @xcite . the summation rule is implemented by introducing a set of  @xmath153 sampling atoms collected in an index set  @xmath154 .
the  @xmath155 interactions connected to these sampling atoms are stored in an index set  @xmath156 .
then , the summation step yields the following approximation , @xmath157 where  @xmath158 are the weight factors associated with atom sites , and  @xmath159 are corresponding weight factors associated with interactions .    because the approximate incremental energy in eq .   only depends on the internal variables associated with  @xmath155 interactions , one can introduce a reduced state variable of the qc system in the form  @xmath160 , where  @xmath161 is a column that stores the internal variables of all sampling interactions .
the reduced state space then reads  @xmath162 .      if the enrichment of section  [ subsect : enrichment ] is adopted and a coarse triangle is cut by a crack , eq
.   still holds , but the selection of the individual sampling atoms  @xmath163 and their weights  @xmath158 changes .    in the central summation rule ,
all crack wake atoms are treated as additional triangle edges .
therefore , a cut triangle is split into two regions , bounded by triangle edges and one crack face . for each of the two subregions ,
the standard central summation rule is used . hence ,
if no internal atoms exist , all edge atoms are added as discrete sampling atoms .
otherwise , a central sampling atom is selected with the weight corresponding to the sum of the number of internal atoms plus one half of the number of edge atoms .
atoms at edge intersections are added as discrete sampling atoms . for the detailed algorithmic description and pictorial representation see alg .
[ subsect : xsumrule : alg:2 ] and fig .
[ subsect : xsumrule : fig:1 ] .
note that the above - described summation rule , which in its current form does not take into account crack closure ( recall eq .   and
the discussion therein ) , can be easily generalized to deal also with non - monotonous loading and crack closure merely by adding all crack wake atoms  @xmath70 as discrete sampling atoms .
special summation rule that samples crack faces can also be devised , but this is omitted from our further considerations for the sake of simplicity . because only monotonous loading and crack opening are observed in examples section  [ sect : examples ] , the summation rule presented in alg .
[ subsect : xsumrule : alg:2 ] is adequate yet efficient .    1 .
identify crack wake atoms  @xmath164 .
2 .   * for each triangle  @xmath165 * 1 .   *
if * triangle is not cut by a crack , use standard central summation rule  @xcite .
* else * 1 .   identify triangle s internal  @xmath166 , edge  @xmath167 , vertex  @xmath168 , and all  @xmath169 , atoms
. 2 .   select vertex atoms  @xmath170 as discrete sampling atoms ( @xmath171 ) .
if * no internal atoms exist , select all atoms  @xmath169 as discrete sampling atoms  ( @xmath171 ) .
* else * 1 .   split  @xmath169 in sets  @xmath172 and  @xmath173 according to the sign of  @xmath77 .
2 .   define  @xmath174 . split  @xmath175 into two disjoint sets  @xmath176 and  @xmath177 according to the sign of  @xmath77 and sort both with respect to the number of neighbours inside triangle ( in decreasing order ) .
3 .   * if * @xmath178 is empty , all atoms in  @xmath179 are discrete sampling atoms
* else * * choose  @xmath180 as central sampling atom . * for each  @xmath181 , @xmath182 . *
add all  @xmath183 as discrete sampling atoms ( @xmath171 ) .
* identify discrete sampling atoms on edges resulting from neighbouring triangles  @xmath184 .
* for each  @xmath185 , @xmath186 .
* end *     for the central summation rule .
large blue dots represent the sampling atoms , and the crack points  @xmath68 from   are shown as red dots . ]      in all previous sections  [ subsect : crack ]   [ subsect : summation ] , we have assumed an existing crack and triangulation  @xmath105 .
we are , however , interested in evolving crack tips and fully resolved regions .
therefore , we summarize the mesh refinement , coarsening , and additions of enrichments of the total x - qc framework in alg .
[ subsect : xqc : alg:1 ] .    1 .
initialize the system : apply initial condition  @xmath187 and construct initial ( coarse ) mesh  @xmath132 with required information , e.g.  @xmath188 , @xmath189 , @xmath190 , @xmath191 .
2 .   * for  @xmath192
apply the boundary conditions at time  @xmath8 , @xmath193 , @xmath194 , @xmath195 , @xmath196 , @xmath197 , etc .
2 .   equilibrate the unbalanced system , i.e. solve for  @xmath198 using  @xmath199 of eq .  .
update crack description  @xmath77 , @xmath68 , etc .
3 .   for each element
@xmath134 evaluate its refinement indicator in   and construct  @xmath200 . if  @xmath201 refine and update  @xmath105
4 .   protect newly added repatoms from all previous refinements during the current time step  @xmath8 ; if necessary , protect also other repatoms , e.g. those near the crack tip .
use   to identify  @xmath202 . if possible , coarsen and update  @xmath105 .
if *  @xmath105 changed in ( iii )   ( iv ) , update the system information  @xmath203 , @xmath204 , @xmath205 , @xmath206 , etc .
, and return to  ( ii ) since the system is unbalanced . +
* else if * the mesh has converged ; proceed to  ( vi ) .
store relevant outputs : @xmath207 , @xmath105 , @xmath203 , @xmath204 , @xmath205 , @xmath206 , @xmath68 , etc .
* end *      this section discusses the implications of the adaptive scheme summarized in alg .
[ subsect : xqc : alg:1 ] from an energetic point of view , and extends the discussion presented in  @xcite , section  3.5 .
recall that the motivations for energy considerations are threefold .
first , the entire theory presented so far is based on energy minimization and the adaptive procedure should be consistent with these principles .
second , validity of the solution  @xmath152 , as well as of  @xmath207 , is decided based on the energy balance  , which must hold during the entire evolution of the system .
finally , upon taking into account energy implications , the accuracy of the x - qc framework can be assessed from the energetic point of view .
when refinement and coarsening occur . only the  @xmath208 and  @xmath209 components
are shown for clarity . ]
the issue focused on here is illustrated in fig .
[ subsect : enimplications : fig:1 ] . at the end of the previous time step  @xmath210 ,
the system is relaxed for some mesh  @xmath211 .
after the next increment is applied for mesh  @xmath212 and the system is equilibrated ( stage  ( ii ) of alg .
[ subsect : xqc : alg:1 ] ) , refinement and coarsening may be required ( ( ii )   ( iv ) ) . upon convergence , this yields a triangulation  @xmath213 . in order to accurately construct the energy increments , one should distinguish physical increments , computed on mesh  @xmath214 and denoted with subscripts  `` p '' , and changes in energy which are due to the mesh adaptation from  @xmath215 to  @xmath216 , denoted with subscripts  `` a '' .
the artificial energy increments due to mesh adaptation are computed by projecting converged solutions onto an artificial mesh  @xmath217 , which establishes communication between the two systems ( note that these two systems have different numbers of dofs and internal variables ) .
the artificial mesh therefore contains the union of repatoms and sampling interactions present in both systems and hence , it is locally the finest triangulation of the two , cf .
[ subsect : enimplications : fig:1 ] . using the above - described procedure and evaluating energies at time instants
@xmath218 provides energy profiles that we will call _
reconstructed_. for further details see  @xcite , section  3.5 .
let us note that we omit dissipative coarsening ( i.e. the interpolation of the internal variables when the mesh is coarsened ) , and assume that the coarsening occurs only in elastic regions .
this assumption is reasonable because damage is localized .
consequently , all damaged bonds are captured accurately by the special function enrichments and associated summation rule .
in this section , we discuss the minimization of the reduced incremental energy  @xmath219 in  , considered in step  ( ii ) of alg .
[ subsect : xqc : alg:1 ] .
a fixed triangulation  @xmath105 is again assumed .
for the solution of the dns , the interested reader is referred to  @xcite , section  4.1 .
we first reduce the incremental energy @xmath220 see e.g.  @xcite , section  3.1.2 , or  @xcite , and rewrite   as @xmath221 second , is solved by a standard newton algorithm because  @xmath222 is assumed to be smooth enough in  @xmath89 . using the chain rule in the taylor expansion of  @xmath223
, the following stationarity condition is obtained :    @xmath224    where    [ subsect : sol : eq:2 ] @xmath225    the partial derivatives in eqs .
read    [ subsect : sol : eq:3 ] @xmath226    where the reduced site energies , @xmath227 , are introduced in analogy to   by condensing out the internal variables .
.   presents a system of linear equations for increments  @xmath228 , and iterating eqs .   and   until convergence of  @xmath229 yields  @xmath230 .
note that explicit expressions for the gradients and hessians can be found in  @xcite , appendix  a. here , the quantities are derived for lattice structures with isotropic hardening plasticity and hence they have to be multiplied by  @xmath231 if a bond is loaded in tension .
moreover , the dependencies on all plastic internal variables have to be dropped .    before the gradient and hessian in eqs .
are computed , the approximate incremental energy must be minimized with respect to  @xmath232 according to eq .  . to this end , the interaction - based formulation of the incremental energy specified in eq .
the minimization then results in  @xmath155 independent one - dimensional box - constrained minimization problems @xmath233 for further details see  @xcite , section  4.1 and appendix  a.
in this section the framework is applied to two examples discussed already in  @xcite , where coarsening was not considered .
consequently , we can make a direct comparison in terms of accuracy and efficiency between the adaptive qc frameworks with and without coarsening .
the same material model is used for each lattice spring ( the superscripts  @xmath55 are dropped for brevity ) : @xmath234 where  @xmath235 is the young s modulus , @xmath236 the cross - sectional area , and  @xmath237 the bond strain .
the dissipation function  @xmath238 in eq .   is adopted from  @xcite , appendix  a , providing an exponential softening rule @xmath239 in eq .
, @xmath240 describes the softening branch of the associated stress - strain diagram , @xmath241 is the bond stress , @xmath242 measures the initial slope of the softening branch , and  @xmath243 is the limit elastic strain , cf .
[ sect : examples : fig:1 ] .
the employed parameters are specified in tab .
[ sect : examples : tab:1 ] .     and  @xmath244 , cf .
the elastically stored energy , @xmath245 , corresponds to the area of the dashed triangle , and the dissipated energy , @xmath246 , to the area of the upper triangle with the curved edge . ]
.dimensionless material and geometric parameters for both examples .
[ cols="<,>,<,>,<",options="header " , ]     [ sect : examples : tab:1 ]       denotes the applied vertical displacement . ]
the first example focuses on an l - shaped plate test , see fig .
[ subsect : simpleex : fig:0 ] .
the lattice properties are homogeneous except for the vicinity of  @xmath247 , where a vertical displacement is applied . here , in order to prevent any damage , the young s modulus is  @xmath248 times larger than elsewhere and the limit elastic strain  @xmath243 is infinite .
the reference domain comprises  @xmath249 atoms , @xmath250 interactions , and is fixed at the bottom part of the boundary , @xmath251 .
the evolution of the system is controlled by the crack mouth opening displacement ( cmod ) through a pair of crack tip atoms highlighted in fig .
[ subsect : simpleex : fig:0 ] .
the numerical simulation is performed for two x - qc systems , each with a different safety margin  @xmath109 , cf .
the first system is referred to as the _ moderate x - qc _ ( @xmath119 ) and the second one to as the _ progressive x - qc _ ( @xmath252 ) ; in both cases , @xmath121 , recall eq .  .
relatively low safety margins are necessary due to the highly fluctuating interaction energies near the crack tip ( recall fig .
[ subsect : nodeprotection : fig:1 ] and the discussion therein ) .
the results are compared to those of the adaptive version of the variational qc , reported in  @xcite , for the same choices of the refinement safety margins , i.e.  @xmath119 and  @xmath252 , and for the central summation rule .
these systems are referred to as the _ moderate adaptive qc _ and the _ progressive adaptive qc_. for completeness , results of the dns are also shown .    the deformed configuration computed by the dns at  @xmath253
is depicted in fig .
[ subsect : simpleex : fig:1a ] .
the crack path initiates at the inner corner ( @xmath254 ) and propagates to the left .
this is also predicted by the qc frameworks . in fig .
[ subsect : simpleex : fig:1b ] , the reaction force  @xmath255 is presented as a function of  @xmath256 .
the responses predicted by the x - qc are practically identical to those of the adaptive qc without coarsening .
( ldeformed ) at ( 0,0 ) ; ( lreact ) at ( 6.5,0 ) ; ( lzoom ) at ( 11.5,-0.72 )   predicted by the dns , ( b )  the force - displacement diagram for the reaction force  @xmath255 as a function of  @xmath256 . the instants at which mesh refinement and coarsening occur for the first time are indicated by dashed vertical lines of the respective colours.,title="fig : " ] ;    ( 8,-1.1 ) rectangle ( 8.7,-0.6 ) ;    ( 10.3,-1.3 ) rectangle ( 12.7,0 ) ;    ( 8.35,-0.6 )  ( 10.3,0 ) ; ( 8.35,-1.1 )  ( 10.3,-1.3 ) ;    ( text ) at ( 11.5,0.2 ) zoom in ;    the energy profiles appear in fig .
[ subsect : simpleex : fig:2a ] as functions of  @xmath256 .
the energy balance   is clearly satisfied along the entire loading path , i.e. the thin dotted lines presenting the work performed by external forces  @xmath257 lie on top of the thick dashed lines representing  @xmath258 .
second , similarly to the reaction force profiles , we notice that the energy paths are almost indistinguishable from those of adaptive qc , demonstrating good accuracy of the x - qc approach . in fig .
[ subsect : simpleex : fig:2b ] , we present the energy components that are exchanged ( the artificial energies ) in the moderate x - qc simulation . their high magnitudes compared to the curves of fig .
[ subsect : simpleex : fig:2a ] indicate the importance of the energy reconstruction procedure , discussed in section  [ subsect : enimplications ] .
[ subsect : simpleex : fig:3a ] presents the relative number of repatoms as a function of  @xmath256 .
note that the number of enriched repatoms  @xmath141 , defined in eq .
, is negligible ( in particular , @xmath259 ) .
both x - qc systems initially follow the trends of the corresponding adaptive qc systems ( the refinement stage and crack localization ) . after the crack localizes and starts to propagate however , the crack wake unloads and the mesh coarsens , resulting in a substantial drop in the number of repatoms of the progressive x - qc . at the end of the simulation , the specimen is almost stress - free , allowing for a very coarse representation .
note that the number of dofs associated with the x - qc approach is less than one half compared to the adaptive qc in the case of the moderate approach , and approximately a quarter in the case of the progressive approach . in both cases ,
the final relative number of dofs as well as repatoms remain below  @xmath260 , whereas they never exceed  @xmath261 during the simulation .
a similar behaviour can be observed also for the relative number of sampling atoms , presented in fig .
[ subsect : simpleex : fig:3b ] .
finally , figs .
[ subsect : simpleex : fig:4 ] and  [ subsect : simpleex : fig:5 ] present some of the meshes of the x - qc and the adaptive qc frameworks associated with fig .
[ subsect : simpleex : fig:3 ] .
although at the initial stages the progressive and moderate x - qc meshes differ substantially , they are almost identical when the crack localizes .
the effect of coarsening in both cases is evident .
note that in the progressive case coarsening occurs not only along the crack path , but also along the right edge of the lower part of the plate .
this explains the particularly large gain in number of repatoms in this case , i.e. the factor of four visible in fig .
[ subsect : simpleex : fig:3a ] .     as a function of  @xmath256 , and ( b )  relative number of sampling atoms.,title="fig : " ] +     +     +          the second example is a four - point bending test , presented e.g. in  @xcite , cf . also fig .
[ subsect : complexex : fig:1 ] .
the rectangular reference domain  @xmath95 is pre - notched at the top to initiate a crack .
the otherwise homogeneous x - braced lattice is again locally stiffened where the dirichlet and neumann boundary conditions are applied ( the young s modulus is  @xmath248 times larger than elsewhere and the limit elastic strain  @xmath243 is infinite to prevent damage ) .
the domain comprises  @xmath262 atoms connected by  @xmath263 interactions and is loaded by a pair of vertical forces  @xmath264 and  @xmath265 , given by @xmath266 where  @xmath267 represents the indirect displacement control parameter , see e.g.  @xcite , section  22.2.3 . unlike the l - shaped plate example , the sum of cmod and cmsd ( crack mouth sliding displacement ) is now used to control the simulation . in fig .
[ subsect : complexex : fig:2a ] the loading program is specified .    again two x - qc systems are studied : the moderate x - qc with  @xmath119 , and the progressive x - qc with  @xmath252 .
the results are again compared to those of the adaptive qc ( with the same values of  @xmath109 ) and to those of the dns .
the deformed configuration predicted by the dns at  @xmath268 is depicted in fig .
[ subsect : complexex : fig:2b ] .
the x - qc crack paths , shown in fig .
[ subsect : complexex : fig:3a ] , are practically identical to those of the adaptive qc ( the maximum distance amounts to one lattice spacing ) .
they also match the crack path predicted by the dns , where the distance does not exceed four lattice spacings in all cases . the same quality of approximation is achieved also for the total applied force , @xmath269 , plotted in fig .
[ subsect : complexex : fig:3b ] against  @xmath270 .
some minor differences can be observed between the responses of the x - qc and the adaptive qc around the peak load ( when the crack localizes ) . at this stage ,
the x - qc meshes start to coarsen . at later stages ,
all curves become practically indistinguishable as the specimen becomes stress - free .
( bcrack ) at ( 0,0 ) ; ( breact ) at ( 9,0.2 ) ; ( breactzoom ) at ( 4.2,0.7 )    versus  @xmath270 .
note that the colour scheme applies to both diagrams .
the instants at which mesh refinement and coarsening occur for the first time are indicated by vertical dashed lines of the respective colours.,title="fig : " ] ;    ( 7,2.3 ) rectangle ( 7.7,1.5 ) ;    ( 2.8,-0.8 ) rectangle ( 5.6,2 ) ;    ( 7.35,2.3 )  ( 5.6,2 ) ; ( 7.35,1.5 )  ( 5.6,-0.8 ) ;    ( text ) at ( 4.2,2.2 ) zoom in ;    the energy profiles for all employed approaches are presented in fig .
[ subsect : complexex : fig:4a ] .
the energy balance   clearly holds .
one can also notice that the x - qc energies match the adaptive qc and dns energies well . in fig .
[ subsect : complexex : fig:4b ] the energy exchanges of the moderate x - qc approach are shown . again , like in the previous example , their contributions are substantial .    in fig .
[ subsect : complexex : fig:5a ] , the relative number of repatoms of the qc schemes are shown as a function of  @xmath270 . the number of enriched repatoms is again negligible compared to the number of repatoms ( @xmath271 ) . before the crack localizes , the number of repatoms rapidly increases .
but it drops significantly after crack localization , and converges to approximately  @xmath272 for both x - qc approaches .
similar behaviour can be observed for the relative number of sampling atoms shown in fig .
[ subsect : complexex : fig:5b ] .
this result , in combination with the accuracy captured by figs .
[ subsect : complexex : fig:3 ] and  [ subsect : complexex : fig:4a ] , clearly demonstrates the strength of the x - qc framework . finally , in figs .
[ subsect : complexex : fig:6 ] and  [ subsect : complexex : fig:7 ] several snapshots capturing the evolution of the triangulations  @xmath105 are presented .     as a function of  @xmath270 , and ( b )  relative number of sampling atoms.,title="fig : " ] +     +   +   +     +   +   +
in this contribution , an adaptive mesh algorithm for the variational qc for lattice networks with localized damage is developed , which refines the mesh where necessary , but also coarsens it where possible . to achieve the latter , xfem - like
enrichment functions are required .
it is shown that the efficiency of the qc framework applied to lattices with localised damage increases significantly .
the main results can be summarized as follows :    1 .
the two qc steps , interpolation and summation , were generalized to incorporate special enrichment functions , reflecting cracks within the framework of the fully - nonlocal variational qc .
2 .   to determine regions
to be coarsened , a heuristic marking strategy was proposed with a variable parameter controlling the mesh coarseness .
the mesh refinement and coarsening were discussed from an energetic point of view .
4 .   for the investigated examples
the extended qc methodology , including coarsening , reduces the number of repatoms and sampling atoms with at least  @xmath273 of that of a simpler adaptive qc method with refinement only .
the accuracy remains the same , however .
one can note that the employed coarsening procedure neglected the diffusive dissipation .
this was justified for the application to localized damage , but if more ductile behaviour is modelled ( due to elastic - plastic models with damage ) , this must be incorporated .
therefore , a further generalization is required .
moreover , instead of the heuristic marking strategies used here , goal - oriented error estimators may be used to further improve the accuracy and efficiency of the qc methodology .
the variational foundation of the method provides an ideal platform for this , cf .
e.g.  @xcite .
the same x - qc framework can also be employed for the description of heterogeneities in otherwise homogeneous lattice networks .
these aspects are , however , outside the scope of this contribution and will be reported separately .
this work was supported by the czech science foundation ( gar ) , through project no .  14 - 00420s .
75 natexlab#1#1[1]`#1 ` [ 2]#2 [ 1]#1 [ 1]http://dx.doi.org/#1 [ ] [ 1]pmid:#1 [ ] [ 2]#2 , , , a. . , .
http://www.sciencedirect.com/science/article/pii/s0045782516301050 , .
, , , , b. . , .
http://dx.doi.org/10.1002/nme.4982 , . .
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s0022509615000630 , .
, , , . . ,
http://dx.doi.org/10.1002/nme.2542 , .
, , , . . ,
http://dx.doi.org/10.1002/nme.2675 , .
http://www.sciencedirect.com/science/article/pii/s0045782511003082 , .
http://dx.doi.org/10.1002/(sici)1097-0207(19970228)40:4<727::aid-nme86>3.0.co;2-n[http://dx.doi.org/10.1002/(sici)1097-0207(19970228)40:4<727::aid-nme86>3.0.co;2-n ] , .
http://journals.cambridge.org/article_s0962492901000010 , .
, , , , a. . , .
http://www.sciencedirect.com/science/article/pii/s0045782514002047 , .
http://dx.doi.org/10.1002/nme.3134 , .
, , , b. . , .
http://www.sciencedirect.com/science/article/pii/s0022509614001100 , .
, , , c. . , .
http://www.sciencedirect.com/science/article/pii/s0022509613002445 , .
, , , d. . , .
http://www.sciencedirect.com/science/article/pii/s004578251300279x , .
, , , , a. . , .
http://www.sciencedirect.com/science/article/pii/s0167663614001495 , .
, , , , b. . , .
http://dx.doi.org/10.1002/gamm.201510018 , .
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s1359835x13000134 , .
http://dx.doi.org/10.1002/(sici)1097-0207(19990620)45:5<601::aid-nme598>3.0.co;2-s[http://dx.doi.org/10.1002/(sici)1097-0207(19990620)45:5<601::aid-nme598>3.0.co;2-s ] , .
, , , , . . ,
http://dx.doi.org/10.1002/1097-0207(20010210)50:4<993::aid-nme164>3.0.co;2-m[http://dx.doi.org/10.1002/1097-0207(20010210)50:4<993::aid-nme164>3.0.co;2-m ] , .
, , , a. . , .
http://www.sciencedirect.com/science/article/pii/s016766361500157x , .
, , , b. . , .
http://www.sciencedirect.com/science/article/pii/s0020768314004600 , . , .
. , . . , , ,
http://www.sciencedirect.com/science/article/pii/s0022509699000289 , .
, , , . . ,
http://epubs.siam.org/doi/abs/10.1137/080741033 , .
, , , . . ,
, , . . , . , ,
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s0045782505003956 , .
, , , , . . ,
http://www.sciencedirect.com/science/article/pii/s0013794415000053 , . , .
http://dx.doi.org/10.1002/nme.2259 , .
http://dx.doi.org/10.1002/nme.3299 , . , , .
http://dx.doi.org/10.1002/nme.2914 , .
, , , . . ,
http://dx.doi.org/10.1002/nme.2488 , .
http://www.sciencedirect.com/science/article/pii/s0020768309004752 , .
http://dx.doi.org/10.1137/s1064827599355840 , , http://arxiv.org/abs/http://dx.doi.org/10.1137/s1064827599355840 [ arxiv : http://dx.doi.org/10.1137/s1064827599355840 ] . , , . .
http://dx.doi.org/10.1007/s10704-012-9753-8 , . , , .
http://www.sciencedirect.com/science/article/pii/s0022509610002425 , .
http://www.sciencedirect.com/science/article/pii/s0020768315002620 , .
http://www.sciencedirect.com/science/article/pii/s0167663612000683 , .
, , , , . . ,
http://dx.doi.org/10.1002/nme.1370 , .
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s016784421000039x , .
http://journals.cambridge.org/article_s0962492913000068 , .
http://www.sciencedirect.com/science/article/pii/s0045782596010870 , .
, , , . . ,
http://dx.doi.org/10.1007/s00466-015-1127-4 , .
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s004578251400423x , . , .
http://dx.doi.org/10.1007/s00161-003-0120-x , .
http://www.springer.com/us/book/9781493927050 , . , , , .
http://dx.doi.org/10.1007/s002050200194 , .
http://dx.doi.org/10.1023/a%3a1026098010127 , .
, , , . . ,
http://dx.doi.org/10.1002/(sici)1097-0207(19990910)46:1<131::aid-nme726>3.0.co;2-j[http://dx.doi.org/10.1002/(sici)1097-0207(19990910)46:1<131::aid-nme726>3.0.co;2-j ] , .
http://www.sciencedirect.com/science/article/pii/002076839190004y , .
http://www.sciencedirect.com/science/article/pii/s0021999102971834 , .
http://www.sciencedirect.com/science/article/pii/s1359835x04002593 , .
http://www.sciencedirect.com/science/article/pii/s002250961100055x , .
, , , . . ,
http://www.sciencedirect.com/science/article/pii/s0045782598002308 , . , , , .
http://www.sciencedirect.com/science/article/pii/s0022509610001341 , . , .
. , . . , , , , a.
. http://arxiv.org/abs/1601.00625[arxiv:1601.00625 ] .
, , , , b. . http://arxiv.org/abs/1601.00625[arxiv:1601.00625 ] . , . .
ph.d . thesis .
technische universiteit delft . .
http://www.sciencedirect.com/science/article/pii/095894659290004f , .
, , , , , , . . ,
http://www.sciencedirect.com/science/article/pii/s0022509698000519 , .
http://www.sciencedirect.com/science/article/pii/s0045782599000729 , .
, , , b. . , .
http://dx.doi.org/10.1002/(sici)1097-0207(20000320)47:8<1401::aid-nme835>3.0.co;2-8[http://dx.doi.org/10.1002/(sici)1097-0207(20000320)47:8<1401::aid-nme835>3.0.co;2-8 ] , .
, , , . . ,
http://dx.doi.org/10.1080/01418619608243000 , , http://arxiv.org/abs/http://dx.doi.org/10.1080/01418619608243000 [ arxiv : http://dx.doi.org/10.1080/01418619608243000 ] . , , , , , . .
http://dx.doi.org/10.1002/nme.2387 , .
http://dx.doi.org/10.1002/nme.143 , .
http://www.sciencedirect.com/science/article/pii/s0045782514003545 , .