For safe and economical spent fuel management, assessing the integrity of the cladding, which is the first barrier to the escape of radioactive material, is very important. For the sake of risk assessment, it is essential to calculate the probability of failure of the spent fuel rods loaded inside the cask during the transportation or storage. However, due to the large amounts of calculations required, it is not practical to analyze every detail of the spent fuel rods and assemblies. This study presents a methodology to perform a cask-level analysis by sequentially simplifying the fuel rods and spent fuel assemblies for the calculation of fuel rod failure probability. A simplified single fuel rod model was generated by considering the material properties of a high burnup fuel rod stored in dry storage for approximately 5 years and the interfacial bonding conditions of the cladding tube. The simplified model produces the same deflection as the detailed model at the critical moment that produces a fracture plastic strain of 1%. The developed single fuel rod simplified model is assembled in a CE 16×16 configuration, and a methodology is presented in which the CE 16×16 assembly model is once again replaced by a simplified model with a cuboidal shape. Compression analyses were performed on each part of the CE 16×16 model to obtain isotropic property data, and a simplified model was created based on those data and the cross-sectional second moment values of the parts. A cask drop analysis was performed to validate the similarity of the CE 16×16 model and the simplified model by comparing important structural responses such as impact acceleration. The 20 simplified fuel assembly models and one detailed model were loaded into a cask to perform the drop analysis. For the detailed model, the impact acceleration was extracted for different loading positions and the corresponding impact load and pinch load were derived. The spring force and contact force corresponding to the pinch load were extracted by applying a Python script technique to extract the maximum value of them exerted on each fuel rod. The vulnerability of spent fuel rods to bending loads and the failure criteria were considered during the simplification process of a single fuel rod. From the extracted impact and pinch loads, the probability of failure of the spent fuel rods as a function of impact acceleration can be calculated.
Laser welding is used in various industries due to fast welding speed, and prediction with FEM is needed to design. The multi-layered heat source model that can cover various models has been recently proposed to analyze laser welding, but it takes a lot of time because more than 2,000 cases are compared to derive the results. In order to reduce the time, the simplified model should be suggested. In order to derive a simplified model, laser welding heat transfer analysis was performed using the welding direction length and the convection coefficient of the welding direction surface as parameters. As a result, the model whose difference is less than 0.1% compared to reference model is deducted. The analysis time was reduced by about 90% from 48 hours to 5 hours. Also difference of convection coefficient on the welding direction surface does not affect the temperature distribution much.
As laser welding has low thermal distortion and fast welding speed, the needs in the field is increasing in varous industries. Prediction of distortion and design with that prediction are very important for using welding process in the field. However many types of heat source models for laser welding are suggested, field engineers feels difffucult for using a proper model for a specific case. Thus I, author, suggested a representative model which can cover most of existing models with multi-layered heat source model. This method is very powerful, but there are much time consumption with analysis and comparision among more than 1,000 candidates. To solve these shortcomings, this study focused on to find the simplified model. In order to construct the simplified model, the length of the orthogonal direction to the welding direction was reduced from 300mm to 35mm, which makes time reduction of 75% with sustaining the quality of the original model. That can cover the disadvantage of multi-layered heat source model for laser welding, enhance the prediction of welding distortion after laser welding
For economic and safe management of Spent Nuclear Fuel (SNF), it is very important to maintain the structural integrity of SNF and to keep the fuel undamaged and handleable. The cladding surrounding nuclear fuel must be protected from physical and mechanical deterioration. The structural evaluation of SNF is very complicated and numerically demanding and it is essential to develop a simplified model for the fuel rod. In this study, a simplified model was developed using a new cladding failure criterion. The simplified model was developed considering only the horizontal or lateral static load utilizing the cladding material properties of irradiated Zirclaoy-4, and applicability in horizontal and vertical drop impacts was investigated. When a fuel rod is subject to bending, a very complicated 3D stress state is generated within the vicinity of the pellet–pellet interface. A very localized stress concentration is observed in the area where the edges of the pellets contact the cladding. If the failure strain criteria obtained from the uniaxial tension test or biaxial tube test is applied, failure is predicted at the beginning stage of loading with premature through-thickness stress or strain development. The localized contact stress or strain is self-limiting and is not a good candidate for the cladding failure criteria. In this work, a new cladding failure criterion is proposed, which can account for the localized stress concentration and the through-thickness stress development. The failure of the cladding is determined by the membrane plus bending stress generated through the thickness of the cladding, which can be calculated by a process called stress linearization along the stress classification line. The failure criterion for SNF was selected as the membrane plus bending stress through stress linearization in the cross-sections through the thickness of the cladding. Because the stress concentration in the cladding around the vicinity of the pellet–pellet interface cannot be simulated in a simplified beam model, a stress correction factor is derived through a comparison of the simplified model and detailed model. The applicability of the developed simplified model is checked through horizontal and vertical drop impact simulations. It is shown that the stress correction factor derived considering static bending loading can be effectively applied to the dynamic impact analyses in both horizontal and vertical orientations.
The estimation of heat source model is very important for heat transfer analysis with finite element method. Part I of this study used adaptive simulated annealing which is one of the global optimization algorithm for anticipating the parameters of the Goldak model. Although the analysis with 3D model which depicted the real situation produced the correct answer, that took too much time with moving heat source model based on Fortran and Abaqus. This research suggests the procedure which can reduce time with maintaining quality of analysis. The lead time with 2D model is reduced by 90% comparing that of 3D model, the temperature distribution is similar to each other. That is based on the saturation of heat transfer among the direction of heat source movement. Adaptive simulated annealing with 2D model can be used to estimate more proper heat source model and which could enhance to reduce the resources and time for experiments.
Anticipation of welding deformation with finite element method is a very interested topic in the industries, adequate heat source model is essential for concluding reasonable results. This study is related to estimate the parameters of Goldak heat source model, and global optimization algorithm is applied to this research. The heat affected zone (HAZ) boundary line of bead on plate (BOP) welding is used as the target, parameters of heat sources are used as the variables. Adaptive simulated annealing is applied and the optimal result is obtained out of 1,000 candidates. The convergence of finite element method and the global optimization is meaningful for estimation of welding deformation, which could enhance to reduce the resources and time for experiments.
In this paper, for a seismic analysis of an offshore subsea manifold, Response Spectrum Analysis(RSA) and Time History Analysis(THA) were conducted under a various analysis conditions. Response spectrum and seismic design procedure have followed ISO19901-2 code. In case of THA, The response spectrum were converted into artificial earthquake history and both of Explicit and Implicit solvers were used to examine the characteristics of seismic analysis. For the verification, Various seismic analysis methods were applied on a single degree of freedom beam model and a simplified model of the actual manifold. The difference between the results of RSA and THA on the simplified manyfold model evaluated for the analysis of the actual manifold. Because THA is impossible in case of real complex structure such as a manifold, Safety of the actual manifold structure was accessed by using the RSA and the difference between the results of RSA and THA from the simplified model.
The polyurethane spring, which provides the restoring force for a sliding bridge bearing, governs the overall behavior of the bearing. An analytical model, therefore, which reflects actual geometrical shape and can be utilized easily in engineering practice, is needed to be developed. In this paper, a simplified model for analyzing the polyurethane spring utilized in the sliding bridge bearing is presented. The model is developed for a cylindrical tube based on Gent’s model, and can be utilized for various height-diameter ratios. The accuracy of the proposed model are compared to existing models via experimental data as well as numerical analysis results. The comparison result shows that the proposed model can be used in estimating mechanical behavior of the cylindrical tube shaped polyurethane spring.
본 논문에서는 앵커볼트의 체결력을 고려한 유한요소 모델을 제안함으로써 앵커볼트로 연결된 부재의 해석적 평가방법을 제안하였다. LS-DYNA를 사용한 유한요소 모델링은 복잡한 3차원 상세모델보다는 단순화된 앵커모델들을 적용함으로써 해석 효율성을 고려하였다. 앵커볼트는 Beam 혹은 Solid 요소로 토크 조임에 따른 앵커볼트 긴장상태를 반영하였고, 토크 조임에 따른 체결력을 고려하기 위해 너트면에 압축력을 도입하거나 너트를 Shell 혹은 Solid로 고려하여 작용 토크값으로 산정되는 체결력과 등가의 하중을 도입하였다. 외력 작용 시 체결력과 마찰력에 의한 앵커볼트의 하중전달은 nodal rigid 혹은 contact 조건으로 고려하였다. 체결력을 고려한 세 종류의 앵커모델을 적용한 해석결과, Model I과 Model III는 볼트 축력과 전단력이 매우 유사한 값으로 계산되었고, Model II의 경우, 볼트 축력과 전단력이 다소 과소평가되는 것으로 나타났다. Model I은 다른 두 모델에 비해 수치해석적으로 효율적인 것으로 분석되었다.
2009년말 기준으로 11,811 다발의 경수로 사용후핵연료가 방출되었으며, 지금까지 각 사용후핵연료에 대해 방사선원항을 가중하여 설계에 반영하기는 사실상 불가능하여, 원자력 관련시설 설계시 보수성을 갖는 기준 사용후핵연료를 선정하고 이를 바탕으로 시스템 설계를 수행하여 왔다. 방사선원항에 대한 단순모델을 적용하면 각 사용후핵연료에 대한 방사선원항을 가중함으로써 이와 같은 보수성을 배제할 수 있으므로 본 연구에서 웨스팅하우스형 원전에 사용된 사용후핵연료를 대상으로 방사선원항, 즉, 붕괴열, 방사능세기, 섭취위해도 등을 예측하기 위한 회귀모형을 개발하였다. 개발된 회귀식을 통해 예측된 방사선원항값은 ORIGEN-ARP 코드로 계산된 값과 약 5% 이내에서 잘 일치함을 확인하였으며, 이의 유용성을 검토한 결과 각각의 사용후핵연료에 대한 방사선원항을 가중하여 설계에 반영하면 보수성을 줄일 수 있음을 확인하였다. 따라서 본 연구에서 개발된 회귀식은 사용후핵연료의 저장 및 처분과 관련한 원자력시설 설계시 개념설계 단계에서 유용하게 사용될 수 있을 것으로 판단된다.
특정 구획으로 유입된 오염물질이 해당 구획 내부에 순간적으로 균일하게 분포한다는 구획모델에 대한 기본가정의 한계로 인해, 전통적인 단일구획모델로는 불포화대에서 오염물질의 이동현상을 적절하게 예측할 수 없다. 한편 물리적으로 동일한 불포화대를 여러 개의 구획으로 구분한 다중구획모델링 기법은 실제 불포화대에서의 오염물질 이동 지연효과를 적절하게 설명할 수 있으나, 지금까지 일반적인 해석해가 보고된 바 없으며 고려하는 구획의 개수가 증가할수록 모델링에 많은 시간이 소요되는 등의 한계가 있다. 이러한 문제점을 해결하기 위하여, 이류가 지배적인 조건 하에서 불포화대의 오염물질 이동현상에 대한 다중구획모델을 수립하고 이를 해석적인 방법으로 계산할 수 있는 일반해를 유도한 후, 다중구획모델을 단일구획모델로 근사할 수 있는 수학적 제약조건을 도출하였다. 단순화된 근사방법론의 유효성은 가상적인 조건 하에서 간단한 수치해석적 방법을 통해 검증하였다. 물리적으로 동일한 특성을 갖는 불포화대를 단일 구획으로 가정할 경우, 불포화대로부터 포화대로 유입되는 오염물질의 전이율은 상수가 아닌 시간 종속적인 명목전이율로 표현할 수 있음을 증명하였다. 또한 명목전이율은 불포화대 구획간 전이율에 민감하며 오염층으로부터의 전이율에 대한 민감도는 미미한 것으로 나타났다. 이 연구에서 개발된 단순화된 근사방법론은 많은 시간이 요구되는 다중구획 모델링을 통하지 않고 불포화대 오염물질 이동현상을 신속하고 합리적으로 예측하기 위한 목적으로 활용될 수 있을 것으로 기대된다.
A steel frame is one of the most commonly used structural systems due to its resistance to various types of applied loads. Many studies have been conducted to investigate the effects of connection flexibility, support conditions, and beam-to-column stiffness ratio on the story drift of a frame. Based on the results of these studies, several design guides have been proposed. This research has been conducted to predict the actual behavior of a double angle connection, and to establish its effect on the story drift and the maximum allowable load of a steel frame. For these purposes, several experimental tests were conducted and a simplified analytical model was proposed. This simplified analytical model consists of four spring elements as well as a column member. In addition, a point bracing system was proposed to control the excessive story drift of an unbraced steel frame.
본 연구에서는 사각형 모듈의 인발성형된 복합재료 바닥판의 휨 거동에 대한 해석 모델을 개발하였다. FRP 바닥판의 해석 모델은 FSDT 평판 이론을 기반으로 임의 적층각을 지닌 FRP 바닥판의 처짐을 예측할 수 있었다. 수치적 예제에서는 네 변이 단순 지지된 등분포 하중을 받는 사각형 모듈의 FRP 바닥판을 2차원 평판 유한 요소해석을 적용하여 수행하였고, 해석 결과에 대해서는 바닥판 길이-높이의 비와 화이버 각도의 변화에 따른 효과에 대해 역점을 두고 다루었다. 연구 결과, 본 연구에서 제안한 해석 모델이 FRP 바닥판의 휨 거동을 해석하고 예측하는데 효과적이고 정확하다는 것이 입증되었다. 또한, FRP 바닥판의 높이가 높아질수록 plate 해석 이론에 있어서 일차전단변형이론(First order Shear Deformable laminated plate Theory : FSDT)이 아닌 고차전단변형(Higher order Shear Deformable plate Theory : HSDT)의 필요성에 대해 언급하였다.