• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.737-752
• /
• 1992

Analysis of elasto-dynamic deformation of flexible linkage mechanism is conducted using the finite element method. The equations of motion of the system are derived from the static structural problem in which dynamic inertia, gravitational and driving forces are treated as external loads. Linear spring model is included in the formulation of equation of motions to represent the effects of deformation of elastic bearings of revolute joints on the system behavior. A computer program is constructed and applied to analyze a specific crank-lever 4-bar mechanism. The algorithm of the program is as follows. First, the natural frequencies and the mode shapes of the system are calculated by solving the eigenproblem of the mechanism system which can be considered as a static structure by assuming the input shaft (crank shaft) to be fixed at any given configuration of mechanism. And finally, the elasto-dynamic deformation of the whole system is obtained using mode superposition method for the case of constant input speed. The effect of geometric stiffness on the mechamism is included in the program with the axial forces of links obtained through the quasi-static displacement analysis. It is found that the geometric stiffness exerts an important effect upon the elasto-dynamic behavior of the flexible linkage mechanism. Elastic deformation of bearing lowers the natural frequencies of the system, resulting smaller elastic displacement at the mid-point of the links and bigger elestic displacement at the ends of the links than rigid bearing. The above investigation of flexible linkage mechanism shows that the effects of the elastic deformation of bearing on the mechanism should be considered to design the mechanism which satisfies more preciously the purpose and the condition of design.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.463-474
• /
• 2009

Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

• Journal of the Korean GEO-environmental Society
• /
• v.24 no.9
• /
• pp.15-24
• /
• 2023

The behavior of buried pipes is directly influenced by the nonlinearity and complex characteristics of the surrounding soil. However, the simplified beam-spring model, which ignores the nonlinearity and complex behavior of soil, is commonly used in practice. In response, several studies have employed continuum analysis methods to account for the nonlinear and complex behavior of the soil. This paper presents various numerical continuum analysis techniques and verifies their comparison with full-scale tests. The study found that reaction force results close to the full-scale test could be obtained by applying contact surface characteristics that take into account the interaction between the ground and the buried pipe. In the case of sharing pipe and soil node method and ignoring the interaction between pipe and soil, excessive reaction force was derived, and the failure shapes were different. In addition, this study applied the dynamic explicit analysis method, ALE method, and CEL method. It was confirmed that the displacement-reaction relationship and failure shape are similar to those of the static analysis.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.3
• /
• pp.39-55
• /
• 1992

In this research, a Three Level Decomposition technique has been developed for configuration design optimization of truss structures. In the first level, as design variables, behavior variables are used and the strain energy has been treated as the cost function to be maximized so that the truss structure can absorb maximum energy. For design constraint of the optimal design problem, allowable stress, buckling stress, and displacement under multi-loading conditions are considered. In the second level, design problem is formulated using the cross-sectional area as the design variable and the weight of the truss structure as the cost function. As for the design constraint, the equilibrium equation with the optimal displacement obtained in the first level is used. In the third level, the nodal point coordinates of the truss structure are used as coordinating variable and the weight has been taken as the cost function. An advantage of the Three Level Decomposition technique is that the first and second level design problems are simple because they are linear programming problems. Moreover, the method is efficient because it is not necessary to carry out time consuming structural analysis and techniques for sensitivity analysis during the design optimization process. By treating the nodal point coordinates as design variables, the third level becomes unconstrained optimal design problems which is easier to solve. Moreover, by using different convergence criteria at each level of design problem, improved convergence can be obtained. The proposed technique has been tested using four different truss structures to yield almost identical optimum designs in the literature with efficient convergence rate regardless of constraint types and configuration of truss structures.

• The korean journal of orthodontics
• /
• v.31 no.4 s.87
• /
• pp.425-438
• /
• 2001

The delivery of optimal orthodontic treatment is greatly influenced by clinician's ability to predict and control tooth movement by applying well-known force system to dentition. It is very important to determine the location of the centers of resistance of a tooth or teeth in order to have better understanding the nature of displacement characteristics under various force levels. In this study, three dimensional finite element analysis was used to measure the initial displacement of the consolidated teeth under loading. The purpose of this study was to define the location of the centers of resistance at the upper six anterior segment. To observe the changes of six anterior segment, 200gm, 250gm, 300gm, and 350gm forces at right and left hand side each were imposed toward lingual direction. For this study, two cases, six anterior teeth and six anterior teeth after corticotomy, were reviewed. In addition, it was reviewed the effects of changes on the location of the center of resistance in both cases based on different degree of forces aforementioned. The results were that : 1. The instantaneous center of resistance for the six anterior teeth was vertically located between level 4 and level 5, which is, at 6.76mm, $44.32\%$ apical to the cementoenamel junction level. 2. The instantaneous center of resistance for the six anterior teeth after corticotomy was located vertically between level 4 and level 5, that is, at 7.09mm $46.38\%$ apical to the cementoenamel junction level. 3. Changes of force showed little effect on the location of the center of resistance in each case. 4. It was observed that the location of the instantaneous center of resistance for the six anterior teeth after corticotomy was changed more than the six anterior teeth without corticotomy to the apical part, and the displacement of the consolidated anterior teeth moved further in case of the consolidated teeth after corticotomy.

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

In this research, configuration design optimization of plane truss structure has been tested by using decomposition technique. In the first level, the problem of transferring the nonlinear programming problem to linear programming problem has been effectively solved and the number of the structural analysis necessary for doing the sensitivity analysis can be decreased by developing stress constraint into member stress approximation according to the design space approach which has been proved to be efficient to the sensitivity analysis. And the weight function has been adopted as cost function in order to minimize structures. For the design constraint, allowable stress, buckling stress, displacement constraint under multi-condition and upper and lower constraints of the design variable are considered. In the second level, the nodal point coordinates of the truss structure are used as coordinating variable and the objective function has been taken as the weight function. By treating the nodal point coordinates as design variable, unconstrained optimal design problems are easy to solve. The decomposition method which optimize the section areas in the first level and optimize configuration variables in the second level was applied to the plane truss structures. The numerical comparisons with results which are obtained from numerical test for several truss structures with various shapes and any design criteria show that convergence rate is very fast regardless of constraint types and configuration of truss structures. And the optimal configuration of the truss structures obtained in this study is almost the identical one from other results. The total weight couldbe decreased by 5.4% - 15.4% when optimal configuration was accomplished, though there is some difference.

• Journal of Korean Association for Spatial Structures
• /
• v.8 no.4
• /
• pp.47-55
• /
• 2008

This paper deals with basic theories and some numerical results on structural optimization for geodesic dome. First of all, the space efficiency of geodesic dome is investigated by using the ratio of icosahedron's surface area to the internal volume enclosed by it. The procedure how to create the geodesic dome is also provided in systematic way and implemented and utilized into the design optimization code ISADO-OPT. The mathematical programming technique is introduced to find out the optimum pattern of member size of geodesic dome against a point load. In this study, total weight of structure is considered as the objective function to be minimized and the displacement occurred at loading point and member stresses of geodesic dome are used as the constraint functions. The finite difference method is used to calculate the design sensitivity of objective function with respect to design variables. The SLP, SQP and MFDM available in the optimizer DoT is used to search optimum member size patterns of geodesic dome. It is found to be that the optimum member size pattern can be efficiently obtained by using the proposed design optimization technique and numerical results can be used as benchmark test as a basic reference solution for design optimization of dome structures.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.6
• /
• pp.19-30
• /
• 2005

A coupled three-dimensional pile group analysis method was developed by considering complex behavior of sub-structures (pile-soil-cap) which included soil nonlinearity and the behavior of super-structure (pier). As an intermediate analysis method between FBPier 3.0 and Group 0.0, it took advantages of each method. Among the components of a pile group, individual piles were modeled with stiffness matrices of pile heads and soils with nonlinear load-transfer curves (t-z, q-z and p-y curves). A pile cap was modeled with modified four-node flat shell elements and a pier with three-dimensional beam element, so that a unified analysis could be possible. A nonlinear analysis method was proposed in this study with a mixed incremental and iteration techniques. The proposed method for a pile group subjected to axial and lateral loads was compared with othe. analytical methods (i.e., Group 6.0 and FBPier 3.0). It was found that the proposed method could predict the complex behavior of a pile group well, even though piles were modelled simply in this study by using pile head stiffness matrices which were different from the method introduced in FBPier 3.0.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.2 no.3
• /
• pp.1-13
• /
• 1982

Transition process of plastic region. displacements, stresses and strains ahead the flaw tips were analysed by the finite element method on the steel plate with the circular hole and the one with the elliptical hole under completely alternating load (repetition of tensile loading, unloading and compressive loading). As the results, the followings were obtained. Transition process of elastic failure (yielding) region was estimated. From this the tendency was confirmed that the fracture would be initiated from ahead the flaw tip, and propagated along the $45^{\circ}$ direction. The fundamental data available in estimating the stress intensity factor that was considered as the core in analysing the fracture mechanism of steel plates were obtained. It was indicated that when unloading after tension the effect of compressive loading, and even the compressive reyield, was occured ahead the flaw tip. Similarly it was indicated that when unloading after compression the effect of tensile loading, and even the tensile reyield, was occured ahead the flaw tip. It was considered that these phenomena were occured because the unloading effect was constrained by the residual strains when unloading. It was considered that the fatigue phenomenon was occured ahead, the flaw tip by repetition of tensile yield, the above compressive reyield, compressive yeild and the above tensile reyield. In addition, the tendency was confirmed that the fracture ahead the flaw tip was occured as the flaw was changed from the circular hole to the elliptical hole and became to be the crack lastly.

• Composites Research
• /
• v.13 no.5
• /
• pp.50-59
• /
• 2000

In this paper, two methods to suppress flutter of the composite panel are examined. First, in the active control method, a controller based on the linear optimal control theory is designed and control input voltage is applied on the actuators and a PZT is used as actuator. Second, a new technique, passive suppression scheme, is suggested for suppression of the nonlinear panel flutter. In the passive suppression scheme, a shunt circuit which consists of inductor-resistor is used to increase damping of the system and as a result the flutter can be attenuated. A passive damping technology, which is believed to be more robust suppression system in practical operation, requires very little or no electrical power and additional apparatuses such as sensor system and controller are not needed. To achieve the great actuating force/damping effect, the optimal shape and location of the actuators are determined by using genetic algorithms. The governing equations are derived by using extended Hamilton's principle. They are based on the nonlinear von Karman strain-displacement relationship for the panel structure and quasi-steady first-order piston theory for the supersonic airflow. The discretized finite element equations are obtained by using 4-node conforming plate element. A modal reduction is performed to the finite element equations in order to suppress the panel flutter effectively and nonlinear-coupled modal equations are obtained. Numerical suppression results, which are based on the reduced nonlinear modal equations, are presented in time domain by using Newmark nonlinear time integration method.