• Earthquakes and Structures
• /
• 제24권6호
• /
• pp.393-404
• /
• 2023

This paper presents the application of hybrid-simulation-based adapter elements for the non-linear two-scale analysis of reinforced concrete frames with masonry infills under seismic-like demands. The approach provides communication and distribution of the computations carried out by two or more remote or locally distributed numerical models connected through the OpenFresco Framework. The modeling consists of a global analysis formed by macro-elements to represent frames and walls, and to reduce global degrees of freedom, portions of the structure that require advanced analysis are substituted by experimental elements and dimensional couplings acting as interfaces with their respective sub-assemblies. The local sub-assemblies are modeled by solid finite elements where the non-linear behavior of concrete matrix and masonry infill adopt a continuum damage representation and the reinforcement steel a discrete one, the conditions at interfaces between concrete and masonry are considered through a contact model. The methodology is illustrated through the analysis of a frame-wall system subjected to lateral loads comparing the results of using macro-elements, finite element model and experimental observations. Finally, to further assess and validate the methodology proposed, the paper presents the pushover analysis of two more complex structures applying both modeling scales to obtain their corresponding capacity curves.

• Earthquakes and Structures
• /
• 제9권1호
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제3권2호
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

• Coupled systems mechanics
• /
• 제13권3호
• /
• pp.231-245
• /
• 2024

In the present work, a microelogated thermoelastic model based on Lord-Shulman (1967) and Green-Lindsay (1972) theories of thermoelasticity has been constructed. The governing equations for the simulated model are converted into two-dimensional case and made dimensionless for further simplification. Laplace and Hankel transforms followed by eigen value approach has been employed to solve the problem. The use of eigen value approach hasthe advantage of finding the solution of governing equationsin matrix form notations. This approach is straight forward and convenient for numerical computation and avoids the complicate nature of the problem. The components of displacement,stress and temperature distribution are obtained in the transformed domain. Numerical inversion techniques have been used to invert the resulting quantities in the physical domain. Graphical representation of the resulting quantities for describing the effect of microelongation are presented. A special case is also deduced from the present investigation. The problem find application in many engineering problems like thick-walled pressure vesselsuch as a nuclear containment vessel, a cylindricalroller etc.

• 한국항만학회지
• /
• 제14권4호
• /
• pp.451-461
• /
• 2000

The design and maintenance of navigation channel and water facilities of an harbor which is located at the mouth of river or at the estuary area are difficult due to the complexity of estuarial water and sediment circulation. Effects of deepening navigable waterways, of changing coastline configurations, or of discharging dredged material to the open sea are necessary to be investigated and predicted in terms of water quality and possible physical changes to the coastal environment. A borad analysis of the transport mechanism in the estuary area was made in terms of sediment property, falling velocity, concentration and flow characteristics. In order to simulate the transport processes, a two-dimensional finite element model is developed, which includes erosion, transport and deposition mechanism of suspended sediments. Galerkin’s weighted residual method is used to solve the transient convection-diffusion equation. The fluid domain is subdivided into a series of triangular elements in which a quadratic approximation is made for suspended sediment concentration. Model could deal with a continuous aggregation by stipulating the settling velocity of the flocs in each element. The model provides suspended sediment concentration, bed shear stress, erosion versus deposition rate and bed profile at the given time step.

• 한국지능시스템학회:학술대회논문집
• /
• 한국퍼지및지능시스템학회 2003년도 춘계 학술대회 학술발표 논문집
• /
• pp.303-306
• /
• 2003

In radar tracking, since the sensor measures range, azimuth and elevation angle of a target, the measurement equation is nonlinear and the extended Kalman filter (EKF) is applied to nonlinear estimation. The conventional EKF has been widely used as a nonlinear filter for radar tracking, but the considerably large measurement error due to the linearization of nonlinear function in highly nonlinear situations may deteriorate the performance of the EKF To solve this problem, a fuzzy-model-based Kalman filter (FMBKF) is proposed for radar tracking. The FMBKF uses a local model approximation based on a TS fuzzy model instead of a Jacobian matrix to linearize nonlinear measurement equation. The hybrid GA and RLS method is used to identify the premise and the consequent parameters and the rule numbers of this TS fuzzy model. In two-dimensional radar tracking problem, the proposed method is compared with the conventional EKF.

• Steel and Composite Structures
• /
• 제32권5호
• /
• pp.657-670
• /
• 2019

This paper proposes a one-dimensional fiber beam element model taking account of materially non-linear behavior, benefiting the highly efficient elastic-plastic analysis of girders with shear-lag effects. Based on the displacement-based fiber beam-column element, two additional degrees of freedom (DOFs) are added into the proposed model to consider the shear-lag warping deformations of the slabs. The new finite element (FE) formulations of the tangent stiffness matrix and resisting force vector are deduced with the variational principle of the minimum potential energy. Then the proposed element is implemented in the OpenSees computational framework as a newly developed element, and the full Newton iteration method is adopted for an iterative solution. The typical materially non-linear behaviors, including the cracking and crushing of concrete, as well as the plasticity of the reinforcement and steel girder, are all considered in the model. The proposed model is applied to several test cases under elastic or plastic loading states and compared with the solutions of theoretical models, tests, and shell/solid refined FE models. The results of these comparisons indicate the accuracy and applicability of the proposed model for the analysis of both concrete box girders and steel-concrete composite girders, under either elastic or plastic states.

• 한국의학물리학회:학술대회논문집
• /
• 한국의학물리학회 2002년도 Proceedings
• /
• pp.439-441
• /
• 2002

The aim of this study is to reconstruct the 3D target volume from multiple projection images. It was assumed that we were already aware of the target position exactly, and all processes were performed in Target Coordinates whose origin was the center of the target. We used six projections: two projections were used to make a Reconstruction Box and four projections were for image acquisition. Reconstruction Box was made up of voxels of 3D matrix. Projection images were transformed into 3D volume in this virtual box using geometrical based back-projection method. Algorithm was applied to an ellipsoid model and horse-shoe shaped model. Projection images were created using C program language by geometrical method and reconstruction was also accomplished using C program language and Matlab(The Mathwork Inc., USA). For ellipsoid model, reconstructed volume was slightly overestimated but target shape and position was proved to be correct. For horse-shoe shaped model, reconstructed volume was somewhat different from original target model but there was a considerable improvement in target volume determination.

• 국제학술발표논문집
• /
• The 6th International Conference on Construction Engineering and Project Management
• /
• pp.496-500
• /
• 2015

International construction projects are inherently more risky than domestic projects with multi-dimensional uncertainties that require complementary risk management at both the country and project levels. However, despite a growing need for systematic country evaluations, most studies have focused on project-level decisions and lack country-based approaches for firms in the construction industry. Accordingly, this study suggests data-driven approaches for evaluating countries using two quantitative models. The first is a two-stage country segmentation model that not only screens negative countries based on country attractiveness (macro-segmentation) but also identifies promising countries based on the level of past project performance in a given country (micro-segmentation). The second is a multi-criteria country segmentation model that combines a firm's business objective with the country evaluation process based on Kraljic's matrix and fuzzy preference relations (FPR). These models utilize not only secondary data from internationally reputable institutions but also performance data on Korean firms from 1990 to 2014 to evaluate 29 countries. The proposed approaches enable firms to enhance their decision-making capacity for evaluating and selecting countries at the early stage of corporate strategy development.

• Nuclear Engineering and Technology
• /
• 제41권8호
• /
• pp.1109-1114
• /
• 2009

The evolution of a high burnup structure (HBS) in a light water reactor (LWR) $UO_2$ fuel was simulated using the Potts model. A simulation system for the Potts model was defined as a two-dimensional triangular lattice, for which the stored energy was calculated from both the irradiation damage of the $UO_2$ matrix and the formation of a grain boundary in the newly recrystallized small HBS grains. In the simulation, the evolution probability of the HBS is calculated by the system energy difference between before and after the Monte Carlo simulation step. The simulated local threshold burnup for the HBS formation was 62 MWd/kgU, consistent with the observed threshold burnup range of 60-80 MWd/kgU. The simulation revealed that the HBS was heterogeneously nucleated on the intergranular bubbles in the proximity of the threshold burnup and then additionally on the intragranular bubbles for a burnup above 86 MWd/kgU. In addition, the simulation carried out under a condition of no bubbles indicated that the bubbles played an important role in lowering the threshold burnup for the HBS formation, thereby enabling the HBS to be observed in the burnup range of conventional high burnup fuels.