IJEM International Journal of Engineering Mathematics 2314-6109 2356-7007 Hindawi 10.1155/2017/7159738 7159738 Research Article Modelling of the Deformation Diffusion Areas on a Para-Aramid Fabric with B-Spline Curves http://orcid.org/0000-0001-6685-236X Kuşak Samancı Hatice 1 Yağcı Filiz 2 Çalişkan Ali 3 Carbone Giuseppe 1 Department of Mathematics Bitlis Eren University Bitlis Turkey bitliseren.edu.tr 2 Department of Mathematics Uludağ University Bursa Turkey uludag.edu.tr 3 Department of Mathematics Ege University İzmir Turkey ege.edu.tr 2017 20 02 2017 2017 06 06 2016 29 08 2016 15 09 2016 20 02 2017 2017 Copyright © 2017 Hatice Kuşak Samancı et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The geometrical modelling of the planar energy diffusion behaviors of the deformations on a para-aramid fabric has been performed. In the application process of the study, in the experimental period, drop test with bullets of different weights has been applied. The B-spline curve-generating technique has been used in the study. This is an efficient method for geometrical modelling of the deformation diffusion areas formed after the drop test. Proper control points have been chosen to be able to draw the borders of the diffusion areas on the fabric which is deformed, and then the De Casteljau and De Boor algorithms have been used. The Holditch area calculation according to the beams taken at certain fixed lengths has been performed for the B-spline border curve obtained as a closed form. After the calculations, it has been determined that the diffusion area where the bullet with pointed end was dropped on a para-aramid fabric is bigger and the diffusion area where the bullet with rounded end was dropped is smaller when compared with the areas where other bullets with different ends were dropped.

Uludag Üniversitesi UAP(F)-2012/21
1. Introduction

The exploration of the use of parametric curves and surfaces can be viewed as the origin of Computer Aided Geometric Design (CAGD). The major breakthroughs in CAGD were undoubtedly the theory of Bezier surfaces and Coons patches, later combined with B-spline methods. Bezier curves and surfaces were independently developed by P. De Casteljau at Citröen and by P. Bezier and Renault . Gordon and Riesenfeld proposed curves and surfaces which use basis splines as blending functions. These are called B-spline curves and surfaces. B-spline curves, like Bezier curves, are defined by polygon vertices and have properties similar to those of Bezier curves. That is, B-spline curves are expressed as a convex combination of polygon vertex position vectors and also have the variation diminishing property. Curve shapes are smoothed versions of the polygon shapes and can be roughly predicted from the polygon shapes . In plane geometry Holditch Theorem states if a chord of fixed length is allowed to rotate inside a convex closed curve then the locus of a point on the chord a distance p from one end and a distance q from the other is a closed curve whose area is less than that of the original curve by π p q . The theorem was published in 1858 by Rev. Amnet Holditch. While not mentioned by Holditch, the proof of the theorem requires an assumption that the chord be short enough that the traced locus is a simple closed curve .

In this research the Classic Holditch Theorem is given in . Also in the paper  the authors used Holditch Theorem to calculate the areas between two quadratic uniform B-spline curves. In this study, the mathematical modelling of the deformation areas on para-aramid fabric has been investigated using the Holditch Theorem which calculated the areas between two uniform B-spline curves.

On the other hand, para-aramid fabrics produced from high-tenacity fibres are used in ballistic protection due to their high energy absorption ability and low tenacity/weight ratio. Ballistic behavior of textile fibres and fabrics has been investigated experimentally and ballistic behavior of textile fabric systems has been estimated in . In our paper, we have focused on a type of para-aramid fabric and we have investigated geometrical modelling of the planar energy diffusion behaviors of the deformations on the ballistic fabric.

2. Material

Para-aramid fabrics which are a class of heat-resistant and strong synthetic fabrics are used in aerospace and military applications for ballistic protection. In our research, Twaron CT type 710 fabric which is a type of para-aramid fabrics is used. Twaron fabric is one of ballistic fabrics used in steel vest production and developed in the early 1970s; see Figure 1.

Structure of para-aramid fabric.

With the help of textile engineers, the warp and wept information of the para-aramid fabric, which is in CT710 type, and other parameters of these fabrics as density and weight are shown in Table 1. Fabrics were plied 1, 2, 4, 6, and 8 ply numbers and sewn 2,5 cm from edges with plain stich with high twist para-aramid sewing yarn. Fabric dimensions are 31 × 31  cm and parameters of these fabrics are given in Table 1.

Properties of Twaron CT710 type fabrics used in this research, .

Fabric Type Count Warp/Weft Meterial Type Warp/Weft Weave Density Fabric Weight Treatment
CT 710 930/930 2040/2000 Plain 117/117 220 Scouring/Water repellent treatment

In our experiment step which was performed in the textile laboratory, the drop test was used under the N i j standards. The first material consists of three different bullets 45 mm in diameter, 1 kg in weight, consisting of a round tip shown with A, a middle butt-end tip shown with B, and a sharp tip shown with C in Figure 2.

Bullets tip types.

The second material is the clay which has the characteristics of the human skin. Then the third material is the panel arm which is 1 meter in size and with which the bullet shots were made. The fourth material is the 50 mm pipe used for guiding the weight. Without using cloth for three bullet types and using 1 fold, 2 folds, 6 folds, and 8 folds, the number of the shots was 72, and 72 photos which were 2-dimensional were taken in Figure 3.

Test apparatus used in drop tests [7, 8].

The dropping tests are applied onto the fabric layers which are put on a clay base. The clay used is a special one which is in agreement with the characteristics of the human tissue. After the tests, the deformation on the clay is accepted as the possible deformation on the human tissue. The computation proceeded in four steps. Firstly, landmarks are identified from photographs and then landmark coordinates are determined. The coordinates of homologous landmarks are obtained from two shapes. These two shapes may represent either individual specimens or the means of two sets of shapes corresponding target landmarks. A transformation from a source shape to a target shape involves the displacement of the source landmarks to the corresponding target landmarks.

Note. Type Roma Plastilina (clay) which is given in N i j standards was used as clay material in this research. Clay was filled in a box as pointed out in the standard.

3. Methods 3.1. Classical Holditch Theorem

Let an AB chord with a constant length of a + b on a circle (C) with a radius r in Euclidean plane E 2 be divided by a point D into two segments with lengths of a and b , respectively. When the end points A and B of the chord draw the circle in full, then geometric location of D forms an inner circle. The difference between the areas of two closed curves can be calculated from the formula by π · a · b as following Figure 4 in .

Classical Holditch Theorem .

3.2. Using B-Splines for Geometric Modelling

In the early 1960s, J. C. Ferguson from the Boeing Airplane Company in the USA developed a method which defined the curves as vectors. A Ferguson Curve is a cubic vector function obtained with the parameters after determining the starting and ending points of a curve. Later in 1964, S. A. Coons introduced a surface definition method which gave the mathematical definition that provided the border conditions. In this method, the border curves and the position vectors in the 4 corners of the surface patch were taken into consideration. In 1974, P. Bezier from the Renault Company in France defined a curve representation by giving a polygon and called it as the Bezier curve. The Bezier curves were applied in practice to Renault automobile body designs. Gordon and Riesenfeld introduced the curves that used the base splines as combined functions and these curves were called B-spline curves. Like the Bezier curves, the B-spline curves are also defined by polygon corners and have similar properties. However, the B-spline curves differ from the Bezier curves in these points: The control points of the Bezier curves ensure only global control; however, the control points of the B-spline curves ensure local control. In other words, if a point is moved, the whole of the curve changes in the Bezier curve, but in the B-spline curve, only some parts are affected, and the other parts remain unchanged. For this reason, because of their ease in practice, the B-spline curves will be used in our study because they ensure local changes and the other parts of the curve do not change.

3.3. De Casteljau Algorithm

Given: b 0 , b 1 , , b n E 3 and t I R

Set: b i r t = 1 - t b i r - 1 t + t b i + 1 r - 1 t , where: r = 1 , , n , i = 0 , , n - r

and b i 0 ( t ) = b i . Then b 0 n ( t ) is the point with parameter value t on the Bezier curve b n .

The polygon P formed by b 0 , b 1 , , b n is called the Bezier polygon or control polygon of the curve b n . Similarly the polygon vertices b i are called control points or Bezier points. Sometimes we also write b n ( t ) = B [ b 0 , , b n ; t ] = B [ P , t ] or, shorter, b n = B [ b 0 , , b n ] = B P . This notation defines B to be the linear operator that associates the Bezier curve with its control polygon. We say that the curve B [ b 0 , , b n ] is the Bernstein-Bezier approximation to the control polygon, a terminology borrowed from approximation theory. In the Bezier curves, the De Casteljau algorithm is used, because it ensures that the control points are taken in a certain order.

3.4. De Boor Algorithm

The De Boor algorithm is used in the design of the B-spline curves, because it ensures a certain method in taking the control points. Let u [ u 1 , u I + 1 ) [ u n - 1 , u L + n - 1 ] . Define (1) d i k u = u i + n - k - u u i + n - k - u i - 1 d i - 1 k - 1 u + u - u i - 1 u i + n - k - u i - 1 d i k - 1 u f o r k = 1 , , n - r and i = I - n + k + 1 , , I - r + 1 . Then s ( u ) = [ B n P ] ( u ) = d I - r + 1 n - r ( u ) is the value of the B-spline curve at parameter value u . Here, r denotes the multiplicity of u if it was already one of the knots. If it was not, set r = 0 . As usual, we set d i 0 ( u ) = d i , [1, 2].

3.5. Uniform Closed B-Spline

A class C M - 2 closed curve (of degree M - 1 ) consisting of n + 1 curve segments P i ( t ) , i = 0,1 , , n , with a curve defining polygon determined by the n + 1 position vectors Q 0 , Q 1 , , Q n in Figure 5 can be expressed by the following function: (2) P i t = j = i - 1 i + M + 2 N j - i + 1 , M t Q j m o d n + 1 , 0 t 1 , where the initial and finish points (3) P i 0 = P i , i = 0,1 , , n , are equal to each other, in [1, 2]. We planned to use the uniform closed B-spline curves in our study because the photos of the deformation traces were taken in 2D during the experiment.

Uniform closed B-spline.

3.6. Location of Joints Uniform Quadratic B-Spline Segments

In , specifically each part is expressed as parabolas consisting of control points with three points because local control is ensured in B-spline curves in Figure 6. The matrix formats of these quadratic B-spline curves are as shown in the following equation: (4) P i t = 1 2 · t 2 t 1 · 1 - 2 1 - 2 2 0 1 1 0 · Q i - 1 Q i Q i + 1 . In this case, only 3 control vertices are used in Figure 7. By joining these curves which correspond to these matrices, the closed curve is formed. In this situation, the midpoints of the two control points from the starting point are considered as being like tangents to the curves. Since the joining places of these curves are the tangents, the curves are connected to each other continuously with the C1 Continuity in Figure 8.

Location of joints uniform quadratic B-spline segments.

The curves with the C1 Continuity.

When we expand this method by taking the n point, the matrices in the following formulas are formed for each part: (5) P 0 t = 1 2 · t 2 t 1 · 1 - 2 1 - 2 2 0 1 1 0 · Q 0 Q 1 Q 2 , P 1 t = 1 2 · t 2 t 1 · 1 - 2 1 - 2 2 0 1 1 0 · Q 1 Q 2 Q 3 , P n t = 1 2 · t 2 t 1 · 1 - 2 1 - 2 2 0 1 1 0 · Q n Q 0 Q 1 , β : P 1 P 2 P n = β t . Uniform B-splines are convenient to represent closed curves. The only thing needed is a change in the number of the curve segments. Figure 9 shows a closed quadratic B-spline curve produced by 6 control points. Also we can find location of joints between uniform quadratic B-spline segments. If we take the vertices of this polygon as control points, quadratic uniform B-spline curve will pass from midpoints of sides. Equations of the curve segments will give us a closed curve for n control points.

Different way of modelling uniform closed B-spline .

3.7. The Holditch Theorem by Using B-Splines

The a and b chords taken in the Classic Holditch Theorem are taken especially in equal lengths in [3, 4]. We can find the area between the B-spline closed curves as (6) A r e a = π · a 2 by the help of Holditch Classical Theorem in Figure 10.

Holditch Theorem with B-spline curves .

4. Main Results

We applied these methods to the deformation photos at the end of the experiment. Our purpose is to calculate the area between the projection of the bullet and its effective area border curve with the help of the Holditch Theorem after drawing the B-spline curves. By doing so and by observing the deformation area in the cloth, we aim to determine which weight-end was used; see Figures 11 and 12.

Deformation area with ball type.

Deformation areas with sharp type.

4.1. Modelling of Deformation Splines for Ball Tip

The coordinates of big polygon for the ball tip in Figure 13 are given in the following:

Modelling deformation splines with ball type.

 Q 0 L Q 1 L Q 2 L Q 3 L Q 4 L Q 5 L X −12,7 12,15 15,21 5,04 −12,09 −15,62 Y 12,11 23,35 −4,55 −15,19 −11,56 2,79 P 0 L P 1 L P 2 L P 3 L P 4 L P 5 L X −0,275 13,68 10,125 −3,525 −13,855 −14,16 Y 17,73 9,4 −9,87 −13,375 −4,385 7,45

Spline Interpolation of Big Polygon for Ball Tip (7) P 0 L = - 10.835 t 2 + 24.85 t - 0.275 , - 19.57 t 2 + 11.24 t + 17.73 , P 1 L = - 6.615 t 2 + 3.06 t + 13.68,8.63 t 2 - 27.9 t + 9.4 , P 2 L = - 3.48 t 2 - 10.17 t + 10.126,7.135 t 2 - 10.64 t - 9.87 , P 3 L = 6.8 t 2 - 17.13 t - 3.525,5.36 t 2 + 3.63 t - 13.375 , P 4 L = 3.225 t 2 - 3.53 t - 13.855 , - 2.515 t 2 + 14.35 t - 4.385 .

Small Polygon for the Ball Tip

 Q 0 S Q 1 S Q 2 S Q 3 S Q 4 S Q 5 S X −3,62 −0,09 3,54 3,54 0,01 −3,66 Y 2,22 4,25 2,08 −2,29 −4,32 −2,15 P 0 S P 1 S P 2 S P 3 S P 4 S P 5 S X −1,855 1,725 3,54 1,775 −1,825 −3,64 Y 3,235 3,165 −0,105 −3,305 −3,235 0,035

Spline Interpolation of Small Polygon for Ball Tip (8) P 0 S = 0.05 t 2 + 3.53 t - 1.855 , - 2.1 t 2 + 2.03 t + 3.235 , P 1 S = - 1.81 t 2 + 3.63 t + 1.725 , - 1.1 t 2 - 2.17 t + 3.165 , P 2 S = - 1.765 t 2 + 3.54,1.17 t 2 - 4.37 t - 0.105 , P 3 S = - 0.07 t 2 - 3.53 t + 1.775,2.1 t 2 - 2.03 t - 3.305 , P 4 S = 1.855 t 2 - 3.67 t - 1.825,1.1 t 2 + 2.17 t - 3.235 .

4.2. Modelling of Deformation Splines for Sharp Tip

Big Polygon for the Sharp Tip in Figure 14

Modelling deformation splines with sharp type.

 Q 0 L Q 1 L Q 2 L Q 3 L Q 4 L Q 5 L X −9,79 −1,09 10,19 10,55 1,97 −11,71 Y 3,91 12,47 4,13 −3,62 −11,12 −3,86 P 0 L P 1 L P 2 L P 3 L P 4 L P 5 L X −5,44 4,55 10,37 6,26 −4,87 −10,75 Y 8,19 8,3 0,255 −7,37 −7,49 0,025

Spline Interpolation of Big Polygon for Sharp Tip (9) P 0 L = 1.29 t 2 + 8.7 t - 5.44 , - 8.45 t 2 + 8.56 t + 8.19 , P 1 L = - 5.46 t 2 + 11.28 t + 4.55,0.295 t 2 - 8.34 t + 8.3 , P 2 L = - 4.47 t 2 + 0.36 t + 10.37,0.125 t 2 - 7.75 t + 0.255 , P 3 L = - 2.55 t 2 - 8.58 t + 6.26,7.38 t 2 - 7.5 t - 7.37 , P 4 L = 7.80 t 2 - 13.68 t - 4.87,0.255 t 2 + 7.26 t - 7.49 .

Small Polygon for the Sharp Tip

 Q 0 S Q 1 S Q 2 S Q 3 S Q 4 S Q 5 S X −3,59 0 3,66 3,66 0,06 −3,59 Y 2,19 4,37 2,17 −2,17 −4,36 −2,18 P 0 S P 1 S P 2 S P 3 S P 4 S P 5 S X −1,795 1,83 3,66 1,86 −1,765 −3,59 Y 3,28 3,27 0 −3,265 −3,27 0,005

Spline Interpolation of Small Polygon for Sharp Tip (10) P 0 S = 0.035 t 2 + 3.59 t - 1.795 , - 2.19 t 2 + 2.18 t + 3.28 , P 1 S = - 1.83 t 2 + 3.66 t + 1.83 , - 1.07 t 2 - 2.2 t + 3.27 , P 2 S = - 1.8 t 2 + 3.66,1.075 t 2 - 4.34 t , P 3 S = - 0.025 t 2 - 3.6 t + 1.86,2.185 t 2 - 2.19 t - 3.265 , P 4 S = 1.825 t 2 - 3.65 t - 1.765,1.095 t 2 + 2.18 t - 3.27 . These splines can be modelling in the MATLAB program; see Figure 15. For example, the spline interpolation of ball tip (11) P 0 L = - 10.835 t 2 + 24.85 t - 0.275 , - 19.57 t 2 + 11.24 t + 17.73 can be represented in the MATLAB program as shown in Figure 15.

An example of geometric modelling by MATLAB.

Also, the area between the B-spline closed curves can be found by using (6); thus Table 2 shows the areas for all tips.

Calculating the expansion areas of tips.

8 Ply of fabric
Type-A, ball type 971.6465
Type-B, blunt type 598.7112
Type-C, sharp type 529.1348
5. Conclusion

In our paper, the diffusion areas were modelled by a geometrical spline method. Some special points around the expansion area were taken on the deformation photos. Then the geometric spline modelling was done. As was mentioned above, the data of the Holditch areas that were found by using spline method after the geometric modelling study for each tip are given. The geometric spline method is very useful for the estimation of the expansion areas of tips. When we consider the tables, in estimating the bullet tip by observing the spread area, we can say that the spread area that is left by the bullet with round tip is more than the area which is left by the bullet with sharp tip. In our research we can see that using the Holditch area theorem with geometric spline method is very useful method in determining the tips from deformation photographs.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This study has received support from the Uludağ University Scientific Research Unit with the Project Reference no. UAP(F)-2012/21. The authors would like to thank Professor Ali Çalışkan for his valuable comments and suggestions in improving the quality of the paper.

Farin G. Curves and Surfaces for Computer Aided Geometric Design, A Practical Guide 1990 2nd Academic Press MR1058011 Yamaguchi F. Curves and Surfaces in Computer Aided Geometric Design 1988 Berlin, Germany Springer 10.1007/978-3-642-48952-5 MR971750 http://en.wikipedia.org/wiki/Holditch%27s_theorem Yucekaya Gülay K. Hacısalihoğlu H. H. Holditch's Theorem for Circles in 2-Dimensional Euclidean Space April 2009 18 Dumlupınar Universitesi Fen Bilimleri Dergisi Çalışkan A. Tantay B. Application of Holditch theorem to bezier and quadratic uniform B-Spline curves Mathematical and Computational Applications 2001 6 3 169 176 10.3390/mca6030169 2-s2.0-0034899708 Karahan M. Kuş A. Eren R. An investigation into ballistic performance and energy absorption capabilities of woven aramid fabrics International Journal of Impact Engineering 2008 35 6 499 510 10.1016/j.ijimpeng.2007.04.003 2-s2.0-38949139053 Karahan M. Gülsoy F. Gündoğan S. The determination of energy propagating behaviour of woven para-aramid fabrics by 2-D thin plate spline method Proceedings of the SAMPE Symposium Proceeding CD June 2007 Baltimore, Md, USA Gulsoy F. Kuşak Samancı H. Çalışkan A. Investigation of the bending energy and the expansion areas of the energy on K-flex fabrics Proceedings of the 14th World Textile Conference (AUTEX '14) 2014 Karahan M. The effect of fibre volume fraction on damage initiation and propagation of woven carbon-epoxy multi-layer composites Textile Research Journal 2012 82 1 45 61 10.1177/0040517511416282 2-s2.0-83455214106 Karahan M. Comparison of ballistic performance and energy absorption capabilities of woven and unidirectional aramid fabrics Textile Research Journal 2008 78 8 718 730 10.1177/0040517508090487 2-s2.0-51649085409 Fahool M. Sabet A. R. Parametric study of energy absorption mechanism in Twaron fabric impregnated with a shear thickening fluid International Journal of Impact Engineering 2016 90 61 71 10.1016/j.ijimpeng.2015.11.016 2-s2.0-84951309685 Dimeski D. Srebrenkoska V. Mirceska N. Ballistic impact resistance mechanism of woven fabrics and their composites Journal of Engineering Research & Technology 2015 4 12 107 111 Shim V. P. W. Tan V. B. C. Tay T. E. Modelling deformation and damage characteristics of woven fabric under small projectile impact International Journal of Impact Engineering 1995 16 4 585 605 10.1016/0734-743X(94)00063-3 2-s2.0-0029356559 Tabiei A. Nilakantan G. Ballistic impact of dry woven fabric composites: a review Applied Mechanics Reviews 2008 61 1 13 010801 10.1115/1.2821711 2-s2.0-59449101487 Cunniff P. M. An analysis of the system effects in woven fabrics under ballistic impact Textile Research Journal 1992 62 9 495 509 10.1177/004051759206200902 2-s2.0-0026918439