• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.560-565
• /
• 2007

Motivated by needs such as those in the aerospace industry, this paper demonstrates ability to significantly increase buckling loads of perforated composite laminated plates by synergizing FEM and a genetic optimization algorithm (GA). Plate geometry is discretized into specially-developed 3D degenerated eight-node shell isoparametric layered composite elements. General shell theory, involving incremental nonlinear finite element equilibrium equation, is employed. Fiber orientation within individual plies of each element is controlled independently by the genetic algorithm. Eigen buckling analysis is performed using the subspace iteration method. Available results demonstrate the approach is superior to more conventional methodologies such as modifying ply thickness or the stacking sequence of individual rectilinear plies having common fiber orientation through the plate.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1408-1416
• /
• 1990

A powerful and efficient method for finding approximate solutions to initial-boundary-value problems in the transient coupled thermoelasticity is formulated in time domain using the finite element technique with time-marching strategy. The final system equations can be derived by the Guritin's variational principle using the definition of convolution integral. But, the finite element formulation for the equations of motion is modified by differentiating in time. Numerical results to some test problems are compared with analytical and other sophisticated approximate solutions. Stable responces are observed in all the given examples irrespective of incremental time steps and mesh shapes. In addition, it is shown that good numerical results are obtained even in coarser mesh or larger time step comparing to other numerical methods.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.5
• /
• pp.153-162
• /
• 2007

For the reinforced earth wall constructed on a soft ground in parallel with replacing soft soils, the behavior of the wall according to variations of thickness and stiffness of soft layer, replacement depth, and wall height is investigated using a finite element method, in which incremental construction steps including consolidation of soft soil layer are considered. The behavior of wall is characterized by investigating displacements and settlements developing at the wall, and shear strains developing in a soil deposit. The stability of wall is, then, evaluated by comparing these values with the safety criteria determined on the basis of the literature. Based on the investigation, it is shown that the behavior of wall is influenced naturally from soft soil thickness(t), replacement depth(d) and wall height(h), but more significantly from d and h. In addition, it is also shown that the normalized replacement depth, d/h, required for the safety of wall is not influenced significantly by the variations of t and h. Consequently, it can be concluded that the proper replacement depth can be suggested in an equivalent value in terms of d/h, even for the cases where the wall height is varying with stations, but the variation is not significant.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.746-752
• /
• 2001

Design sensitivity analysis for nonlinear structural problems has been emerged in the last decade as a glowing area of engineering research. As a result, theoretical formulations and computational algorithms have already developed for design sensitivity of nonlinear structural problems. There is not enough research for practical nonlinear problems using multi-element, due to difficulties of implementation into FEA. Therefore, nonlinear response analysis for stiffened shell which consists of Mindlin plate and Timoshenko beam, was considered. Specially, it presents the backward-Euler method which is adopted to describe an exact yield state in the stress computation procedure. Then, design sensitivity analysis of nonlinear structures, particularly elasto-perfectly-plastic structure, is developed using direct differentiation method. The accuracy of the developed sensitivity analysis was compared with the central finite difference method. Finally, on the basis of above results, design improvement for stiffened shell is suggested.

• Journal of Korean Society of Steel Construction
• /
• v.15 no.3
• /
• pp.299-311
• /
• 2003

This paper presents a finite element procedure to estimate the ultimate strength of space frames considering spread of plasticity. The improved displacement field is introduced based on the inclusion of second-order terms of finite rotations. All the non-linear terms due to bending moment, torsional moment, and axial force are precisely considered. The concept of plastic hinges is introduced and the incremental load/displacement method is applied for elasto-plastic analyses. The initial yield surface is defined based on the residual stress, and the full plastification surface is considered under the combined action of axial forces, bending and torsional moments. The elasto-plastic stiffness matrices are derived using the flow rule and the normality condition of the limit function. Finite element solutions for the ultimate strength of space frames are compared with available solutions and experimental results.

• Restorative Dentistry and Endodontics
• /
• v.25 no.1
• /
• pp.91-108
• /
• 2000

Most of cervical abrasion and erosion lesions show gingival margin where the cavosurface angle is on cementum or dentin. Composite resin restoration of cervical lesion shrink toward enamel margin due to polymerization contraction. This shrinkage has clinical problem such as microleakage and secondary caries. Several methods to diminish contraction stress of composite resin restoration, such as modifying cavity form and building up restorations in several increments have been attempted. The purpose of this study was to compare polymerization contraction stress of composite resin in Class V cavity subjected to cavity forms and placement methods. In this study, finite element model of 5 types of Class V cavity was developed on computer tomogram of maxillary central incisor. The types are : 1) Box cavity 2) Box cavity with incisal bevel 3) V shape cavity 4) V shape cavity with incisal bevel 5) Saucer shape cavity. The placement methods are 1) Incisal first oblique incremental curing 2) Bulk curing. An FEM based program for light activated polymerization is not available. For simulation of curing dynamics, time dependent transient thermal conduction analysis was conducted on each cavity and each placement method. For simulation of polymerization shrinkage, thermal stress analysis was performed with each cavity and each placement method. The time-temperature dependent volume shrinkage rate, elastic modulus, and Poisson's ratio were determined in thermal conduction data. The results were as follows : 1. With all five Class V cavifies, the highest Von Mises stress at the composite-tooth interface occurred at gingival margin. 2. With box cavity, V shape cavity and saucer cavity, Von Mises stress at gingival margin of V shape cavity was lower than the others. And that of box cavity was lower than that of saucer cavity. 3. Preparing bevel at incisal cavosurface margin decreased the rate of stress development in early polymerization stage. 4. Preparing bevel at incisal cavosurface margin of V shape cavity increased the Von Mises stress at gingival margin, but decreased at incisal margin. 5. At incisal margin, stress development by bulk curing method was rapid at early stage. Stress development by first increment of incremental curing method was also rapid but lower than that by bulk curing method, however after second increment curing final stress was the same for two placement methods. 6. At gingival margin, stress development by incremental curing method was suddenly rapid at early stage of second increment curing, but final stress was the same for two placement methods.

• Structural Engineering and Mechanics
• /
• v.58 no.2
• /
• pp.243-257
• /
• 2016

The aim of this study is to investigate the seismic resistance of dry stone arches under in-plane seismic loading. For that purpose, several numerical analyses were performed using the combined finite-discrete element method (FDEM). Twelve types of arches with different ratios of a rise at the mid-span to the span, different thicknesses of stone blocks and different numbers of stone blocks in the arch were subjected to an incremental dynamic analysis based on excitation from three real horizontal and vertical ground motions. The minimum value of the failure peak ground acceleration that caused the collapse of the arch was adopted as a measure of the seismic resistance. In this study, the collapse mechanisms of each type of stone arch, as well as the influence of the geometry of stone blocks and stone arches on the seismic resistance of structures were observed. The conclusions obtained on the basis of the performed numerical analyses can be used as guidelines for the design of dry stone arches.

• Structural Engineering and Mechanics
• /
• v.77 no.1
• /
• pp.19-35
• /
• 2021

Existing methods to estimate the probability of seismic pounding occurrence of adjacent buildings do not account for nonlinear behavior or only apply to simple lumped mass systems. The present study proposes an efficient method based on subset simulation for fragility and risk assessment of seismic pounding occurrence between nonlinear adjacent buildings neglecting pounding effects with application to finite element models. The proposed method is first applied to adjacent buildings modeled as elastoplastic systems with substantially different dynamic properties for different structural parameters. Seismic pounding fragility and risk of adjacent frame structures with different floor levels is then assessed, paying special attention to modeling the non-linear material behavior in finite element models. Difference in natural periods and impact location are identified to affect the pounding fragility simultaneously. The reliability levels of the minimum code-specified separation distances are also determined. In addition, the incremental dynamic analysis method is extended to assess seismic pounding fragility of the adjacent frame structures, resulting in higher fragility estimates for separation distances larger than the minimum code-specified ones in comparison with the proposed method.

• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1570-1578
• /
• 2022

During a severe accident, steam generator (SG) tubes undergo rapid changes in the pressure and temperature. Therefore, an appropriate creep model to predict a short term creep damage is essential. In this paper, a novel creep model for Alloy 690 SG tube material was proposed. It is based on the theta (θ) projection method that can represent all three stages of the creep process. The original θ projection method poses a limitation owing to its inability to represent experimental creep curves for SG tube materials for a large strain rate in the tertiary creep region. Therefore, a new modified θ projection method is proposed; subsequently, a master curve for Alloy 690 SG material is also proposed to optimize the creep model parameters, θ_i (i = 1-5). To adapt the implicit creep scheme to the finite element code, a partial derivative of incremental creep with respect to the stress is necessary. Accordingly, creep model parameters with a strictly linear relationship with the stress and temperature were proposed. The effectiveness of the model was validated using a commercial finite element analysis software. The creep model can be applied to evaluate the creep rupture behavior of SG tubes in nuclear power plants.

• Structural Engineering and Mechanics
• /
• v.89 no.1
• /
• pp.23-31
• /
• 2024

Ratcheting is the cyclic buildup of inelastic strain on a structure resulting from a combination of primary and secondary cyclic stress. It can lead to excessive plastic deformation, incremental collapse, or fatigue. Ratcheting has been numerically investigated on a cantilever beam, considering the current study's primary and secondary bending loads. In addition, the effect of input frequency on the onset of ratcheting has been investigated. The non-linear dynamic elastic-plastic approach has been utilized. Analogous to Yamashita's bending-bending ratchet diagram, a non-dimensional ratchet diagram with a frequency effect is proposed. The result presents that the secondary stress values fall sequentially with the increase of primary stress values. Moreover, a displacement amplification factor graph is also established to explain the effect of frequency on ratchet occurrence conditions. In terms of frequency effect, it has been observed that the lower frequency (0.25 times the natural frequency) was more detrimental for ratchet occurrence conditions than the higher frequency (2 times the natural frequency) due to the effect of dynamic displacement. Finally, the effect of material modeling of ratcheting behavior on a beam is shown using different hardening coefficients of kinematic hardening material modeling.