SciELO - Scientific Electronic Library Online

 
vol.85 número207Study the spatial variability of the noise levels inside two commercial poultry housing with different adiabatic evaporative cooling systemsStudy of the structural and magnetic properties of the system Bi1-xYxFeO3 x = 0 and 0.07 using Mössbauer spectroscopy índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google

Compartilhar


DYNA

versão impressa ISSN 0012-7353

Dyna rev.fac.nac.minas vol.85 no.207 Medellín out./dez. 2018

https://doi.org/10.15446/dyna.v85n207.67537 

Artículos

System of cutting force data acquisition in mechanical lathes

Sistema de adquisición de datos de la fuerza de corte en tornos mecánicos

Felipe Sineider Martinsa 

Rodrigo Reina-Muñoza 

Valdemir Martins Liraa 

a Universidade Federal do ABC, Santo André, Brazil. felipe.sineider@ufabc.edu.br, rodrigo.munoz@ufabc.edu.br, valdemir.lira@ufabc.edu.br


Abstract

This report concerns a data acquisition system with the final idea of measuring the cutting force in the threading process done on a conventional lathe. Analytical models are used to find the depth of cut of the machining tool. This is important because the machining force can be responsible for the tool collapse by plastic edge deformation. The cutting force was obtained through direct measurements in selected samples. This allows obtaining more realistic data compared with data obtained by analytical estimation from theory. The data acquisition system was implemented in LabView®, and it allowed analyzing the several variables involved in the machining process on conventional lathes. In such a way, it is possible to find out, in real time, the depth of cut variable. Specifically, the authors present results obtained from analysis of Fredericks, Sandvik and Lira models. This analysis allows drawing conclusions about the cutting force from each model.

Keywords: machining force; machining; data acquisition; depth of cut

Resumen

Este artículo presenta un sistema de adquisición de datos para medir la fuerza de corte en el proceso de roscado hecho en un torno convencional. São usados modelos analíticos para encontrar la profundidad del corte de la herramienta de fabricación. La fuerza de la herramienta de corte puede ser responsable por el colapso de la herramienta debido a la deformación plástica de la arista. La fuerza de corte fue obtenida a través de medidas directas en muestras seleccionadas permitiendo obtener datos confiables comparados con datos obtenidos del análisis teórico. El sistema fue implementado en LabView® analizando diversas variables que entran en juego en el proceso de manufactura en tornos convencionales. Es posible así encontrar en tiempo real la profundidad del corte. Específicamente, los autores presentan resultados usando los modelos de Fredericks, Sandvik y Lira. El análisis permite obtener conclusiones sobre la fuerza de corte de cada modelo.

Palabras clave: fuerza mecánica; mecanizado; adquisición de datos; profundidad del corte

1. Introduction

In the machining process, several variables are influenced by the process conditions [1-2], but only some variables are useful for the monitoring process and can be measured with sensors installed close to the machining tool zone [3]. The variable that can be considered having more relevance is the cutting force, because it is directly related to several performance characteristics and optimization [4]. For example, a determined cutting condition can be observed and the cutting parameters can be adjusted. In such a way, better cutting force efficiency can be obtained and therefore allows extending the tools life time, improving the quality of the machined parts, and contributing to a reduction of the machining costs [5-7].

The tool condition is an important factor to the precision of the machining process [8-9]. The tool wear is highly influenced by the machining force [10-11], so the higher the actuating forces on the cutting tool, the higher the tool wear [12]. The cutting force measurement can be used as an indication of change of the tool related to the progressive wear and prevention of the tool´s disruption in a conventional turning process [9]. Experimental results reported in [13], shows that the gradual growth of flank wear, the rapid cutting edge deterioration and catastrophic failure of the cutting tool, can be detected by using cutting force monitoring.

It is well known that during the machining process, the machining tool wears out, and in turn can increase the cutting force of the cutting tool. Therfore, there is important to look for a mechanism that will better be able to monitor the relevant variables in the machining process in real time, specifically the cutting force variable. Real time monitoring of the cutting force is the best solution [14].

In a real situation, there does not exist any theory related to workpiece machining that can be considered well suited enough, so that the cutting force can be easily predicted. What we have normally is a set of empiric formulations involving a set of constants which are not rapidly available, and also a dependency on several factors which demand a huge amount of data [15].

The turning process shown in Fig. 1 is a machining operation where a mechanical workpiece is fixed in the plate of a lathe and is rotated while the turning tool, fixed in a tool holder moves in a (x, z) plane that ideally contains the rotation axis of the machined workpiece [16]. This motion produces the machining force which is defined as the applied force by the mechanical workpiece to the cutting tool. This is the variable we are interested in to do the monitoring. The turning process described in Fig. 1 is characterized by three types of movements and five types of forces [16-18]:

Source: [6,17-18].

Figure 1 Turning operation of a conventional lathe. 

  • Cut movement: It is characterized as the motion which creates only one chip removal.

  • Advance motion: It is characterized as the motion that promotes a repetitive or continuous lifting chip during several revolutions.

  • Effective cutting motion: It is characterized as the resulting motion due to cutting and advancing movements.

  • Fu: Machining force.

  • Ft (Active force): component of Fu localized in the working plane.

  • Fc (cutting force): projection of Fu on the cutting direction (given by the cutting speed).

  • Ff (advance force): Projection of Fu on the advance direction.

  • Fp (passive force): Component of Fu localized perpendicularly to the working plane.

Taking into account that few reports can be found about the cutting force in the threading process in a conventional lathe, this paper has the objective of doing the monitoring of the cutting force (Fc) in this kind of process, in real time. To do so, it was necessary to develop a data acquisition system based on LabView® programming.

In the research reported by Fredericks [19] and Teets [20], in the tools datasheets [21], and in many CNC programing manuals about the Thread cutting (TC) concerning the advance analyses methods of the tools into mechanical part direction, there are basically two strategies: One of them foresees constant cutting depth (CCD) (with the consequent cutting cross sections areas being variables) as can be seen in Fig. 2a, and the second strategy considering variable cutting depths (VCD) with cutting cross sections constants. This can be seen in Fig. 2b.

Source: The Authors.

Figure 2 Strategies of the depth of cutting: (a) constant (b) variable. 

Usually, developers follow the tools manufacture’s recommendations that guide the user through tables obtained empirically to the VCD strategy, giving the necessary geometrical input parameters to the CNC machine program.

In the proposed methodology there will be a description of the analytical models that consider the advance strategy with variable depth of cutting.

2. Materials and methodology

2.1. Materials

In the setup of the experiments there were used three proof samples that were configured as shown in Fig. 3. In this figure there are shown the dimensions used in these samples. Samples consist of ABNT steel 1045.

Source: The Authors.

Figure 3 Dimensions of the proof samples used in the experiment. 

Fig. 4 shows the setup used. The support for the insert used was R166 OG-16VM01-002 manufactured by Sandvik company (Fig. 4(a)). Also used in the experiment was a mechanical lathe by Romi, model TORMAX 30 A, a tool holder adapted with the load cell, a weight indicator (type Siciliano SDS), and a computer with the LabView® software installed.

Source: The Authors.

Figure 4 Experiment setup used for data acquisition of the cutting force: (a) Disk support and load cell; (b) hardware to signal processing, analog to digital conversion and exhibition of the measured cutting force; (c) Computer to data acquisition and storage. 

2.2. Methodology

Fig. 5 shows variables and parameters considered in the cutting force monitoring. The depth of the cutting parameter will be the only input variable. The rest of the parameters will be held constant. This parameter was chosen because it promotes a significant influence on the cutting force [22].

Source: The Authors.

Figure 5 Input and output parameters of the experiment. 

This parameter is calculated using different methods. Markedly, the Fredericks´s model [19], the Sandvik´s model [21], and Lira´s model [23]. All of them were used to find the radial advance of the tool to the machining of the cylindrical threads in a conventional lathe.

The input parameters that were held constant were defined according to: The total number of passes (each pass is related with the infeed of a sample into the lathe) that are made by the tools (six passes in this experiment), the total thread depth is 0.94 mm and the fillet tip radius of the insert 0,2 mm. Also considered, was a rotation of 900 rpm. The input parameter that is considered variable is defined by the models:

- Sandvik´s model [21]: In this model is used the depth of cutting and the number of passes, for determination of absolute radial advance according with eq.(1):

Where: Δapx is absolute radial advance (mm), ap is the depth of cutting (mm), nap is the effective number of passes and j is a parameter that has a 0,3 value in the first trial, and for the rest of the passes it is calculated accordingly with the next equation:

- Lira´s Model [23]: In this model, in the calculus of the absolute radial depth is considered the tool tip radius. Calculus is made accordingly with eq.(3):

Where R f is the tool tip radius, i is the current trial, N is the total number of passes, and h is the depth of the thread.

Fredericks´s model [19]: In this model, the radial depth is calculated according to eq.(4):

Where: n is the current pass and N is the number of passes.

Using the data obtained from the eqs. (1), (3) and (4) there was a graph constructed as shown in Fig. 6. It is inferred that the Fredericks algorithm is the one that presents the greatest initial radial advance as well as the greater variation of the incremental radial advance. The model proposed by Sandvik shows the smallest incremental radial advance variation. The mathematical model of Lira indicates an intermediate behavior.

Source: The Authors.

Figure 6 Curves representing the data obtained for the absolute radial advance for each model considered in the respective pass. 

With the data obtained by the eqs. (1), (3) and (4), there were obtained the values of the effective radial advance in each pass of each mathematical model as shown in Table 1.

Table 1 Radial effective progress in each pass (mm). 

Source: The Authors.

2.3. LabView® data acquisition

The experimental setup is illustrated in Fig. 4. In that figure (part A), the load cell detects the analog signal of the cutting force and a weight indicator executes the conversion of the measured signal to the digital representation and exhibits the measured signal (part B of Fig. 4). This signal, in turn, is transmitted to the computer where data acquired from the weight indicator is processed and stored. Processing is done using an interface (program) developed in LabView® (part C in Fig. 4).

Fig. 7 shows a schematic of the setup used in the experiment to collect the data. This schematic is a simplified representation of Fig. 4. Part A in the figure corresponds with the load cell to acquire the analog signal. This signal goes into the weight indicator (part B) where there is a signal conditioning to obtain a digital signal. The digital signal goes via USB (part C) to the computer (part D). Finally, in the computer the data obtained is processed and stored using an application developed in LabView®.

Source: The Authors.

Figure 7 Schematic of the setup used to collect data of the cutting force. 

Fig. 8(a) shows the main user interface in the LabView® Control panel. Using this interface, we developed the rest of the programming as shown in Fig. 8(b). As can be seen in Fig. 8(b), the programming was divided in six parts:

  • - A: In this part of the LabView® program the reading of the serial port is specified by the user and the preparation to recognize the termination character. This ending character also can be modified as necessary.

  • - B: In this block there is a processing of the reading of a specific number of bytes and the returning of the read data as a string.

  • - C: In this part of the programming is done the conversion of the read string into a numeric format. Also, completed in this block the counting and exhibition of the execution time.

  • - D: The number found in the previous item is converted into dynamic data as well as is done the association of this information with time.

  • - E: In this block it is implemented a graph and the data tabulation. As soon as the stop button is pressed the data in the table is sent to an Excel file.

  • - F: Correspond with the program completion or then causes the stopping of the program in case of an error in the serial port can be detected.

Source: The Authors.

Figure 8(a) Main interface of the program in the LabView control panel. 

Source: The Authors.

Figure 8(b) LabView® programming to acquisition and data storage. 

3. Results and discussions

In this experiment there were used three different mathematical models to determine the radial advancement of the tool to each pass. There were completed six passes for each model. The average data of the force that were obtained with the data acquisition system are shown in Fig. 9. This figure shows the average cutting force of each pass for each model considered. It can be seen that for each model, there is a different force profile with a different response as a function of the progression of passes.

Source: The Authors.

Figure 9 Average cutting force for each model and each trial 

In a machining process the cutting force increases proportionally with the incrementation of the cutting depth. But in the specific case of threading, the force variation during the passes is caused by the contact area variation of the cutting tool with the machined part. The greater the contact area, the greater is the resistance to the cutting tool advance, and therefore the greater is the force required to be applied to the machined part.

It can be inferred from Fig. 9:

  • - With Sandvik’s model there was a progressive increment of the cutting tool until the fifth pass. In the sixth pass, there was a minor force in comparison with fifth pass. This occurs due to the fact that in the fifth pass the thread was almost finished requiring a minor force. It was observed that the first and fifth passes are those having the minor and greater strength respectively between the models considered in these experiments. This model shows that the cutting force increases progressively and decreases in the last pass. The main problem observed in this model was that it exhibited a greater cutting force increment between the passes. A greater cutting force variation between the passes can cause a greater likelihood of catastrophic failure of the cutting tool.

  • - Lira’s model. Also was observed a progressive increment of the cutting force but only up to the fourth pass. For the remaining passes the cutting force decreased. This is due to the fact that after the fourth pass the thread was almost finished, and therefore, in the remaining two passes the cutting force required is minor. This model is very interesting compared to Sandvik’s model due to the cutting force increment between the passes was the minor observed with the models chosen to run the experiments.

  • - Fredericks’s model. Results with this model were somewhat different compared with previous models, showing the cutting force decreasing progressively along the passes. The model shows the minor maximum cutting force achieved in the first pass as well as the minor average cutting force of all models studied. The behavior observed with this model indicates that the cutting tool is exposed to a minor average cutting force when compared to the other models. It means that the cutting tool can have minor wear when this model is considered.

The knowledge about the threading process is very important as it allows us to minimize cost of the machining process. It is possible because allows predicting the cutting tools wear or failure probability and therefore avoiding the loss coming from unnecessary maintenance. As it is well known, usually it is adopted by the industry as a criterion of end of life of the cutting tool, the number of machined parts and safety factors guidelines. This approach can lead to additional machining costs. Also it would be possible avoid loss due to bad sizing of machines and tools, as it would otherwise be without a cutting tool force study, helping to estimate with enhanced precision the power required to specify the machines (lathes) and tools selection, as wells as a better sizing of fixing systems.

4. Conclusion

We presented here a system for data acquisition of the cutting force in mechanical lathes. The data acquisition system was implemented in LabView®. The results are as highlighted: The passes that have the higher average cutting force are: In the Sandvik´s model it occurs in the fifth trial (48.41 kgf). In Lira´s model, it occurs in the fourth trial (42.55 kgf). Finally, in Fredericks´s model, the higher cutting force occurs in the first trial (41.22 kgf). As can be expected, the passes that exhibited the highest cutting force values are the more critical passes, and are the ones with higher probability of catastrophic failure of the cutting tool.

Considering the average cutting force obtained with the three models (Sandvik, Lira and Fredericks models) the conclusion that can be drawn is that Sandvik´s model is the one having the highest cutting force variation between each trial (8.56 kgf) and a higher average cutting force of (36.45 kgf), and also having the passes with higher (48.41 kgf) and lower average cutting force (18.39 kgf). Therefore, such observed behavior can promote a bigger stress on the cutting tool which means a higher probability of failure of the cutting tool. In the case of Fredericks´s model, it was observed that it exhibit the lower average cutting force (33.99 kgf), and lower maximum peak of cutting force (41.22 kgf) taking place in the first trial. Lira´s model has an intermediate average cutting force (35.34 kgf) when compared with Sandvik and Fredericks models. This model exhibits the lower cutting force variations between passes (5.20 kgf).

This implemented data acquisition system allowed us to know, in real time, the intensity of the cutting force in conventional lathes. It allowed us to introduce machining strategies in such a way to preserve the cutting tool during the machining utilization, increasing its lifespan and allowing for obtaining products with the best quality. These experimental results can be useful as a manufacturing tool to predict the probability of failure of the cutting tool, and helping to reduce the cost of maintenance. Also important in this context is the possibility of having a mechanism that can help manufacturers to improve sizing of the machines and tools.

Acknowledgements

Authors would like to thank the Federal University of ABC - UFABC by supporting this applied research.

References

[1] Paul, P.S., Varadarajan, A.S. and Lawrance, G., Effect of material on damping characteristics of impact mass during hard turning. Emerging Materials Research, 4(1), pp. 81-88, 2015, DOI: 10.1680/emr.15.00007. [ Links ]

[2] Das, S.R., Mohapatra, D.K., Routray, P.C. and Rout, B.A., Effects and optimization of machining parameters in hard turning process: A review. Journal of Engineering Innovation and Research [online]. 5(1), pp. 1-9, 2015. [date of reference: August 20th of 2015]. Available at: Available at: https://jeir.org/paperfile/8786544_EFFECTS_AND_OPTIMIZATION_OF_MACHINING_PARAMETERS_IN_HARD_TURNING_PROCESS-A_REVIEW.pdfLinks ]

[3] Teti, R., Advanced IT methods of signal processing and decision making for zero defect manufacturing in machining, Proceedings of the 3rd CIRP Global Web Conference - Production Engineering Research: Advancement beyond state of the art, 28, pp. 3-15, 2015. DOI: 10.1016/j.procir.2015.04.003. [ Links ]

[4] Janković, P., Madić, M. and Lazarević, A., Robust conditions for cutting force minimization in polyamide turning process. FME Transactions, 43(2), pp. 114-118, 2015. DOI: 10.5937/fmet1502114L. [ Links ]

[5] Campocasso, S., Costes, J.P., Fromentin, G., Bissey-Breton, S. and Poulachon, G., A generalized geometrical model of turning operations for cutting force modelling using edge discretization. Applied Mathematical Modelling, 39(21), 2015. DOI: 10.1016/j.apm.2015.02.008. [ Links ]

[6] Chen, L.H., Xu, Y.W., Zhang, L.H. and Man, J., A new cutting force model for predicting stability of interrupted turning. Materials Research Innovations, 19(6), pp. S6-108-S6-110, 2015. DOI: 10.1179/1432891715Z.0000000001460. [ Links ]

[7] Khechana, M., Djamaa, M.C., Djebala, A. and Boumazza, Z., Identification of structural damage in the turning process of a disk based on the analysis of cutting force signals. The International Journal of Advanced Manufacturing Technology, 80(5-8), pp. 1363-1368, 2015. DOI: 10.1007/s00170-015-7110-6. [ Links ]

[8] Chen, X., Cheng, K. and Wang, C., Design of a smart turning tool with application to in-process cutting force measurement in ultraprecision and micro cutting. Manufacturing Letters, 2(4), pp. 112-117, 2014. DOI: 10.1016/j.mfglet.2014.07.001. [ Links ]

[9] Wang, C., Cheng, K., Nelson, N., Worapong, S. and Rakowski, R., Cutting force-based analysis and correlative observations on the tool wear in diamond turning of single-crystal silicon. Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 229(100), pp. 1867-1873, 2014. DOI: 10.1177/0954405414543316. [ Links ]

[10] Gao, D., Zhirong, L., Zekun, L.V. and Yong, L.U., Multi-scale statistical signal processing of cutting force in cutting tool condition monitoring. The International Journal of Advanced Manufacturing Technology, Springer London, 80(9-12), 2015. DOI: 10.1007/s00170-015-7116-0. [ Links ]

[11] Sheng, J., A modeling method for turning parameters coupling based on minimum cutting tool wear. The International Journal of Advanced Manufacturing Technology, 76(1-4), pp. 705-712, 2014. DOI: 10.1007/s00170-014-6318-1. [ Links ]

[12] Banerjee, S., Deka, A., Sarmah, G.D. and Bhardwaj, N.A., Artificial neural network modeling of the effect of cutting conditions on cutting force components during orthogonal turning. International Journal of Current Engineering and Technology, Special Issue 2, pp. 127-130, 2014. DOI: 10.14741/ijcet/spl.2.2014.23. [ Links ]

[13] Chinchanikar, S. and Choudhuryb, S.K., Characteristic of wear, force and their Inter-relationship: In-process monitoring of tool within different phases of the tool life, Proceedings of the International Conference on Advances in Manufacturing and Materials Engineering, ICAMME 2014 [Online]. pp. 1424-1433, 2014. [date of reference: June 20th of 2015]. Available at: https://www.sciencedirect.com/journal/procedia-materials-science/vol/5?page-size=100&page=2Links ]

[14] Khechana, M., Djamaa, M.C., Djebala, A. and Boumazza, Z., Identification of structural damage in the turning process of a disk based on the analysis of cutting force signals. The International Journal of Advanced Manufacturing Technology, 80(5-8), pp. 1363-1368, 2015. DOI: 10.1007/s00170-015-7110-6. [ Links ]

[15] Rao, C.J., Sreeamulu, D. and Mathew, A.T., Analysis of tool life during turning operation by determining optimal process parameters, Proceedings of the 12th Global Congress on Manufacturing and Management, GCMM, 97, pp. 241-250, 2014. DOI: 10.1016/j.proeng.2014.12.247. [ Links ]

[16] Trent, M.C. and Wright, P.K., Metal cutting principles. USA: Butterworth-Heinemann, 4th Edition, 2000. [ Links ]

[17] Zhao, Y., Zhao, Y., Liang, S. and Zhou, G. A high performance sensor for triaxial cutting force measurement in turning. Sensors, 15(4), pp. 7969-7984, Elsevier, 2015. DOI: 10.3390/s150407969. [ Links ]

[18] Abainia, S. and Ouelaa, N., Experimental study of the combined influence of the tool geometry parameters on the cutting forces and tool vibrations. The International Journal of Advanced Manufacturing Technology, 79(5-8), pp. 1127-1138, 2015. DOI: 10.1007/s00170-015-6885-9. [ Links ]

[19] Fredericks, G., Threading algorithm suitable for computer numerically controlled turning center. International Journal of Machine Tools Design & Research, 18(1), pp. 29-31, Oxford, 1978. DOI: 10.1016/0020-7357(78)90017-3. [ Links ]

[20] Teets, B.R., Single-point threading with today’s tooling options. Modern Machine Shop [Online], 62(8), 1990. [date of reference: May 15th of 2014]. Available at: Available at: http://www.bradleyteets.com/work_mms_1/Bradley%20Teets%20Article%20on%20Single%20Point%20Threading.htmLinks ]

[21] Sandvik, Threading: formulas and definitions [Online]. [date of reference: February 10th of 2015]. Available at: https://www.sandvik.coromant.com/en-gb/knowledge/machining-formulas-definitions/pages/threading.aspxLinks ]

[22] Rao, C.J., Rao, D.N. and Srihari, P. Influence of cutting parameters on cutting force and surface finish in turning operation. Procedia Engineering International Conference on Design and Manufacturing, 64, pp. 1405-1415, 2013. DOI: 10.1016/j.proeng.2013.09.222. [ Links ]

[23] Lira, V.M., Cut parameters optimization in the threading machining on CNC Lathes, MSc. Dissertation, Mechanical Engineering Postgraduate Program of the Universidade Estadual de Campinas, Campinas, Brazil, 2000. [ Links ]

How to cite: Martins, F.S, Reina-Muñoz, R. and Lira, V.M, System of cutting force data acquisition in mechanical lathes. DYNA, 85(207), pp. 16-21, Octubre - Diciembre, 2019.

F.S. Martins, received the BSc. in Science and Technology in 2011 and the BSc. Eng in Instrumentation, Automation and Robotics Engineering in 2017, all of them from the Federal University of ABC, Brazil. His research interests include simulation of metal cutting processes; cutting tool condition monitoring; and optimization of manufacturing process. ORCID: 0000-0002-3940-7938

V.M. Lira, is a full assistant professor IV of the Federal University of the ABC (UFABC), Brazil, with the MSc. in Mechanical Engineering from the Faculty of Mechanical Engineering of the State University of Campinas (Unicamp - Brazil) in 2000 and PhD. degree in Mechanical Engineering by Polytechnic School of the University of São Paulo (USP-Brazil) in 2008. He is a researcher and author of several articles in magazines mechanical engineering. ORCID: 0000-0003-3258-4664

R. Reina-Muñoz, was born in Cali, Colombia. He received the BSc. degree from Valle University, Cali, Colombia, in 1989 and the MSc. and PhD. degrees from the University of São Paulo, São Paulo, Brazil, in 1993 and 1998, respectively, all in electrical engineering (microelectronics). He is currently a professor with the Federal University of ABC (UFABC), Santo André, Brazil. His current research interests include analog integrated circuits design, power management, FPGA, VHDL and Verilog, radio−frequency CMOS integrated circuits, and pressure sensors and pressure transmitters applications. ORCID: 0000-0002-5305-1595

Received: September 06, 2017; Revised: April 15, 2018; Accepted: July 15, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License