• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.151-158
• /
• 2005

For the safety analysis of large structures such as nuclear containment buildings, we conventionally prefer to use analytical approach using finite element method rather than empirical test. Therefor, this paper is mainly focused to develop low-order solid finite element model with the elasto-plastic material model for the safety analysis of nuclear containment building. Drucker-Prager failure criteria in uncracked concrete and maximum tensile stress criteria in cracked concrete are used to model the constitutive behavior of concrete. The concrete material model takes into account the aspects of tensile strain, compression strength reduction of concrete and shear transfer to improve the accuracy of the finite element analysis. Finally, numerical simulation to compare the performance of the developed model with experimental results is employed. The numerical results in this study agree very well with the experimental data.

• Transactions of Materials Processing
• /
• v.25 no.5
• /
• pp.325-331
• /
• 2016

In this paper, a three-dimensional part model is constructed for the finite element analysis of hemming models where hemming defects frequently occur. The roller path is modeled as the boundary condition with the one-dimensional beam element and the revolute joint model. With the constructed part model and the roller movement, a finite element analysis has been pursued in order to identify the hemming load and hemming defects such as wrinkling in the flange region. The analysis result shows that the maximum hemming load occurs in the intake situation while oscillatory behavior of the load is found especially when hemming the curved model because of wrinkle generation. This paper compares the amplitude and the period of wrinkle between the analysis result and the experiment, which shows good agreement with each other.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2003.04a
• /
• pp.352-361
• /
• 2003

A fatigue damage accumulation model based on the continuum damage mechanics theory was develope(1 where modules decay ratios in tension and shear on used as indicators for damage variables D . In the model, the damage variables are considered to be second-order tensors. Then the maximum principal damage variable, $D^*$ is introduced According to the similarity to the Principal stress, $D^*$ is obtained as the maximum eigen value of damage tensor [D']. Under proportional tension and torsion loadings, fatigue lives were satisfactorily predicted at any combined stress ratios using the present model in which the fatigue characteristics only under uniaxial tension and pure torsion loadings on needed. Fatigue life prediction under uniaxial tension and pure torsion loadings, was performed based on the damage mechanics using boundary element method.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.13 no.6
• /
• pp.64-73
• /
• 2004

A fatigue damage accumulation model based on the continuum damage mechanics theory was developed where modulus decay ratios in tension and shear were used as indicators for damage variables D. In the model, the damage variables are considered to be second-order tensors. Then, the maximum principal damage variable, $D^*$ is introduced. According to the similarity to the principal stress, $D^*$ is obtained as the maximum eigen value of damage tensor [D]. Under proportional tension and torsion loadings, fatigue lives were satisfactorily predicted at any combined stress ratios using the present model in which the Fatigue characteristics only under uniaxial tension and pure torsion loadings were needed. Fatigue life prediction under uniaxial tension and pure torsion loadings, was performed based on the damage mechanics using boundary element method.

• Journal of the Korea Concrete Institute
• /
• v.11 no.3
• /
• pp.3-12
• /
• 1999

This study presents a suggestion of regulation of skewed slab bridge. In order to find the characteristic behavior of skew bridge, many cases of skew bridges were analyzed with changed angle of skew. The comparison of design methods for cantilever part in pier was also made. It was found that : (1) The lower the skew angle was, the higher the maximum support reaction forces at the end point were. (2) The higher the ratio of L/B was, the higher the maximum support reaction force at the point was. (3) The effect of skew may be neglected for skew angles of $70^{\circ}$or more. (4) If elastic springs are applied to the boundary conditions to simulate the rubber pad bearings, the results will be more reasonable. (5) The shear deformation effect must be considered in the analysis of cantilever part of substructure. (6) Using strut and tie model to design cantilever part of pier, it will be more simple than finite element method with same accuracy and more accurate than using frame element.

• Journal of applied mathematics & informatics
• /
• v.37 no.5_6
• /
• pp.357-382
• /
• 2019

In this paper, we consider boundary value problem for a weakly coupled system of two singularly perturbed differential equations of convection diffusion type with discontinuous source term. In general, solution of this type of problems exhibits interior and boundary layers. A numerical method based on streamline diffusiom finite element and Shishkin meshes is presented. We derive an error estimate of order $O(N^{-2}\;{\ln}^2\;N$) in the maximum norm with respect to the perturbation parameters. Numerical experiments are also presented to support our theoritical results.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.2
• /
• pp.67-75
• /
• 2004

This study aims to investigate the dynamic behaviour of a parabolic arch with initial deflection by using the elasto-plastic finite element model where the von-Mises yield criteria have been adopted. The initial deflection of arch was assumed by the high order polynomial of ${\omega}_i$ = ${\omega}_o$${(1-{(2x/L)}^m)}^n$) and the sinusoidal profile of ${\omega}_i$ = ${\omega}_o$$\sin$(n$\pi$x/L). Several numerical examples were tested considering symmetric initial deflection modes when the maximum initial deflection of an arch is fixed as L/500, L/1000, L/2000 or L/5000. The effects of polynomials order on the dynamic behavior of arch were not conspicuous. The most unfavorite dynamic response occurs when the maximum initial deflection varies from L/1000 to L/4000 if the initial deflection mode is represented by high order polynomials.

• Computational Structural Engineering : An International Journal
• /
• v.3 no.1
• /
• pp.19-29
• /
• 2003

In the present study, finite element linear buckling analysis is performed for grid-stiffened composite plates. A hybrid element with drilling degrees of freedom is employed to reduce the effect of the sensitivity of mesh distortion and to match the degrees of freedom between skins and stiffeners. The preliminary static stress distribution is analyzed for the determination of accurate load distribution. Parametric study of grid structures is performed and three types of buckling modes are observed. The maximum limit of buckling load was found at the local skin-buckling mode. In order to maximize buckling loads, stiffened panels need to be designed to be buckled in skin-buckling mode.

• Structural Engineering and Mechanics
• /
• v.59 no.2
• /
• pp.243-259
• /
• 2016

This research develops a finite element code for the transient dynamic analysis of tapered fiber reinforced polymer (FRP) poles with hollow circular cross-section and flexible joints used in power transmission lines. The FRP poles are modeled by tapered beam elements and their flexible joints by a rotational spring. To solve the time equations of transient dynamic analysis, precise time integration method is utilized. In order to verify the utilized formulations, a typical jointed FRP pole under step, triangular and sine pulses is analyzed by the developed finite element code and also ANSYS commercial finite element software for comparison. Thereafter, the effect of joint flexibility on its dynamic behavior is investigated. It is observed that by increasing the joint stiffness, the amplitude of the pole tip deflection history decreases, and the time of occurrence of the maximum deflection is earlier.

• Nuclear Engineering and Technology
• /
• v.52 no.6
• /
• pp.1137-1147
• /
• 2020

The discrete ordinates method (S_N) is one of the major shielding calculation method, which is suitable for solving deep-penetration transport problems. Our objective is to explore the available quadrature sets and to improve the accuracy in shielding problems involving strong anisotropy. The linear discontinuous finite-element (LDFE) quadrature sets based on the icosahedron (in short, ICLDFE quadrature sets) are developed by defining projected points on the surfaces of the icosahedron. Weights are then introduced in the integration of the discontinuous finite-element basis functions in the relevant angular regions. The multivariate secant method is used to optimize the discrete directions and their corresponding weights. The numerical integration of polynomials in the direction cosines and the Kobayashi benchmark are used to analyze and verify the properties of these new quadrature sets. Results show that the ICLDFE quadrature sets can exactly integrate the zero-order and first-order of the spherical harmonic functions over one-twentieth of the spherical surface. As for the Kobayashi benchmark problem, the maximum relative error between the fifth-order ICLDFE quadrature sets and references is only -0.55%. The ICLDFE quadrature sets provide better integration precision of the spherical harmonic functions in local discrete angle domains and higher accuracy for simple shielding problems.