• Computational Structural Engineering
• /
• v.11 no.4
• /
• pp.155-168
• /
• 1998

보강된 판 및 쉘구조의 안정성 및 후좌굴을 포함하는 기하학적 비선형 해석을 수행하기 위하여, total Lagrangian formulation에 근거한 연속체의 증분평형방정식으로부터 변형된 쉘요소인 유한요소이론을 제시하였다. 쉘구조의 곡률이 불연속적으로 변하거나 쉘부재들이 유한한 각도로 만나는 보강된 판 및 쉘구조의 비선형 해석이 가능하도록 주부재와 보강재 간의 연결점에 대한 일반적인 변환관계를 제시하였으며 좌굴해석 및 기하학적 비선형해석의 경우에 해의 정확성 및 수렴성을 개선시키기 위하여 접선강도행렬 산정시 회전각의 2차항을 포함시켰다. 또한, shear locking 현상을 극복하기 위하여 감차적분을 적용하였고 쉘구조의 좌굴해석에서는 power method를 적용하여 해석의 효율을 높였으며, 후좌굴해석에서는 변위 및 하중증분법을 적절히 결합시켜 보강된 쉘구조의 후좌굴 거동추적이 용이하였다. 또한, 입력자료를 손쉽게 준비하고 좌굴모드 및 후좌굴거동을 효율적으로 분석하기 위하여 전, 후 처리 프로그램을 개발하였고 다양한 해석예제를 통하여 다른 문헌의 해석결과를 비교함으로써 본 연구에서 개발된 유한요소 해석프로그램의 타당성 및 정확성을 입증하였다.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.161-170
• /
• 1998

A general formulation of the long cylinder tolerance design for the joint structure is here presented. The aim of this paper is to calculate the tolerance of joint by defining tolerance as a kind of uncertainty and to obtain the robustness of the joint structure. It is formulated on the bases of the minimum sensitivity theorem. The objective function is the tolerance sensitivity for the Von-Mises stress. It also took into full account the stress, displacement and weight constraints. PLBA(Pshenichny-Lim-Belegundu-Arora) algorithm is used to solve the constrained nonlinear optimization problem. The finite element analysis is performed with CST(Constant-Strain-Triangle) axisymmetric element. Sensitivities for design variables are calculated by the direct differentiation method. The numerical result is presented for the cylindrical structure where the joint tolerance is treated as random variables.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.4
• /
• pp.19-36
• /
• 1987

In this research the coupling problem between the elastic structure and the fluid, specially the hydroelastic harmonic vibration problem, is studied. In order to couple the domains, i.e., the structural domain and the fluid domain, the boundary integral method(direct boundary integral formulation) is used in the fluid domain in combination with the finite element method for the structure. The boundary integral method has been widely developed to apply it to the hydroelastic vibration problem. The hybrid boundary integral method using eigenfunctions on the radiation boundaries and the boundary integral method using the series form image-functions to replace the even bottom and free surface boundaries in case of high frequencies have been developed and tested. According to the boundary conditions and the frequency ranges the different boundary integral methods with the different idealizations of the fluid boundaries have been studied. Using the same interpolation functions for the pressure distribution and the displacement the two domains have been coupled and using Hamilton principle the solution of the hydroelastic have been obtained through the direct minimizing process. It has become evident that the finite-boundary element method combining with the eigenfunction or the image-function method give good results in comparison with the experimental ones and the other numerical results by the finite element method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.3-17
• /
• 2003

This paper deals with the problems and their possible solutions in the development of finite element for analysis of shell. Based on these solution schemes, a series of flat shell elements are established which show no signs of membrane locking and other defects even though the coarse meshes are used. In the element formulation, non-conforming displacement modes are extensively used for improvement of element behaviors. A number of numerical tests are performed to prove the validity of the solutions to the problems involved in establishing a series of high performance flat shell elements. The test results reveal among others that the high accuracy and fast convergence characteristics of the elements are obtainable by the use of various non-conforming modes and that the ‘Direct Modification Method’ is a very useful tool for non-conforming elements to pass the patch tests. Furthermore, hierarchical and higher order non-conforming modes are proved to be very efficient not only to make an element insensitive to the mesh distortion but also to remove the membrane locking. Some numerical examples are solved to demonstrate the validity and applicability of the presented elements to practical engineering shell problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.937-956
• /
• 2014

The MEMS structures usually are made from silicon; consideration of the viscoelastic effect in microbeams duo to the phenomena of silicon creep is necessary. Application of the fractional model of microbeams made from viscoelastic materials is studied in this paper. Quasi-static and dynamical responses of an electrically actuated viscoelastic microbeam are investigated. For this purpose, a nonlinear finite element formulation of viscoelastic beams in combination with the fractional derivative constitutive equations is elucidated. The four-parameter fractional derivative model is used to describe the constitutive equations. The electric force acting on the microbeam is introduced and numerical methods for solving the nonlinear algebraic equation of quasi-static response and nonlinear equation of motion of dynamical response are described. The deflected configurations of a microbeam for different purely DC voltages and the tip displacement of the microbeam under a combined DC and AC voltages are presented. The validity of the present analysis is confirmed by comparing the results with those of the corresponding cases available in the literature.

• Computational Structural Engineering
• /
• v.3 no.1
• /
• pp.97-107
• /
• 1990

The development of an improved degenerated shell element is presented in this paper. In the formulation of this element, an enhanced interpolation of transverse shear strains in the natural coordinate system is used to overcome the shear locking problem ; the reduced integration technique in in-plane strains is applied to avoid the membrane locking behavior ; and selective addition of the nonconforming displacement modes improve the element performances. This element is free of serious locking problems and undesirable compatible or commutable spurious kinematic deformation modes, and passes the patch tests. To illustrate the performance of this improved degenerated shell element, some benchmark problems are presented. Numerical results indicate that the new element shows fast convergence and reliable solutions.

• Structural Engineering and Mechanics
• /
• v.54 no.6
• /
• pp.1153-1174
• /
• 2015

Deployable structures have gained more and more applications in space and civil structures, while it takes a large amount of computational resources to analyze this kind of multibody systems using common analysis methods. This paper presents a new approach for dynamic analysis of multibody systems consisting of both rigid bars and arbitrarily shaped rigid bodies. The bars and rigid bodies are connected through their nodes by ideal pin joints, which are usually fundamental components of deployable structures. Utilizing the Moore-Penrose generalized inverse matrix, equations of motion and constraint equations of the bars and rigid bodies are formulated with nodal Cartesian coordinates as unknowns. Based on the constraint equations, the nodal displacements are expressed as linear combination of the independent modes of the rigid body displacements, i.e., the null space orthogonal basis of the constraint matrix. The proposed method has less unknowns and a simple formulation compared with common multibody dynamic methods. An analysis program for the proposed method is developed, and its validity and efficiency are investigated by analyses of several representative numerical examples, where good accuracy and efficiency are demonstrated through comparison with commercial software package ADAMS.

• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Structural Engineering and Mechanics
• /
• v.33 no.3
• /
• pp.373-385
• /
• 2009

The purpose of this paper is to study shear locking-free parametric earthquake analysis of thick and thin plates using Mindlin's theory, to determine the effects of the thickness/span ratio, the aspect ratio and the boundary conditions on the linear responses of thick and thin plates subjected to earthquake excitations. In the analysis, finite element method is used for spatial integration and the Newmark-${\beta}$ method is used for the time integration. Finite element formulation of the equations of the thick plate theory is derived by using higher order displacement shape functions. A computer program using finite element method is coded in C++ to analyze the plates clamped or simply supported along all four edges. In the analysis, 17-noded finite element is used. Graphs are presented that should help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that 17-noded finite element can be effectively used in the earthquake analysis of thick and thin plates. It is also concluded that, in general, the changes in the thickness/span ratio are more effective on the maximum responses considered in this study than the changes in the aspect ratio.