• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.1
• /
• pp.1-8
• /
• 1992

A method analizing contact pressure between plate and elastic half space is presented by using F.E.M. With the method, the pressure intensities at surface nodes of half space cae be directly calculated by using flexibility matrix of half space. The method is originally presented by Y.K. Cheung et al.(3) Insted of Y.K. Cheung's method, which use a conception of equi-contact pressure area around each surface nodes of half space in the noded rectanqular element area. We use the equi-contact pressure area around the Gaussian integration points of half space surface in the noded isoparametric element area. Numarical examples are presented and compared with other's studies.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.99-107
• /
• 1993

In this study, material nonlinear finite element program is developed to analyze reinforced concrete shell of arbitrary geometry considering tension stiffening effects. This study is capable of tracing the load-deformation response and crack propagation, as well as determining the internal concrete and steel stresses through the elastic, inelastic and ultimate ranges in one continuous computer analysis. The cracked shear retention factor is introduced to estimate the effective shear modulus including aggregate interlock and dowel action. The concrete is assumed to be brittle in tension and elasto-plastic in compression. The Drucker-Prager yield criterion and the associated flow rule are adopted to govern the plastic behavior of the concrete. The reinforcing bars are considered as a steel layer of equivalent thickness. A layered isoparametric flat finite element considering the coupling effect between the in-plane and the bending action was developed. Mindlin plate theory taking account of transverse shear deformation was used. An incremental tangential stiffness method is used to obtain a numerical solution. Numerical examples about reinforced concrete shell are presented. Validity of this method is studied by comparing with the experimential results of Hedgren and the numerical analysis of Lin.

• Restorative Dentistry and Endodontics
• /
• v.19 no.2
• /
• pp.384-408
• /
• 1994

The purpose of this study was to evaluate the fracture resistance of tooth restored with gold inlay. A profound understanding of the isthmus width factor, which is one of the several parameters of cavity designs, would facilitate the appropriate cavity preparation in a specific clinical situation. In this study, the cavities for gold inlay were prepared in maxillary left first premolar. A three-dimensional model was designed using I-DEAS program. The model was composed of 2515- nodes and 2172 isoparametric brick elements. In the model isthmus width was varied into 1/4, 1/3 and 1/2 of intercuspal width respectively, and numeric values of the material properties of enamel, dentin and gold was set. Three types of load : concentrated load, divided load and distributed load was 500N. The empty cavities in the model were also examined using divided load and distributed load. The three - dimensional Finite Element Method was used to analysis the displacement and stress distribution. The results were as follows : 1. All of the experimental models which were filled with gold inlay revealed similar direction of displacement to that of the natural tooth model under the same load type. But in the models with empty cavities, as the isthmus width increased, the degree of displacement increased in the case of divided load type. 2. All experimental models which were filled with gold inlay showed stress concentration at load points, but in the models with empty cavities at divided load type, as isthmus width increased, stress was concentrated at the comer of the pulpal floor. 3. In the models with empty cavities at divided load type, tooth fracture was expected regardless of isthmus width, but all experimental models which were filled with gold inlay after cavity preparation were not susceptible to fracture. 4. In all experimental models which were filled with gold inlay after cavity preparation, displacement patterns were similar under both concentrated and divided load types. In the models with empty cavities, a divided load resulted in a bucco-lingual cuspal displacemenat in both sides, but a distributed load resulted in a lingual displacement of the tooth.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.619-628
• /
• 2002

This paper deals with the FE analysis for the free vibration of sloshing in horizontal cylindrical tank with baffles. We use Laplace equation based on potential theory as governing equation. This problem is solved by FEM using lineal isoparametric elements. We assume that the tank as well as baffles is rigid body and by separating nodes into two at the baffle location, baffle effect is obtained by separating nodes into two at the baffle location. For the calculation of natural frequencies and mode shapes, we introduce Lanczos transformation and Jacobi iteration methods. Numerical results of the first longitudinal and transverse modes, while comparing with literature cited, are very good. In order for the baffle effects on the free vibration of sloshing, various combinations of baffle parameters, which are location, inner diameter and number, are examined.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.82-91
• /
• 1996

This is a basic study for the static and dynamic analysis on the elasto-plastic and elasto-viscoplastic of an axi-symmetric shell. The objective of this study was to investigate the mechanical characteristics of a nuclear reinforced concrete containment structure, which was selected as a model, by a numerical analysis using a finite element method. The structure was modeled with discrete ring elements of 8-noded isoparametric element rotating against the symmetrical axis, and the interaction between the foundation and the structure was modeled by Winkler's model. Also, the meridional tendon was modeled with 2-node truss elements, and the hoop tendon was done with point elements in two degrees of freedom. The effect of the tendon was considered without the increasement in total degree of freedom as the stiffness matrix of modeled tendon elements was assembled on the stiffness matrix of ring elements linked with the tendon. The results obtained from the analysis of an example were summarized as follows : 1. The stresses in the hoop direction on the interior and exterior surfaces of the structure were shown in changes of similar trend, and high stresses appeared on the structure wall 2. The stresses in the meridional direction on the interior and exterior surfaces were shown in change of different trend. Especially, the stresses at the junctions between the dome and the wall and between the wall and the bottom plate of the structure were very high, compared with those at other parts of the structure. 3. The stress changes in the direction of thickness on the crown of the dome were much linearly distributed. However, as the amount of tendon increased, the stresses in the upper and lower parts of the wall established with the tendon were shown stress concentration. 4. The stress changes in the direction of thickness on the center of the structure wall was linearly distributed in the all cases, and special stress due to the use of the tendon was not shown.

• The Journal of Korean Academy of Prosthodontics
• /
• v.29 no.1
• /
• pp.167-213
• /
• 1991

The successful replacement of missing teeth has been one driving aim behind the emergence of implant dentistry as both a technology and clinical vocation for over four decades. To date, a multitude of dental implant devices had been designed and utilized in the patient population. Most of these devices have been designed without support of the engineering criteria. The long-term success of any dental implant is dependent upon the optimization of stresses which occurs during oral function and parafunction. Although many studies have examined the biologic interactions between dental implants and living tissue, few studies have been reported on the biomechanical aspects of dental implants. The purpose of this study was to analyze the stress distribution of osseointegrated prosthesis on certain conditions, such as amount of load, location of load, length of fixtures, number of fixtures used, arch shape, bone quality, etc. Three dimentional finite element analysis was used for this study. FEM models were created using commercial software(Super SAP. for IBM 16 bit AT computer. All elements were 8-node brick, isoparametric. Mandible and prosthesis was modeled with 780 elements and 1074 nodes. The results were as follows : 1. In case of cantilever extension, there was a compressive stress at the base of the first implant and a tensile stress at the base of the second implant. 2. The stresses were linearly proportional to the amount of load. 3. The stresses were linearly proportional to the length of cantilever. 4. There was a stress concentration at the neck of the implant and bone under horizontal loads.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• Computers and Concrete
• /
• v.11 no.2
• /
• pp.123-148
• /
• 2013

The nonlinear finite element method with eight noded isoparametric quadrilateral element for concrete and two noded element for reinforcement is used for the prediction of the behavior of reinforcement concrete structures. The disturbed state concept (DSC) including the hierarchical single surface (HISS) plasticity model with associated flow rule with modifications is used to characterize the constitutive behavior of concrete both in compression and in tension which is named DSC/HISS-CT. The HISS model is applied to shows the plastic behavior of concrete, and DSC for microcracking, fracture and softening simulations of concrete. It should be noted that the DSC expresses the behavior of a material element as a mixture of two interacting components and can include both softening and stiffening, while the classical damage approach assumes that cracks (damage) induced in a material treated acts as a void, with no strength. The DSC/HISS-CT is a unified model with different mechanism, which expresses the observed behavior in terms of interacting behavior of components; thus the mechanism in the DSC is much different than that of the damage model, which is based on physical cracks which has no strength and interaction with the undamaged part. This is the first time the DSC/HISS-CT model, with the capacity to account for both compression and tension yields, is applied for concrete materials. The DSC model allows also for the characterization of non-associative behavior through the use of disturbance. Elastic perfectly plastic behavior is assumed for modeling of steel reinforcement. The DSC model is validated at two levels: (1) specimen and (2) practical boundary value problem. For the specimen level, the predictions are obtained by the integration of the incremental constitutive relations. The FE procedure with DSC/HISS-CT model is used to obtain predictions for practical boundary value problems. Based on the comparisons between DSC/HISS-CT predictions, test data and ANSYS software predictions, it is found that the model provides highly satisfactory predictions. The model allows computation of microcracking during deformation leading to the fracture and failure; in the model, the critical disturbance, Dc, identifies fracture and failure.

• Restorative Dentistry and Endodontics
• /
• v.19 no.2
• /
• pp.345-371
• /
• 1994

Restorative procedures can lead to weakening tooth due to reduction and alteraton of tooth structure. It is essential to prevent fractures to conserve tooth. Among the several parameters in cavity designs, cavity isthmus and depth are very important. In this study, MO amalgam cavity was prepared on maxillary first premolar. Three dimensional. finite element models were made by serial photographic method and cavity depth(1.7mm, 2.4mm) and isthmus (11 4, 1/3, 1/2 of intercuspal distance) were varied. linear, eight and six-nodal, isoparametric brick elements were used for the three dimensional finite element model. The periodontal ligament and alveolar bone surrounding the tooth were excluded in these models. Three types model(B, G and R model) were developed. B model was assumed perfect bonding between the restoration and cavity wall. Both compressive and tensile forces were distributed directly to the adjacent regions. G model(Gap Distance: 0.000001mm) was assumed the possibility of play at the interface simulated the lack of real bonding between the amalgam and cavity wall (enamel and dentin). When compression occurred along the interface, the forces were transferred to the adjacent regions. However, tensile forces perpendicular to the interface were excluded. R model was assumed non-connection between the restoration and cavity wall. No force was transferred to the adjacent regions. A load of 500N was applied vertically at the first node from the lingual slope of the buccal cusp tip. This study analysed the displacement, von Mises stress, 1 and 2 direction normal stress and strain with FEM software ABAQUS Version 5.2 and hardware IRIS 4D/310 VGX Work-station. The results were as follows: 1. G model showed stress and strain patterns between Band R model. 2. B model and G model showed the bending phenomenon in the displacement. 3. R model showed the greatest amount of the displacement of the buccal cusp followed by G and B model in descending order. G model showed the greatest amount of the displacement of the lingual cusp followed by B and R model in descending order. 4. B model showed no change of the displacement as increasing depth and width of the cavity. G and R model showed greater displacement of the buccal cusp as increasing depth and width of the cavity, but no change in the displacement of the lingual cusp. 5. As increasing of the width of the cavity, stress and strain were not changed in B model. Stress and strain were increased on the distal marginal ridge and buccopulpal line angle in G and R model. The possibility of the tooth fracture was increased. 6. As increasing of the depth of the cavity, stress and strain were not changed in B and G model. Stress and strain were increased on the distal marginal ridge and buccopulpal line angle in R model. The possibility of the tooth fracture was increased.

• Restorative Dentistry and Endodontics
• /
• v.20 no.2
• /
• pp.432-461
• /
• 1995

The basic principles in the design of Class II amalgam cavity preparations have been modified but not changed in essence over the last 90 years. The early essential principle was "extension for prevention". Most of the modifications have served to reduce the extent of preparation and, thus, increase the conservation of sound tooth structure. A more recent concept relating to conservative Class II cavity preparations involves elimination of occlusal preparation if no carious lesion exists in this area. To evaluate the ideal ClassII cavity preparation design, if carious lesion exists only in the interproximal area, three cavity design conditions were studied: Rodda's conventional cavity, simple proximal box cavity and proximal box cavity with retention grooves. In this study, MO amalgam cavity was prepared on maxillary first premolar. Three dimensional finite element models were made by serial photographic method. Linear, eight and six-nodal, isoparametric brick elements were used for the three dimensional finite element model. The periodontal ligament and alveolar bone surrounding the tooth were excluded in these models. Three types model(B option, Gap option and R option model) were developed. B option model was assumed perfect bonding between the restoration and cavty wall. Gap option model(Gap distance: $2{\mu}m$) was assumed the possibility of play at the interface simulated the lack of real bonding between the amalgam and cavity wall (enamel and dentin). R option model was assumed non-connection between the restoration and cavty wall. A load of 500N was applied vertically at the first node from the lingual slope of the buccal cusp tip. This study analysed the displacement, 1 and 2 direction normal stress and strain with FEM software ABAQUS Version 5.2 and hardware IRIS 4D/310 VGX Work-station. The results were as followed. 1. Rodda's cavity form model showed greater amount of displacement with other two models. 2. The stress and strain were increased on the distal marginal ridge and buccopulpal line angle in Rodda's cavity form model. 3. The stress and strain were increased on the central groove and a part of distal marginal ridge in simple proximal box model and proximal box model with retention grooves. 4. With Gap option, Rodda's cavity form model showed the greatest amount of the stress on distal marginal ridge followed by proximal box model with retention grooves and simple proximal box model in descending order. 5. With Gap option, simple proximal box model showed greater amount of stress on the central groove with proximal box model with retention grooves. 6. Retention grooves in the proximal box played the role of supporting the restorations opposing to loads.