• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.1 s.244
• /
• pp.66-75
• /
• 2006

Nonlinear response structural optimization using equivalent static loads (NROESL) has been proposed. Nonlinear response optimization is solved by sequential linear response optimization with equivalent static loads which are generated from the nonlinear responses and linear stiffness matrix. The linear stiffness matrix should be obtained in NROESL, and this process can be fairly difficult for some applications. Proportional transformation of loads (PTL) is proposed to overcome the difficulties. Equivalent static loads are obtained by PTL. It is the same as NROESL except for the process of calculating equivalent static loads. PTL is developed for large-scale probems. First, linear and nonlinear responses are evaluated from linear and nonlinear analyses, respectively. At a DOF of the finite element method, the ratio of the two responses is calculated and an equivalent static load is made by multiplying the ratio and the loads for linear analysis. Therefore, the mumber of the equivalent static loads is as many as that of DOF's and an equivalent static load is used with the reponse for the corresponding DOF in the optimization process. All the equivalent static loads are used as multiple loading conditions during linear response optimization. The process iterates until it converges. Examples are solved by using the proposed method and the results are compared with conventional methods.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.609-620
• /
• 1992

A numerical solution of temperature and thermally induced stress in a wafer during rapid thermal processing (R.T.P) is obtained, and an analysis of onset and propagation of slip is performed and compared with experiment. In order to calculate temperature distribution of a wafer in R.T.P system, heat conduction equation that incorporated with radiative and convective heat transfer model is solved, and the solution of the equation is calculated numerically using alternating direction implicit (A.D.I) method. In dealing with radiative heat transfer, a partially transparent body that absorbs the radiation energy is assumed and this transparent body undergoes multiple internal reflections and absorptions. Two dimensional (assuming plane stress) thermoelastic constitutive equation is used to calculate thermal stress induced in a wafer and finite element method is employed to solve the equation numerically. The stress resolved in the slip directions on the slip planes of silicon is compared with the yield stress of silicon in order to predict the slip. The result of the analysis shows that the wafer temperature at which slip occurs is affected by the heating rate of the R.T.P system. It is observed that once slip occurs in the wafer, the slip grows.

• Smart Structures and Systems
• /
• v.32 no.3
• /
• pp.153-166
• /
• 2023

Health monitoring of civil infrastructures, in particular, old bridges that are still in service, has become more than necessary, given the risk that a possible degradation or failure of these infrastructures can induce on the safety of users in addition to the resulting commercial and economic impact. Bridge integrity assessment has attracted significant research efforts over the past forty years with the aim of developing new damage identification methods applicable to real structures. The bridge of Ouled Mimoun (Tlemcen, Algeria) is one of the oldest railway structure in the country. It was built in 1889. This bridge, which is too low with respect to the level of the road, has suffered multiple shocks from various machines that caused considerable damage to its central part. The present work aims to analyze the stability of this bridge by identifying damages and evaluating the damage rate in different parts of the structure on the basis of a finite element model. The applied method is based on an inverse analysis of the normal stress responses that were calculated from the corresponding recorded strains, during the passage of a real train, by means of a set of strain gauges placed on certain elements of the bridge. The results obtained from the inverse analysis made it possible to successfully locate areas that were really damaged and to estimate the damage rate. These results were also used to detect an excessive rigidity in certain elements due to the presence of plates, which were neglected in the numerical reference model. In the case of the continuous bridge monitoring, this developed method will be a very powerful tool as a smart health monitoring system, allowing engineers to take in time decisions in the event of bridge damage.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2008.06a
• /
• pp.480-480
• /
• 2008

Vacuum circuit breaker(VCB) is now emerging as an alternative of gas circuit breaker(GCB) which uses SF6 gas as insulating material whose dielectric strength is outstanding. But we have to reduce SF6 gas because SF6 gas is one of greenhouse gas and efforts to reduce greenhouse gas are now trend of the world. Therefore, we can say VCB is the optimal alternative of GCB because vacuum is environmentally friendly. The vacuum interrupter is the core part of VCB to interrupt arcing current. There are mainly two methods to extinguish arc. One is radial magnetic field (RMF) method and the other is axial magnetic field (AMF) method. We deals with RMF method in this paper. Compared with AMP, RMF arc quenching method has different principle to extinguish arc. In case of RMF method, pinch effect is much larger than AMF method. Because of pinch effect RMF type contact electrodes have the single large spot which is severly damaged and melted while AMF type contact electrodes have small and multiple spots which are slightly damaged and melted. To prevent contact electrode being damaged and melted from high temperature-arc, RMF method uses Lorentz force to move arc. In this paper we calculated and compared the arc driving force of two cases and we analyzed the force acting on each part of arc by means of commercial finite element method software Maxwell 3D. They have 3petals and we considered two cases. One is the case when fixed(upper) and movable(lower) contacts are in mirror arrangement (Case 1). The other is the case when one of two contacts (movable contact) is revolved at maximum angle as possible as it can be (Case 2). And at each case above, we analyzed arc driving force at two positions, position 1 is the closest to the center of contact and position 2 is near the edge of petal on fixed contact. As a result we could find that Case 2 generated stronger arc driving force than Case 1 at position 1. But at position 2 Case 1 generated stronger arc driving force than Case 2. This simulation method can contribute to optimizing spiral-type electrode designs in a view of arc driving force.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.288-293
• /
• 2008

The standard prestressed concrete I-girder bridge (PSC I-girder bridge) is one of the most prevalent types for small and medium bridges in Korea. When determining the member forces in a section to assess the safety of girder in this type of bridge, the general practice is to use the simplified practical equations or the live load distribution factors proposed in design standards rather than the precise analysis through the finite element method or so. Meanwhile, the live load distribution factors currently used in Korean design practice are just a reflection of overseas research results or design standards without alterations. Therefore, it is necessary to develop an equation of the live load distribution factors fit for the design conditions of Korea, considering the standardized section of standard PSC I-girder bridges and the design strength of concrete. In this study, to develop an equation of the live load distribution factors, a parametric analysis and sensitivity analysis were carried out on the parameters such as width of bridge, span length, girder spacing, width of traffic lane, etc. Then, an equation of live load distribution factors was developed through the multiple linear regression analysis on the results of parametric analysis. When the actual practice engineers design a bridge with the equation of live load distribution factors developed here, they will determine the design of member forces ensuring the appropriate safety rate more easily. Moreover, in the preliminary design, this model is expected to save much time for the repetitive design to improve the structural efficiency of PSC I-girder bridges.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.20 no.5
• /
• pp.608-619
• /
• 2017

In this paper, the effects of parameters related to the residual stress induced due to laser shock peening process to determine mitigation of the initial tensile residual stresses are discussed, such as the maximum pressure, pressure pulse duration, laser spot size and number of laser shots. In order to estimate the influence of the initial tensile residual stresses, which is generated by welding in 35CD4 50HRC steel alloy, the initial condition option was employed in the finite element code. It is found that $2{\times}HEL$ maximum pressure and a certain range of the pressure pulse duration time can produce maximum mitigation effects near the surface and depth, regardless of the magnitudes of tensile residual stess. But plastically affected depth increase with increasing maximum pressure and pressure pulse duration time. For the laser spot size, maximum compressive residual stresses have almost constant values. But LSP is more effective with increasing the magnitudes of tensile residual stress. For the multiple LSP, magnitudes of compressive residual stresses and plastically affected depths are found to increase with increasing number of laser shots, but the effect is less pronounced for more laser shots. And to conclude, even though the initial tensile residual stresses such as weld residual stress field are existed, LSP is enough to make the surface and depth reinforcement effects.

• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Tribology and Lubricants
• /
• v.32 no.2
• /
• pp.61-66
• /
• 2016

Total thickness variation (TTV), BOW, and surface roughness are essential characteristics for high quality sapphire substrates. Many researchers have attempted to increase removal rate by controlling the key process parameters like pressure and velocity owing to the high cost of consumables in sapphire chemical mechanical polishing (CMP). In case of the pressure approach, increased pressure owing to higher deviation of pressure over the wafer leads to significant degradation of the TTV. In this study, the authors focused on reducing TTV under the high-pressure conditions. When the production equipment polishes multiple wafers attached on a carrier, higher loads seem to be concentrated around the leading edge of the head; this occurs because of frictional force generated by the combination of table rotation and the height of the gimbal of the polishing head. We believe the skewed pressure distribution during polishing to be the main reason of within-wafer non-uniformity (WIWNU). The insertion of a hub ring between the polishing head and substrate carrier helped reduce the pressure deviation. Adjusting the location of the hub ring enables tuning of the pressure distribution. The results indicated that the position of the hub ring strongly affected the removal profile, which confirmed that the position of the hub ring changes the pressure distribution. Furthermore, we analyzed the deformation of the head via finite element method (FEM) to verify the pressure non-uniformity over the contact area Based on experiment and FEM results, we determined the optimal position of hub ring for achieving uniform polishing of the substrate.

• Structural Monitoring and Maintenance
• /
• v.5 no.2
• /
• pp.205-229
• /
• 2018

Energy induced by minor earthquake and micro vibration cannot be dissipated by traditional buckling-restrained braces (BRBs). To solve this problem, a new type of hybrid passive control device, named as VE-BRB, which is configured by a BRB with high-damping viscoelastic (VE) layers, is developed and studied. Theoretical analysis, performance tests, numerical simulation and case analysis are conducted to study the seismic behavior of VE-BRBs. The results indicate that the combination of hysteretic and damping devices lead to a multi-phased nature and good performance. VE-BRB's working state can be divided into three phases: before yielding of the steel core, VE layers provide sufficient damping ratio to mitigate minor vibrations; after yielding of the steel core, the steel's hysteretic deformations provide supplemental dissipative capacity for structures; after rupture of the steel core, VE layers are still able to work normally and provide multiple security assurance for structures. The simulation results agreed well with the experimental results, validating the finite element analysis method, constitutive models and the identified parameters. The comparison of the time history analysis on a 6-story frame with VE-BRBs and BRBs verified the advantages of VE-BRB for seismic protection of structures compared with traditional BRB. In general, VE-BRB had the potential to provide better control effect on structural displacement and shear in all stages than BRB as expected.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.7
• /
• pp.16-25
• /
• 2002

This paper presents a methodology and application for extracting parasitic inductances in a multi-level interconnect semiconductor structure by a numerical technique. In order to calculate the parasitic inductances, the distrubution of electric potential and current density in the metal lines are calculated by finite element method (FEM). Thereafter, the magneto-static energy caused by the current density in metal lines was calculated. The result of simulation is compared with the result of Grover equation about analytic simple structures, and 4 bit ROM array with a dimension of $13{\times}10.25{\times}8.25{\mu}m^3$ was simulated to extract the parasitic inductnaces. In this calculation, 6,358 nodes with 31,941 tetrahedra were used in ULTRA 10 workstation. The total CPU time for the simulation was about 150 seconds, while the memory size of 20 MB was required.