• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.7 s.166
• /
• pp.1205-1214
• /
• 1999

A new quadrilateral plane stress element is developed for efficient and accurate analysis of thermal stress problems. It is convenient to use the same mesh and the same shape functions for thermal analysis and stress analysis. But, because of the inconsistency between deformation related strain field and thermal strain field, oscillatory responses and considerable errors in stresses are resulted in. To avoid undesired oscillations, strain approximation is enhanced by supplementing several assumed strain terms based on the variational principle. Thermal deformation is incorporated into the generalized mixed variational principle for displacement, strain and stress fields, and basic equations for the modified enhanced assumed strain method are derived. For the stress approximation of bilinear elements, the $5{\beta}$ version of Pian and Sumihara is adopted. The numerical results for several problems show that the present element behaves well and reduces oscillatory responses. it also results in almost the same magnitude of error as compared with the quadratic element.

• Proceedings of the KIEE Conference
• /
• 1990.07a
• /
• pp.433-436
• /
• 1990

An optimal Stabilization technique for a bilinear in duction model is introduced. This technique includes to o parts; the one is an stabilization control using Lyap unov Function which has the form of a sum of linear and quadratic function of the state variables, and the other is an optimal control using the performance index which depends on the choice of the elements of the Ly apunov matrices concerning both the state variables and the input variables. Therefore, induction motor is drived with the shorter transient time of the state variables and with the smaller overshoot of the ones, simulation results are obtained from a digital computer. Experimental ones are obtained from implementation of the optimizing controller using 8086 microprocessor kits and analog circuits are compared.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.291-300
• /
• 1998

This paper demonstrates how ambient vibration measurements at a limited number of locations can be effectively utilized to estimate parameters of a finite element model of a large-scale structural system involving a large number of elements. System identification using ambient vibration measurements presents a challenge requiring the use of special identification techniques, which ran deal with very small magnitudes of ambient vibration contaminated by noise without the knowledge of input farces. In the present study, the modal parameters such as natural frequencies, damping ratios, and mode shapes of the structural system were estimated by means of appropriate system identification techniques including the random decrement method. Moreover, estimation of parameters such as the stiffness matrix of the finite element model from the system response measured by a limited number of sensors is another challenge. In this study, the system stiffness matrix was estimated by using the quadratic optimization involving the computed and measured modal strain energy of the system, with the aid of a sensitivity relationship between each element stiffness and the modal parameters established by the second order inverse modal perturbation theory. The finite element models thus identified represent the actual structural system very well, as their calculated dynamic characteristics satisfactorily matched the observed ones from the ambient vibration test performed on a large-scale structural system subjected primarily to ambient wind excitations. The dynamic models identified by this study will be used for design of an active mass damper system to be installed on this structure fer suppressing its wind vibration.

• Structural Engineering and Mechanics
• /
• v.23 no.1
• /
• pp.15-27
• /
• 2006

A locally varying temperature field or a mixture of two or more different materials can cause local variation of elasticity properties of a beam. In this paper, a new Euler-Bernoulli beam element with varying Young's modulus along its longitudinal axis is presented. The influence of axial forces according to the linearized 2nd order beam theory is considered, as well. The stiffness matrix of this element contains the transfer constants which depend on Young's modulus variation and on axial forces. Occurrence of the polynomial variation of Young's modulus has been assumed. Such approach can be also used for smooth local variation of Young's modulus. The critical loads of the straight slender columns were studied using the new beam element. The influence of position of the local Young's modulus variation and its type (such as linear, quadratic, etc.) on the critical load value and rate of convergence was investigated. The obtained results based on the new beam element were compared with ANSYS solutions, where the number of elements gradually increased. Our results show significant influence of the locally varying Young's modulus on the critical load value and the convergence rate.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.1
• /
• pp.75-82
• /
• 1999

일반 해양 구조물이나 해저면에 설치되는 해저 관로는 외력에 의한 변형이 발생된다. 구조물 형상이 복잡하거나, 구성 요소의 개수가 많을 경우 응력해석 시 많은 초기값이 필요하고 해석 시간 또는 장 시간 소요된다. 해양 구조물에 작용하는 대표적인 외력은 파도, 조류, 바람이고 이런 외력은 구조물의 사용 기간(operation life)동안 계속적으로 작용하기 때문에 구조물의 변형율은 항상 허용치 안에서 발생되도록 설계되어야 한다. 허용 변형은 탄성범위 내에 존재해야 하며, 비교적 큰 변형을 일으키는 구조물이나 해저파이프라인의 응력해석을 수치적으로 접근하는 방법을 고찰하였다. 평행상태의 하중 벡터값만 직각 좌표계에서 인트린직(intrinsic) 좌료로 변환시킬 때 변형이 발생함으로, 본 논문에서 소개하는 이차 요소(quadratic element)방법을 사용할 경우 수치해석 시 많은 장점이 있다는 것을 보여준다. 본 방법을 도입함으로써 비교적 큰 변형이 발생되는 구조물 해석 시 일반 수치해석 방법과 그 결과는 같으나 해석 시간을 단축시킬 수 있다는 장점이 있다. 응력 해석 시 국부 강도 행열(element stiffness matrix)은 방향과 무관하며 이차요소 방법을 사용하여 각 요소 벡터를 발생시켰다. 해저관로 일점 상승 시 관로에 작용하는 변형과 상승력에 따른 휨 모멘트를 산출하여 일반적으로 사용되는 선형이론과 비교하였다.

• Structural Engineering and Mechanics
• /
• v.45 no.2
• /
• pp.233-257
• /
• 2013

Given the yield shear of a single-story inelastic structure with simple eccentricity, the problem of strength distribution among the resisting elements is investigated, with respect to minimize its torsional response during a ground motion. Making the hypothesis that the peak accelerations, of both modes of vibration, are determined from the inelastic acceleration spectrum, and assuming further that a peak response quantity is obtained by an appropriate combination rule (square root of sum of squares-SRSS or complete quadratic combination-CQC), the first aim of this study is to present an interaction relationship between the yield shear and the maximum torque that may be developed in such systems. It is shown that this torque may be developed, with equal probability, in both directions (clockwise and anticlockwise), but as it is not concurrent with the yield shear, a rational design should be based on a combination of the yield shear with a fraction of the peak torque. The second aim is to examine the response of such model structures under characteristic ground motions. These models provide a rather small peak rotation and code provisions that are based on such principles (NBCC-1995, UBC-1994, EAK-2000, NZS-1992) are superiors to EC8 (1993) and to systems with a stiffness proportional strength distribution.

• Journal of Korean Port Research
• /
• v.14 no.4
• /
• pp.451-461
• /
• 2000

The design and maintenance of navigation channel and water facilities of an harbor which is located at the mouth of river or at the estuary area are difficult due to the complexity of estuarial water and sediment circulation. Effects of deepening navigable waterways, of changing coastline configurations, or of discharging dredged material to the open sea are necessary to be investigated and predicted in terms of water quality and possible physical changes to the coastal environment. A borad analysis of the transport mechanism in the estuary area was made in terms of sediment property, falling velocity, concentration and flow characteristics. In order to simulate the transport processes, a two-dimensional finite element model is developed, which includes erosion, transport and deposition mechanism of suspended sediments. Galerkin’s weighted residual method is used to solve the transient convection-diffusion equation. The fluid domain is subdivided into a series of triangular elements in which a quadratic approximation is made for suspended sediment concentration. Model could deal with a continuous aggregation by stipulating the settling velocity of the flocs in each element. The model provides suspended sediment concentration, bed shear stress, erosion versus deposition rate and bed profile at the given time step.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.8
• /
• pp.3295-3300
• /
• 2012

Three Dimensional finite element method (FEM) was used to obtain the stress intensity factor for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model. Nodes are generated by the bubble packing, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. To examine accuracy and efficiency of the present system, the stress intensity factor for a semi-elliptical surface crack in a plate subjected to uniform tension is calculated, and compared with Raju-Newman's solutions. Then the system is applied to analyze interaction effects of two dissimilar semi-elliptical cracks in a plate subjected to uniform tension.

• Structural Engineering and Mechanics
• /
• v.76 no.6
• /
• pp.799-821
• /
• 2020

The optimum design of truss structures is one of the significant categories in structural optimization that has widely been applied by researchers. In the present study, new mathematical programming called Consistent Approximation (CONAP) method is utilized for the simultaneous optimization of the size and shape of truss structures. The CONAP algorithm has already been introduced to optimize some structures and functions. In the CONAP algorithm, some important parameters are designed by employing design sensitivities to enhance the capability of the method and its consistency in various optimum design problems, especially structural optimization. The cross-sectional area of the bar elements and the nodal coordinates of the truss are assumed to be the size and shape design variables, respectively. The displacement, allowable stress and the Euler buckling stress are taken as the design constraints for the problem. In the proposed method, the primary optimization problem is replaced with a sequence of explicit sub-problems. Each sub-problem is efficiently solved using the sequential quadratic programming (SQP) algorithm. Several truss structures are designed by employing the CONAP method to illustrate the efficiency of the algorithm for simultaneous shape and size optimization. The optimal solutions are compared with some of the mathematical programming algorithms, the approximation methods and metaheuristic algorithms those reported in the literature. Results demonstrate that the accuracy of the optimization is improved and the convergence rate speeds up.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.1-13
• /
• 1998

In the present paper, the hydrodynamic characteristics of three dimensional hydrofoils moving with a constant speed below the free surface using a higher-order boundary element method based on 9-node Lagrangian curvilinear elements are investigated. A bi-quadratic spline scheme is employed to improve the numerical results on the free surface. To validate the present scheme, the calculated results are compared with the analytic solutions for a submerged sphere and a spheroid showing a good agreement. For the validation of the hydrofoil study, the computed lift and drag of a hydrofoil having $NACA64_{1}A412$ section with aspect ratio(A.R.) of 4 are compared with the experimental data by Wadlin et al.[28]. The comparison covers a number of variations of angle of attack and submergence depth. Then, using an A.R. hydrofoil with NACA0012 section, the free surface on the lift and drag are investigated and these are compared with the previous results. The wave elevations and patterns created by the aforementioned submerged bodies are also investigated with Froude numbers and submergences.