Historic monuments, drywall structures, and graphite blocks in AGR nuclear power plants are block-like structures that have to withstand rocking when subject to seismic excitation of their base, which can lead to over...
详细信息
Historic monuments, drywall structures, and graphite blocks in AGR nuclear power plants are block-like structures that have to withstand rocking when subject to seismic excitation of their base, which can lead to overturning of some of their components and results in the collapse of the whole structure. We revisit the known nonlinear equations of motion for a dual-block stack and present the conditions for transition between the eight possible rocking configurations (due to initiation of rocking, opening of new contacts, and collisions between blocks). An algorithm for the numerical simulation of rocking of the dual-block stack is developed using the Newmark integration method, the Newton-Raphson iteration method, and a novel contact detection and resolution procedure. The algorithm is used to evaluate rocking stability of five dual-block stacks, one of which is compared to the results available in the literature. In parallel, a novel experimental program is designed and implemented, to validate the numerically obtained results using a shaking table. While most of the excitation conditions leading to stable rocking and limit values leading to overturning have been successfully validated, some discrepancies between the numerically and experimentally obtained results still exist and point to the need for improvement of the algorithm used, possibly through a more realistic energy-loss mechanism. Most importantly, we have confirmed the known theoretical prediction that splitting a single block into two half-size blocks benefits its rocking stability.
This paper establishes a general computational framework to solve the muti-scale contact problem by integrating the statistical contact model with the finite element format. Compared to existing models, the proposed m...
详细信息
This paper establishes a general computational framework to solve the muti-scale contact problem by integrating the statistical contact model with the finite element format. Compared to existing models, the proposed method is applicable to most geometric configurations and can effectively evaluate the pressure distribution. In this work, a modified Karush-Kuhn-Tucker (KKT) condition is proposed by the assumption that asperity height obeys the Gaussian distribution. Therefore, in the variational formula, the contact contribution is decomposed into body contribution and asperity contribution, corresponding to the nominal smooth surface and roughness, respectively. Then the linearization and constraint enforcement of these two components are derived, followed by a nonlinear Newton-Raphson-based iterative algorithm. The contact patch test and Hertz contact test are conducted, and the predicted results are consistent with the theoretical and experimental values, confirming the effectiveness and accuracy of the proposed approach. It is worth noting that in the Hertz contact test, the contact pressure distribution varies progressively with the roughness level and external force, tending to the Hertz limit or Gaussian limit. This means that the proposed method can be applied to any roughness and load. Finally, the contact behaviors of the transmission interface of a piezoelectric actuator, i.e., a typical multi-scale contact problem, are studied as an engineering application case.
Future space exploration technology requires a long-life and reliable power source that is not reliant on solar energy. Space micro-reactors are able to meet this need, with Heat Pipe Cooled Reactors (HPR) emerging as...
详细信息
Future space exploration technology requires a long-life and reliable power source that is not reliant on solar energy. Space micro-reactors are able to meet this need, with Heat Pipe Cooled Reactors (HPR) emerging as a notable type of space micro-reactor that has attracted widespread attention in recent years. The HPR utilizes high-temperature alkali metal heat pipes for heat transfer, which presents certain complexities due to the solid state of the alkali metals working medium at room temperature. This results in a three-phase transition during the high-temperature heat pipes start-up process, which significantly impacts the heat transfer characteristics and dynamic behavior of the HPR start-up process. Consequently, thorough research is necessary in this area. numerical simulation is a crucial tool that can effectively analyze, predict, and guide experiments. This article utilizes the Finite Volume Method (FVM) to develop a simulation code for high-temperature heat pipe frozen start-up. Various physical models are integrated to describe different components of the heat pipe: the container is represented by a two-dimensional axisymmetric heat conduction equation, the wick region utilizes a Fixed Grid Method (FGM) to depict the melting process of the medium, and the vapor channel is described through a two-dimensional axisymmetric compressible laminar flow. The wick region and vapor channel are coupled through the evaporation and condensation of the medium. For the vapor channel, numerical methods such as SIMPLE and PISO are used for solving. Adaptive time step and OpenMP acceleration are employed in the code to enhance computational efficiency. Finally, by comparing the calculated results with experimental data, the feasibility and accuracy of the code are assessed, highlighting special phenomena during the start-up process. The findings confirm that the developed code accurately predicts parameter changes during start-up, and can serve as a heat pipe analysis mo
It is known that the generating function of a sequence of Toeplitz matrices may not describe the asymptotic distribution of the eigenvalues of the considered matrix sequence in the non-Hermitian setting. In a recent w...
详细信息
It is known that the generating function of a sequence of Toeplitz matrices may not describe the asymptotic distribution of the eigenvalues of the considered matrix sequence in the non-Hermitian setting. In a recent work, under the assumption that the eigenvalues are real, admitting an asymptotic expansion whose first term is the distribution function, fast algorithms computing all the spectra were proposed in different settings. In the current work, we extend this idea to non-Hermitian Toeplitz matrices with complex eigenvalues, in the case where the range of the generating function does not disconnect the complex field or the limiting set of the spectra, as the matrix-size tends to infinity, has one nonclosed analytic arc. For a generating function having a power singularity, we prove the existence of an asymptotic expansion, that can be used as a theoretical base for the respective numerical algorithm. Different generating functions are explored, highlighting different numerical and theoretical aspects;for example, non-Hermitian and complex symmetric matrix sequences, the reconstruction of the generating function, a consistent eigenvalue ordering, the requirements of high-precision data types. Several numerical experiments are reported and critically discussed, and avenues of possible future research are presented.
We study the numerical solutions for a class of quasilinear Schr & ouml;dinger equations arising from the self-channeling of high-power ultra short lasers in matter, which are associated with energy functionals wi...
详细信息
We study the numerical solutions for a class of quasilinear Schr & ouml;dinger equations arising from the self-channeling of high-power ultra short lasers in matter, which are associated with energy functionals with nonlinear principle parts so that the classical algorithm cannot directly be used. By the method of variable replacement to transform the quasilinear Schr & ouml;dinger equation into an semilinear elliptic equation, the numerical mountain pass algorithm is then applied. Some numerical experiments are also performed, including zero potential, nonzero constant potential and singular problem.
In this article, we deal with a numerical method for the approximation of a class of coupled shape optimization problems, which consist in minimizing an appropriate general volume cost functional subjected to coupled ...
详细信息
In this article, we deal with a numerical method for the approximation of a class of coupled shape optimization problems, which consist in minimizing an appropriate general volume cost functional subjected to coupled boundary value problems by means of a Neumann boundary transmission condition. We show the existence of the shape derivative of the cost functional and express it by means of support functions, using a new formula of shape derivative on a family of convex domains. This allows us to avoid the disadvantages related to the classical shape derivative method using vectors field. Then the numerical discretization is performed using the dual reciprocity boundary element method in order to avert the remeshing task required for the finite element method. Finally, we give some numerical results, based on the gradient method, showing the efficiency of the proposed approach.
In certain pressure and temperature ranges, high-precision thermal equations of state possess multiple roots, with only one being the desired physical solution to the problem. This paper considers methods for setting ...
详细信息
In certain pressure and temperature ranges, high-precision thermal equations of state possess multiple roots, with only one being the desired physical solution to the problem. This paper considers methods for setting initial values to calculate the desired root of the equation. A method for verifying the found root has been developed, enabling the exclusion of nonphysical solutions. The proposed approach allows for the construction of a robust algorithm for numerically solving transcendental equations of state for practical applications.
For automated fiber placement (AFP), it is crucial to research the contact characteristics between tow and roller during the tow feeding process. However, it has not been widely explored. In this paper, a theoretical ...
详细信息
For automated fiber placement (AFP), it is crucial to research the contact characteristics between tow and roller during the tow feeding process. However, it has not been widely explored. In this paper, a theoretical model of contact characteristics for the tow feeding process is established by analyzing the distribution of normal stress, tangential stress, adhesion zone, and slip zone, respectively. The model proposed a numerical algorithm for calculating the contact characteristic. Then the contact characteristic of the tow feeding system can be acquired and validated by the experiment and finite element method. The model that provides theoretical support for AFP optimization. Furthermore, proper pressure (P = 0.6Mpa) and tension (T = 3 N) for towpreg feeding under experimental conditions were obtained.
We give a study to the algorithm for semi-linear parabolic PDEs in Henry-Labordere (2012) and then generalize it to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating the branching process, the...
详细信息
We give a study to the algorithm for semi-linear parabolic PDEs in Henry-Labordere (2012) and then generalize it to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating the branching process, the algorithm does not need any backward regression. To prove that the numerical algorithm converges to the solution of BSDEs, we use the notion of viscosity solution of path dependent PDEs introduced by Ekren et al. (to appear) [5] and extended in Ekren et al. (2012) [6,7]. (C) 2013 Elsevier B.V. All rights reserved.
In this paper, we consider the numerical solutions of three-dimensional axisymmetric nonlinear boundary integral equations with logarithmic kernel. A numerical algorithm with using extrapolation twice is developed to ...
详细信息
In this paper, we consider the numerical solutions of three-dimensional axisymmetric nonlinear boundary integral equations with logarithmic kernel. A numerical algorithm with using extrapolation twice is developed to solve the equations, which possesses the low computing complexities and high accuracy. The asymptotic compact operator theory is used to prove the convergence of the algorithm. The efficiency of the algorithm is illustrated by numerical examples.
暂无评论