Symmetry of in-shoe force signatures in athletes under field conditions : an application for injury prevention

The aetiology behind overuse injuries such as stress fractures is complex and multi-factorial. In sporting events where the loading is likely to be uneven (e.g. hurdling and jumps), research has suggested that the frequency of stress fractures seems to favour the athlete’s dominant limb. The tendency for an individual to have a preferred limb for voluntary motor acts makes limb selection a possible factor behind the development of unilateral overuse injuries, particularly when repeatedly used during high loading activities. The event of sprint hurdling is well suited for the study of loading asymmetry as the hurdling technique is repetitive and the limb movement asymmetrical. Of relevance to this study is the high incidence of Navicular Stress Fractures (NSF) in hurdlers, with suggestions there is a tendency for the fracture to develop in the trail leg foot, although this is not fully accepted. The Ground Reaction Force (GRF) with each foot contact is influenced by the hurdle action, with research finding step-to-step loading variations. However, it is unknown if this loading asymmetry extends to individual forefoot joints, thereby influencing stress fracture development. The first part of the study involved a series of investigations using a commercially available matrix style in-shoe sensor system (FscanTM, Tekscan Inc.). The suitability of insole sensor systems and custom made discrete sensors for use in hurdling-related training activities was assessed. The methodology used to analyse foot loading with each technology was investigated. The insole and discrete sensors systems tested proved to be unsuitable for use during full pace hurdling. Instead, a running barrier task designed to replicate the four repetitive foot contacts present during hurdling was assessed. This involved the clearance of a series of 6 barriers (low training hurdles), place in a straight line, using 4 strides between each. The second part of the study involved the analysis of "inter-limb" and "within foot loading asymmetries" using stance duration as well as vertical GRF under the Hallux (T1), the first metatarsal head (M1) and the central forefoot peak pressure site (M2), during walking, running, and running with barrier clearances. The contribution to loading asymmetry that each of the four repetitive foot contacts made during a series of barrier clearances was also assessed. Inter-limb asymmetry, in forefoot loading, occurred at discrete forefoot sites in a non-uniform manner across the three gait conditions. When the individual barrier foot contacts were compared, the stance duration was asymmetrical and the proportion of total forefoot load at M2 was asymmetrical. There were no significant differences between the proportion of forefoot load at M1, compared to M2; for any of the steps involved in the barrier clearance. A case study testing experimental (discrete) sensors during full pace sprinting and hurdling found that during both gait conditions, the trail limb experienced the greater vertical GRF at M1 and M2. During full pace hurdling, increased stance duration and vertical loading was a characteristic of the trail limb hurdle foot contacts. Commercially available in-shoe systems are not suitable for on field assessment of full pace hurdling. For the use of discrete sensor technology to become commonplace in the field, more robust sensors need to be developed.

Application of high voltage, high frequency pulsed electromagnetic field on cortical bone tissue

Over the last few decades, electric and electromagnetic fields have achieved important role as stimulator and therapeutic facility in biology and medicine. In particular, low magnitude, low frequency, pulsed electromagnetic field has shown significant positive effect on bone fracture healing and some bone diseases treatment. Nevertheless, to date, little attention has been paid to investigate the possible effect of high frequency, high magnitude pulsed electromagnetic field (pulse power) on functional behaviour and biomechanical properties of bone tissue. Bone is a dynamic, complex organ, which is made of bone materials (consisting of organic components, inorganic mineral and water) known as extracellular matrix, and bone cells (live part). The cells give the bone the capability of self-repairing by adapting itself to its mechanical environment. The specific bone material composite comprising of collagen matrix reinforced with mineral apatite provides the bone with particular biomechanical properties in an anisotropic, inhomogeneous structure. This project hypothesized to investigate the possible effect of pulse power signals on cortical bone characteristics through evaluating the fundamental mechanical properties of bone material. A positive buck-boost converter was applied to generate adjustable high voltage, high frequency pulses up to 500 V and 10 kHz. Bone shows distinctive characteristics in different loading mode. Thus, functional behaviour of bone in response to pulse power excitation were elucidated by using three different conventional mechanical tests applying three-point bending load in elastic region, tensile and compressive loading until failure. Flexural stiffness, tensile and compressive strength, hysteresis and total fracture energy were determined as measure of main bone characteristics. To assess bone structure variation due to pulse power excitation in deeper aspect, a supplementary fractographic study was also conducted using scanning electron micrograph from tensile fracture surfaces. Furthermore, a non-destructive ultrasonic technique was applied for determination and comparison of bone elasticity before and after pulse power stimulation. This method provided the ability to evaluate the stiffness of millimetre-sized bone samples in three orthogonal directions. According to the results of non-destructive bending test, the flexural elasticity of cortical bone samples appeared to remain unchanged due to pulse power excitation. Similar results were observed in the bone stiffness for all three orthogonal directions obtained from ultrasonic technique and in the bone stiffness from the compression test. From tensile tests, no significant changes were found in tensile strength and total strain energy absorption of the bone samples exposed to pulse power compared with those of the control samples. Also, the apparent microstructure of the fracture surfaces of PP-exposed samples (including porosity and microcracks diffusion) showed no significant variation due to pulse power stimulation. Nevertheless, the compressive strength and toughness of millimetre-sized samples appeared to increase when the samples were exposed to 66 hours high power pulsed electromagnetic field through screws with small contact cross-section (increasing the pulsed electric field intensity) compare to the control samples. This can show the different load-bearing characteristics of cortical bone tissue in response to pulse power excitation and effectiveness of this type of stimulation on smaller-sized samples. These overall results may address that although, the pulse power stimulation can influence the arrangement or the quality of the collagen network causing the bone strength and toughness augmentation, it apparently did not affect the mineral phase of the cortical bone material. The results also confirmed that the indirect application of high power pulsed electromagnetic field at 500 V and 10 kHz through capacitive coupling method, was athermal and did not damage the bone tissue construction.

Assessing occupational stress in psychiatric nurses using the full job strain model the value of social support to nurses

The present research examined the effects of occupational stress in psychiatric nursing on employee well!being using the full Job Strain Model.The Job Strain Model was assessed for its ability to predict employee well!being in terms of job satisfaction and mental health. The original Job Strain Model was expanded to include social support based on previous research concerning the impact of social support on well!being[ In the present study\ both work support and non-work were assessed for their contribution to wellbeing.The results of this study indicate that the full Job Strain Model can be used to significantly predict job satisfaction and mental health in this sample of Australian psychiatric nurses. Furthermore social support was shown to be an important component of the Job Strain Model.

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.

Evaluation of strain energy based damage indices for flexural crack detection of RC beams

Vibration Based Damage Identification Techniques which use modal data or their functions, have received significant research interest in recent years due to their ability to detect damage in structures and hence contribute towards the safety of the structures. In this context, Strain Energy Based Damage Indices (SEDIs), based on modal strain energy, have been successful in localising damage in structuers made of homogeneous materials such as steel. However, their application to reinforced concrete (RC) structures needs further investigation due to the significant difference in the prominent damage type, the flexural crack. The work reported in this paper is an integral part of a comprehensive research program to develop and apply effective strain energy based damage indices to assess damage in reinforced concrete flexural members. This research program established (i) a suitable flexural crack simulation technique, (ii) four improved SEDI's and (iii) programmable sequentional steps to minimise effects of noise. This paper evaluates and ranks the four newly developed SEDIs and existing seven SEDIs for their ability to detect and localise flexural cracks in RC beams. Based on the results of the evaluations, it recommends the SEDIs for use with single and multiple vibration modes.

Damage assessment in reinforced concrete flexural members using modal strain energy based method

Damage assessment (damage detection, localization and quantification) in structures and appropriate retrofitting will enable the safe and efficient function of the structures. In this context, many Vibration Based Damage Identification Techniques (VBDIT) have emerged with potential for accurate damage assessment. VBDITs have achieved significant research interest in recent years, mainly due to their non-destructive nature and ability to assess inaccessible and invisible damage locations. Damage Index (DI) methods are also vibration based, but they are not based on the structural model. DI methods are fast and inexpensive compared to the model-based methods and have the ability to automate the damage detection process. DI method analyses the change in vibration response of the structure between two states so that the damage can be identified. Extensive research has been carried out to apply the DI method to assess damage in steel structures. Comparatively, there has been very little research interest in the use of DI methods to assess damage in Reinforced Concrete (RC) structures due to the complexity of simulating the predominant damage type, the flexural crack. Flexural cracks in RC beams distribute non- linearly and propagate along all directions. Secondary cracks extend more rapidly along the longitudinal and transverse directions of a RC structure than propagation of existing cracks in the depth direction due to stress distribution caused by the tensile reinforcement. Simplified damage simulation techniques (such as reductions in the modulus or section depth or use of rotational spring elements) that have been extensively used with research on steel structures, cannot be applied to simulate flexural cracks in RC elements. This highlights a big gap in knowledge and as a consequence VBDITs have not been successfully applied to damage assessment in RC structures. This research will address the above gap in knowledge and will develop and apply a modal strain energy based DI method to assess damage in RC flexural members. Firstly, this research evaluated different damage simulation techniques and recommended an appropriate technique to simulate the post cracking behaviour of RC structures. The ABAQUS finite element package was used throughout the study with properly validated material models. The damaged plasticity model was recommended as the method which can correctly simulate the post cracking behaviour of RC structures and was used in the rest of this study. Four different forms of Modal Strain Energy based Damage Indices (MSEDIs) were proposed to improve the damage assessment capability by minimising the numbers and intensities of false alarms. The developed MSEDIs were then used to automate the damage detection process by incorporating programmable algorithms. The developed algorithms have the ability to identify common issues associated with the vibration properties such as mode shifting and phase change. To minimise the effect of noise on the DI calculation process, this research proposed a sequential order of curve fitting technique. Finally, a statistical based damage assessment scheme was proposed to enhance the reliability of the damage assessment results. The proposed techniques were applied to locate damage in RC beams and slabs on girder bridge model to demonstrate their accuracy and efficiency. The outcomes of this research will make a significant contribution to the technical knowledge of VBDIT and will enhance the accuracy of damage assessment in RC structures. The application of the research findings to RC flexural members will enable their safe and efficient performance.

The time course of in vivo recovery of transverse strain in high-stress tendons following exercise

Objective To evaluate the time course of the recovery of transverse strain in the Achilles and patellar tendon following a bout of resistance exercise. Methods Seventeen healthy adults underwent sonographic examination of the right patellar (n=9) and Achilles (n=8) tendons immediately prior to and following 90 repetitions of weight-bearing quadriceps and gastrocnemius-resistance exercise performed against an effective resistance of 175% and 250% body weight, respectively. Sagittal tendon thickness was determined 20 mm from the enthesis and transverse strain, as defined by the stretch ratio, was repeatedly monitored over a 24 h recovery period. Results Resistance exercise resulted in an immediate decrease in Achilles (t7=10.6, p<0.01) and patellar (t8=8.9, p<0.01) tendon thickness, resulting in an average transverse stretch ratio of 0.86±0.04 and 0.82±0.05, which was not significantly different between tendons. The magnitude of the immediate transverse strain response, however, was reduced with advancing age (r=0.63, p<0.01). Recovery in transverse strain was prolonged compared with the duration of loading and exponential in nature. The average primary recovery time was not significantly different between the Achilles (6.5±3.2 h) and patellar (7.1±3.2 h) tendons. Body weight accounted for 62% and 64% of the variation in recovery time, respectively. Conclusions Despite structural and biochemical differences between the Achilles and patellar tendon, the mechanisms underlying transverse creep recovery in vivo appear similar and are highly time dependent. These novel findings have important implications concerning the time required for the mechanical recovery of high-stress tendons following an acute bout of exercise.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

Fracture analysis of wind turbine main shaft

Shaft fracture at an early stage of operation is a common problem for a certain type of wind turbine. To determine the cause of shaft failure a series of experimental tests were conducted to evaluate the chemical composition and mechanical properties. A detail analysis involving macroscopic feature and microstructure analysis of the material of the shaft was also performed to have an in depth knowledge of the cause of fracture. The experimental tests and analysis results show that there are no significant differences in the material property of the main shaft when comparing it with the Standard, EN10083-3:2006. The results show that stress concentration on the shaft surface close to the critical section of the shaft due to rubbing of the annular ring and coupled with high stress concentration caused by the change of inner diameter of the main shaft are the main reasons that result in fracture of the main shaft. In addition, inhomogeneity of the main shaft micro-structure also accelerates up the fracture process of the main shaft. In addition, the theoretical calculation of equivalent stress at the end of the shaft was performed, which demonstrate that cracks can easily occur under the action of impact loads. The contribution of this paper is to provide a reference in fracture analysis of similar main shaft of wind turbines.

Acute hamstring strain injury in track and field athletes : a 3-year observational study at the Penn Relay Carnival

Purpose To observe the incidence rates of hamstring strain injuries (HSIs) across different competition levels and ages during the Penn Relays Carnival. Methods Over a 3-year period all injuries treated by the medical staff were recorded. The type of injury, anatomic location, event in which the injury occurred, competition level and demographic data were documented. Absolute and relative HSI (per 1000 participants) were determined and odds ratios (OR) were calculated between genders, competition levels and events. Results Throughout the study period 48,473 athletes registered to participate in the Penn Relays Carnival, with 118 HSIs treated by the medical team. High school females displayed lesser risk of HSI than high school males (OR = 0.55, p = 0.021), and masters athletes were more likely than high school (OR = 4.26, p < 0.001) and college (OR = 3.55, p = 0.001) level athletes to suffer a HSI. The 4x400m relay displayed a greater likelihood of HSI compared to the 4x100m relay (OR = 1.77, p = 0.008). Conclusions High school males and masters levels athletes are most likely to suffer HSI, and there is higher risk in 400m events compared to 100m events.

Predicting crack initiation due to ratchetting in rail heads using critical element analysis

This paper presents a strategy to predict the lifetime of rails subjected to large rolling contact loads that induce ratchetting strains in the rail head. A critical element concept is used to calculate the number of loading cycles needed for crack initiation to occur in the rail head surface. In this technique the finite element method (FEM) is used to determine the maximum equivalent ratchetting strain per load cycle, which is calculated by combining longitudinal and shear stains in the critical element. This technique builds on a previously developed critical plane concept that has been used to calculate the number of cycles to crack initiation in rolling contact fatigue under ratchetting failure conditions. The critical element concept simplifies the analytical difficulties of critical plane analysis. Finite element analysis (FEA) is used to identify the critical element in the mesh, and then the strain values of the critical element are used to calculate the ratchetting rate analytically. Finally, a ratchetting criterion is used to calculate the number of cycles to crack initiation from the ratchetting rate calculated.

Heat strain and hydration status of surface mine blast crew workers

Objective Dehydration and symptoms of heat illness are common among the surface mining workforce. This investigation aimed to determine whether heat strain and hydration status exceeded recommended limits. Methods Fifteen blast crew personnel operating in the tropics were monitored across a 12-hour shift. Heart rate, core body temperature, and urine-specific gravity were continuously recorded. Participants self-reported fluid consumption and completed a heat illness symptom inventory. Results Core body temperature averaged 37.46 +/- 0.13[degrees]C, with the group maximum 37.98 +/- 0.19[degrees]C. Mean urine-specific gravity was 1.024 +/- 0.007, with 78.6% of samples 1.020 or more. Seventy-three percent of workers reported at least one symptom of heat illness during the shift. Conclusions Core body temperature remained within the recommended limits; however, more than 80% of workers were dehydrated before commencing the shift, and tended to remain so for the duration.

Hybrid genetic algorithm for elimination of severe stress concentration in railhead ends

This paper presents a computational method for eliminating severe stress concentration at the unsupported railhead ends in rail joints through innovative shape optimization of the contact zone, which is complex due to near field nonlinear contact. With a view to minimizing the computational efforts, hybrid genetic algorithm method coupled with parametric finite element has been developed and compared with the traditional genetic algorithm (GA). The shape of railhead top surface where the wheel contacts nonlinearly was optimized using the hybridized GA method. Comparative study of the optimal result and the search efficiency between the traditional and hybrid GA methods has shown that the hybridized GA provides the optimal shape in fewer computational cycles without losing accuracy. The method will be beneficial to solving complex engineering problems involving contact nonlinearity.

Influence of tilt and surface roughness on the outflow wind field of an impinging jet

A physical and numerical steady flow impinging jet has been used to simulate the bulk characteristics of a downburst-like wind field. The influence of downdraft tilt and surface roughness on the ensuing wall jet flow has been investigated. It was found that a simulated downdraft impinging the surface at a non-normal angle has the potential for causing larger structural loads than the normal impingement case. It was also found that for the current impinging jet simulations, surface roughness played a minor role in determining the storm maximum wind structure, but this influence increased as the wall jet diverged. However, through comparison with previous research it was found that the influence of surface roughness is Reynolds number dependent and therefore may differ from that reported herein for full-scale downburst cases. Using the current experimental results an empirical model has been developed for laboratory-scale impinging jet velocity structure that includes the influence of both jet tilt and surface roughness.

The effect of high voltage, high frequency pulsed electric field on slain ovine cortical bone

High power, high frequency pulsed electric fields known as pulsed power (PP) has been applied recently in biology and medicine. However, little attention has been paid to investigate the application of pulse power in musculoskeletal system and its possible effect on functional behavior and biomechanical properties of bone tissue. This paper presents the first research investigating whether or not PP can be applied safely on bone tissue as a stimuli and what will be the possible effect of these signals on the characteristics of cortical bone by comparing the mechanical properties of this type of bone pre and post expose to PP and in comparison with the control samples. A positive buck‑boost converter was applied to generate adjustable high voltage, high frequency pulses (up to 500 V and 10 kHz). The functional behavior of bone in response to pulse power excitation was elucidated by applying compressive loading until failure. The stiffness, failure stress (strength) and the total fracture energy (bone toughness) were determined as a measure of the main bone characteristics. Furthermore, an ultrasonic technique was applied to determine and comprise bone elasticity before and after pulse power stimulation. The elastic property of cortical bone samples appeared to remain unchanged following exposure to pulse power excitation for all three orthogonal directions obtained from ultrasonic technique and similarly from the compression test. Nevertheless, the compressive strength and toughness of bone samples were increased when they were exposed to 66 h of high power pulsed electromagnetic field compared to the control samples. As the toughness and the strength of the cortical bone tissue are directly associated with the quality and integrity of the collagen matrix whereas its stiffness is primarily related to bone mineral content these overall results may address that although, the pulse power stimulation can influence the arrangement or the quality of the collagen network causing the bone strength and toughness augmentation, it apparently did not affect the mineral phase of the cortical bone material. The results also confirmed that the indirect application of high power pulsed electric field at 500 V and 10 kHz through capacitive coupling method was safe and did not destroy the bone tissue construction.