• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.615-626
• /
• 2010

The discrete-layer model is proposed to analyze the edge-effect problem of laminates under extension and flexure. Based on three-dimensional elasticity theory, the displacement fields of each layer in a laminate have been treated discretely in terms of three displacement components across the thickness. The displacement fields at bottom and top surfaces within a layer are approximated by two-dimensional shape functions. Then two surfaces are connected by one-dimensional high order shape functions. Thus the p-convergent refinement on approximated one- and two-dimensional shape functions can be implemented independently of each other. The quality of present model is mostly determined by polynomial degrees of shape functions for given displacement fields. For nodal modes with physical meaning, the linear Lagrangian polynomials are considered. Additional modes without physical meaning, which are created by increasing nodeless degrees of shape functions, are derived from integrals of Legendre polynomials which have an orthogonality property. Also, it is assumed that mapping functions are linear in the light of shape of laminated plates. The results obtained by this proposed model are compared with those available in literatures. Especially, three-dimensional out-of-plane stresses in the interior and near the free edges are evaluated and convergence performance of the present model is established with the stress results.

• Journal of Biomedical Engineering Research
• /
• v.14 no.2
• /
• pp.155-162
• /
• 1993

The human's biomechanical structure keeps an optimal state by adapting the original biomechanical structure according to a change in the physical environment. This phenomenon is believed to be the main cause of loosening of the total hip replacement which is used widely in these days. In this study the bone density change due to artificial hip joint, which is generally believed as bone-remodeling, was investigated by the finite element method. For this, 2-D FEM models with 4 nodal point elements were constructed for intact and implanted cases. The density was calculated by comparing the relative amounts of effective stress for these two cases. In this way, calculated new density values were used in the next step as input values and this procedure repeated until convergence was obtained. Severe density change was detected at the femoral cortex of the proximal-medial side as expected. Moreover, following surprising result was found from this analysis. Titanium alloy prosthesis showed less density change compared to stainless steel prosthesis at earlier stage, however, almost same amount of the density change was detected at final stage. It was also found that other design parameters could not significantly affect its density change.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.10a
• /
• pp.133-138
• /
• 1991

The aim of the present paper is to develop a computer program predicting ultimate fracture strength of initially cracked structure under monotonically increasing external loads. For this purpose, two kinds of 3-D isoparametric solid elements, one 6-node wedge element and another 8-node brick element are formulated along the small deformation theory. Plasticity in the element is checked using von Mises' yield criterion. Elasto-plastic stiffness matrix of the element is calculated taking account of strain hardening effect. If the principal strain at crack tip which is one nodal point exceeds the critical strain dependin on the material property, crack tip is supposed to be opened and the crack tip node which was previously constrained in the direction perpendicular to the crack line is released. After that, the crack lay be propagated to the adjacent node. Once a crack tip node is fractured, the energy of the newly fractured node should be released which is to be absorbed by the remaining part. The accumulated reaction force which was carried by the newly fractured node so far is then applied in the opposite direction. During the action of crack tip relief force, since unloading may be occured in the plastic element, unloading check should be made. If a plastic element unloads, elastic stress-strain equation is used in the calculation of the stiffness matrix of the element, while for a loading element, elasto-plastic stress-strain equation is continuously used. Verification of the computer program is made comparing with the experimental results for center cracked panel subjected to uniform tensile load. Also some factors affecting ultimate fracture strength of initially cracked plate are investigated. It is concluded that the computer program developed here gives an accurate solution and becomes useful tool for predicting ultimate fracture load of initially cracked structural system under monotonically increasing external loads.

• Composites Research
• /
• v.19 no.4
• /
• pp.31-38
• /
• 2006

In this study, effects of curing temperature and autofrettage pressure on a Type 3 cryogenic propellant tank, which is composed of composite hoop/helical layers and a metal liner, were investigated by thermo elastic analysis and composite/aluminum ring specimen tests. Temperature field of a Type 3 tank was obtained from solving the heat transfer problem and, in turn, was used as nodal temperature boundary conditions during the elastic analyses for curing temperature and autofrettage pressure effects. As a result, it was shown that the higher curing temperature was, the more residual compressive stress and tensile stress were induced in composites and metal liner, respectively. On the contrary, autofrettage pressure brought the reduction of these residual thermal stresses caused by cryogenic environments to the tank structure. This tradeoff for curing temperature and autofrettage pressure must be considered in the design and manufacturing stages for a Type 3 cryogenic tank.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.95-109
• /
• 1993

In this study, the optimal configuration of arch structure has been tested by a decomposition technique. The object of this study is to provide the method of optimizing the shapes of both two hinged and fixed arches. The problem of optimal configuration of arch structures includes the interaction formulas, the working stress, and the buckling stress constraints on the assumption that arch ribs can be approximated by a finite number of straight members. On the first level, buckling loads are calculated from the relation of the stiffness matrix and the geometric stiffness matrix by using Rayleigh-Ritz method, and the number of the structural analyses can be decreased by approximating member forces through sensitivity analysis using the design space approach. The objective function is formulated as the total weight of the structures, and the constraints are derived by including the working stress, the buckling stress, and the side limit. On the second level, the nodal point coordinates of the arch structures are used as design variables and the objective function has been taken as the weight function. By treating the nodal point coordinates as design variable, the problem of optimization can be reduced to unconstrained optimal design problem which is easy to solve. Numerical comparisons with results which are obtained from numerical tests for several arch structures with various shapes and constraints show that convergence rate is very fast regardless of constraint types and configuration of arch structures. And the optimal configuration or the arch structures obtained in this study is almost the identical one from other results. The total weight could be decreased by 17.7%-91.7% when an optimal configuration is accomplished.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.1
• /
• pp.9-16
• /
• 2018

When designing a reinforced concrete member using nonlinear finite element analysis results, the bending moment at the critical section should be calculated. In this paper, a bending moment calculation method using the results of reinforced concrete finite element analysis(FEA) using continuum elements is presented and the optimum element size according to the order of the displacement function of the finite element is proposed. The bending moments calculated by integrating the stresses from the FEA are compared with the bending moments calculated using the static equilibrium conditions. In the method of integrating the stress, both the stress due to the reinforcing bar and the stress of the concrete are considered. In addition, various factors affecting the accuracy of the stresses calculated by the FEA were analyzed and the influence of the displacement function and the element size was verified. If the purpose of the analysis is to roughly observe the behavior of the members, it is appropriate to use the first order displacement function and the element size should be about 25% of the section height of the analytical model. When the bending moment of a member with high accuracy is required, it is suggested that the secondary displacement function be used and the element size be 12.5%.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• The Journal of Korean Academy of Prosthodontics
• /
• v.47 no.4
• /
• pp.394-405
• /
• 2009

Statement of problem: Crestal bone loss, a common problem associated with dental implant, has been attributed to excessive bone stresses. Design of implant's transgingival (TG) part may affect the crestal bone stresses. Purpose: To investigate if concavely designed geometry at a dental implant's TG part reduces peri-implant bone stresses. Material and methods: A total of five differently configured TG parts were compared. Base model was the ITI one piece implant (Straumann, Waldenburg, Switzerland) characterized by straight TG part. Other 4 experimental models, i.e. Model-1 to Model-4, were designed to have concave TG part. Finite element analyses were carried out using an axisymmetric assumption. A vertical load of 50 N or an oblique load of 50 N acting at $30^{\circ}$ with the implant's long axis was applied. For a systematic stress comparison, a total of 19 reference points were defined on nodal points around the implant. The peak crestal bone stress acting at the intersection of implant and crestal bone was estimated using regression analysis from the stress results obtained at 5 reference points defined along the mid plane of the crestal bone. Results: Base Model with straight configuration at the transgingival part created highest stresses on the crestal bone. Stress level was reduced when concavity was imposed. The greater the concavity and the closer the concavity to the crestal bone level, the less the crestal stresses. Conclusion: The transgingival part of dental implant affect the crestal bone stress. And that concavely designed one may be used to reduce bone stress.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.4
• /
• pp.221-234
• /
• 1999

This paper deals with the application of joint element in the finite element modeling of discontinuities encountered during rock tunneling. A nodal displacement joint element was implemented in a two dimensional finite element program GEOFE2D. The applicability of the joint element for modeling of discontinuities and the numerical stability of the implemented algorithm were examined by comparing the results of reduced small scale model tests as well as commercially available FEM program. The GEOFE2D was then used to analyze a tunnel crossed by a major discontinuity for the purpose of understanding the effect of discontinuity on the tunnel behavior. In addition, a modeling technique for the junction of discontinuity and shotcrete lining was presented. The results of analysis indicated that the stress-strain field around the tunnel is significantly altered by the presence of discontinuity, and that the stresses in the shotcrete lining considerably increase at the junction of the shotcrete lining and the discontinuity. It is therefore concluded that the major discontinuities must be carefully modeled in the finite element analysis of a tunneling problem in order to obtain more reliable results close to actual tunnel behavior.