• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.3
• /
• pp.7-17
• /
• 1990

The objective of this paper is to provide a method of optimizing are as of the members as well as shape of both truss and arch structures. The design process includes satisfaction of stress and Euler buckling stress constraints for truss and combined stress constraints for arch structures. In order to reduce the number of detailed finite element analysis, the Force Approximation Method is used. A finite element analysis of the initial structure is performed and the gradients of the member end forces are calculated with respect to the areas and nodal coordinates. The gradients are used to form an approximate structural analysis based on first order Taylor series expansions of the member end forces. Using move limits, a numerical optimizer minimizes the volume of the structure with information from the approximate structural analysis. Numerical examples are performed and compared with other methods to demonstrate the efficiency and reliability of the Force Approximation Method for shape optimization. It is shown that the number of finite element analysis is greatly reduced and that it leads to a highly efficient method of shape optimization of structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.1-15
• /
• 1986

The algorithm Proposed utilizes the two-levels technique. In the first level which consists of two phases, the cross-sectional area of the truss member is optimized by transforming the nonlinear problem into SUMT, and solving SUMT utilizing the modified Newton-Rahson method. In the second level, the geometric shape is optimized utilizing the unindirectional search technique of the Powell method which make it possible to minimize only the objective function. The algorithm Proposed in this study is numerically tested for several truss structures with various shapes, loading conditions and design criteria, and compared with the results of the other algorithms to examine its applicability and stability. The numerical comparisons show that the two-Levels algorithm Proposed in this study is safely applicable to any design criteria, and the convergency rate is relathely fast and stable compared with other iteration methods for the geometric optimization of truss structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• The Journal of Engineering Geology
• /
• v.16 no.3 s.49
• /
• pp.255-263
• /
• 2006

In an electrical tomographic survey using an inclined borehole with a pole-dipole array, we must consider several factors: a singular point associated with zero potential difference, a spatial discrepancy between electrode and nodal point in a model due to a inclined borehole, and a variation of geometric factors in connection with a irregular topography. Singular points which are represented by the normal distance from current source to the ground surface can be represented by serveral regions due to a irregular topography of ground surface. The method of element division can be applied to the region in which the borehole is curved, inclined or the distance between the electrodes is shorter than that of nodal points, because the coordinate of each electrode cannot be assigned directly to the nodal point if several electrodes are in an element. Test on a three-dimensional (3-D) synthetic model produces good images of conductive target and shoves stable convergence.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.132-139
• /
• 1992

This paper presents geometrically nonlinear formulation of shell problems using the three-dimensional curved shell element, which includs large displacements and large rotations. Formulations of the geometrically nonlinear problems can be derived in a variety of ways, but most of them have been obtained by assuming that nodal rotations are small. Hence, the tangent stiffness matrix is derived under the assumptions that rotational increments are infinitesimal and the effect of finite rotational increments have to be considered during the equilibrium iterations. To study the large displacement and large rotation problems, the restrictions are removed and the formulations of the curved shell element including the effect of large rotational increments are developed in this paper. The displacement based finite element method using this improved formulation are applied to the analyses of the geometrically nonlinear behaviors of the single and double curved shells, which are compared with the results by others.

• 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.

• 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.

• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.21-30
• /
• 2004

This study deals with free vibrations of laminated composite structures with a channel section using finite element method. In this paper, the mixed finite element method using Lagrangian and Hermite interpolation functions is adopted and a high-order plate theory is used to analyze laminated composite non-prismatic folded plates with a channel section more accurately for free vibration. The theory accounts for parabolic distribution of the transverse shear stress and requires no shear correction factors supposed in the first-order plate theory. An 32×32 matrix is assembled to transform the system element matrices from the local to global coordinates using a coordinate transformation matrix, in which an eighth drilling degree of freedom (DOF) per node is appended to the existing 7-DOF system. The results in this study are compared with those of available literatures for the conventional and first-order plate theory. Sample studies are carried out for various layup configurations and length-thickness ratio, and geometric shapes of plates. The significance of the high-order plate theory in analyzing complex composite structures with a channel section is enunciated in this paper.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.25-30
• /
• 1992

최근 전자계산기를 이용한 진동해석 방법이 눈부시게 발달하여, 일반 구조물 이나 기계 구조물 등의 동특성을 설계 단계에서 정도 높게 예측하는 것이 가능하게 되었다. 그러나 종래의 구조해석은 주어진 시스템의 동특성을 위한 것으로 얻어진 동특성으로부터 질량, 관성제원 및 스프링상수값 등의 설계상 수값을 규명하는 연구는 미미한 실정이다. 이것에 대한 해결방법으로 크게 해석적인 방법과 실험적인 방법으로의 접근이 있어 왔다. 해석적인 방법으로 유한요소해석에서 얻은 모드좌표를 물리좌표로 변환하는 방법으로 Guyan의 정축소와 같은 절점축소를 행하는 방법이 고찰되었다. 실험적인 방법으로 가 진실험에서 얻은 전달함수나 모드파라미터로부터 [M], [K] 행렬을 결정하는 연구가 있었지만 어떤것도 질량, 스프링상수 등의 설계상수를 완전히 규명하 지는 못하였다. 또한, 설계 단계에서 필요한 질량, 관성제원 또는 스프링상수 등의 최적한 값이나, 원하는 시스템특성을 얻을 수 있는 설계상수의 적정한 폭을 구하는 연구는 설계자의 경험과 반복된 시행착오에 의존하는 실정이다. 감도해석은 이러한 문제점을 개선하는 수단으로 설계변수에 대한 동특성의 변화율을 구하는 것이다. 감도해석을 수행하는 것은 어느 설계변수를 수정하 는 것이 주어진 동특성에 부합되는 지를 알려주고, 어느 것을 수정하는 것이 원하는 방향의 동특성변화에 가장 효과적인지를 알려주는 것이다. 따라서 감 도해석을 이용하여 설계의 최적화 프로그램을 만들수 있고, 이것은 설계자가 요구하는 동특성을 목적함수로 하여 주어진 구조물을 최적화하는 설계상수 값을 얻을 수 있게 한다. 본 논문에서는 강체모델의 동특성으로부터 모델의 설계 상수를 규명하고, 동특성의 개선을 위하여 설계변수의 변경량을 물리좌 표계에서 얻는것을 목적으로 한다. 강체 마운트계의 관성제원 및 마운트강성 의 규명을 위하여 임으로 주어진 설계상수를 모델데이타로 하여 관성제원과 스프링 강성을 구하였다. 관성제원의 규명은 주어진 모델의 관성값을 모르는 것으로 하여 임의의 초기 관성값으로 감도해석에 의해 주어진 계의 관성값 을 물리 좌표계에서 규명하였다. 마운트 강성의 규명도 관성제원의 규명과 같은 방법으로 임의의 강성값으로 감도해석을 하여 강성값을 규명하였다. 또 한 감도해석에 의한 동특성 변경은 특정한 고유진동 수의 변경이 필요할 때, 고유진동수의 이동을 위한 관성제원의 변경 및 마운트 강성변경값을 예측할 수 있다. 본 연구수행의 기본적인 흐름도는 Fig.1.1과 같다. 위와 같은 작업 으로 엔진 마운트와 같은 강체 모델의 시스템 규명을 행하는 경우에 유한요 소해석 및 가진 실험으로 얻은 고유진동수의 정보 또는 원하는 고유진동수 의 특성을 기본으로 실제 설계에서 사용이 가능하도록 물리 좌표계에서 관 성 제원 및 스프링상수를 구할 수 있을 것이다.