• Magazine of the Korea Concrete Institute
• /
• v.7 no.5
• /
• pp.172-178
• /
• 1995

Concrete contains numerous microcracks at initially poured. The growth and propagation of nicrockacsk are believed tc finally incur the faiure of concrete. These processings are understood as a damage. Damage IS represented as a second-order tensor and crack is treated as a con tinuum phenomenon. In this paper, damage is characterized through the effective stress concept together with the hypothesis of elastic energy equivalence, and damage evolution law and constitutive equation of a damage model are derived by using the Helmholtz frte eriergy and the dissipation potential by means of the thermodynamic principles. The constitutive equation of the model includes the effects of elasticity, anisotropic damage and plasticity of concrete. There are two effective tangent stiffness tensors in this model : one is for elastic-darnage and the other for plastic damage. For the verification of the model, finite element analysis was performed for the analysis of concrete subjec:t to uniaxial and biaxial loading and the results obtained were compared with test results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.257-270
• /
• 2000

The objective of this study is to develop a damage model based on damage mechanics that can be used to analyze the mechanical behavior of structures with defects and the global behavior of damaged structures. A modified second order damage tensor that can be applied to finite element analysis is used to reflect the effect of damage. The damage stress computed from the effective stress is considered as an additional loading term acting on nodes and can represent the effect of crack surface. The accuracy of the proposed algorithm is verified by comparing the analysis results with the experimental data from other studies and the analysis results based on transverse isotropic theory. The developed damage model is applied to the analyses of structures with cracks under linear elastic condition. The comparisons confirmed that the quantitative analysis of the structural behavior due to crack orientation and multiple sets of cracks is possible. Also, the damage caused by rock excavation and fault zone is analyzed. The results also showed that the developed model can effectively analyze the global behavior of damaged structures.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Structural Engineering and Mechanics
• /
• v.8 no.6
• /
• pp.607-622
• /
• 1999

Being a significant mode of deformation, shear effect in addition to the other modes of stretching and bending have been considered to develop two finite element models for the analysis of beams on elastic foundation. The first beam model is developed utilizing the differential-equation approach; in which the complex variables obtained from the solution of the differential equations are used as interpolation functions for the displacement field in this beam element. A single element is sufficient to exactly represent a continuous part of a beam on Winkler foundation for cases involving end-loadings, thus providing a benchmark solution to validate the other model developed. The second beam model is developed utilizing the hybrid-mixed formulation, i.e., Hellinger-Reissner variational principle; in which both displacement and stress fields for the beam as well as the foundation are approxmated separately in order to eliminate the well-known phenomenon of shear locking, as well as the newly-identified problem of "foundation-locking" that can arise in cases involving foundations with extreme rigidities. This latter model is versatile and indented for utilization in general applications; i.e., for thin-thick beams, general loadings, and a wide variation of the underlying foundation rigidity with respect to beam stiffness. A set of numerical examples are given to demonstrate and assess the performance of the developed beam models in practical applications involving shear deformation effect.

• Journal of Biomedical Engineering Research
• /
• v.15 no.2
• /
• pp.217-224
• /
• 1994

In order to study the shape and dimensions of heart, a procedure to reconstruct a three dimensional left ventricular geometry from two dimensional echocardiographic images was studied including the coordinate transformation, curve fitting and interpolation utilizing three dimensional position registration arm. Nonlinear material property of the left ventricular myocardium was obtained by finite element method performed on the reconstructed geometry and by optimization techniques which compared the computer predicted 3D deformation with the experimentally determined deformation. Elastic modulus ranged from 3.5g/$cm^2$ at early diastole to l53g/$cm^2$ at around end diastole showing slightly nonlinear relationship between the modulus and the pressure. Afterwards using the obtained nonlinear material propertry the stress distribution related with oxyzen consumption rate was analyzed. The maximum and minimum of ${\sigma}_1$ (max. principal stress) occurred at nodes on the second level intersection points of x-axis with endocardium and with epicardium, respectively. And the tendency of the interventricular septum to be flattened was observed from the compressive ${\sigma}_1$ on the anterior, posterior nodes of left ventricle and from the most significant change of dimension in $D_{RL}$ (septal-lateral dimension of right ventricle).

• Coupled systems mechanics
• /
• v.8 no.4
• /
• pp.351-375
• /
• 2019

The present work investigates the use of a rigid rocking block as a tool to reduce vibrations in a frame structure. The study is based on a simplified model composed by a 2-DOF linear system, meant to represent a general M-DOF frame structure, coupled with a rocking rigid block through a linear visco-elastic device, which connects only the lower part of the 2-DOF system. The possibility to restrain the block directly to the ground, by means of a second visco-elastic device, is investigated as well. The dynamic response of the model under an harmonic base excitation is then analysed in order to evaluate the effectiveness of the coupling in reducing the displacements and the drift of the 2-DOF system. The nonlinear equations of motion of the coupled assemblage 2-DOF-block are obtained by a Lagrangian approach and then numerically integrated considering some reference mechanical and geometrical quantities as variable parameters. It follows an extensive parametric analysis, whose results are summarized through behaviour maps, which portray the ratio between the maximum displacements and drifts of the system, with and without the coupling with the rigid block, for several combinations of system's parameters. When the ratio of the displacements is less than unity, the coupling is considered effective. Results show that the presence of the rocking rigid block improves the dynamics of the system in large ranges of the characterizing parameters.

• Geotechnical Engineering
• /
• v.4 no.3
• /
• pp.19-26
• /
• 1988

A stochastic numerical model for predictions of differential settlement of foundation Eoils is developed in this Paper. The differential settlement is highly dependent on the spatial variability of elastic modulus of soil. The Kriging method is used to account for the spatial variability of the elastic modulus. This technique provides the best linear unbiased estimator of a parameter and its minimum variance from a limited number of measured data. The stochastic finite element method, employing the first-order second-moment analysis for computations of error Propagation, is used to obtain the means, ariances, and covariances of nodal displacements. Finally, a reliability model of differential settlement is proposed by using the results of the stochastic FEM analysis. It is found that maximum differential settlement occurs when the distance between two foundations is approximately same It with the scale of fluctuation in horizontal direction, and the probability that differential settlement exceeds the allot.able vague might be significant.

• Computers and Concrete
• /
• v.14 no.1
• /
• pp.41-57
• /
• 2014

Modern trend in structural design is to use smaller elements in order to ensure several purposes such as economy, functionality and aesthetic in appearance. However, because of decreasing rigidity of the structural elements, the system displacements increases and displacements become an important subject in this kind of structures takes into account both geometrical changes and the carrying capacity of the material after linear-elastic boundary. In this study, a method is proposed to calculate the failure loads and to analyse the reinforced concrete space frame systems. The numerical examples gathered from the literature survey are solved with this method utilising the prepared computer program and the comparable results are presented. The results show that the method is sufficiently accurate.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.239-246
• /
• 1997

The general formulation of free vibration and stability analysis of unsymmetric thin-walled space frames and beam-columns is presented. The kinetic and total potential energy is derived by applying the extended virtual work principle, introducing displacement parameters defined at the arbitrarily chosen axis and including second order terms of finite semitangential rotations. In formulating the finite element procedure, cubic Hermitian polynomials are utilized as shape functions of the two node space frame element. Mass, elastic stiffness, and geometric stiffness matrices for the unsymmetric thin-walled section are evaluated. In order to illustrate the accuracy and practical usefulness of this formulation, finite element solutions for the free vibration and stability problems of thin-walled beam-columns and space frames are presented and compared with available solutions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.