• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.3
• /
• pp.313-321
• /
• 2012

Vibration analysis of a thin-walled structure can be performed with a consistent mass matrix determined by the shape functions of all degrees of freedom (d.o.f.) used for construction of conventional stiffness matrix, or with a lumped mass matrix. In similar way stability of a structure can be analysed with consistent geometric stiffness matrix or geometric stiffness matrix with lumped buckling load, related only to the rotational d.o.f. Recently, the simplified mass matrix is constructed employing shape functions of in-plane displacements for plate deflection. In this paper the same approach is used for construction of simplified geometric stiffness matrix. Beam element, and triangular and rectangular plate element are considered. Application of the new geometric stiffness is illustrated in the case of simply supported beam and square plate. The same problems are solved with consistent and lumped geometric stiffness matrix, and the obtained results are compared with the analytical solution. Also, a combination of simplified and lumped geometric stiffness matrix is analysed in order to increase accuracy of stability analysis.

• 제어로봇시스템학회:학술대회논문집
• /
• 1996.10b
• /
• pp.565-568
• /
• 1996

In this paper, a solution method is proposed to calculate the optimum solution to discrete optimal H$_{.inf}$ control problem for feedback of linear time-invariant system states and disturbance variable. From the results of this study, condition of existence and uniqueness of its solution is that transfer matrix of controlled variable to input variable is left invertible and has no invariant zeros on the unit circle of the z-domain as well as extra geometric conditions given in this paper. Through a numerical example, the noniterative solution method proposed in this paper is illustrated.

• Architectural research
• /
• v.13 no.1
• /
• pp.17-24
• /
• 2011

Isogeometric solution of Poisson equation is provided. NURBS (NonUniform B-spline Surface) is introduced to express both geometry of structure and unknown field of governing equation. The terms of stiffness matrix and load vector are consistently derived with very accurate geometric definition. The validity of the isogeometric formulation is demonstrated by using two numerical examples such as square plate and L-shape plate. From numerical results, the present solutions have a good agreement with analytical and finite element (FE) solutions with the use of a few cells in isogeometric analysis.

• Structural Engineering and Mechanics
• /
• v.24 no.6
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.

• Steel and Composite Structures
• /
• v.27 no.1
• /
• pp.13-25
• /
• 2018

The main objective of this study is to develop a dual approach for geometrically nonlinear finite element analysis of plane truss structures. The geometric nonlinearity is considered using the Total Lagrangian formulation. The nonlinear solution is obtained by introducing and minimizing an objective function subjected to displacement-type constraints. The proposed method can fully trace the whole equilibrium path of geometrically nonlinear plane truss structures not only before the limit point but also after it. No stiffness matrix is used in the main approach and the solution is acquired only based on the direct classical stress-strain formulations. As a result, produced errors caused by linearization and approximation of the main equilibrium equation will be eliminated. The suggested algorithm can predict both pre- and post-buckling behavior of the steel plane truss structures as well as any arbitrary point of equilibrium path. In addition, an equilibrium path with multiple limit points and snap-back phenomenon can be followed in this approach. To demonstrate the accuracy, efficiency and robustness of the proposed procedure, numerical results of the suggested approach are compared with theoretical solution, modified arc-length method, and those of reported in the literature.

• Communications of the Korean Mathematical Society
• /
• v.12 no.1
• /
• pp.177-192
• /
• 1997

The basic theory of the quasi-birth-and-death process is extended to a process which is inhomogeous in levels. Several key results in the standard homogeneous theory hold in a more general context than that usually stated, in particular not requiring positive recurrence. Theser results are subsumed under our development. The treatment is entirely probabilistic.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.631-641
• /
• 1988

Computational procedures associated with finite element nonlinear analysis of plane frame structures were examined and new solution schemes were suggested. Element stiffness matrix was derived from the principle of virtual displacements. Geometric and material nonlinearities were considered in the formulation. Solution method was based upon the constant displacement length method in conjunction with the Newton-Raphson method. New solution schemes were introduced in determining the initial load increment and the sign of load increments and predicting the length of displacement increment to improve user convenience, efficiency and stability. Numerical experiments were performed for several typical problems and suggested schemes were found efficient and convenient for analyzing nonlinear frame structures.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.28 no.1
• /
• pp.79-92
• /
• 1992

This is analyzed using the finite element method which is appling excellent isoparametric curve element in the aspect of large usages of dynamic responses in which is regarding geometric and material nonlinear of a large scale shell structure of an airplane, a submarine, a ship, and an ocean structure. The solution of dynamic equations is got by direct integration method using time-stepping procedure and regarding Central Difference Method of the both solutions. But because formal matrix factorization is not necessary in each time step and it does not take less time to compute relatively, this method must be regarded very few time steps on the condition. Axisymmatric shell problems are inspected using 8 node Isoparametric element in this paper. Partial axisymmatric spherical shell is used as a model to analyze axisymmatric nonlinear dynamic behavior regarding. Total Lagrangian formulation in geometric nonlinear behavior and elastio-viscoplastic in material nonlinear behavior.

• Journal of Korean Institute of Industrial Engineers
• /
• v.21 no.3
• /
• pp.345-356
• /
• 1995

The aim of this paper is to analyze the batch service queueing model M/M(a, b ; ${\mu}_k/1$) under general bulk service rule with mean service rate ${\mu}_k$ for a batch of k units, where $a{\leq}k{\leq}b$. This queueing model consists of the two-dimensional state space so that it is characterized by two-dimensional state Markov process. The steady-state solution and performane measure of this process are derived by using Matrix Geometric method. Meanwhile, a new approach is suggested to calculate the two-dimensional traffic density R which is used to obtain the steady-state solution. In addition, to determine the optimal service initiation threshold a, a decision model of this queueing system is developed evaluating cost of service per batch and cost of waiting per customer. In a job order production system, the decision-making procedure presented in this paper can be applicable to determining when production should be started.

• Structural Engineering and Mechanics
• /
• v.39 no.1
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.