• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.79-87
• /
• 2013

An isogeometric topological shape optimization method is developed using the level sets and Heaviside enrichments. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set functions, which facilitates to handle complicated topological shape changes. The Heaviside enrichment improves the isogeometric analysis by adding some enrichment functions to model the internal boundaries. The proposed topological shape optimization method has several benefits: exact geometric models can be obtained using the isogeometric approach and the limitation of tensor-product patches can be overcome using the Heaviside enrichments to represent the internal voids. Even in a single patch, discontinuous displacement fields as well as smooth stress field can be obtained. Since the level sets offer the implicit moving boundary inside the domain, it is easy to represent the topological shape variations in the isogeometric analysis using Heaviside enrichments.

• Computers and Concrete
• /
• v.2 no.1
• /
• pp.55-77
• /
• 2005

The closed form solution of the equilibrium equations in the ultimate design of reinforced concrete sections under biaxial bending is presented. The stresses in the materials are described by the Model Code 1990 equations. Computation of the integral equations is performed generally in terms of all variables. The deformed shape of the section in the ultimate conditions is defined by Heaviside functions. The procedure is convenient for the use of mathematical manipulation programs and the results are easily included into nonlinear analysis codes. The equations developed for rectangular sections can be applied for other sections, such as T, L, I for instance, by decomposition into rectangles. Numerical examples of the developed model for rectangular sections and composed sections are included.

• Computers and Concrete
• /
• v.10 no.5
• /
• pp.541-555
• /
• 2012

This paper presents an algorithm for the evaluation of stresses in reinforced concrete sections under service loads. The algorithm is applicable to any section defined by polygonal contours and is based on an analytical integration of the stresses. The nonlinear behaviour of concrete is represented by the parabola-rectangle law used in the Eurocode-2 for the ultimate concrete design. An integrated definition of the strains in concrete and steel is possible by the use of Heaviside functions, similarly to what is done for ultimate section design in Barros et al. (2004). Other constitutive equations for the definition of the stresses in the concrete or steel can be easily incorporated into the code. The examples presented consist in the evaluation of resulting axial load and bending moment in an irregular section and in a section in L shape. The results, for service stresses, can also be plotted in terms of design abacus; a rectangular doubly reinforced section is presented as example.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.41-45
• /
• 2002

The zone-conditioned CMC equations are derived by taking an unconditional average of the generic conservation equations multiplied by delta and Heaviside functions in terms of mixture fraction and reaction progress variable. The resulting equations are essentially in the same form as the single zone CMC equations except for separate flow fields for burned and unburned gas. The zone-conditioned two-fluid equations are applied to a stagnating turbulent premixed flame brush of Cheng and Shepherd[5l. It is shown that the flame stretch factor is of crucial importance to accurately reproduce the measured mean reaction progress variable and conditional velocities. Further work is in progress for the relationship between surface and volume averages and extension to partially premixed combustion on the basis of a triple flame structure, e. g. in a lifted turbulent diffusion flame.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.4 s.97
• /
• pp.395-403
• /
• 2005

The dynamic characteristics of delaminated smart composite laminates are studied using animproved layerwise laminate theory. The theory is capable of capturing interlaminar shear stresses that are critical to delamination. The presence of discrete delamination is modeled through the use of Heaviside unit step functions. Stress free boundary conditions are enforced at all free surfaces. Continuity in displacement field and transverse shear stresses are enforced at each ply level. In modeling piezoelectric composite plates, a coupled piezoelectric-mechanical formulation is used in the development of the constitutive equations. Numerical analysis is conducted to investigate the effect of nonlinearity in the transient vibration of bimodular behavior caused by the contact impact of delaminated interfaces. Composite plates with delamination, subject to external loads and voltage history from surface bonded sensors, are investigated and the results are compared with corresponding experimental results and plates without delamination.

• Computers and Concrete
• /
• v.19 no.4
• /
• pp.405-411
• /
• 2017

Carbon Fiber Reinforced Plastic (CFRP) has commonly been used to strengthen existing RC structures. Wrapping the whole component with CFRP is an effective method and simple to execute. Besides, specific configuration of CFRP sheets (L, X and T shape) has also been considered in some experiments to examine CFRP effects in advance. This study aimed to provide an optimal CFRP configuration to effectively retrofit the beam-column connection using continuous material topology optimization procedure. In addition, Moved and Regularized Heaviside Functions and penalization factors were also considered. Furthermore, a multi-material procedure was also used to compare with the results from the single material procedure.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.11a
• /
• pp.469-475
• /
• 2004

The transient analysis of delaminated smart composite laminates is studied using an improved layerwise laminate theory. The theory is capable of capturing interlaminar shear stresses that are critical to delamination. The presence of discrete delamination is modeled through the use of Heaviside unit step functions. Stress free boundary conditions are enforced at all fee surfaces. Continuity in displacement field and transverse shear stresses are enforced at each ply level. In modeling piezoelectric composite plates, a coupled piezoelectric-mechanical formulation is used in the development of the constitutive equations. Numerical analysis is conducted to investigate the effect of nonlinearity in the transient vibration of bimodular behavior caused by the contact impact of delaminated interfaces. Composite plates with delamination, subject to external loads and voltage history from surface bonded sensors, are investigated and the results are compared with corresponding experimental results and plates without delamination.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.60 no.11
• /
• pp.2097-2102
• /
• 2011

Discharge analysis technique for dielectric liquid was presented by using the Finite Element Analysis (FEA) under a lightning impulse incorporating two-phase flow phenomena which described gas and liquid phases in discharge space. Until now, the response of step voltage has been extensively explored, but that of lightning impulse voltage was rarely viewed in the literature. We, therefore, developed an analyzing technique for dielectric liquid in a tip-sphere electrode stressed by a high electric field. To capture the bubble phase, the Heaviside function was introduced mathematically and the material functions for the ionization, dissociation, recombination, and attachment were defined in liquid and bubble, respectively. By using this numerical setup, the molecular dissociation and ionization mechanisms were tested under low and high electric fields resulted from the lightning impulse voltage of 1.2/50 ${\mu}s$. To verify our numerical results, the velocity of electric field wave was measured and compared to the previous experimental results which can be viewed in many papers. Those results had good agreement with each other.

• Steel and Composite Structures
• /
• v.29 no.5
• /
• pp.635-645
• /
• 2018

Beam-column joints play a significant role in static and dynamic performances of reinforced concrete frame structures. This study contributes a numerical approach of topologically optimal design of carbon fiber reinforced plastics (CFRP) to retrofit existing beam-column connections with crack patterns. In recent, CFRP is used commonly in the rehabilitation and strengthening of concrete members due to the remarkable properties, such as lightweight, anti-corrosion and simplicity to execute construction. With the target to provide an optimal CFRP configuration to effectively retrofit the beam-column connection under semi-failure situation such as given cracks, extended finite element method (X-FEM) is used by combining with multi-material topology optimization (MTO) as a mechanical description approach for strong discontinuity state to mechanically model cracked structures. The well founded mathematical formulation of topology optimization problem for cracked structures by using multiple materials is described in detail in this study. In addition, moved and regularized Heaviside functions (MRHF), that have the role of a filter in multiple materials case, is also considered. The numerical example results illustrated in two cases of beam-column joints with stationary cracks verify the validity, benefit and supremacy of the proposed method.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.807-823
• /
• 2008

This paper presents a technique for estimating discharge rating curve parameters. In typical practical applications, the original non-linear rating curve is transformed into a simple linear regression model by log-transforming the measurement without examining the effect of log transformation. The model of pseudo-likelihood estimation is developed in this study to deal with heteroscedasticity of residuals in the original non-linear model. The parameters of rating curves and variance functions of errors are simultaneously estimated by the pseudo-likelihood estimation(P-LE) method. Simulated annealing, a global optimization technique, is adapted to minimize the log likelihood of the weighted residuals. The P-LE model was then applied to a hypothetical site where stage-discharge data were generated by incorporating various errors. Results of the P-LE model show reduced error values and narrower confidence intervals than those of the common log-transform linear least squares(LT-LR) model. Also, the limit of water levels for segmentation of discharge rating curve is estimated in the process of P-LE using the Heaviside function. Finally, model performance of the conventional log-transformed linear regression and the developed model, P-LE are computed and compared. After statistical simulation, the developed method is then applied to the real data sets from 5 gauge stations in the Geum River basin. It can be suggested that this developed strategy is applied to real sites to successfully determine weights taking into account error distributions from the observed discharge data.