• Structural Engineering and Mechanics
• /
• 제55권2호
• /
• pp.315-334
• /
• 2015

This study's primary aim is to check the existence of a representative volume element for granular materials and determine the link between the properties (responses) of macro structures and the size of the discrete particle assembly used to represent a constitutive relation in a two-scale model. In our two-scale method the boundary value problem on the macro level was solved using finite element method, based on the Cosserat continuum; the macro stresses and modulus were obtained using a solution of discrete particle assemblies at certain element integration points. Meanwhile, discrete particle assemblies were solved using discrete element method under boundary conditions provided by the macro deformation. Our investigations focused largely on the size effects of the discrete particle assembly and the radius of the particle on macro properties, such as deformation stiffness, bearing capacity and the residual strength of the granular structure. According to the numerical results, we suggest fitting formulas linking the values of different macro properties (responses) and size of discrete particle assemblies. In addition, this study also concerns the configuration and displacement fluctuation of discrete particle assemblies on the micro level, accompanied with the evolution of bearing capacity and deformation on the macro level.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness of materials. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of hierarchical void element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2003년도 춘계학술대회논문집
• /
• pp.1038-1043
• /
• 2003

There are many industrial applications including thin-body structures such as fins. For the numerical modeling of radiation of sound from thin bodies, the conventional boundary element method (BEM) using the Helmholtz integral equation fails to yield a reliable solution. Therefore, many researchers have tried to solve the thin-body acoustic problems. In the area of the design sensitivity analysis (DSA) and optimization methods, however, there has been just a few study reported. Especially fur the thin-body acoustics, however, no further study in the DSA and optimization fields has been reported. In this research, the normal derivative integral equation is adopted as an analysis formulation in the thin-body acoustics, and then used for the sizing DSA and optimization. Since the gradient-based method is used for the optimization, it is important to have accurate gradients (design sensitivities) of the objective function and constraints with respect to the design variables. The DSA formulations are derived through chain-ruled derivatives using the finite element method (FEM) and BEM by using the direct differentiation and continuum variation concepts. The proposed approaches are implemented and validated using a numerical example.

• 한국해양공학회지
• /
• 제27권2호
• /
• pp.53-58
• /
• 2013

Recently, some fluid-structure interaction (FSI) problems involving the fluid impact loads interacting with structures, such as sloshing, slamming, green-water, etc., have been considered, especially in the ocean engineering field. The governing equations for both an elastic solid model and flow model were originally derived from similar continuum mechanics principles. In this study, an elastic model based on a particle method, the MPS method, was developed for simulating the FSI problems. The developed model was first applied to a simple cantilever deflection problem for verification. Then, the model was coupled with the fluid flow model, the PNU (Pusan National University modified)-MPS method, and applied to the numerical investigation of the coupling effects between a cantilever and a mass of water, which has variable density, free-falling to the end of the cantilever.

• Bulletin of the Korean Chemical Society
• /
• 제35권4호
• /
• pp.1029-1035
• /
• 2014

The protolytic dissociation process of hydrochloric acid (HCl) and hydrofluoric acid (HF) is studied using the B3LYP and MP2 methods with the 6-311+G(d,p) basis set in the gas phase and in aqueous solution. To study the phenomena in detail, discrete and discrete/continuum models were applied by placing water molecules in various positions around the acid. The dissociation process was studied using the thermodynamic cycle involving the structures optimized both in the gas phase and in aqueous solution and was analyzed with two key energy factors, relaxation free energy (${\Delta}G_{Rex(g)}$) and solvation free energy (${\Delta}G_s$). Based on the results, we could understand the dissociation mechanism and wish to propose the best way to study acid dissociation process using the CPCM methodology in aqueous solution.

• Wind and Structures
• /
• 제21권5호
• /
• pp.461-487
• /
• 2015

In recent years, the Graphics Processing Unit (GPU) has become a competitive computing technology in comparison with the standard Central Processing Unit (CPU) technology due to reduced unit cost, energy and computing time. This paper describes the derivation and implementation of GPU-based algorithms for the analysis of wind loading uncertainty on high-rise systems, in line with the research field of probability-based wind engineering. The study begins by presenting an application of the GPU technology to basic linear algebra problems to demonstrate advantages and limitations. Subsequently, Monte-Carlo integration and synthetic generation of wind turbulence are examined. Finally, the GPU architecture is used for the dynamic analysis of three high-rise structural systems under uncertain wind loads. In the first example the fragility analysis of a single degree-of-freedom structure is illustrated. Since fragility analysis employs sampling-based Monte Carlo simulation, it is feasible to distribute the evaluation of different random parameters among different GPU threads and to compute the results in parallel. In the second case the fragility analysis is carried out on a continuum structure, i.e., a tall building, in which double integration is required to evaluate the generalized turbulent wind load and the dynamic response in the frequency domain. The third example examines the computation of the generalized coupled wind load and response on a tall building in both along-wind and cross-wind directions. It is concluded that the GPU can perform computational tasks on average 10 times faster than the CPU.

• Earthquakes and Structures
• /
• 제3권5호
• /
• pp.675-694
• /
• 2012

The 2008 Iwate-Miyagi Nairiku earthquake ($M_w$ 6.9, $M_{jma}$ 7.2) occurred on 14 June 2008 in Japan. The amplification and asymmetric waveform of the vertical acceleration at the ground surface recorded by accelerometers at station IWTH25, situated 3 km from the source, were remarkable in two ways. First, the vertical acceleration was extremely large (PGA = 38.66 $m/s^2$ for the vertical component, PGA = 42.78 $m/s^2$ for the sum of the three components). Second, an unusual asymmetric waveform, which is too far above the zero acceleration axis, as well as large upward spikes were observed. Using a multidegree-of-freedom (MDF) system consisting of a one-dimensional continuum subjected to vertical acceleration recorded at a depth of 260 m below ground level, the present paper clarifies numerically that these singular phenomena in the surface vertical acceleration records occurred as a result of the jumping and collision of a layer in vertical motion. We herein propose a new mechanism for such jumping and collision of ground layers. The unexpected extensive landslides that occurred in the area around the epicenter are believed to have been produced by such jumping under the influence of vertical acceleration.

• 천문학회지
• /
• 제36권3호
• /
• pp.177-187
• /
• 2003

Active galactic nuclei (AGNs) are distant, powerful sources of radiation over the entire electromagnetic spectrum, from radio waves to gamma-rays. There is much evidence that they are driven by gravitational accretion of stars, dust, and gas, onto central massive black holes (MBHs) imprisoning anywhere from $\~$1 to $\~$10,000 million solar masses; such objects may naturally form in the centers of galaxies during their normal dynamical evolution. A small fraction of AGNs, of the radio-loud type (RLAGNs), are somehow able to generate powerful synchrotron-emitting structures (cores, jets, lobes) with sizes ranging from pc to Mpc. A brief summary of AGN observations and theories is given, with an emphasis on RLAGNs. Preliminary results from the imaging of 10000 extragalactic radio sources observed in the MITVLA snapshot survey, and from a new analytic theory of the time-variable power output from Kerr black hole magnetospheres, are presented. To better understand the complex physical processes within the central engines of AGNs, it is important to confront the observations with theories, from the viewpoint of analyzing the time-variable behaviours of AGNs - which have been recorded over both 'short' human ($10^0-10^9\;s$) and 'long' cosmic ($10^{13} - 10^{17}\;s$) timescales. Some key ingredients of a basic mathematical formalism are outlined, which may help in building detailed Monte-Carlo models of evolving AGN populations; such numerical calculations should be potentially important tools for useful interpretation of the large amounts of statistical data now publicly available for both AGNs and RLAGNs.

• Structural Engineering and Mechanics
• /
• 제19권5호
• /
• pp.581-602
• /
• 2005

In this paper the accumulation of local damage during the cyclic loading in reinforced high-strength concrete columns is experimentally investigated. Two identical column specimens with annular cross-section and spiral reinforcement were designed and two tests, up to failure, under the action of a constant vertical concentrated force and a time-dependent concentrated horizontal force, were carried out at the LNEC shaking tables facility. Sine type signals, controlled in amplitude, frequency and time duration were used for these experiments. The concept of local damage based on local stiffness degradation is considered in detail and illustrated by experimental results. The specimens were designed and reinforced in such a way that the accumulation of damage was predicted by dominating deformations (cracking and crushing of the concrete) while the increasing of the loading values was a dominating factor of damage. It was observed that the local damage of HSC columns has exposed their anisotropic local behaviour. The damage accumulation was slightly different from the expected in accordance with the continuum damage concept, and a partial random character was observed.

• 대한기계학회논문집A
• /
• 제30권12호
• /
• pp.1518-1525
• /
• 2006

Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.