• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.113-125
• /
• 2007

Improved numerical method to obtain the exact stiffness matrices is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric and open/closed thin-walled beam on elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column This numerical technique is accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Next polynomial expressions as trial solutions are assumed for displacement parameters corresponding to zero eigenvalues and the eigenmodes containing undetermined parameters equal to the number of zero eigenvalues are determined by invoking the identity condition. And then the exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions. In order to illustrate the accuracy and the practical usefulness of this study, the numerical solutions are compared with results obtained from the thin-walled beam and shell elements.

• Journal of the Korean Geosynthetics Society
• /
• v.9 no.4
• /
• pp.9-15
• /
• 2010

In this study, a non-open cut method using a round or square pipe which has been commonly employed in Korea was experimentally evaluated and the behavioral features resulting from the friction while the pipe is penetrated into the ground was identified through the scale model test. To that end, a test device was fabricated by type of penetration pipe, by which the surface displacement caused by surrounding friction resistance was monitored. To simulate the settlement and heaving by excavation stage, the test was conducted based on generalized friction condition and surface displacement and the result therof was compared and analyzed, considering the type of penetration pipe.

• Structural Engineering and Mechanics
• /
• v.48 no.6
• /
• pp.879-914
• /
• 2013

In part I of the article, formulation and characteristics of the several well-known structural geometrical nonlinear solution techniques were studied. In the present paper, the efficiencies and capabilities of residual load minimization, normal plane, updated normal plane, cylindrical arc length, work control, residual displacement minimization, generalized displacement control and modified normal flow will be evaluated. To achieve this goal, a comprehensive comparison of these solution methods will be performed. Due to limit page of the article, only the findings of 17 numerical problems, including 2-D and 3-D trusses, 2-D and 3-D frames, and shells, will be presented. Performance of the solution strategies will be considered by doing more than 12500 nonlinear analyses, and conclusions will be drawn based on the outcomes. Most of the mentioned structures have complex nonlinear behavior, including load limit and snap-back points. In this investigation, criteria like number of diverged and complete analyses, the ability of passing load limit and snap-back points, the total number of steps and analysis iterations, the analysis running time and divergence points will be examined. Numerical properties of each problem, like, maximum allowed iteration, divergence tolerance, maximum and minimum size of the load factor, load increment changes and the target point will be selected in such a way that comparison result to be highly reliable. Following this, capabilities and deficiencies of each solution technique will be surveyed in comparison with the other ones, and superior solution schemes will be introduced.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.463-469
• /
• 2002

The governing equation and force-displacement rotations of a beam-column element on elastic foundation we derived based on variational approach of total potential energy. An exact static and dynamic 4×4 element stiffness matrix of the beam-column element is established via a generalized lineal-eigenvalue problem by introducing 4 displacement parameters and a system of linear algebraic equations with complex matrices. The structure stiffness matrix is established by the conventional direct stiffness method. In addition the F. E. procedure is presented by using Hermitian polynomials as shape function and evaluating the corresponding elastic and geometric stiffness and the mass matrix. In order to verify the efficiency and accuracy of the beam-column element using exact dynamic stiffness matrix, buckling loads and natural frequencies are calculated for the continuous beam structures and the results are compared with F E. solutions.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.1
• /
• pp.81-92
• /
• 2008

Generally the connection of members is defined as hinge or rigid. But, real joints on structure have to be considered semi-rigid connections because this permits relative rotation for members on joints. The purpose of this study is to derive a generalized tangential stiffness matrix of frames with semi-rigid connections and to develop a buckling analysis program. For the exact stiffness matrix, an accurate displacement field is introduced using an equilibrium equation for beam-columns under the bending and axial forces. Also, stability functions that consider sway deformation and force-displacement relations with rotational spring on ends were defined. In order to illustrate the accuracy of this study and the characteristics of semi-rigid for system buckling load, samples of angle-, portal- and 3-story frames with semi-rigid connections are presented, where the proposed approach is found to be in excellent agreement with other research results. Meanwhile, the application of codes such as Eurocode 3 and LRFD led to significant inaccuracies.

• Structural Engineering and Mechanics
• /
• v.16 no.6
• /
• pp.675-693
• /
• 2003

In order to consider the modified seismic response of framed structures in the presence of masonry infills, proper models have to be formulated. Because of the complexity of the problem, a careful definition of a diagonal pin-jointed strut, able to represent the horizontal force-interstorey displacement cyclic law of the actual infill, may be a solution. In this connection the present paper shows a generalized criterion for the determination of the ideal cross-section of the strut mentioned before. The procedure is based on the equivalence between the lateral stiffness of the actual infilled frame scheme during the conventional elastic stage of the response and the lateral stiffness of the same frame stiffened by a strut at the same stage. Unlike the usual empirical approaches available in the literature, the proposed technique involves the axial stiffness of the columns of the frame more than their flexural stiffness. Further, the influence of the bidimensional behaviour of the infill is stressed and, consequently, the dependence of the dimensions of the equivalent pin-jointed strut on the Poisson ratio of the material constituting the infill is also shown. The proposed approach is extended to the case of infills with openings, which is very common in practical applications.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.11 s.116
• /
• pp.1115-1123
• /
• 2006

Nuclear fuel rods are exposed to axial flow in a reactor, and flow-induced-vibration due to the flow usually causes damage in the fuel rods. Thus a prior knowledge about dynamic behavior of a fuel rod exposed to the flow condition should be provided. This paper shows that dynamic characteristics of a nuclear fuel rod depend on axial flow velocity. Assuming small lateral displacement, the effects of uniform axial flow are investigated. The analytic results show that axial flow generally reduces fuel rod stiffness and raises its damping in normal condition. Also, the critical axial velocities which make the fuel rod behavior unstable were found. That is, solving generalized eigenvalue equation of the fuel rod dynamic system, the eigenvalues with positive real part are detected. Based on the simulation results, on the other hand, it turns out that the ordinary axial flow in nuclear reactors does not affect to stability of a nuclear fuel rod even in the conservative condition.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.1
• /
• pp.69-81
• /
• 2015

The present paper deals with the determination of temperature, displacement and thermal stresses in a semi-infinite hollow circular disk due to internal heat generation within it. Initially the disk is kept at arbitrary temperature F(r, z). For times t > 0 heat is generated within the circular disk at a rate of g(r, z, t) $Btu/hr.ft^3$. The heat flux is applied on the inner circular boundary (r = a) and the outer circular boundary (r = b). Also, the lower surface (z = 0) is kept at temperature $Q_3(r,t)$ and the upper surface ($Z={\infty}$) is kept at zero temperature. Hollow circular disk extends in the z-direction from z = 0 to infinity. The governing heat conduction equation has been solved by using finite Hankel transform and the generalized finite Fourier transform. As a special case mathematical model is constructed for different metallic disk have been considered. The results are obtained in series form in terms of Bessel's functions. These have been computed numerically and illustrated graphically.

• Geomechanics and Engineering
• /
• v.22 no.6
• /
• pp.553-564
• /
• 2020

In this study, a simple numerical approach for a circular tunnel opening in strain-softening surrounding rock is proposed considering out-of-plane stress and seepage force based on Biot's effective stress principle. The plastic region of strain-softening surrounding rock was divided into a finite number of concentric rings, of which the thickness was determined by the internal equilibrium equation. The increments of stress and strain for each ring, starting from the elastic-plastic interface, were obtained by successively incorporating the effect of out-of-plane stress and Biot's effective stress principle. The initial value of the outmost ring was determined using equilibrium and compatibility equations. Based on the Mohr-Coulomb (M-C) and generalized Hoek-Brown (H-B) failure criteria, the stress-increment approach for solving stress, displacement, and plastic radius was improved by considering the effects of Biot's effective stress principle and the nonlinear degradation of strength and deformation parameters in plastic zone incorporating out-of-plane stress. The correctness of the proposed approach is validated by numerical simulation.

• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.