A Time-Domain Numerical Simulation for Free Motion Responses of Two Ships Advancing in Head Waves

The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.

Time-Domain Higher-Order Boundary Element Method for Simulating High Forward-Speed Ship Motions in Waves

The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.

Inversion of Seabed Geotechnical Properties in the Arctic Chukchi Deep Sea Basin Based on Time Domain Adaptive Search Matching Algorithm

The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.

Finite Time Domain Dynamics of Dirac Fields

In the present paper we consider the case of a Dirac field in a finite time domain and coupled to an external field. We decompose the field and its Hamiltonian in terms of creation and annihilation operators and path integrate it via Grassmannian variables techniques. In that way we obtain its finite time domain Green function. We use it in the perturbative study of the interaction of Dirac particles with classical electromagnetic waves.

Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.

Simulation of pavement roughness based on time domain model

In order to describe pavement roughness more intuitively and effectively, a method of pavement roughness simulation, i.e., the stochastic sinusoidal wave, is introduced. The method is based on the primary idea that pavement roughness is denoted as the sum of numerous sines or cosines with stochastic phases, and uses the discrete spectrum to approach the target stochastic process. It is a discrete numerical method used to simulate pavement roughness. According to a given pavement power spectral density （PSD） coefficient, under the condition that the character of displacement frequency based on the time domain model is in accordance with the given pavement surface spectrum, the pavement roughness is optimized to stochastic equivalent vibrations by computer simulation, and the curves that describe pavement roughness under each grade are obtained. The results show that the stochastic sinusoidal wave is suitable for simulation of measured pavement surface spectra based on the time domain model. The method of the stochastic sinusoidal wave is important to the research on vehicle ride comfort due to its rigorous mathematical derivation, extensive application range and intuitive simulation curve. Finally, a roughness index defined as the nominal roughness index （NRI） is introduced, and it has correlation with the PSD coefficient.

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.

A Method in Mixed Frequency-time domain for Fatigue Life Estimation of Steel Girders of Cable-stayed Bridges Due to Buffeting

发展了斜拉桥钢主梁抖振疲劳寿命的混合时－频域方法．在本方法中，首先建立了桥面处的风速风向概率密度函数，导出的应力谱因其包括背景分量而不再是窄带谱．据应力谱模拟出了疲劳应力的时间历程，并用雨流计数法统计了其幅值特性．基于修正的Ｍｉｎｅｒ定理导出了考虑风速风向作用的桥梁主梁抖振疲劳寿命的估算公式，并用此分析了杨浦大桥钢主梁的抖振疲劳寿命．结果表明，风向对抖振疲劳寿命有很大影响；杨浦大桥主梁的疲劳寿命远大于设计寿命．

Full-Wave Analysis of Slotline Using Time-Domain Finite-Difference Method

The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.

Fully Coupled Effects of Hull,Mooring and Risers Model in Time Domain Based on An Innovative Deep Draft Multi-Spar

The coupled hull, mooring and riser analysis techniques in time domain are widely recognized as the unique approach to predict the accurate global motions. However, these complex issues have not been perfectly solved due to a large number of nonlinear factors, e.g. forces nonlinearity, mooring nonlinearity, motion nonlinearity and so on. This paper investigates the coupled effects through the numerical uncoupled model, mooring coupled model and fully coupled model accounting mooring and risers based on a novel deep draft multi-spar which is especially designed for deepwater in 2009. The numerical static-offset, free-decay, wind-action tests are executed, and finally three hours simulations are conducted under 100-year return period of GOM conditions involving wave, wind and current actions. The damping contributions, response characteristics and mooring line tensions are emphatically studied.

Statistical moment-based structural damage detection method in time domain

A novel structural damage detection method with a new damage index,i.e.,the statistical moment-based damage detection（SMBDD） method in the frequency domain,has been recently proposed.The aim of this study is to extend the SMBDD method in the frequency domain to the time domain for building structures subjected to non-Gaussian and non-stationary excitations.The applicability and effectiveness of the SMBDD method in the time domainis verified both numerically and experimentally.Shear buildings with various damage scenarios are first numerically investigated in the time domain taking into account the effect of measurement noise.The applicability of the proposed method in the time domain to building structures subjected to non-Gaussian and non-stationary excitations is then experimentally investigated through a series of shaking table tests,in which two three-story shear building models with four damage scenarios aretested.The identified damage locations and severities are then compared with the preset values.The comparative results are found to be satisfactory,and the SMBDD method is shown to be feasible and effective for building structures subjected to non-Gaussian and non-stationary excitations.

3D viscous-spring artificial boundary in time domain

After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.

Study on generalized magneto-thermoelastic problems by FEM in time domain

This paper presents an investigation of temperature, displacement, stress, and induced magnetic field in a half space perfectly-conductive plate. Finite element equations regarding generalized magneto-thermoelasticity problems with two relaxation times （i.e., the G-L theory） are derived using the principle of virtual work. For avoiding numerical complication involved in inverse Laplace and Fourier transformation and low precision thereof, the equations are solved directly in time-domain. As a numerical example, the derived equation is used to investigate the generalized magneto-thermoelastic behavior of a semi-infinite plate under magnetic field and subjecting to a thermal shock loading. The results demonstrate that FEM can faithfully predict the deformation of the plate and the induced magnetic field, and most importantly can reveal the sophisticated second sound effect of heat conduction in two-dimensional generalized thermoelastic solids, which is usually difficult to model by routine transformation methods. A peak can be observed in the distribution of stress and induced front and the magnitude of magnetic field at the heat wave the peak decreases with time, which can not be obtained by transformation methods. The new method can also be used to study generalized piezo-thermoelastic problems.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.

Time-and Frequency-Domain Analysis of Ice-Induced Displacement Responses Offshore Fixed Platforms

On the basis of ice- induced forced vibration model, ice- induced displacement responses of offshore fixed platforms are investigated in both time domain and frequency domain. The relationships of ice-induced displacement responses with ice breaking modes, ice acting directions and platform structures are analyzed and determined. The results lead to an important conclusion obtained for the first time that ice breaking frequency and the natural frequency of the first mode of the platform are the two main factors that dominate the degree of vibration. The present work provides a firm basis for both design and operation of fixed platforms against ice loading.

Digital coherent detection research on Brillouin optical time domain reflectometry with simplex pulse codes

The digital coherent detection technique has been investigated without any frequency-scanning device in the Brillouin optical time domain reflectometry （BOTDR）, where the simplex pulse codes are applied in the sensing system. The time domain signal of every code sequence is collected by the data acquisition card （DAQ）. A shift-averaging technique is applied in the frequency domain for the reason that the local oscillator （LO） in the coherent detection is fix-frequency deviated from the primary source. With the 31-bit simplex code, the signal-to-noise ratio （SNR） has 3.5-dB enhancement with the same single pulse traces, accordant with the theoretical analysis. The frequency fluctuation for simplex codes is 14.01 MHz less than that for a single pulse as to 4-m spatial resolution. The results are believed to be beneficial for the BOTDR performance improvement.

Influence of laser linewidth on performance of Brillouin optical time domain reflectometry

The effects of optical sources with different laser linewidths on Brillouin optical time domain reflectometry （BOTDR） are investigated numerically and experimentally. Simulation results show that the spectral linewidth of spontaneous Brillouin scattering remains almost constant when the laser linewidth is less than 1 MHz at the same pulse width; otherwise, it increases sharply. A comparison between a fiber laser （FL） with 4-kHz linewidth at 3 dB and a distributed feedback （DFB） laser with 3-MHz linewidth is made experimentally. When a constant laser power is launched into the sensing fiber, the fitting linewidths of the beat signals （backscattered Brillouin light and local oscillator （LO）） is about 5 MHz wider for the DFB laser than for the FL and the intensity of the beat signal is about a half. Furthermore, the frequency fluctuation in the long sensing fiber is lower for the FL source, yielding about 2 MHz less than that of the DFB laser, indicating higher temperature/strain resolution. The experimental results are in good agreement with the numerical simulations.

Time Domain Simulation of Transient Responses of Very Large Floating Structures Under Unsteady External Loads

A time domain finite element method （FEM） for the analysis of transient elastic response of a very large floating structure （VLFS） subjected to arbitrary time-dependent external loads is presented. This method is developed directly in time domain and the hydrodynamic problem is formulated based on linear, inviscid and slightly compressible fluid theory and the structural response is analyzed on the thin plate assumption. The time domain finite element procedure herein is validated by comparing numerical results with available experimental data. Finally, the transient elastic response of a pontoon-type VLFS under the landing of an airplane is computed by the proposed time domain FEM. The time histories of the applied force and the position and velocity of an airplane during landing are modeled with data from a Boeing 747-400 jumbo jet.

Dynamic Load Identification in Time Domain

In the time domain method the dynamic load is successfully identified when the accelerations, velocities and displacements or velocities and displacements of the structure are known. But in engineering practice or experiments usually only the acceleration response is recorded. In this paper an improved time domain method is proposed for dynamic load identification. In this method by using of Duhamel integral, only the acceleration response is required for load identification. As an application of the present method, the dynamic ice load on a Bohai offshore platform is identified based on some measured acceleration. The identified values of ice load are in good agreement with the measured ones.

Domain Decomposition and Matching for Time-Domain Analysis of Motions of Ships Advancing in Head Sea

A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare.