Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

Time-Domain Analysis of Body Freedom Flutter Based on 6DOF Equation

The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.

An Improved Time-Domain Inverse Technique for Localization and Quantifcation of Rotating Sound Sources

The time-domain inverse technique based on the time-domain rotating equivalent source method has been proposed to localize and quantify rotating sound sources. However, this technique encounters two problems to be addressed: one is the time-consuming process of solving the transcendental equation at each time step, and the other is the difculty of controlling the instability problem due to the time-varying transfer matrix. In view of that, an improved technique is proposed in this paper to resolve these two problems. In the improved technique, a de-Dopplerization method in the time-domain rotating reference frame is frst applied to eliminate the Doppler efect caused by the source rotation in the measured pressure signals, and then the restored pressure signals without the Doppler efect are used as the inputs of the time-domain stationary equivalent source method to locate and quantify sound sources. Compared with the original technique, the improved technique can avoid solving the transcendental equation at each time step, and facilitate the treatment of the instability problem because the transfer matrix does not change with time. Numerical simulation and experimental results show that the improved technique can eliminate the Doppler efect efectively, and then localize and quantify the rotating nonstationary or broadband sources accurately. The results also demonstrate that the improved technique can guarantee a more stable reconstruction and compute more efciently than the original one.

Integrated system of traditional THz time-domain spectroscopy and asynchronous optical sampling

Terahertz time-domain spectroscopy(THz-TDS)system,as a new means of spectral analysis and detection,plays an increasingly pivotal role in basic scientific research.However,owing to the long scanning time of the traditional THz-TDS system and the complex control of the asynchronous optical scanning(ASOPS)system,which requires frequent calibration,we combine traditional THz-TDS and ASOPS systems to form a composite system and propose an all-fiber trigger signal generation method based on the time overlapping interference signal generated by the collinear motion of two laser pulses.Finally,the time-domain and frequency-domain spectra are obtained by using two independent systems in the integrated systems.It is found that the full width at half maximum(FWHM)of the time-domain spectra and the spectral width of the frequency-domain spectra are almost the same,but the sampling speed of the ASOPS system is significantly faster than that of the traditional THz-TDS system,which conduces to the study of the transient characteristics of substances.

PC-based artif icial neural network inversion for airborne time-domain electromagnetic data

Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.

Solution and Analysis of Chatter Stability for End Milling in the Time-domain

In this paper, the instantaneous undeformed chip thickness is modeled to include the dynamic modulation caused by the tool vibration while the dynamic regenerative effects are taken into account. The numerical method is used to solve the differential equations goveming the dynamics of the milling system. Several chatter detection criteria are applied synthetically to the simulated signals and the stability diagram is obtained in time-domain. The simulation results in time-domain show a good agreement with the analytical prediction, which is validated by the cutting experiments. By simulating the chatter stability lobes in the time-domain and analyzing the influences of different spindle speeds on the vibration amplitudes of the tool under a Fixed chip-load condition, conclusions could be drawn as follows： In rough milling, higher machining efficiency can be achieved by selecting a spindle speed corresponding to the axial depth of cut in accordance with the simulated chatter stability lobes, and in Fmish milling, lower surface roughness can be achieved by selecting a spindle speed well beyond the resonant frequency of machining system.

Comparison of the time-domain electromagnetic field from an infinitesimal point charge and dipole source

An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.

Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation

To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.

Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.

Three-dimensional arbitrarily anisotropic modeling for time-domain airborne electromagnetic surveys

Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Investigation of copper sulfate pentahydrate dehydration by terahertz time-domain spectroscopy

Copper sulfate pentahydrate is investigated by terahertz time-domain spectroscopy. It is shown that the terahertz absorption coefficients are correlated with the particle size of the samples, as well as the heating rates of the ambient temperature. Furthermore, the water molecules of copper sulfate pentahydrate can be quantitatively characterized due to the high sensitivity of the terahertz wave to water molecules. Based on such results, the status of water incorporated in mineral opal is also characterized using terahertz time-domain spectroscopy. It indicates that terahertz technology can be considered as an efficient method to detect the dehydration of minerals.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

Investigating path loss characteristics of UWB signals in vacancy indoor environment based on time-domain technique

The path loss analysis model based on 5 rays in vacancy indoor environment is proposed. The relationship between multipath overlapping and the path loss is analyzed mathematically. Time-domain technique is introduced to compute reflection coefficient in a very short time interval. A 5 rays path loss calculation method, which is satisfactory accurate, is developed. 5 typical environments are involved to analyze and generalize the common path loss characteristics in vacancy indoor environment. The simulation result shows that the path loss can be characterized as 3 zones with different path loss exponent as distance between transmitter and receiver increasing.

Broadband terahertz time-domain spectroscopy and fast FMCW imaging: Principle and applications

We report a broadband terahertz time-domain spectroscopy(THz-TDS)which enables twenty vibrational modes of adenosine nucleoside to be resolved in a wide frequency range of 1-20 THz.The observed spectroscopic features of adenosine are in good agreement with the published spectra obtained using Fourier transform infrared spectroscopy(FTIR)and Raman spectroscopy.This much extended bandwidth leads to enhanced material characterization capability as it provides spectroscopic information on both intra-and inter-molecular vibrations.In addition,we also report a low-cost frequency modulation continuous wave(FMCW)imaging system which has a fast measurement speed of 40000 waveforms per second.Cross-sectional imaging capability through cardboard has also been demonstrated using its excellent penetration capability at a frequency range of 76-81 GHz.We anticipate that the integration of these two complementary imaging technologies would be highly desirable for many real-world applications because it provides both spectroscopic discrimination and penetration capabilities in a single instrument.

Time-Domain Simulation for Coupled Motions of Three Barges Moored Side-by-Side in Floatover Operation

Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.

Fast Evaluation of Time-Domain Green Function for Finite Water Depth

For computation of large amplitude motions of ships fastened to a dock, a fast evaluation scheme is implemented for computation of the time-domain Green function for finite water depth. Based on accurate evaluation of the Green function directly, a fast approximation method for the Green function is developed by use of Chebyshev polynomials. Examinations are carried out of the accuracy of the Green function and its derivatives from the scheme. It is shown that when an appropriate number of polynomial terms are used, very accurate approximation can be obtained.

Time-Domain Analysis of Underground Station-Layered Soil Interaction Based on High-Order Doubly Asymptotic Transmitting Boundary

Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.

Brillouin optical time-domain analysis for geotechnical monitoring

In this paper, we show some recent experimental applications of Brillouin optical time-domain analysis （BOTDA） based sensors for geotechnical monitoring. In particular, how these sensors can be applied to detecting early movements of soil slopes by the direct embedding of suitable fiber cables in the ground is presented. Furthermore, the same technology can be used to realize innovative inclinometers, as well as smart foundation anchors.

Terahertz Time-Domain Spectroscopy method for optical parameter extraction of plastic materials

A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.