Research Article  Open Access
Modelling of the Deformation Diffusion Areas on a ParaAramid Fabric with BSpline Curves
Abstract
The geometrical modelling of the planar energy diffusion behaviors of the deformations on a paraaramid 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 Bspline curvegenerating 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 Bspline 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 paraaramid 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.
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 Bspline methods. Bezier curves and surfaces were independently developed by P. De Casteljau at Citröen and by P. Bezier and Renault [1]. Gordon and Riesenfeld proposed curves and surfaces which use basis splines as blending functions. These are called Bspline curves and surfaces. Bspline curves, like Bezier curves, are defined by polygon vertices and have properties similar to those of Bezier curves. That is, Bspline 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 [2]. 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 from one end and a distance from the other is a closed curve whose area is less than that of the original curve by . 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 [3].
In this research the Classic Holditch Theorem is given in [4]. Also in the paper [5] the authors used Holditch Theorem to calculate the areas between two quadratic uniform Bspline curves. In this study, the mathematical modelling of the deformation areas on paraaramid fabric has been investigated using the Holditch Theorem which calculated the areas between two uniform Bspline curves.
On the other hand, paraaramid fabrics produced from hightenacity 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 [6–15]. In our paper, we have focused on a type of paraaramid fabric and we have investigated geometrical modelling of the planar energy diffusion behaviors of the deformations on the ballistic fabric.
2. Material
Paraaramid fabrics which are a class of heatresistant 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 paraaramid fabrics is used. Twaron fabric is one of ballistic fabrics used in steel vest production and developed in the early 1970s; see Figure 1.
With the help of textile engineers, the warp and wept information of the paraaramid 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 paraaramid sewing yarn. Fabric dimensions are cm and parameters of these fabrics are given in Table 1.

In our experiment step which was performed in the textile laboratory, the drop test was used under the 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 buttend tip shown with B, and a sharp tip shown with C in Figure 2.
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 2dimensional were taken in Figure 3.
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 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 on a circle (C) with a radius in Euclidean plane be divided by a point D into two segments with lengths of and , 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 as following Figure 4 in [4].
3.2. Using BSplines 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 Bspline curves. Like the Bezier curves, the Bspline curves are also defined by polygon corners and have similar properties. However, the Bspline 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 Bspline 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 Bspline curve, only some parts are affected, and the other parts remain unchanged. For this reason, because of their ease in practice, the Bspline 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: and Set: and Then is the point with parameter value on the Bezier curve .
The polygon formed by is called the Bezier polygon or control polygon of the curve . Similarly the polygon vertices are called control points or Bezier points. Sometimes we also write or, shorter, . This notation defines to be the linear operator that associates the Bezier curve with its control polygon. We say that the curve is the BernsteinBezier 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 Bspline curves, because it ensures a certain method in taking the control points. Let . Defineand Then is the value of the Bspline curve at parameter value . Here, denotes the multiplicity of if it was already one of the knots. If it was not, set As usual, we set , [1, 2].
3.5. Uniform Closed BSpline
A class closed curve (of degree ) consisting of curve segments , , with a curve defining polygon determined by the position vectors in Figure 5 can be expressed by the following function: where the initial and finish points are equal to each other, in [1, 2]. We planned to use the uniform closed Bspline curves in our study because the photos of the deformation traces were taken in 2D during the experiment.
3.6. Location of Joints Uniform Quadratic BSpline Segments
In [5], specifically each part is expressed as parabolas consisting of control points with three points because local control is ensured in Bspline curves in Figure 6. The matrix formats of these quadratic Bspline curves are as shown in the following equation: 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.
When we expand this method by taking the point, the matrices in the following formulas are formed for each part:Uniform Bsplines 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 Bspline curve produced by 6 control points. Also we can find location of joints between uniform quadratic Bspline segments. If we take the vertices of this polygon as control points, quadratic uniform Bspline curve will pass from midpoints of sides. Equations of the curve segments will give us a closed curve for control points.
3.7. The Holditch Theorem by Using BSplines
The and chords taken in the Classic Holditch Theorem are taken especially in equal lengths in [3, 4]. We can find the area between the Bspline closed curves asby the help of Holditch Classical Theorem in Figure 10.
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 Bspline curves. By doing so and by observing the deformation area in the cloth, we aim to determine which weightend was used; see Figures 11 and 12.
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:
−12,7  12,15  15,21  5,04  −12,09  −15,62  
12,11  23,35  −4,55  −15,19  −11,56  2,79  
−0,275  13,68  10,125  −3,525  −13,855  −14,16  
17,73  9,4  −9,87  −13,375  −4,385  7,45 
Spline Interpolation of Big Polygon for Ball Tip
Small Polygon for the Ball Tip
−3,62  −0,09  3,54  3,54  0,01  −3,66  
2,22  4,25  2,08  −2,29  −4,32  −2,15  
−1,855  1,725  3,54  1,775  −1,825  −3,64  
3,235  3,165  −0,105  −3,305  −3,235  0,035 
Spline Interpolation of Small Polygon for Ball Tip
4.2. Modelling of Deformation Splines for Sharp Tip
Big Polygon for the Sharp Tip in Figure 14
−9,79  −1,09  10,19  10,55  1,97  −11,71  
3,91  12,47  4,13  −3,62  −11,12  −3,86  
−5,44  4,55  10,37  6,26  −4,87  −10,75  
8,19  8,3  0,255  −7,37  −7,49  0,025 
Spline Interpolation of Big Polygon for Sharp Tip
Small Polygon for the Sharp Tip
−3,59  0  3,66  3,66  0,06  −3,59  
2,19  4,37  2,17  −2,17  −4,36  −2,18  
−1,795  1,83  3,66  1,86  −1,765  −3,59  
3,28  3,27  0  −3,265  −3,27  0,005 
Spline Interpolation of Small Polygon for Sharp TipThese splines can be modelling in the MATLAB program; see Figure 15. For example, the spline interpolation of ball tip can be represented in the MATLAB program as shown in Figure 15.
Also, the area between the Bspline closed curves can be found by using (6); thus Table 2 shows the areas for all tips.

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.
References
 G. Farin, Curves and Surfaces for Computer Aided Geometric Design, A Practical Guide, Academic Press, 2nd edition, 1990. View at: MathSciNet
 F. Yamaguchi, Curves and Surfaces in Computer Aided Geometric Design, Springer, Berlin, Germany, 1988. View at: Publisher Site  MathSciNet
 http://en.wikipedia.org/wiki/Holditch%27s_theorem.
 K. Yucekaya Gülay and H. H. Hacısalihoğlu, Holditch's Theorem for Circles in 2Dimensional Euclidean Space, vol. 18, Dumlupınar Universitesi Fen Bilimleri Dergisi, April 2009.
 A. Çalışkan and B. Tantay, “Application of Holditch theorem to bezier and quadratic uniform BSpline curves,” Mathematical and Computational Applications, vol. 6, no. 3, pp. 169–176, 2001. View at: Publisher Site  Google Scholar
 M. Karahan, A. Kuş, and R. Eren, “An investigation into ballistic performance and energy absorption capabilities of woven aramid fabrics,” International Journal of Impact Engineering, vol. 35, no. 6, pp. 499–510, 2008. View at: Publisher Site  Google Scholar
 M. Karahan, F. Gülsoy, and S. Gündoğan, “The determination of energy propagating behaviour of woven paraaramid fabrics by 2D thin plate spline method,” in Proceedings of the SAMPE Symposium Proceeding CD, Baltimore, Md, USA, June 2007. View at: Google Scholar
 F. Gulsoy, H. Kuşak Samancı, and A. Çalışkan, “Investigation of the bending energy and the expansion areas of the energy on Kflex fabrics,” in Proceedings of the 14th World Textile Conference (AUTEX '14), 2014. View at: Google Scholar
 M. Karahan, “The effect of fibre volume fraction on damage initiation and propagation of woven carbonepoxy multilayer composites,” Textile Research Journal, vol. 82, no. 1, pp. 45–61, 2012. View at: Publisher Site  Google Scholar
 M. Karahan, “Comparison of ballistic performance and energy absorption capabilities of woven and unidirectional aramid fabrics,” Textile Research Journal, vol. 78, no. 8, pp. 718–730, 2008. View at: Publisher Site  Google Scholar
 M. Fahool and A. R. Sabet, “Parametric study of energy absorption mechanism in Twaron fabric impregnated with a shear thickening fluid,” International Journal of Impact Engineering, vol. 90, pp. 61–71, 2016. View at: Publisher Site  Google Scholar
 D. Dimeski, V. Srebrenkoska, and N. Mirceska, “Ballistic impact resistance mechanism of woven fabrics and their composites,” Journal of Engineering Research & Technology, vol. 4, no. 12, pp. 107–111, 2015. View at: Google Scholar
 V. P. W. Shim, V. B. C. Tan, and T. E. Tay, “Modelling deformation and damage characteristics of woven fabric under small projectile impact,” International Journal of Impact Engineering, vol. 16, no. 4, pp. 585–605, 1995. View at: Publisher Site  Google Scholar
 A. Tabiei and G. Nilakantan, “Ballistic impact of dry woven fabric composites: a review,” Applied Mechanics Reviews, vol. 61, no. 1, Article ID 010801, 13 pages, 2008. View at: Publisher Site  Google Scholar
 P. M. Cunniff, “An analysis of the system effects in woven fabrics under ballistic impact,” Textile Research Journal, vol. 62, no. 9, pp. 495–509, 1992. View at: Publisher Site  Google Scholar
Copyright
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.