27.12.2012 Views

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

ARUP; ISBN: 978-0-9562121-5-3 - CMBBE 2012 - Cardiff University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

used the PAC constitutive constants (obtained from the study by Jin et al. [15]) that<br />

correspond to the tissue relaxed state: 1) Matrix Young’s modulus E=26.4 kPa; 2) Fibre<br />

constants (defined as in [14] and [15]), C3=56 kPa, C4=11, C5=8.45 MPa,<br />

C6=-8.39 MPa ; and 3) The stretch value where the fibers are straightened * =1.32.<br />

Edge A Edge B<br />

4. RESULTS<br />

Fig. 4 Meshless model applied when<br />

modelling the experiments on elongation of<br />

samples of pia-arachnoid-complex tissue<br />

conducted by Jin et al. [13]. Dimensions:<br />

height (along Y axis) of 12.7 mm, length<br />

(along X axis) of 20.0 mm, thickness of<br />

0.0246 mm. Edge A was rigidly constrained,<br />

and Edge B was subjected to elongation in X<br />

(horizontal) axis direction.<br />

4.1 Verification of meshless algorithm for soft tissue fragmentation modelling<br />

The results obtained using our meshless algorithm with visibility criterion for<br />

fragmentation modelling were very close to those calculated by the ABAQUS finite<br />

element solver. The maximum difference in the nodal displacements along X axis was<br />

0.37 mm (i.e. less than 2% of the total model elongation) and occurred in the vicinity of<br />

the notch tip where the finite element mesh exhibited large distortion (Fig. 5). The<br />

average absolute difference in the nodal displacements was order of magnitude smaller<br />

than the maximum difference and equaled 0.029 mm (which is only 0.15% of the total<br />

model elongation).<br />

Fig. 5 Results of verification of the<br />

meshless algorithm in modelling of<br />

discontinuities (crack) due to tissue<br />

fragmentation. The model set-up used in<br />

this verification is shown in Fig. 3.<br />

Distribution of the absolute differences (in<br />

metres) in nodal displacements in X<br />

direction between the meshless algorithm<br />

and ABAQUS finite element solver overimposed<br />

on deformed meshless grid for<br />

elongation of 20% of the initial length of<br />

the model. Model dimensions are in metres.<br />

4.2 Evaluation of meshless algorithm performance: Modelling of pia-arachnoid complex<br />

rupture<br />

The nominal stress–strain relationship of the PAC sample predicted by the model<br />

implemented using our meshless algorithm with visibility criterion was within the<br />

envelope of the experimental responses (Fig. 6). Overall deformation of the sample<br />

predicted by the model agreed well with that observed in the experiments (Fig. 7).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!