Tool Condition Monitoring Based on Nonlinear Output Frequency Response Functions and Multivariate Control Chart

Tool condition monitoring(TCM)is a key technology for intelligent manufacturing.The objective is to monitor the tool operation status and detect tool breakage so that the tool can be changed in time to avoid significant damage to workpieces and reduce manufacturing costs.Recently,an innovative TCM approach based on sensor data modelling and model frequency analysis has been proposed.Different from traditional signal feature-based monitoring,the data from sensors are utilized to build a dynamic process model.Then,the nonlinear output frequency response functions,a concept which extends the linear system frequency response function to the nonlinear case,over the frequency range of the tooth passing frequency of the machining process are extracted to reveal tool health conditions.In order to extend the novel sensor data modelling and model frequency analysis to unsupervised condition monitoring of cutting tools,in the present study,a multivariate control chart is proposed for TCM based on the frequency domain properties of machining processes derived from the innovative sensor data modelling and model frequency analysis.The feature dimension is reduced by principal component analysis first.Then the moving average strategy is exploited to generate monitoring variables and overcome the effects of noises.The milling experiments of titanium alloys are conducted to verify the effectiveness of the proposed approach in detecting excessive flank wear of solid carbide end mills.The results demonstrate the advantages of the new approach over conventional TCM techniques and its potential in industrial applications.

Nutrient Response Functions of Sorghum for Miesso District Central Rift Valley of Ethiopia

This study was executed to offer the basis for optimized profit from fertilizer use for sorghum yield and to determine robust crop nutrient response function and economic rate for the production of sorghum at Miesso Central Rift Valley of Ethiopia. Trails were conducted at six experimental sites, sorghum yield response to N and P fertilizers application and economically optimum rates of nitrogen (EONR) and phosphorus (EOPR) were evaluated on a vertisols within the semi-arid Miesso districts west Hararge zone of Oromia region. The nutrient rates in 2014 cropping season four levels of Nitrogen (N) alone, these levels with 20 kg·ha−1 Phosphorus (P) and without N, 69 kg·ha−1 N with three levels of P treatments including the zero control were evaluated. In 2015, cropping season similar rates of N alone, the same rate N with 20 kg·ha−1 P, 92 kg·ha−1 N with three rates of P including the zero control were evaluated. The treatments were arranged in a randomized complete block with three replications in factorial design. Nutrient responses of sorghum were determined using asymptotic quadratic plateau functions. The significantly highest nitrogen rate was 46 kg·ha−1 alone in 2014 season, which gave grain yield of 2.56 Mg·ha−1 with a maximum yield advantage of 43%. P rates in both seasons and combined (sites + seasons) were not significantly influenced sorghum yield. Nitrogen agronomic and partial factor productivity peaked at 23 kg N ha−1 but declined with increasing N rate. The EONR combined (sites + seasons) were 37, 45, 52 and 60 kg·ha−1 and for the profit to cost ratio (PCR) were 2.43, 3.65, 4.86 and 5.79 at difference cost to grain price ratios (CP) = 3.6, 2.3, 1.6 and 1.2 respectively at Miesso Ethiopia. Nitrogen application had economically profitable than P. The study concluded that the application of N at 37 or 60 kg N ha−1 to sorghum production could be economically profitable for those economically constrained farmers or economically not constrained farmers. Validation should be farther conducted on farmers’ fields for refining the results obtained.

Monte-Carlo Simulation of Response Functions for Natural Gamma-Rays in LaBr_3 Detector System with Complex Borehole Configurations

Usually, there are several methods, e.g. experiment, interpolation experiment-based, analytic function, and Monte-Carlo simulation, to calculate the response functions in LaBr3（Ce） detectors. In logging applications, the experiment-based methods cannot be adopted because of their limitations. Analytic function has the advantage of fast calculating speed, but it is very difficult to take into account many effects that occur in practical applications. On the contrary, Monte-Carlo simulation can deal with physical and geometric configurations very tactfully. It has a distinct advantage for calculating the functions with complex configurations in borehole. A new application of LaBr3（Ce） detector is in natural gamma-rays borehole spectrometer for uranium well logging. Calculation of response functions must consider a series of physical and geometric factors under complex logging conditions, including earth formations and its relevant parameters, different energies, material and thickness of the casings, the fluid between the two tubes, and relative position of the LaBr3（Ce） crystal to steel ingot at the front of logging tube. The present work establishes Monte-Carlo simulation models for the above-mentioned situations, and then performs calculations for main gamma-rays from natural radio-elements series. The response functions can offer experimental directions for the design of borehole detection system, and provide technique basis and basic data for spectral analysis of natural gamma-rays, and for sonrceless calibration in uranium quantitative interpretation.

Pilot study on the Non-Invasive Detectability of Femoral Neck Fractures with Frequency Response Functions

A suspicion of a femoral neck fracture is a frequently recurring situation, especially in nursing homes. For the clarification of such a suspicion normally imaging techniques are used. Such equipment is expensive and therefore is located in hospitals. In addition to the costs, a transport causes stress for the patient. This pilot study is devoted to the question whether the detection of a femoral neck fracture with vibration measurements is possible in principal. In such a case, the clarification could be done on-site by an ordinary person using much cheaper equipment. This would reduce the stress for the patient and save money. For this purpose vibration measurements on a dead body with intact, with partially fractured and with complete cut femoral neck have been performed. Two different methods for the vibration initiation have been investigated, the so called impact testing and the shaker testing. The frequency response function has been determined for all combinations on both sides of the body. It turned out that there is a clear difference in the frequency response functions of the fractured bone with respect to the intact bone when shaker testing is used. This indicates that the method could have the potential to be a cost-saving alternative to imaging techniques. However, in a next step a statistically reliable clinical survey on living persons needs to be done.

Bifurcation and Turing Pattern Formation in a Diffusion Modified Leslie-Gower Predator-Prey Model with Crowley-Martin Functional Response

In this paper, we study a modified Leslie-Gower predator-prey model with Smith growth subject to homogeneous Neumann boundary condition, in which the functional response is the Crowley-Martin functional response term. Firstly, for ODE model, the local stability of equilibrium point is given. And by using bifurcation theory and selecting suitable bifurcation parameters, we find many kinds of bifurcation phenomena, including Transcritical bifurcation and Hopf bifurcation. For the reaction-diffusion model, we find that Turing instability occurs. Besides, it is proved that Hopf bifurcation exists in the model. Finally, numerical simulations are presented to verify and illustrate the theoretical results.

Dynamics Analysis of Active Variable Stiffness Vibration Isolator for Whole-Spacecraft Systems Based on Nonlinear Output Frequency Response Functions

In order to improve the harsh dynamic environment experienced by heavy rockets during different external excitations,this study presents a novel active variable stiffness vibration isolator(AVS-VI)used as the vibration isolation device to reduce excessive vibration of the whole-spacecraft isolation system.The AVS-VI is composed of horizontal stiffness spring,positive stiffness spring,parallelogram linkage mechanism,piezoelectric actuator,acceleration sensor,viscoelastic damping,and PID active controller.Based on the AVS-VI,the generalized vibration transmissibility determined by the nonlinear output frequency response functions and the energy absorption rate is applied to analyze the isolation performance of the whole-spacecraft system with AVS-VI.The AVS-VI can conduct adaptive vibration suppression with variable stiffness to the whole-spacecraft system,and the analysis results indicate that the AVS-VI is efTective in reducing the extravagant vibration of the whole-spacecraft system,where the vibration isolation is decreased up to above 65%under different acceleration excitations.Finally,different parameters of AVS-VI are considered to optimize the whole-spacecraft system based on the generalized vibration transmissibility and the energy absorption rate.

TLR signaling that induces weak inflammatory response and SHIP1 enhances osteogenic functions

Toll-like receptor （TLR）-mediated inflammatory response could negatively affect bone metabolism. In this study, we determined how osteogenesis is regulated during inflammatory responses that are downstream of TLR signaling. Human primary osteoblasts were cultured in collagen gels. Pam3CSK4 （P3C） and Escherichia coli lipopolysaccharide （EcLPS） were used as TLR2 and TLR4 ligand respectively. Porphyromonas gingivalis LPS having TLR2 activity with either TLR4 agonism （Pg1690） or TLR4 antagonism （Pg1449） and mutant E. coli LPS （LPxE/LPxF/WSK） were used. IL-lp, SH2-containing inositol phosphatase-1 （SHIP1） that has regulatory roles in osteogenesis, alkaline phosphatase and mineralization were analyzed. 3α-Aminocholestane （3AC） was used to inhibit SHIP1. Our results suggest that osteoblasts stimulated by P3C, poorly induced IL-1β but strongly upregulated SHIP1 and enhanced osteogenic mediators. On the contrary, EcLPS significantly induced IL-1β and osteogenic mediators were not induced. While Pg1690 downmodulated osteogenic mediators, Pg1449 enhanced osteogenic responses, suggesting that TLR4 signaling annuls osteogenesis even with TLR2 activity. Interestingly, mutant E. coli LPS that induces weak inflammation upregulated osteogenesis, but SHIP1 was not induced. Moreover, inhibiting SHIP1 significantly upregulated TLR2-mediated inflammatory response and downmodulated osteogenesis. In conclusion, these results suggest that induction of weak inflammatory response through TLR2 （with SHIP1 activity） and mutant TLR4 ligands could enhance osteogenesis.

Improvement of a temperature response function for divertor heat flux monitoring in fusion devices

Temperature response functions have been developed to investigate sensor design and divertor heat flux estimation in magnetically confined plasmas. The time-dependent heat flux can be derived by fitting the response function to experimental thermocouple(TC) data. Because the TC signals have a time delay to transit events such as discharge start or confinement transition, the time delay is taken into account in a temperature response function. Such a function accurately describes the signal from each TC channel with time delay in a sensor test using a neutral beam injection. Measurement for commercial TCs shows that the time delay is caused by the finite heat capacity of TC wire and contact heat resistance between TC and target surface.

Stability Analysis of a Self-Memory Prey-Predator Diffusion Model Based on Bazykin Functional Response

To investigate the effects of self-memory diffusion on predator-prey models, we consider a predator-prey model with Bazykin functional response of self- memory diffusion. The uniqueness, boundedness, positivity, existence and stability of equilibrium point of the model are studied. In this paper, the uniqueness of the solution is discussed under the non-negative initial function and Neumann boundary conditions satisfying a specific space. The boundness of the solution is proved by the comparison principle of parabolic equations, and the positivity of the solution is proved by the strong maximum principle of parabolic equations. Hurwitz criterion and Lyapunov function construction are used to analyze the local stability and global stability of feasible equilibrium points. The results show that the system solution is unique non-negative and bounded. The model is unstable at the trivial equilibrium point E0 and the boundary equilibrium point E1, and the condition of whether the positive equilibrium point E2 is stable under certain conditions is given.

Dynamic Analysis of a Predator-Prey Model with Holling-II Functional Response

A predator-prey model with linear capture term Holling-II functional response was studied by using differential equation theory. The existence and the stabilities of non-negative equilibrium points of the model were discussed. The results show that under certain limited conditions, these two groups can maintain a balanced position, which provides a theoretical reference for relevant departments to make decisions on ecological protection.

A novel algorithm for evaluating cement azimuthal density based on perturbation theory in horizontal well

Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution and casing eccentricity in horizontal wells often complicates the accurate evaluation of cement azimuthal density.In this regard,this paper proposes an algorithm to calculate the cement azimuthal density in horizontal wells using a multi-detector gamma-ray detection system.The spatial dynamic response functions are simulated to obtain the influence of cement density on gamma-ray counts by the perturbation theory,and the contribution of cement density in six sectors to the gamma-ray recorded by different detectors is obtained by integrating the spatial dynamic response functions.Combined with the relationship between gamma-ray counts and cement density,a multi-parameter calculation equation system is established,and the regularized Newton iteration method is employed to invert casing eccentricity and cement azimuthal density.This approach ensures the stability of the inversion process while simultaneously achieving an accuracy of 0.05 g/cm^(3) for the cement azimuthal density.This accuracy level is ten times higher compared to density accuracy calculated using calibration equations.Overall,this algorithm enhances the accuracy of cement azimuthal density evaluation,provides valuable technical support for the monitoring of cement azimuthal density in the oil and gas industry.

Parametric characteristic of the random vibration response of nonlinear systems

Volterra series is a powerful mathematical tool for nonlinear system analysis,and there is a wide range of nonlinear engineering systems and structures that can be represented by a Volterra series model.In the present study,the random vibration of nonlinear systems is investigated using Volterra series.Analytical expressions were derived for the calculation of the output power spectral density（PSD） and input-output cross-PSD for nonlinear systems subjected to Gaussian excitation.Based on these expressions,it was revealed that both the output PSD and the input-output crossPSD can be expressed as polynomial functions of the nonlinear characteristic parameters or the input intensity.Numerical studies were carried out to verify the theoretical analysis result and to demonstrate the effectiveness of the derived relationship.The results reached in this study are of significance to the analysis and design of the nonlinear engineering systems and structures which can be represented by a Volterra series model.

Examining approaches for modeling individual tree growth response to thinning in Norway spruce

Using periodic measurements from permanent plots in non-thinned and thinned Norway spruce(Picea abies(L.)H.Karst.)stands in Norway,individual-tree growth models were developed to predict annual diameter increment,height increment,and height to crown base increment.Based on long-term data across a range of thinning regimes and stand conditions,alternative approaches for modeling response to treatment were assessed.Dynamic thinning response functions in the form of multiplicative modifiers that predict no effect at the time of thinning,a rapid increase followed by an early maximum before the effect gradually declines to zero could not be fitted to initially derived baseline models without thinning related predictors.However,alternative approaches were used and found to perform well.Specifically,indicator variables representing varying time periods after thinning were statistically significant and behaved in a robust manner as well as consistent with general expectations.In addition,they improved overall prediction accuracy when incorporated as fixed effects into the baseline models for diameter and height to crown base increment.Further,more simply,including exponentially decreasing multiplicative thinning response functions improved prediction accuracy for height increment and height to crown base increment.Irrespective of studied attribute and modelling approach,improvement in performance of these extended models was relatively limited when compared to the corresponding baseline models and more pronounced in trees from thinned stands.We conclude that the largely varying and often multi-year measurement intervals of the periodic data used in this study likely prevented the development of more sophisticated thinning response functions.However,based on the evaluation of the final models’overall performance such complex response functions may not to be necessary to reliably predict individual tree growth after thinning for certain conditions or species,which should be further considered in future analyses of similar nature.

Electromagnetic response function for the semi-annual variations estimated by the single-station method using the night-time values of the geomagnetic fields

The geomagnetic night-time values were used to estimate the electromagnetic response function Q1 for half-year period. If the spatial structure of the source field can be described by the approximation, one can estimate the Q1 value using the single-station Z/H method. This technique enables us to carry out regional deep gcomagnetic sounding by the method. The data used for analyses are geomagnetic night-time values for about, typically, 26 years from 5 good-quality stations and for several years from 34 stations distributed over the globe. The results indicate that the night-time values yield more reliable response estimates for half-year period compared to the usual estimates obtained from daily means. It implies that the approximation for the night-time fields holds good for the half-year period, but the daily means are not suitable for estimating the response function of the semi-annual variations by using the single-station method. Source field analyses for daily means data and night-time means data have also been carried out in this paper.

EXISTENCE OF POSITIVE SOLUTION FOR A TWO-PATCHES COMPETITION SYSTEM WITH DIFFUSION AND TIME DELAY AND FUNCTIONAL RESPONSE

By using the continuation theorem of coincidence theory, the existence of a positive periodic solution for a two patches competition system with diffusion and time delay and functional responsex [FK(W1*3/4。*2/3]′ 1 (t)=x 1(t)a 1(t)-b 1(t)x 1(t)-c 1(t)y(t)1+m(t)x 1(t)+D 1(t)[x 2(t)-x 1(t)], x [FK(W1*3/4。*2/3]′ 2 (t)=x 2(t)a 2(t)-b 2(t)x 2(t)-c 2(t)∫ 0 -τ k(s)x 2(t+s) d s+D 2(t)[x 1(t)-x 2(t)], y′(t)=y(t)a 3(t)-b 3(t)y(t)-c 3(t)x 1(t)1+m(t)x 1(t)is established, where a i(t),b i(t),c i(t)(i=1,2,3),m(t) and D i(t)(i=1,2) are all positive periodic continuous functions with period w >0, τ is a nonnegative constant and k(s) is a continuous nonnegative function on [- τ ,0].

Mechanism of Functional Responses to Loading of Carbon Fiber Reinforced Cement-based Composites

Single fiber pull-out testing was conducted to study the origin of the functional responses to loading of carbon fiber reinforced cement-based composites. The variation of electrical resistance with the bonding force on the fiber-matrix interface was measured. Single fiber electromechanical testing was also conducted by measuring the electrical resistance under static tension. Comparison of the results shows that the resistance increasing during single fiber pull-out is mainly due to the changes at the interface. The conduction mechanism of the composite can be explained by the tunneling model. The interfacial stress causes the deformation of interfacial structure and the interfacial debonding, which have influences on the tunneling effect and result in the change of resistance.

Uncertainty analysis of strain modal parameters by Bayesian method using frequency response function

Structural strain modes are able to detect changes in local structural performance, but errors are inevitably intermixed in the measured data. In this paper, strain modal parameters are considered as random variables, and their uncertainty is analyzed by a Bayesian method based on the structural frequency response function （FRF）. The estimates of strain modal parameters with maximal posterior probability are determined. Several independent measurements of the FRF of a four-story reinforced concrete flame structural model were performed in the laboratory. The ability to identify the stiffness change in a concrete column using the strain mode was verified. It is shown that the uncertainty of the natural frequency is very small. Compared with the displacement mode shape, the variations of strain mode shapes at each point are quite different. The damping ratios are more affected by the types of test systems. Except for the case where a high order strain mode does not identify local damage, the first order strain mode can provide an exact indication of the damage location.

Research on the iterative method for model updating based on the frequency response function

Model reduction technique is usually employed in model updating process. In this paper, a new model updat- ing method named as cross-model cross-frequency response function （CMCF） method is proposed and a new iterative method associating the model updating method with the mo- del reduction technique is investigated. The new model up- dating method utilizes the frequency response function to avoid the modal analysis process and it does not need to pair or scale the measured and the analytical frequency re- sponse function, which could greatly increase the number of the equations and the updating parameters. Based on the traditional iterative method, a correction term related to the errors resulting from the replacement of the reduction ma- trix of the experimental model with that of the finite element model is added in the new iterative method. Comparisons be- tween the traditional iterative method and the proposed itera- tive method are shown by model updating examples of solar panels, and both of these two iterative methods combine the CMCF method and the succession-level approximate reduc- tion technique. Results show the effectiveness of the CMCF method and the proposed iterative method .

Effects of Temperature on Functional Response of Anagrus nilaparvatae Pang et Wang (Hymenoptera:Mymaridae) on the Eggs of Whitebacked Planthopper,Sogatella furcifera Horvth and Brown Planthopper,Nilaparvata lugens Stl

Understanding the temperature affecting parasitic efficiency is critical to succeed in utilizing parasitoid as natural enemy in pest management. Laboratory studies were carried out to determine the effects of temperature on parasitoid preference of female Anagrus nilaparvatae Pang et Wang （Hymenoptera：Mymaridae） to the eggs of whitebacked planthopper （WBPH）, Sogatella furcifera Horváth and brown planthopper （BPH）, Nilaparvata lugens Stl to build a composite model describing changes in parasitic response along a temperature gradient （18, 22, 26, 30, 34°C）. The results showed that attack responses of A. nilaparvatae on WBPH and BPH were the best described by a Type II functional response. The two parameters, attack rates （a） and handling times （Th）, of A. nilaparvatae to both eggs were influenced by the temperature. The maximum attack rates to WBPH （1.235） and BPH （1.049） were at 26 and 34°C, respectively, and the shortest handling times to WBPH （0.063） and BPH （0.057） were at 30 and 26°C, respectively. However, the optimal temperature for parasitic efficiency of A. nilaparvatae to WBPH and BPH eggs was both at 26°C, which showed that the present microclimate temperature of the habitat in the paddyfield was beneficial to A. nilaparvatae and indicated that parasitic efficiency of A. nilaparvatae would be impaired by global warming.

Response Characteristics of Load on Vessels in Waves

Considering the requirement of direct design and fatigue test for ships and floating structures by use of FEM technique, a computational procedure of spectral analysis for wave load on the hull surface is developed in this paper. The response of hydrodynamic pressure on the body surface to a designated sea state for ships and floating structures is calculated by use of the revised strip method with the hull bound perturbation flow concept. The spectral function of wave load for the defined point on the body surface can be determined from the Wiener-Khinhine theorem and the characteristic load value can be also obtained from spectral moment analysis. A container ship is taken as a computational example acid the sample of wave load with a certain probability and corresponding encountered frequency is provided.