• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.2 s.173
• /
• pp.329-337
• /
• 2000

An investigation has been performed to develop an analysis tool based on a nonlinear beam theory, which can be used to predict the long-term behavior of an artificial hip joint. The nonlinear behav ior of the femur arise from the coupled dependence of the bone density and the mechanical properties on each other. The beam theory together with its numerical algorithm is developed to take into account the nonlinear bone remodeling process of the femur that is long enough to be assumed as a beam. A piecewise linear curve for the bone remodeling rate is used in the bone remodeling theory and the surface area density of bone is modeled as the third order polynomial function of bone density. At each section of the beam, a constant curvature is assumed and the longitudinal strains are also assumed to vary linearly across the section. The Newton-Rhapson iteration method is used to solve the nonlinear equations for each cross section of the bone and a backward method is used to march along the time. The density and the remodeling signal ar, calculated along with time for the various time steps, and the developed beam theory has been verified by comparing with the results of finite element analysis of a remodeling bone with an artificial hip joint of titanium prosthesis subjected to uni-axial loads and pure bending moment. It is concluded that the developed beam theory can be used to predict the long-term behavior of the femur and thus to design the artificial hip prosthesis.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.642-651
• /
• 1988

The present paper develops finite element procedures to calculate displacements, strains and stresses in initially stressed elastic solids subjected to static or time-dependent loading conditions. As a point of departure, we employ Hamilton's principle to obtain nonlinear equations of motion characterizing the displacement in a solid. The equations of motion reduce to linear equations of motion if incremental stresses are assumed to be infinitesimal. In the case of linear problem, finite element solutions are obtained by Newmark's direct integration method and by modal analysis. An analytic solution is referred to compare with the linear finite element solution. In the case of nonlinear problem, finite element solutions are obtained by Newton-Raphson iteration method and compared with the linear solution. Finally, the effect of the order of Gauss-Legendre numerical integration on the nonlinear finite element solution, has been investigated.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.4
• /
• pp.163-172
• /
• 2000

The postbuckling behavior and progressive damage of composite laminated cylindrical shell under uniform external pressure were investigated by nonlinear finite element method programming. For the finite element analysis, nine-node 3-D degenerated elements were utilized, and arc-length method including line search was adopted for the iteration and load-increment along postbuckling equilibrium path. As results. buckling load, postbucking behavior, and progressive failure f3r various composite laminated cylindrical shells were discussed.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.955-971
• /
• 2014

Large post-buckling behavior of Timoshenko beams subjected to non-follower axial compression loads are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. Two types of support conditions for the beams are considered. In the case of beams subjected to compression loads, load rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of lower-Carbon Steel. In the study, the relationships between deflections, rotational angles, critical buckling loads, post-buckling configuration, Cauchy stress of the beams and load rising are illustrated in detail in post-buckling case.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1990.10a
• /
• pp.16-21
• /
• 1990

Recently, the finite element method has been sucessfully extended to treat the rather couplet phenomena such as nonlinear buckling problems which are of considerable practical interest. In this study, a finite element program to evaluate the elasto-plastic buckling stress is developed. The Stowell's deformation theory for the plastic buckling of flat plates, which is in good agreement with experimental results, is used to evaluate bending stiffness matrix. A bifurcation analysis is performed to compute the elasto-plastic buckling stress. The subspace iteration method is employed to find the eigenvalues. The results are compared with corresponding enact solutions to the governing equations presented by Stowell and also with experimental data due to Pride. The developed program Is applied to obtain elastic and elasto-plastic buckling stresses for various loafing cases. The effect of different plate aspect ratio is also investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Korea Safety Management & Science
• /
• v.13 no.2
• /
• pp.97-102
• /
• 2011

본 연구에서는 원형 박판 구조물의 안정성에 대하여 해석 하였다. 임계하중은 하중을 점차적으로 증가하여 구조물이 파괴가 발생하여 안정성을 상실 하는 상태에서 가장 작은 하중을 의미한다. 판구조의 안정성을 임계하중의 크기로 기초를 두고 해석 하였다. 원형 박판구조의 차분해석은 일반 판구조와 같으므로 최근에 많은 연구의 대상이 되어왔다. 차분법은 복잡한 구조물에서도 물론, 다양한 경계조건을 포함하는 문제에 이르기까지 효과적인 수치방법이다. 본 연구에서는 기본 박판구조의 지배방정식을 유도하고 차분화 하여 직접적으로 접근하였다. 원 둘레 의 지점이 힌지 받침으로, 등분포 하중을 받고 있는 박판을 기하학적 비선형 해석으로 수행하여 원형 박판의 처짐 및 응력을 해석 하였다.

• Communications for Statistical Applications and Methods
• /
• v.12 no.2
• /
• pp.285-294
• /
• 2005

This paper deals with the comparison of parameter estimation methods in a 3-parameter Kappa distribution which is sometimes used in flood frequency analysis. Method of moment estimation(MME), L-moment estimation(L-ME), and maximum likelihood estimation(MLE) are applied to estimate three parameters. The performance of these methods are compared by Monte-carlo simulations. Especially for computing MME and L-ME, three dimensional nonlinear equations are simplified to one dimensional equation which is calculated by the Newton-Raphson iteration under constraint. Based on the criterion of the mean squared error, L-ME (or MME) is recommended to use for small sample size( n$\le$100) while MLE is good for large sample size.

• Journal of the Society of Disaster Information
• /
• v.17 no.4
• /
• pp.839-847
• /
• 2021

Purpose: This study investigates mechanical behavior of functionally graded (FG) carbon nanotube-reinforced composite (CNTRC) plate in flexure. Isogeometric analysis (IGA) method coupled with shear deformable theory of higher-order (HSDT) to analyze the nonlinear bending response is presented. Method: Shear deformable plate theory into which a polynomial shear shape function and the von Karman type geometric nonlinearity are incorporated is used to derive the nonlinear equations of equilibrium for FG-CNTRC plate in bending. The modified Newton-Raphson iteration is adopted to solve the system equations. Result: The dispersion pattern of carbon nanotubes, plate geometric parameter and boundary condition have significant effects on the nonlinear flexural behavior of FG-CNTRC plate. Conclusion: The proposed IGA method coupled with the HSDT can successfully predict the flexural behavior of FG-CNTRC plate.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.