• Steel and Composite Structures
• /
• v.15 no.3
• /
• pp.247-266
• /
• 2013

This paper presents a comparison of results obtained by a newly developed numerical model for predicting the behaviour of structures under fire with experimental study carried out on heated and simply supported steel beam elements. A newly developed numerical model consists of three submodels: 3D beam model designed for calculating the inner forces in the structure, 2D model designed for calculation of stress and strain distribution over the cross section, including the section stiffness, and 3D transient nonlinear heat transfer model that is capable of calculating the temperature distribution along the structure, and the distribution over the cross section as well. Predictions of the calculated temperatures and vertical deflections obtained by the numerical model are compared with the results of the inhouse experiment in which steel beam element under load was heated for 90 minutes.

• Structural Engineering and Mechanics
• /
• v.61 no.3
• /
• pp.335-345
• /
• 2017

Free vibration analysis for continuous bridge under any number of vehicles is conducted in this paper. Calculation strategy for natural frequency and mode shape is proposed based on Euler-Bernoulli beam theory and numerical assembly method. Firstly, a half-car planar model is adopted; equations of motion and displacement functions for bridge and vehicle are established, respectively. Secondly, the undermined coefficient matrices for wheels, vehicles, intermediate support, left-end support and right-end support are derived. Then, the numerical assembly technique for conventional finite element method is adopted to construct the overall matrix of coefficients for whole system. Finally, natural frequencies and corresponding mode shapes are determined based on iterative method and overall matrix solution. Numerical simulation is presented to verify the effectiveness of the proposed method. The results reveal that the solutions of present method are exact ones. Natural frequencies and associate modal shapes of continuous bridge under different conditions of vehicles are investigated. The influences of vehicle parameters on natural frequencies are also demonstrated.

• Nuclear Engineering and Technology
• /
• v.50 no.5
• /
• pp.758-766
• /
• 2018

We developed numerical models to efficiently simulate the low-cycle fatigue behavior of a pipe elbow. To verify the model, in-plane cyclic bending tests of pipe elbow specimens were conducted, and a through crack occurred in the vicinity of the crown. Numerical models based on the erosion method and tie-break method are developed, and the numerical results are compared with experimental results. The calculated results of both models are in good agreement with experimental results, and the model using the tie-break method possesses two times faster calculation speed. Therefore, the numerical model based on the tie-break method would be beneficial to evaluate the strength of piping systems under seismic loadings.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.1
• /
• pp.49-66
• /
• 2022

The G-Euler process has been proposed to overcome the difficulties of the calculation of the exponential function of the Jacobian. It is an explicit method that uses the exponential function of the scalar skew-symmetric matrix. We define the moving shapes of true solutions and the moving shapes of numerical solutions. It is discussed whether the moving shape of the numerical solution matches the moving shape of the true solution. The match rates of these two kinds of moving shapes are sequentially calculated by the G-Euler process without using the true solution. It is shown that the closer the minimum match rate is to 100%, the more closely the numerical solutions follow the true solutions to the end. The minimum match rate indicates the reliability of the numerical solution calculated by the G-Euler process. The graphs of the Lorenz system in Perko [1] are different from those drawn by the G-Euler process. By the way, there is no basis for claiming that the Perko's graphs are reliable.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.10
• /
• pp.1667-1675
• /
• 2003

This paper has considered the calculation of lateral stiffness of radial tire's sidewall, which consists of cord stiffness and rubber sheet stiffness, by using the material constants of rubber compounds of tire. We have suggested and illustrated how to calculate the rubber sheet lateral stiffness by considering the following aspects. First, the rubber sheet consists of various kinds of rubber compounds with different thickness along the sidewall in the radial direction. Secondly, equivalent Young's modulus of the rubber sheet can be calculated by using available experimental data of rubber compounds. The present method enables us to divide the calculation domain as many as we want, which can reduce numerical error in the calculation of geometrical and mechanical properties. We have illustrated the calculation by using the data of the radial tire for passenger car of P205/60R15.

• Radiation Oncology Journal
• /
• v.2 no.1
• /
• pp.129-137
• /
• 1984

In brachytherapy, it is important to determine the positions of the radiation sources which are inserted into a patient and to estimate the dose resulting from the treatment. Calculation of the dose distribution throughout an implant is so laborious that it is rarely done by manual methods except for model cases. It is possible to calculate isodose distributions and tumor doses for individual patients by the use of a microcomputer. In this program, the dose rate and dose distributions are calculated by numerical integration of point source and the localization of radiation sources are obtained from two radiographs at right angles taken by a simulator developed for the treatment planning. By using microcomputer for brachytherapy, we obtained the result as following 1. Dose calculation and irradiation time for tumor could be calculated under one or five seconds after input data. 2. It was same value under$\pm2\%$ error between dose calculation by computer program and measurement dose. 3. It took about five minutes to reconstruct completely dose distribution for intracavitary irradiation. 4. Calculating by computer made remarkly reduction of dose errors compared with Quimby's calculation in interstitial radiation implantation. 5. It could calculate the biological isoffect dose for high and low dose rate activities.

• Steel and Composite Structures
• /
• v.26 no.5
• /
• pp.595-606
• /
• 2018

A calculation method was presented to calculate the deflection of GFRP-concrete-steel beams with full or partial shear connections. First, the sectional analysis method was improved by considering concrete nonlinearity and shear connection stiffness variation along the beam direction. Then the equivalent slip strain was used to take into consideration of variable cross-sections. Experiments and nonlinear finite element analysis were performed to validate the calculation method. The experimental results showed the deflection of composite beams could be accurately predicted by using the theoretical model or the finite element simulation. Furthermore, more finite element models were established to verify the accuracy of the theoretical model, which included different GFRP plates and different numbers of shear connectors. The theoretical results agreed well with the numerical results. In addition, parametric studies using theoretical method were also performed to find out the effect of parameters on the deflection. Based on the parametric studies, a simplified calculation formula of GFRP-concrete-steel composite beam was exhibited. In general, the calculation method could provide a more accurate theoretical result without complex finite element simulation, and serve for the further study of continuous GFRP-concrete-steel composite beams.

• Transactions of Materials Processing
• /
• v.28 no.3
• /
• pp.159-166
• /
• 2019

Electromagnetic forming (EMF) simulations consider 3-dimensionally coupled electromagnetic-mechanical phenomenon using LS-DYNA, therefore the calculation cost is normally expensive. In this study, a sensitivity analysis in regard to the simulation variables affecting the calculation time was carried out. The EMF experiments were conducted to form an elliptically protruding shape on a high-strength steel sheet, and it was predicted using LS-DYNA simulation. In this particular EMF simulation case, the effect of several simulation variables, viz., element size, contact condition, EM-time step interval, and re-calculation number of the EM matrices, on the shape of elliptical protrusion and the total calculation time was analyzed. As a result, reasonable values of the simulation variables between the simulation precision and calculation time were proposed, and the EMF experiments with respect to the charging voltages were successfully predicted.

• Progress in Superconductivity and Cryogenics
• /
• v.12 no.1
• /
• pp.28-31
• /
• 2010

A magnetic field calculation method for toroidal type winding which has circular section was developed. At first, the equation for magnetic field by single filament coil was extended using numerical integration to estimate the entire interesting region of solenoid, especially winding region itself. And then, the magnetic field by toroidal arrangement of solenoids was computed with a coordinate transformation of vector fields. The superconducting magnet with toroidal arrangement can be made up of several tens of solenoid type double pancake windings for some applications such as superconducting magnetic energy storage system(SMES). In this system, the field calculation on the high-Tc superconducting(HTS) tape itself is very important because the entire system can be reached to a fault by magnetic stress of conductor or the critical current of superconducting tape can be dramatically reduced under its self field condition. To make matters worse, 3-dimensional analysis is indispensable for this type of magnet and the most of commercial programs with finite element method can be taken too much time for analysis and design. In this paper, a magnetic field calculation method for toroidal type winding with circular section was induced.

• Nuclear Engineering and Technology
• /
• v.54 no.12
• /
• pp.4722-4730
• /
• 2022

This paper introduces a new point depletion-based source term calculation code named BESNA (Bateman Equation Solver for Nuclear Applications), which is aimed to estimate nuclide inventories and source terms from spent nuclear fuels. The BESNA code employs a new modified CE/CM (Constant Extrapolation - Constant Midpoint) predictor-corrector scheme in depletion calculations for improving computational efficiency. In this modified CE/CM scheme, the decay components leading to the large norm of the depletion matrix are excluded in the corrector, and hence the corrector calculation involves only the reaction components, which can be efficiently solved with the Talyor Expansion Method (TEM). The numerical test shows that the new scheme substantially reduces computing time without loss of accuracy in comparison with the conventional scheme using CRAM (Chebyshev Rational Approximation Method), especially when the substep calculations are applied. The depletion calculation and source term estimation capability of BESNA are verified and validated through several problems, where results from BESNA are compared with those calculated by other codes as well as measured data. The analysis results show the computational efficiency of the new modified scheme and the reliability of BESNA in both isotopic predictions and source term estimations.