• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4601-4619
• /
• 2022

The classical Probabilistic Safety Analysis (PSA) does not include any time dependence explicitly. However, the success criteria (SC) could evolve during the cycle for some initiating events. In that sense, there is a type of sequence in which this time-dependency is quite important, the family of Anticipated Transient without Scram (ATWS) sequences in Pressurized Water Reactors. Therefore, a new risk-informed approach is proposed in this paper, which makes it possible to obtain the time-dependent SC evolution of the safety functions affected by the Moderator Temperature Coefficient (MTC) value. Then, the evolution of the ATWS conditional core damage probability (CCDP) could be obtained using a PSA model. To quantify the CCDP, the average values of the time-dependent failure probabilities must be computed. Finally, the comparison between the CCDP obtained through the application of the classical PSA approach and the new one makes it possible to quantify the impact of time-dependence on the SC of the headers that this new risk-informed ATWS approach can provide.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1287-1299
• /
• 2023

A new Monte-Carlo-based computer program (RDS-BASIC) is developed to simulate the transport of energetic ions in pure matter. This computer program is utilizing an algorithm that uses detailed numerical solutions for the classical scattering integral for evaluating the outcomes of the binary collision processes. This approach is adopted by several prominent similar simulation programs and is known to provide results with higher accuracy compared to other approaches that use approximations to shorten the simulation time. Furthermore, RDS-BASIC simulation program contains special methods to reduce the displacement energy threshold of surface atoms. This implementation is found essential for accurate simulation results for sputtering yield in the case of very low energy ions irradiation (near sputtering energy threshold) and also successfully solve the problem of simultaneously obtaining an acceptable number of atomic displacements per incident ions. Results of our simulation for several irradiation systems are presented and compared with their respective TRIM (SRIM-2013) and the state-of-the-art SDTrimSP simulation results. Our sputtering simulation results were also compared with available experimental data. The simulation execution time for these different simulation programs has also been compared.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.151-158
• /
• 2005

This research analyzes eigenvalue analysis of stiffened plates on the Pasternak foundations using the finite clement method. For analyzing the stiffened plates, both the Mindlin plate theory and Timoshenko beam-column theory were applied. In application of the finite element method, 8-nodes serendipity clement system and 3-nodes finite element system were used for plate and beam elements, respectively. Elastic foundations were modeled as the Pasternak foundations in which the continuity effect of foundations is considered. In order to verify the theory of this study, solutions obtained by this analysis were compared with the classical solutions in reference, experimental solutions and solutions obtained by SAP 2000. The natural frequency of stiffened plates on Pasternak foundations were determined according to changes or foundation parameters and dimensions of stiffener.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.10
• /
• pp.1427-1435
• /
• 2010

A mixed finite element model was developed using the classical plate theory to analyze the nonlinear bending of a plate. The appropriate weight functions for the constraints integrated over the domain were determined by the Lagrange multiplier method by using the principle of minimum virtual energy; which provides the constitutive relations between force-like variables and strains. All of detail terms of element wise coefficient matrices and associate tangent matrices to be used in the Newton iterative method are presented. Then, the linear solutions of the current model and those of the traditional displacement model under the SS (simple support) boundary conditions were compared with the existing analytical solution. The post-processed images of the nonlinear results of the force-like variables are presented to show the continuity of the solutions at the joint of the element boundaries. Finally, the converged nonlinear finite element solutions of the current model are compared with those of existing traditional displacement model.

• Tribology and Lubricants
• /
• v.31 no.5
• /
• pp.189-194
• /
• 2015

This paper focuses on the conformal contact problem. A typical example of conformal contact is the contact between a pin and hole. In particular, this paper focuses on the condition for assuming a contact stress field to be a Hertzian pressure profile by using well-known classical solutions associated with Hertzian contact. Persson first developed the conformal contact analysis method around half a century ago, but there have been no significant improvements since then. The present research also adopted this method, but developed new solutions from the viewpoint of application to structural design. The analysis began with a comparison between Persson°Øs conformal contact stress and the Hertzian stress fields. The next step was to check the differences in the normalized stress values of both. This study used the tolerance for the difference in the peak stresses of Persson°Øs solution and the Hertz solution to validate the Hertzian assumption. This gave the range for the difference in radii of the pin and hole when the contact force and mechanical properties of the material are specified. The results showed that, at a tolerance of 5%, the Hertzian assumption is valid if half of the contact angle is less than 35°ý. In addition, the Hertzian assumption holds even for a relatively long contact length, in contrast to the general incomplete contact problem. This paper discusses these results along with other aspects of the application to the design.

• Structural Engineering and Mechanics
• /
• v.64 no.1
• /
• pp.81-92
• /
• 2017

Many materials in engineering exhibit different modulus in tension and compression, which are known as bi-modulus materials. Based on the bi-modulus elastic theory, a modified semi-analytical model, by introducing a stress function, is established in this paper to study the mechanical response of a bi-modulus cylinder placed in an axisymmetric temperature field. Meanwhile, a numerical procedure to calculate the temperature stresses in bi-modulus structures is developed. It is proved that the bi-modulus solution can be degenerated to the classical same modulus solution, and is in great accordance with the solutions calculated by the semi-analytical model proposed by Kamiya (1977) and the numerical solutions calculated both by the procedure complied in this paper and by the finite element software ABAQUS, which demonstrates that the semi-analytical model and the numerical procedure are accurate and reliable. The result shows that the modified semi-analytical model simplifies the calculation process and improves the speed of computation. And the numerical procedure simplifies the modeling process and can be extended to study the stress field of bi-modulus structures with complex geometry and boundary conditions. Besides, the necessity to introduce the bi-modulus theory is discussed and some suggestions for the qualitative analysis and the quantitative calculation of such structure are proposed.

• Structural Engineering and Mechanics
• /
• v.16 no.2
• /
• pp.153-176
• /
• 2003

From the equation of motion of a "bare" non-uniform beam (without any concentrated elements), an eigenfunction in term of four unknown integration constants can be obtained. When the last eigenfunction is substituted into the three compatible equations, one force-equilibrium equation, one governing equation for each attaching point of the concentrated element, and the boundary equations for the two ends of the beam, a matrix equation of the form [B]{C} = {0} is obtained. The solution of |B| = 0 (where ${\mid}{\cdot}{\mid}$ denotes a determinant) will give the "exact" natural frequencies of the "constrained" beam (carrying any number of point masses or/and concentrated springs) and the substitution of each corresponding values of {C} into the associated eigenfunction for each attaching point will determine the corresponding mode shapes. Since the order of [B] is 4n + 4, where n is the total number of point masses and concentrated springs, the "explicit" mathematical expression for the existing approach becomes lengthily intractable if n > 2. The "numerical assembly method"(NAM) introduced in this paper aims at improving the last drawback of the existing approach. The "exact"solutions in this paper refer to the numerical results obtained from the "continuum" models for the classical analytical approaches rather than from the "discretized" ones for the conventional finite element methods.

• Journal of Korean Institute of Industrial Engineers
• /
• v.35 no.2
• /
• pp.150-159
• /
• 2009

The heterogeneous fleet vehicle routing problem(HVRP) is a variant of the classical vehicle routing problem in which customers are served by a heterogeneous fleet of vehicles with various capacities, fixed costs and variable costs. We propose a new conceptual HVRPCR(HVRP with customer restriction) model including additional customer restrictions in HVRP. In this paper, we develop hybrid particle swarm optimization(HPSO) algorithm with 2-opt and node exchange technique for HVRP. The solution representation is a n-dimensional particle for HVRP with N customers. The decoding method for this representation starts with the transformation of particle into a priority list of customer to enter route and limit of vehicle to serve each customer. The vehicle routes are then constructed based on the customer priority list and limit of vehicle to serve. The proposed algorithm is tested using 8 benchmark problems and it consistently produces high-quality solutions, including new best solutions. The numerical results show that the proposed algorithm is robust and efficient.

• Steel and Composite Structures
• /
• v.16 no.6
• /
• pp.577-593
• /
• 2014

An exact dynamic stiffness method is introduced for investigating the free vibration characteristics of the steel-concrete composite beams consisting of a reinforced concrete slab and a steel beam which are connected by using the stud connectors. The elementary beam theory is used to define the dynamic behaviors of the two beams and the relative transverse deformation of the connectors is included in the formulation. The dynamic stiffness matrix is formulated from the exact analytical solutions of the governing differential equations of the composite beams in undamped free vibration. The application of the derived dynamic stiffness matrix is illustrated to predict the natural frequencies and mode shapes of the steel-concrete composite beams with seven boundary conditions. The present results are compared to the available solutions in the literature whenever possible.

• Advances in aircraft and spacecraft science
• /
• v.3 no.1
• /
• pp.1-14
• /
• 2016

The Carrera Unified Formulation (CUF) is here extended to perform free-vibrational analyses of rotating structures. CUF is a hierarchical formulation, which enables one to obtain refined structural theories by writing the unknown displacement variables using generic functions of the cross-section coordinates (x, z). In this work, Taylor-like expansions are used. The increase of the theory order leads to three-dimensional solutions while, the classical beam models can be obtained as particular cases of the linear theory. The Finite Element technique is used to solve the weak form of the three-dimensional differential equations of motion in terms of "fundamental nuclei", whose forms do not depend on the adopted approximation. Including both gyroscopic and stiffening contributions, structures rotating about either transversal or longitudinal axis can be considered. In particular, the dynamic characteristics of thin-walled cylinders and composite blades are investigated to predict the frequency variations with the rotational speed. The results reveal that the present one-dimensional approach combines a significant accuracy with a very low computational cost compared with 2D and 3D solutions. The advantages are especially evident when deformable and composite structures are analyzed.