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.展开更多
Relationship between the radiation conditions at infinity and the transmitting boundary for numerical simulation of the near-field wave motion has been studied in this paper. The conclusion is that the transmitting bo...Relationship between the radiation conditions at infinity and the transmitting boundary for numerical simulation of the near-field wave motion has been studied in this paper. The conclusion is that the transmitting boundary is approximately equivalent to the radiation conditions at infinity for a large class of infinite media. And the errors of the approximation are of the same order of magnitude as those of the finite elements or finite differences in numerical simulation of wave motion. This result provides a sound theoretical basis for the transmitting boundary used in the numerical simulation of the near-field wave motion and gives a complete explanation for the major experiences accumulated in applications of the transmitting boundary to the numerical simulation.展开更多
基金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.
基金the Joint Seismological Science Foundation. Prof. Dai Yishan in Harbin Engineering University has read the manuscript and made relevant comments The expansions of 1/r and r using the Legendre polynomials and the proof of eq (29) on the properties of
文摘Relationship between the radiation conditions at infinity and the transmitting boundary for numerical simulation of the near-field wave motion has been studied in this paper. The conclusion is that the transmitting boundary is approximately equivalent to the radiation conditions at infinity for a large class of infinite media. And the errors of the approximation are of the same order of magnitude as those of the finite elements or finite differences in numerical simulation of wave motion. This result provides a sound theoretical basis for the transmitting boundary used in the numerical simulation of the near-field wave motion and gives a complete explanation for the major experiences accumulated in applications of the transmitting boundary to the numerical simulation.
