• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.227-233
• /
• 1992

To study the large displacement and large rotation problems, geometrically nonlinear formulation of eccentric degenerated beam element has been developed, where the restrictions of infinitesimal rotation increments are removed and the incremental equations are derived using the Taylor series expansion of the displacement function at time t+dt. The geometrically nonlinear analyses are carried out for the cases of cantilever, square frame, shallow arch and 45-degree bend beam and all of them are compared with each of the other results published. The element developed in the present research can be efficiently utilized for analysis of the nonlinear behaviours of structures when displacements and rotations are large.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• Structural Engineering and Mechanics
• /
• v.68 no.5
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.4
• /
• pp.163-170
• /
• 2023

In the fluid-structure interaction analysis, the finite element formulation is performed for the wave equation for dynamic fluid pressure, and the dynamic pressure is defined as a degree of freedom at the fluid nodes. Therefore, to connect the fluid to the structure, it is necessary to connect the degree of freedom of fluid dynamic pressure and the degree of freedom of structure displacement through an interface element derived from the relationship between dynamic pressure and displacement. The previously proposed fluid-structure interface elements use conformal finite element meshes in which the fluid and structure match. However, it is challenging to construct conformal meshes when complex models, such as water purification plants and wastewater treatment facilities, are models. Therefore, to increase modeling convenience, a method is required to model the fluid and structure domains by independent finite element meshes and then connect them. In this study, two fluid-structure interface elements, one based on constraints and the other based on the integration of nonsmooth functions, are proposed in nonconformal finite element meshes for structures and fluids, and their accuracy is verified.

• Steel and Composite Structures
• /
• v.48 no.2
• /
• pp.113-130
• /
• 2023

The present manuscript aims to investigate the deviation between the middle surface (MS) and neutral surface (NS) formulations on the static response of bi-directionally functionally graded (BDFG) porous plate. The higher order shear deformation plate theory with a four variable is exploited to define the displacement field of BDFG plate. The displacement field variables based on both NS and on MS are presented in detail. These relations tend to get and derive a new set of boundary conditions (BCs). The porosity distribution is portrayed by cosine function including three different configurations, center, bottom, and top distributions. The elastic foundation including shear and normal stiffnesses by Winkler-Pasternak model is included. The equilibrium equations based on MS and NS are derived by using Hamilton's principles and expressed by variable coefficient partial differential equations. The numerical differential quadrature method (DQM) is adopted to solve the derived partial differential equations with variable coefficient. Rigidities coefficients and stress resultants for both MS and NS formulations are derived. The mathematical formulation is proved with previous published work. Additional numerical and parametric results are developed to present the influences of modified boundary conditions, NS and MS formulations, gradation parameters, elastic foundations coefficients, porosity type and porosity coefficient on the static response of BDFG porous plate. The following model can be used in design and analysis of BDFG structure used in aerospace, vehicle, dental, bio-structure, civil and nuclear structures.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.1
• /
• pp.31-40
• /
• 2017

Absolute nodal coordinate formulation was developed in the mid-1990s, and is used in the flexible dynamic analysis. In the process of deriving the equation of motion, if the order of polynomial referring to the displacement field increases, then the degrees of freedom increase, as well as the analysis time increases. Therefore, in this study, the primary objective was to reduce the analysis time by transforming the dimensional equation of motion to a non-dimensional equation of motion. After the shape function was rearranged to be non-dimensional and the nodal coordinate was rearranged to be in length dimension, the non-dimensional mass matrix, stiffness matrix, and conservative force was derived from the non-dimensional variables. The verification and efficiency of this non-dimensional equation of motion was performed using two examples; cantilever beam which has the exact solution about static deflection and flexible pendulum.

• The Transactions of the Korea Information Processing Society
• /
• v.4 no.9
• /
• pp.2403-2411
• /
• 1997

Direct3D is a recently developed 3D graphic accelerator such that we could make a Windows graphic programs with Direct3D easily. This paper presents the fractal terrain model viewer developed in terms of Direct3D. The object classes and the input dialog box used the model viewer, program initializing process, and the flow of the model viewer are discussed. Finally, the output of terrain formulation in wire-frame, solid, and point type form are presented. To show the results, the data of a terrain model were made by a fractal technique, the midpoint displacement methods with square lattices of points.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.123-132
• /
• 2001

An improved version of the Element-free Galerkin method (EFGM) is presented here for addressing the problem of transverse shear locking in shear-deformable beams with a high length over thickness ratio. Based upon Timoshenko's theory of thick beams, it has been recognized that shear locking will be completely eliminated if the rotation field is constructed to match the field of slope, given by the first derivative of displacement. This criterion is applied directly to the most commonly implemented version of EFGM. However in the numerical process to integrate strain energy, the second derivative of the standard Moving Least Square (MLS) shape functions must be evaluated, thus requiring at least a $C^1$ continuity of MLS shape functions instead of $C^0$ continuity in the conventional EFGM. Yet this hindrance is overcome effortlessly by only using at least a $C^1$ weight function. One-dimensional quartic spline weight function with $C^2$ continuity is therefore adopted for this purpose. Various numerical results in this work indicate that the modified version of the EFGM does not exhibit transverse shear locking, reduces stress oscillations, produces fast convergence, and provides a surprisingly high degree of accuracy even with coarse domain discretizations.

• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.501-518
• /
• 2007

In this work, a new approach is developed for dynamic analysis of a composite beam with an interply crack, based on finite element solution of partial differential equations with the use of the COMSOL Multiphysics package, allowing for fast and simple change of geometric characteristics of the delaminated area. The use of COMSOL Multiphysics package facilitates automatic mesh generation, which is needed if the problem has to be solved many times with different crack lengths. In the model, a physically impossible interpenetration of the crack faces is prevented by imposing a special constraint, leading to taking account of a force of contact interaction of the crack faces and to nonlinearity of the formulated boundary value problem. The model is based on the first order shear deformation theory, i.e., the longitudinal displacement is assumed to vary linearly through the beam's thickness. The shear deformation and rotary inertia terms are included into the formulation, to achieve better accuracy. Nonlinear partial differential equations of motion with boundary conditions are developed and written in the format acceptable by the COMSOL Multiphysics package. An example problem of a clamped-free beam with a piezoelectric actuator is considered, and its finite element solution is obtained. A noticeable difference of forced vibrations of the delaminated and undelaminated beams due to the contact interaction of the crack's faces is predicted by the developed model.

• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.263-278
• /
• 2008

Structural pipe-in-pipe cross-sections have significant potential for application in offshore oil and gas production systems because of their property that combines insulation performance with structural strength in an integrated way. Such cross-sections comprise inner and outer thin walled pipes with the annulus between them fully filled by a selectable thick filler material to impart an appropriate combination of properties. Structural pipe-in-pipe cross-sections can exhibit several different collapse mechanisms and the basis of the preferential occurrence of one over others is of interest. This paper presents an elastic analyses of a structural pipe-in-pipe cross-section when subjected to external hydrostatic pressure. It formulates and solves the static and elastic buckling problem using the variational principle of minimum potential energy. The paper also investigates a simplified formulation of the problem where the outer pipe and its contact with the filler material is considered as a 'pipe on an elastic foundation'. Results are presented to show the variation of elastic buckling pressure with the relative elastic modulus of the filler and pipe materials, the filler thickness and the thicknesses of the inner and outer pipes. The range of applicability of the simplified 'pipe on an elastic foundation' analysis is also presented. A brief review of the types of materials that could be used as the filler is combined with the results of the analysis to draw conclusions about elastic buckling behaviour of structural pipe-in-pipe cross-sections.