This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of...This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of structural collapse.The ET method involves dynamic analysis of a structure under an intensifying record over time.While conventional ET methods can determine the median of the structural response,they lack the ability to calculate its dispersion.To address this limitation,the present study utilizes ET analysis and single-degree-of-freedom(SDOF)systems to develop a method that considers the record-to-record variability for calculating the probability distribution of structural response.The accuracy of this method is evaluated by comparing it with the incremental dynamic analysis(IDA)method using special moment frames.The results demonstrate that the proposed method achieves a reasonably accurate estimation of dispersion while significantly reducing the computational burden(by approximately 95%)compared to the IDA method.展开更多
The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidanc...The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidance system.This paper presents a new method that uses the Finite Difference Time Domain(FDTD)method to calculate the electrical performance of radomes under Thermo-Mechanical-Electrical(TME)coupling.This method can accurately characterize the effects of material dielectric temperature drift and structural deformation on the electrical performance of the radome under flight conditions,enabling high-precision full-wave calculations of the broadband electrical performance of the radome.The method initiates by utilizing a Finite Element Grid Model(FE-GM)of the radome to sequentially acquire the radome's response temperature field and structural deformation field through thermal and mechanical simulations.Subsequently,spatial mapping techniques are developed to accurately incorporate the dielectric temperature drift and structural deformation of the radome into its Yee grid Electromagnetic(EM)simulation model.A verification case was designed to test the proposed method,and the results confirmed its high computational accuracy.Additionally,the effectiveness and necessity of the method were further demonstrated by analyzing the electrical performance of a fused silica ceramic radome used on a high-speed aircraft.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions ...To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions at normal atmospheric temperature. The contents of N, P, K, and organic matter in biogas slurry of straw were determined in different storage times. The results showed that： during the pro-cess of biogas slurry storage, little change occurred in the content of the organic matter while the total content of N, P, K significantly declined; up to 50 days, the total content of N, P, K reduced to nearly 80%-90%. Because the contents of N, P, K in biogas slurry reduced less in airtight storage conditions so that a better re-sult was found on airtight storage methods than open storage methods in fertilizer field of biogas slurry of straw.展开更多
Liquid storage tanks are essential structures that are often located in residential and industrial areas; thus an assessment of their seismic performance is an important engineering issue. In this paper, the seismic r...Liquid storage tanks are essential structures that are often located in residential and industrial areas; thus an assessment of their seismic performance is an important engineering issue. In this paper, the seismic response ofunanchored steel liquid storage tanks is investigated using the endurance time （ET） dynamic analysis procedure and compared to responses obtained for anchored tanks under actual ground motions and intensifying ET records. In most cases, the results from ground motions are properly obtained with negligible differences using ET records. It is observed that uplifting of the tank base, which is closely related to the tank aspect ratio, has the greatest significance in the responses of the tank and can be predicted with reasonable accuracy by using currently available ET records.展开更多
In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum,...In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum, intensifi es with time so they can be used approximately to simulate the average effects of several ground motions scaled to different intensities. In this paper applicability of the ET method for evaluating collapse potential of buildings is investigated. A set of four steel moment frames is used for collapse assessment. The process of using ET method in collapse evaluation is explained and the results are compared with incremental dynamic analysis(IDA) results. It is shown that although the computational effort using the ET method is much less than the IDA analysis, the results of both methods are consistent. Finally collapse fragility curves using ET and IDA methods are produced and it is shown that the probabilities of collapse in different hazard levels are also consistent.展开更多
The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and t...The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and the pH of the solutions were also investigated. A timed complexatian colorimetric analysis method of Al-Fe-ferron in view of the total concentration of {AI + Fe} was then established to determine the species distribution of polymeric Al-Fe. The testing wavelength was recommended at 362 net and the testing pH value was 5. With a comparison of the ratios of n(Al)/n(Fe), the standard adsorption curves of the polymeric Al-Fe solutions were derived from the experimental results. Furthermore, the solutions' composition were carious in both the molar n(Al)/n(Fe) ratios, i.e. 0/0, 5/5, 9/1 and 0/10, and the concentrations associated with the total ( Al + Fe which ranged from 10(-5) to 10(-4) mol/L..展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to ana...In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to stud...In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.展开更多
The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-mov...The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-moving algorithm is introduced to ensure the clouds of points following the movements of bodyboundaries. A dual time method for solving the two-dimenslonal Euler equations in Arbitrary Lagrangian-Eulerian （ALE） formulation is presented. Dual time method allows the real-time step to be chosen on the basis of accuracy rather than stability. It also permits the acceleration techniques, which are commonly used to speed up steady flow calculations, to be used when marching the equations in pseudo time. The spatial derivatives, which are used to estimating the inviscid flux, are directly approximated by using local least-squares curve method. An explicit multistage Runge-Kutta algorithm is used to advance the flow equations in pseudo time. In order to accelerate the solution to convergence, local time stepping technique and residual averaging are employed. The results of NACA0012 airfoil in transonic steady flow are presented to verify the accuracy of the present spatial discretization method. Finally, two AGARD standard test cases in which NACA0012 airfoil and NACA64A010 airfoil oscillate in transonic flow are simulated. The computational results are compared with the experimental data to demonstrate the validity and practicality of the presented method.展开更多
In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments...In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments of Nino3 SST anomalies and Tahiti-Darwin SO index. The results show that the scheme is feasible and ENSO predictable.展开更多
Al Ferron timed complex colorimetric method (AFM) and 27 Al NMR spectroscopy method(ANM) were discussed. For the former, the different colorimetric reagent preparation methods' results indicate that the...Al Ferron timed complex colorimetric method (AFM) and 27 Al NMR spectroscopy method(ANM) were discussed. For the former, the different colorimetric reagent preparation methods' results indicate that there are some differences beteween them, and the combined method can be used as a simplified procedure. For the latter, the small tube method is more accurate. Eventually, the Al 13 (ANM) was compared to the Al b (AFM).展开更多
Probabilistic Seismic Loss Estimation is a methodology used as a quantitative and explicit expression of the performance of buildings using terms that address the interests of both owners and insurance companies. Appl...Probabilistic Seismic Loss Estimation is a methodology used as a quantitative and explicit expression of the performance of buildings using terms that address the interests of both owners and insurance companies. Applying the ATC 58 approach for seismic loss assessment of buildings requires using Incremental Dynamic Analysis （IDA）, which needs hundreds of time-consuming analyses, which in turn hinders its wide application. The Endurance Time Method （ETM） is proposed herein as part of a demand propagation prediction procedure and is shown to be an economical alternative to IDA. Various scenarios were considered to achieve this purpose and their appropriateness has been evaluated using statistical methods. The most precise and efficient scenario was validated through comparison against IDA driven response predictions of 34 code conforming benchmark structures and was proven to be sufficiently precise while offering a great deal of efficiency. The loss values were estimated by replacing IDA with the proposed ETM-based procedure in the ATC 58 procedure and it was fotmd that these values suffer from varying inaccuracies, which were attributed to the discretized nature of damage and loss prediction functions provided by ATC 58.展开更多
The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to ...The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.展开更多
A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme...A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.展开更多
In this work, we will derive numerical schemes for solving 3-coupled nonlinear Schrödinger equations using finite difference method and time splitting method combined with finite difference method. The result...In this work, we will derive numerical schemes for solving 3-coupled nonlinear Schrödinger equations using finite difference method and time splitting method combined with finite difference method. The resulting schemes are highly accurate, unconditionally stable. We use the exact single soliton solution and the conserved quantities to check the accuracy and the efficiency of the proposed schemes. Also, we use these methods to study the interaction dynamics of two solitons. It is found that both elastic and inelastic collision can take place under suitable parametric conditions. We have noticed that the inelastic collision of single solitons occurs in two different manners: enhancement or suppression of the amplitude.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
The endurance time (ET) method is a time history based dynamic analysis in which structures are subjected to gradually intensifying excitations and their performances are judged based on their responses at various exc...The endurance time (ET) method is a time history based dynamic analysis in which structures are subjected to gradually intensifying excitations and their performances are judged based on their responses at various excitation levels. Using this method, the computational effort required for estimating probable seismic demand parameters can be reduced by an order of magnitude. Calculation of the maximum displacement or target displacement is a basic requirement for estimating performance based on structural design. The purpose of this paper is to compare the results of the nonlinear ET method with the nonlinear static pushover (NSP) method of FEMA 356 by evaluating performances and target displacements of steel frames. This study will lead to a deeper insight into the capabilities and limitations of the ET method. The results are further compared with those of the standard nonlinear response history analysis. We conclude that results from the ET analysis are in proper agreement with those from standard procedures.展开更多
文摘This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of structural collapse.The ET method involves dynamic analysis of a structure under an intensifying record over time.While conventional ET methods can determine the median of the structural response,they lack the ability to calculate its dispersion.To address this limitation,the present study utilizes ET analysis and single-degree-of-freedom(SDOF)systems to develop a method that considers the record-to-record variability for calculating the probability distribution of structural response.The accuracy of this method is evaluated by comparing it with the incremental dynamic analysis(IDA)method using special moment frames.The results demonstrate that the proposed method achieves a reasonably accurate estimation of dispersion while significantly reducing the computational burden(by approximately 95%)compared to the IDA method.
文摘The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidance system.This paper presents a new method that uses the Finite Difference Time Domain(FDTD)method to calculate the electrical performance of radomes under Thermo-Mechanical-Electrical(TME)coupling.This method can accurately characterize the effects of material dielectric temperature drift and structural deformation on the electrical performance of the radome under flight conditions,enabling high-precision full-wave calculations of the broadband electrical performance of the radome.The method initiates by utilizing a Finite Element Grid Model(FE-GM)of the radome to sequentially acquire the radome's response temperature field and structural deformation field through thermal and mechanical simulations.Subsequently,spatial mapping techniques are developed to accurately incorporate the dielectric temperature drift and structural deformation of the radome into its Yee grid Electromagnetic(EM)simulation model.A verification case was designed to test the proposed method,and the results confirmed its high computational accuracy.Additionally,the effectiveness and necessity of the method were further demonstrated by analyzing the electrical performance of a fused silica ceramic radome used on a high-speed aircraft.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
基金Supported by Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period(2012BAJ21B04-04)Shandong Environment Bottlenecks Analysis Project(SDHBPJ-ZB-05)~~
文摘To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions at normal atmospheric temperature. The contents of N, P, K, and organic matter in biogas slurry of straw were determined in different storage times. The results showed that： during the pro-cess of biogas slurry storage, little change occurred in the content of the organic matter while the total content of N, P, K significantly declined; up to 50 days, the total content of N, P, K reduced to nearly 80%-90%. Because the contents of N, P, K in biogas slurry reduced less in airtight storage conditions so that a better re-sult was found on airtight storage methods than open storage methods in fertilizer field of biogas slurry of straw.
文摘Liquid storage tanks are essential structures that are often located in residential and industrial areas; thus an assessment of their seismic performance is an important engineering issue. In this paper, the seismic response ofunanchored steel liquid storage tanks is investigated using the endurance time （ET） dynamic analysis procedure and compared to responses obtained for anchored tanks under actual ground motions and intensifying ET records. In most cases, the results from ground motions are properly obtained with negligible differences using ET records. It is observed that uplifting of the tank base, which is closely related to the tank aspect ratio, has the greatest significance in the responses of the tank and can be predicted with reasonable accuracy by using currently available ET records.
文摘In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum, intensifi es with time so they can be used approximately to simulate the average effects of several ground motions scaled to different intensities. In this paper applicability of the ET method for evaluating collapse potential of buildings is investigated. A set of four steel moment frames is used for collapse assessment. The process of using ET method in collapse evaluation is explained and the results are compared with incremental dynamic analysis(IDA) results. It is shown that although the computational effort using the ET method is much less than the IDA analysis, the results of both methods are consistent. Finally collapse fragility curves using ET and IDA methods are produced and it is shown that the probabilities of collapse in different hazard levels are also consistent.
基金TheNationalNaturalScienceFoundationofChina (No .2 96 770 0 4)
文摘The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and the pH of the solutions were also investigated. A timed complexatian colorimetric analysis method of Al-Fe-ferron in view of the total concentration of {AI + Fe} was then established to determine the species distribution of polymeric Al-Fe. The testing wavelength was recommended at 362 net and the testing pH value was 5. With a comparison of the ratios of n(Al)/n(Fe), the standard adsorption curves of the polymeric Al-Fe solutions were derived from the experimental results. Furthermore, the solutions' composition were carious in both the molar n(Al)/n(Fe) ratios, i.e. 0/0, 5/5, 9/1 and 0/10, and the concentrations associated with the total ( Al + Fe which ranged from 10(-5) to 10(-4) mol/L..
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
基金The project partly supported by the Foundation of Zhongshan University Advanced Research Center
文摘In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
基金supported by the Natural Science Foundation of China Government (10902051)the Natural Science Foundation of Jiangsu Province (BK2008046)the German Science Foundation
文摘In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.
文摘The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-moving algorithm is introduced to ensure the clouds of points following the movements of bodyboundaries. A dual time method for solving the two-dimenslonal Euler equations in Arbitrary Lagrangian-Eulerian （ALE） formulation is presented. Dual time method allows the real-time step to be chosen on the basis of accuracy rather than stability. It also permits the acceleration techniques, which are commonly used to speed up steady flow calculations, to be used when marching the equations in pseudo time. The spatial derivatives, which are used to estimating the inviscid flux, are directly approximated by using local least-squares curve method. An explicit multistage Runge-Kutta algorithm is used to advance the flow equations in pseudo time. In order to accelerate the solution to convergence, local time stepping technique and residual averaging are employed. The results of NACA0012 airfoil in transonic steady flow are presented to verify the accuracy of the present spatial discretization method. Finally, two AGARD standard test cases in which NACA0012 airfoil and NACA64A010 airfoil oscillate in transonic flow are simulated. The computational results are compared with the experimental data to demonstrate the validity and practicality of the presented method.
文摘In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments of Nino3 SST anomalies and Tahiti-Darwin SO index. The results show that the scheme is feasible and ENSO predictable.
文摘Al Ferron timed complex colorimetric method (AFM) and 27 Al NMR spectroscopy method(ANM) were discussed. For the former, the different colorimetric reagent preparation methods' results indicate that there are some differences beteween them, and the combined method can be used as a simplified procedure. For the latter, the small tube method is more accurate. Eventually, the Al 13 (ANM) was compared to the Al b (AFM).
文摘Probabilistic Seismic Loss Estimation is a methodology used as a quantitative and explicit expression of the performance of buildings using terms that address the interests of both owners and insurance companies. Applying the ATC 58 approach for seismic loss assessment of buildings requires using Incremental Dynamic Analysis （IDA）, which needs hundreds of time-consuming analyses, which in turn hinders its wide application. The Endurance Time Method （ETM） is proposed herein as part of a demand propagation prediction procedure and is shown to be an economical alternative to IDA. Various scenarios were considered to achieve this purpose and their appropriateness has been evaluated using statistical methods. The most precise and efficient scenario was validated through comparison against IDA driven response predictions of 34 code conforming benchmark structures and was proven to be sufficiently precise while offering a great deal of efficiency. The loss values were estimated by replacing IDA with the proposed ETM-based procedure in the ATC 58 procedure and it was fotmd that these values suffer from varying inaccuracies, which were attributed to the discretized nature of damage and loss prediction functions provided by ATC 58.
基金supported by National Key Research and Development Program of China(No.2021YFA0717100)NationalNatural Science Foundation of China(Nos.12072270,U2013206).
文摘The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.
文摘A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.
文摘In this work, we will derive numerical schemes for solving 3-coupled nonlinear Schrödinger equations using finite difference method and time splitting method combined with finite difference method. The resulting schemes are highly accurate, unconditionally stable. We use the exact single soliton solution and the conserved quantities to check the accuracy and the efficiency of the proposed schemes. Also, we use these methods to study the interaction dynamics of two solitons. It is found that both elastic and inelastic collision can take place under suitable parametric conditions. We have noticed that the inelastic collision of single solitons occurs in two different manners: enhancement or suppression of the amplitude.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
文摘The endurance time (ET) method is a time history based dynamic analysis in which structures are subjected to gradually intensifying excitations and their performances are judged based on their responses at various excitation levels. Using this method, the computational effort required for estimating probable seismic demand parameters can be reduced by an order of magnitude. Calculation of the maximum displacement or target displacement is a basic requirement for estimating performance based on structural design. The purpose of this paper is to compare the results of the nonlinear ET method with the nonlinear static pushover (NSP) method of FEMA 356 by evaluating performances and target displacements of steel frames. This study will lead to a deeper insight into the capabilities and limitations of the ET method. The results are further compared with those of the standard nonlinear response history analysis. We conclude that results from the ET analysis are in proper agreement with those from standard procedures.
