The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag mod...The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag model and homogenization technique to compute the bending wavenumbers via a 6th-order equation.The General Laminate Model(GLM)employs Mindlin’s displacement field to establish displacement-constraint relationships and resolves a quadratic Eigenvalue Problem(EVP)of the dispersion relation.The Wave Finite Element(WFE)scheme formulates the Nonlinear Eigenvalue Problem(NEP)for waves in varying directions and tracks complex wavenumbers using Weighted Wave Assurance Criteria(WWAC).Two approaches are introduced to estimate the Damping Loss Factor(DLF)of HDS,with the average DLF calculated by the modal density at various angles where non-homogeneity is present.Evaluation of robustness and accuracy is made by comparing the wavenumbers and DLF obtained from AHM and GLM with WFE.WFE is finally extended to a sandwich metastructure with a non-homogeneous core,and the Power Input Method(PIM)with Finite Element Method(FEM)data is employed to assess the average DLF,demonstrating an enhanced DLF compared to layered configurations with the same material portion,indicating increased energy dissipation due to the bending-shear coupling effects.展开更多
As a promising numerical tool of structural dynamics in mid- and high frequencies, the wave and finite element method（WFEM） is receiving increasingly attention and applications. In this paper, an enhanced WFEM has b...As a promising numerical tool of structural dynamics in mid- and high frequencies, the wave and finite element method（WFEM） is receiving increasingly attention and applications. In this paper, an enhanced WFEM has been developed with a reduced model and a new eigenvalue scheme. The reduced model is applicable for structures with piezoelectric shunts or local dampers;the new eigenvalue scheme can mitigate the ill-conditioning when the wave basis is calculated. The enhanced WFEM is applied to a thin-wall structure with periodically distributed piezoelectric materials（PZT）. Both free wave characteristics and forced response are analyzed and the influences of the suggested enhancements are presented. It is shown that if the control factors are properly chosen, these enhancements can improve the accuracy while accelerating the calculation. Resulting from the complexity of the application, these enhancements are not optional but imperative.展开更多
A superposing principle, by suitably adding the strain waves from a number of concentrated explosive charges to approximate the waves generated by a cylindrical charge based on the strain wave of a point or small sphe...A superposing principle, by suitably adding the strain waves from a number of concentrated explosive charges to approximate the waves generated by a cylindrical charge based on the strain wave of a point or small spherical explosive charge generated in rock, is used to further study the triggering time of strain gauges installed in radial direction at same distances but different positions surrounding a cylindrical explosive charge in rock. The duration of the first compression phase and peak value of strain wave, and furthermore, their differences are analyzed and some explanations are given. Besides that, the gauge orientation in which the maximum peak value occurs is also discussed. At last, the effect of velocity of detonation(V.O.D.) of a cylindrical explosive charge on the strain waves generated in the surrounding rock is taken as key research and the pattern of peak amplitude of a strain wave varies with the V.O.D. is likely to have been found.展开更多
A numerical wave flume with fully nonlinear free surface boundary conditions is adopted to investigate the temporal characteristics of extreme waves in the presence of wind at various speeds. Incident wave trains are ...A numerical wave flume with fully nonlinear free surface boundary conditions is adopted to investigate the temporal characteristics of extreme waves in the presence of wind at various speeds. Incident wave trains are numerically generated by a piston-type wave maker, and the wind-excited pressure is introduced into dynamic boundary conditions using a pressure distribution over steep crests, as defined by Jeffreys＇ sheltering mechanism.A boundary value problem is solved by a higher-order boundary element method（HOBEM） and a mixed Eulerian-Lagrangian time marching scheme. The proposed model is validated through comparison with published experimental data from a focused wave group. The influence of wind on extreme wave properties,including maximum extreme wave crest, focal position shift, and spectrum evolution, is also studied. To consider the effects of the wind-driven currents on a wave evolution, the simulations assume a uniform current over varying water depth. The results show that wind causes weak increases in the extreme wave crest, and makes the nonlinear energy transfer non-reversible in the focusing and defocusing processes. The numerical results also provide a comparison to demonstrate the shifts at focal points, considering the combined effects of the winds and the wind-driven currents.展开更多
A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressur...A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressure development of the intercalated soft layer under wave cyclic loading. By combining the undrained shear strength with the empirical formula of overconsolidation clay produced by unloading and the development model of pore pressure, the dynamic degradation law that describes the undrained shear strength as a function of cycle number and stress level is derived. Based on the proposed dynamic degradation law and M-C yield criterion, a dynamic finite element method is numerically implemented to predict changes in undrained shear strength of the intercalated soft layer by using the general-purpose FEM software ABAQUS, and the accuracy of the method is verified. The effects of cycle number and amplitude of the wave force on the degradation of the undrained shear strength of the intercalated soft layer and the associated excess pore pressure response are investigated by analyzing an overall distribution and three typical sections underneath the breakwater. By comparing the undrained shear strength distributions obtained by the static method and the quasi-static method with the undrained shear strength distributions obtained by the dynamic finite element method in the three typical sections, the superiority of the dynamic finite element method in predicting changes in undrained shear strength is demonstrated.展开更多
The paper presents Voronoi based micro-structure modeling through elastic wave propagation in a poly-crystalline copper using finite element method. The micro-structural parameters studied here are; the grain size and...The paper presents Voronoi based micro-structure modeling through elastic wave propagation in a poly-crystalline copper using finite element method. The micro-structural parameters studied here are; the grain size and the grain orientation. The poly-crystalline copper is modeled as a randomly oriented Voronoi cells in a fixed 2D computational domain. Tone burst 3-cycle pulse of 1 MHz frequency is used as the line source or point source for testing. Welded contact conditions are used at the interface boundaries of any two mutual cells of the domain. It is reported that wave scattering independent of the shape when the size of the scatterer less than the wavelength. Also, It is concluded that transmission efficiency increases as the cell size decreases.展开更多
In this article, the analytical homogenization method of periodic discrete media（HPDM）and the numerical condensed wave finite element method（CWFEM） are employed to study the longitudinal and transverse vibrations ...In this article, the analytical homogenization method of periodic discrete media（HPDM）and the numerical condensed wave finite element method（CWFEM） are employed to study the longitudinal and transverse vibrations of framed structures. The valid frequency range of the HPDM is re-evaluated using the wave propagation feature identified by the CWFEM. The relative error of the wavenumber by the HPDM compared to that by the CWFEM is illustrated in functions of frequency and scale ratio. A parametric study on the thickness of the structure is carried out where the dispersion relation and the relative error are given for three different thicknesses. The dynamics of a finite structure such as natural frequency and forced response are also investigated using the HPDM and the CWFEM.展开更多
To investigate the dynamics of submersible mussel rafts, the finite element program Aqua-FETM, developed by the University of New Hampshire（UNH）, was applied to rafts moored at the surface and submerged. The submerg...To investigate the dynamics of submersible mussel rafts, the finite element program Aqua-FETM, developed by the University of New Hampshire（UNH）, was applied to rafts moored at the surface and submerged. The submerged configuration is used to reduce wave forcing and to avoid contact with floating ice during winters in northern waters. Each raft consists of three pontoons connected by a grid framework. Rafts are intended to support densely spaced mussel ropes hung from the framework. When submerged, the pontoons are flooded, and the raft is held vertically by floats attached by lines. The computer models were developed in Aqua-FE？ to simulate the effects of waves and current. They were validated by comparison with wave tank results by use of a 1/10 scale raft physical model. Comparisons showed good agreement for the important heave（vertical） and pitch（rotational） motions, though there was a tendency towards conservative results for wave and current drag. Full-scale simulations of surface and submerged single raft and two rafts connected in tandem were performed. Submerged raft wave response was found to be reduced relative to that at the surface for both the single and two-raft configurations. In particular, the vertical motion of mussel rope connection points was significantly reduced by submergence, resulting in reduced potential for mussel drop-off. For example, the maximum vertical velocities of mussel rope attachment points in the submerged two raft case were 7%-20% of the corresponding velocities when at the surface.展开更多
The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat sour...The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat source. Finite element method is applied to analyzed the temperature field in the solder ball, and experi- ment is conducted to test the vibration. the results show that,that, the temperature field flucturates with the same frequency as that of the laser pulse, which in turn causes a forced vibration of the same frequency in the liquid solder ball.展开更多
The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect ...The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.展开更多
If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restric...If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.展开更多
A nonlinear reaction-diffusion equation is studied numerically by a Petrov-Galerkin finite element method, which has been proved to be 2nd-order accurate in time and 4th-order in space. The comparison between the exac...A nonlinear reaction-diffusion equation is studied numerically by a Petrov-Galerkin finite element method, which has been proved to be 2nd-order accurate in time and 4th-order in space. The comparison between the exact and numerical solutions of progressive waves shows that this numerical scheme is quite accurate, stable andefflcient. It is also shown that any local disturbance will spread, have a full growth and finally form two progressive waves propagating in both directions. The shape and the speed of the long term progressive waves are determined by the system itself, and do not depend on the details of the initial values.展开更多
Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground mot...Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from InSAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the Ms7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were con- structed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault, rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.展开更多
The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
This paper presents a novel parallel implementation technology for wave-based structural health monitoring (SHM) in laminated composite plates. The wavelet-based B-spline wavelet on he interval (BSWI) element is cons...This paper presents a novel parallel implementation technology for wave-based structural health monitoring (SHM) in laminated composite plates. The wavelet-based B-spline wavelet on he interval (BSWI) element is constructed according to Hamilton’s principle, and the element by element algorithm is parallelly executed on graphics processing unit (GPU) using compute unified device architecture (CUDA) to get the responses in full wave field accurately. By means of the Fourier spectral analysis method,the Mindlin plate theory is selected for wave modeling of laminated composite plates while the Kirchhoff plate theory predicts unreasonably phase and group velocities. Numerical examples involving wave propagation in laminated composite plates without and with crack are performed and discussed in detail. The parallel implementation on GPU is accelerated 146 times comparing with the same wave motion problem executed on central processing unit (CPU). The validity and accuracy of the proposed parallel implementation are also demonstrated by comparing with conventional finite element method (FEM) and the computation time has been reduced from hours to minutes. The damage size and location have been successfully determined according to wave propagation results based on delay-and-sum (DAS). The results show that the proposed parallel implementation of wavelet finite element method (WFEM) is very appropriate and efficient for wave-based SHM in laminated composite plates.展开更多
A boundary element method(BEM) is presented to compute the transmission spectra of two-dimensional(2-D) phononic crystals of a square lattice which are finite along the x-direction and infinite along the y-direction.T...A boundary element method(BEM) is presented to compute the transmission spectra of two-dimensional(2-D) phononic crystals of a square lattice which are finite along the x-direction and infinite along the y-direction.The cross sections of the scatterers may be circular or square.For a periodic cell,the boundary integral equations of the matrix and the scatterers are formulated.Substituting the periodic boundary conditions and the interface continuity conditions,a linear equation set is formed,from which the elastic wave transmission can be obtained.From the transmission spectra,the band gaps can be identified,which are compared with the band structures of the corresponding infinite systems.It is shown that generally the transmission spectra completely correspond to the band structures.In addition,the accuracy and the efficiency of the boundary element method are analyzed and discussed.展开更多
The boundary element method（BEM） is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required...The boundary element method（BEM） is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required calculations and storage increase rapidly with the increase of the structure scale. Thus, an accelerated method with a low storage is desirable for the wave interaction with a very large structure. A systematic review is given in this paper for the BEM for solving the problem of the wave interaction with a large scale structure. Various integral equations are derived based on different Green functions, the advantages and disadvantages of different discretization schemes of the integral equations by the constant panels, the higher order elements, and the spline functions are discussed. For the higher order element discretization method, the special concerns are given to the numerical calculations of the single-layer potential, the double layer potential and the solid angle coefficients. For a large scale computation problem such as the wave interaction with a very large structure or a large number of bodies, the BEMs with the FMM and p FFT accelerations are discussed, respectively, including the principles of the FMM and the p FFT, and their implementations in various integral equations with different Green functions. Finally, some potential applications of the acceleration methods for problems with large scale computations in the ocean and coastal engineering are introduced.展开更多
In this paper, we first develop the far field asymptotic solutions of the second-order scattering waves for the vertical plane problem taking the second-order Stokes waves as the incident waves. The asymptotic solutio...In this paper, we first develop the far field asymptotic solutions of the second-order scattering waves for the vertical plane problem taking the second-order Stokes waves as the incident waves. The asymptotic solutions satisfy the Laplace equation, the sea bed and free surface boundary conditions and are the out-going waves. Then the radiation conditions of the second-order mattering waves are derived by using the asymptotic solutions. By using the two-dimensinal finite clement method with the radiation conditions imposed on the ar- tificial boundaries, the computer program, known as 'NWF2', for determining nonlinear wave forces on large submerged bodies has been written. As a numerical example, nonlinear wave forces on a semi-circu- lar cylinder lying on the sea bed arc presented.展开更多
基金supported by the Natural Sciences and Engineering Research Council of Canada-Discovery Grant(individual)Program(No.NSEC-DG#355433-2009)funded by the LabEx CeLyA(Centre Lyonnais d’Acoustique,No.ANR-10-LABX-0060)of Universite?de Lyon。
文摘The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag model and homogenization technique to compute the bending wavenumbers via a 6th-order equation.The General Laminate Model(GLM)employs Mindlin’s displacement field to establish displacement-constraint relationships and resolves a quadratic Eigenvalue Problem(EVP)of the dispersion relation.The Wave Finite Element(WFE)scheme formulates the Nonlinear Eigenvalue Problem(NEP)for waves in varying directions and tracks complex wavenumbers using Weighted Wave Assurance Criteria(WWAC).Two approaches are introduced to estimate the Damping Loss Factor(DLF)of HDS,with the average DLF calculated by the modal density at various angles where non-homogeneity is present.Evaluation of robustness and accuracy is made by comparing the wavenumbers and DLF obtained from AHM and GLM with WFE.WFE is finally extended to a sandwich metastructure with a non-homogeneous core,and the Power Input Method(PIM)with Finite Element Method(FEM)data is employed to assess the average DLF,demonstrating an enhanced DLF compared to layered configurations with the same material portion,indicating increased energy dissipation due to the bending-shear coupling effects.
基金the company PSA Peugeot Citroёn for the financial support
文摘As a promising numerical tool of structural dynamics in mid- and high frequencies, the wave and finite element method（WFEM） is receiving increasingly attention and applications. In this paper, an enhanced WFEM has been developed with a reduced model and a new eigenvalue scheme. The reduced model is applicable for structures with piezoelectric shunts or local dampers;the new eigenvalue scheme can mitigate the ill-conditioning when the wave basis is calculated. The enhanced WFEM is applied to a thin-wall structure with periodically distributed piezoelectric materials（PZT）. Both free wave characteristics and forced response are analyzed and the influences of the suggested enhancements are presented. It is shown that if the control factors are properly chosen, these enhancements can improve the accuracy while accelerating the calculation. Resulting from the complexity of the application, these enhancements are not optional but imperative.
基金Projects(51304239,51374243)supported by the National Natural Science Foundation of China
文摘A superposing principle, by suitably adding the strain waves from a number of concentrated explosive charges to approximate the waves generated by a cylindrical charge based on the strain wave of a point or small spherical explosive charge generated in rock, is used to further study the triggering time of strain gauges installed in radial direction at same distances but different positions surrounding a cylindrical explosive charge in rock. The duration of the first compression phase and peak value of strain wave, and furthermore, their differences are analyzed and some explanations are given. Besides that, the gauge orientation in which the maximum peak value occurs is also discussed. At last, the effect of velocity of detonation(V.O.D.) of a cylindrical explosive charge on the strain waves generated in the surrounding rock is taken as key research and the pattern of peak amplitude of a strain wave varies with the V.O.D. is likely to have been found.
基金The National Natural Science Foundation of China under contract Nos 51679036,51490672 and 51709038the Fundamental Research Funds for the Central Universities under contract Nos DUT17GJ202 and DUT16RC(3)113the Open Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering under contract No.2016490111
文摘A numerical wave flume with fully nonlinear free surface boundary conditions is adopted to investigate the temporal characteristics of extreme waves in the presence of wind at various speeds. Incident wave trains are numerically generated by a piston-type wave maker, and the wind-excited pressure is introduced into dynamic boundary conditions using a pressure distribution over steep crests, as defined by Jeffreys＇ sheltering mechanism.A boundary value problem is solved by a higher-order boundary element method（HOBEM） and a mixed Eulerian-Lagrangian time marching scheme. The proposed model is validated through comparison with published experimental data from a focused wave group. The influence of wind on extreme wave properties,including maximum extreme wave crest, focal position shift, and spectrum evolution, is also studied. To consider the effects of the wind-driven currents on a wave evolution, the simulations assume a uniform current over varying water depth. The results show that wind causes weak increases in the extreme wave crest, and makes the nonlinear energy transfer non-reversible in the focusing and defocusing processes. The numerical results also provide a comparison to demonstrate the shifts at focal points, considering the combined effects of the winds and the wind-driven currents.
基金financially supported by the National Natural Science Foundation of China(Grant No.51279128)the National Natural Science Fund for Innovative Research Groups Science Foundation(Grant No.51321065)the Construction Science and Technology Project of Ministry of Transport of the People’s Republic of China(Grant No.2013328224070)
文摘A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressure development of the intercalated soft layer under wave cyclic loading. By combining the undrained shear strength with the empirical formula of overconsolidation clay produced by unloading and the development model of pore pressure, the dynamic degradation law that describes the undrained shear strength as a function of cycle number and stress level is derived. Based on the proposed dynamic degradation law and M-C yield criterion, a dynamic finite element method is numerically implemented to predict changes in undrained shear strength of the intercalated soft layer by using the general-purpose FEM software ABAQUS, and the accuracy of the method is verified. The effects of cycle number and amplitude of the wave force on the degradation of the undrained shear strength of the intercalated soft layer and the associated excess pore pressure response are investigated by analyzing an overall distribution and three typical sections underneath the breakwater. By comparing the undrained shear strength distributions obtained by the static method and the quasi-static method with the undrained shear strength distributions obtained by the dynamic finite element method in the three typical sections, the superiority of the dynamic finite element method in predicting changes in undrained shear strength is demonstrated.
基金BRNS project No.MEE111228BRNSKRIS for their financial support
文摘The paper presents Voronoi based micro-structure modeling through elastic wave propagation in a poly-crystalline copper using finite element method. The micro-structural parameters studied here are; the grain size and the grain orientation. The poly-crystalline copper is modeled as a randomly oriented Voronoi cells in a fixed 2D computational domain. Tone burst 3-cycle pulse of 1 MHz frequency is used as the line source or point source for testing. Welded contact conditions are used at the interface boundaries of any two mutual cells of the domain. It is reported that wave scattering independent of the shape when the size of the scatterer less than the wavelength. Also, It is concluded that transmission efficiency increases as the cell size decreases.
文摘In this article, the analytical homogenization method of periodic discrete media（HPDM）and the numerical condensed wave finite element method（CWFEM） are employed to study the longitudinal and transverse vibrations of framed structures. The valid frequency range of the HPDM is re-evaluated using the wave propagation feature identified by the CWFEM. The relative error of the wavenumber by the HPDM compared to that by the CWFEM is illustrated in functions of frequency and scale ratio. A parametric study on the thickness of the structure is carried out where the dispersion relation and the relative error are given for three different thicknesses. The dynamics of a finite structure such as natural frequency and forced response are also investigated using the HPDM and the CWFEM.
基金financially supported by the Small Business Innovation Research(SBIR)program of the USDA National Institute for Food and Agriculture(NIFA)(Grant No.2013-33610-21190)to Pemaquid Mussel FarmsDuring her time at the University of New Hampshire where this study was completedsupported by a graduate student fellowship funded by the People’s Republic of China
文摘To investigate the dynamics of submersible mussel rafts, the finite element program Aqua-FETM, developed by the University of New Hampshire（UNH）, was applied to rafts moored at the surface and submerged. The submerged configuration is used to reduce wave forcing and to avoid contact with floating ice during winters in northern waters. Each raft consists of three pontoons connected by a grid framework. Rafts are intended to support densely spaced mussel ropes hung from the framework. When submerged, the pontoons are flooded, and the raft is held vertically by floats attached by lines. The computer models were developed in Aqua-FE？ to simulate the effects of waves and current. They were validated by comparison with wave tank results by use of a 1/10 scale raft physical model. Comparisons showed good agreement for the important heave（vertical） and pitch（rotational） motions, though there was a tendency towards conservative results for wave and current drag. Full-scale simulations of surface and submerged single raft and two rafts connected in tandem were performed. Submerged raft wave response was found to be reduced relative to that at the surface for both the single and two-raft configurations. In particular, the vertical motion of mussel rope connection points was significantly reduced by submergence, resulting in reduced potential for mussel drop-off. For example, the maximum vertical velocities of mussel rope attachment points in the submerged two raft case were 7%-20% of the corresponding velocities when at the surface.
文摘The research presented here is focused on the vibration condition of a small volume solder solder ball,which is placed on the surface of a soldering pad and is exerted a pulse modulated continuous wave laser heat source. Finite element method is applied to analyzed the temperature field in the solder ball, and experi- ment is conducted to test the vibration. the results show that,that, the temperature field flucturates with the same frequency as that of the laser pulse, which in turn causes a forced vibration of the same frequency in the liquid solder ball.
文摘The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.
基金National Natural Science Foundation of China (50178065), 973 Program (2002CB412706), National Social Com-monweal Research Foundation (2002DIB30076) and Joint Seismological Science Foundation (101066).
文摘If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.
文摘A nonlinear reaction-diffusion equation is studied numerically by a Petrov-Galerkin finite element method, which has been proved to be 2nd-order accurate in time and 4th-order in space. The comparison between the exact and numerical solutions of progressive waves shows that this numerical scheme is quite accurate, stable andefflcient. It is also shown that any local disturbance will spread, have a full growth and finally form two progressive waves propagating in both directions. The shape and the speed of the long term progressive waves are determined by the system itself, and do not depend on the details of the initial values.
基金supported by the Earthquake Public Welfare Scientific Research Special Project (No.201408014)
文摘Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from InSAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the Ms7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were con- structed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault, rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51421004 & 51405369)the National Key Basic Research Program of China (Grant No. 2015CB057400)+1 种基金the China Postdoctoral Science Foundation (Grant No. 2014M560766)the China Scholarship Council,and the Fundamental Research Funds for the Central Universities(Grant No. xjj2014107)
文摘This paper presents a novel parallel implementation technology for wave-based structural health monitoring (SHM) in laminated composite plates. The wavelet-based B-spline wavelet on he interval (BSWI) element is constructed according to Hamilton’s principle, and the element by element algorithm is parallelly executed on graphics processing unit (GPU) using compute unified device architecture (CUDA) to get the responses in full wave field accurately. By means of the Fourier spectral analysis method,the Mindlin plate theory is selected for wave modeling of laminated composite plates while the Kirchhoff plate theory predicts unreasonably phase and group velocities. Numerical examples involving wave propagation in laminated composite plates without and with crack are performed and discussed in detail. The parallel implementation on GPU is accelerated 146 times comparing with the same wave motion problem executed on central processing unit (CPU). The validity and accuracy of the proposed parallel implementation are also demonstrated by comparing with conventional finite element method (FEM) and the computation time has been reduced from hours to minutes. The damage size and location have been successfully determined according to wave propagation results based on delay-and-sum (DAS). The results show that the proposed parallel implementation of wavelet finite element method (WFEM) is very appropriate and efficient for wave-based SHM in laminated composite plates.
基金supported by the National Natural Science Foundation of China(Grant Nos.11202021,11472249 and 51178037)the Beijing Natural Science Foundation(Grant No.1163008)the Postdoctoral Science Foundation of China(Grant No.2012M510311)
文摘A boundary element method(BEM) is presented to compute the transmission spectra of two-dimensional(2-D) phononic crystals of a square lattice which are finite along the x-direction and infinite along the y-direction.The cross sections of the scatterers may be circular or square.For a periodic cell,the boundary integral equations of the matrix and the scatterers are formulated.Substituting the periodic boundary conditions and the interface continuity conditions,a linear equation set is formed,from which the elastic wave transmission can be obtained.From the transmission spectra,the band gaps can be identified,which are compared with the band structures of the corresponding infinite systems.It is shown that generally the transmission spectra completely correspond to the band structures.In addition,the accuracy and the efficiency of the boundary element method are analyzed and discussed.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51379032,51490672 and 51479026)
文摘The boundary element method（BEM） is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required calculations and storage increase rapidly with the increase of the structure scale. Thus, an accelerated method with a low storage is desirable for the wave interaction with a very large structure. A systematic review is given in this paper for the BEM for solving the problem of the wave interaction with a large scale structure. Various integral equations are derived based on different Green functions, the advantages and disadvantages of different discretization schemes of the integral equations by the constant panels, the higher order elements, and the spline functions are discussed. For the higher order element discretization method, the special concerns are given to the numerical calculations of the single-layer potential, the double layer potential and the solid angle coefficients. For a large scale computation problem such as the wave interaction with a very large structure or a large number of bodies, the BEMs with the FMM and p FFT accelerations are discussed, respectively, including the principles of the FMM and the p FFT, and their implementations in various integral equations with different Green functions. Finally, some potential applications of the acceleration methods for problems with large scale computations in the ocean and coastal engineering are introduced.
文摘In this paper, we first develop the far field asymptotic solutions of the second-order scattering waves for the vertical plane problem taking the second-order Stokes waves as the incident waves. The asymptotic solutions satisfy the Laplace equation, the sea bed and free surface boundary conditions and are the out-going waves. Then the radiation conditions of the second-order mattering waves are derived by using the asymptotic solutions. By using the two-dimensinal finite clement method with the radiation conditions imposed on the ar- tificial boundaries, the computer program, known as 'NWF2', for determining nonlinear wave forces on large submerged bodies has been written. As a numerical example, nonlinear wave forces on a semi-circu- lar cylinder lying on the sea bed arc presented.
