• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.2
• /
• pp.324-345
• /
• 2015

A new procedure for determining properties of thick plate finite elements, based on the modified Mindlin theory for moderately thick plate, is presented. Bending deflection is used as a potential function for the definition of total (bending and shear) deflection and angles of cross-section rotations. As a result of the introduced interdependence among displacements, the shear locking problem, present and solved in known finite element formulations, is avoided. Natural vibration analysis of rectangular plate, utilizing the proposed four-node quadrilateral finite element, shows higher accuracy than the sophisticated finite elements incorporated in some commercial software. In addition, the relation between thick and thin finite element properties is established, and compared with those in relevant literature.

• Structural Engineering and Mechanics
• /
• v.30 no.4
• /
• pp.387-402
• /
• 2008

The Integrated Force Method (IFM) has been developed in recent years for the analysis of civil, mechanical and aerospace engineering structures. In this method all independent or internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. The solution by IFM needs the computation of element equilibrium and flexibility matrices from the assumed displacement, stress-resultant fields and material properties. This paper presents a general purpose code for the automatic generation of element equilibrium and flexibility matrices for plate bending elements using the Integrated Force Method. Kirchhoff and the Mindlin-Reissner plate theories have been employed in the code. Paper illustrates development of element equilibrium and flexibility matrices for the Mindlin-Reissner theory based four node quadrilateral plate bending element using the Integrated Force Method.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.2
• /
• pp.46-53
• /
• 2005

Compare to the large number of curved quadrilateral degenerated-shell elements, there are only a very few curved triangular degenerated-shell elements. Based on the assumed natural strain sampling scheme previously developed for a quadratic degenerated-shell element for linear analysis, this paper devises geometric non-linear six-node degenerated-shell element. The element can be curved and is only equipped with the standard nodal d.o.f.'s. Careful consideration has been exercised to circumvent various locking phenomena that plague degenerated-shell element. Numerical examples are presented to illustrate efficiency.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.139-170
• /
• 2013

We extend a p-hierarchical decomposition of the second degree finite element space of N$\acute{e}$d$\acute{e}$lec for tetrahedral meshes in three dimensions given in [1] to meshes with hexahedral elements, and derive p-hierarchical decompositions of the second degree finite element space of Raviart-Thomas in two dimensions for triangular and quadrilateral meshes. After having proved stability of these subspace decompositions and requiring certain saturation assumptions to hold, we construct a local a posteriori error estimator for fem and bem coupling of a time-harmonic electromagnetic eddy current problem in $\mathbb{R}^3$. We perform some numerical tests to underline reliability and efficiency of the estimator and test its usefulness in an adaptive refinement scheme.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.919-935
• /
• 2014

In this study, the control of the shape of pre-stressed cable structures and the effective control element were examined. The process of deriving the displacement control equations using the force method was explained, and the concurrent control scheme (CCS) and the sequence control scheme (SCS) were proposed. To explain the control scheme process, the quadrilateral cable net model was adopted and classified into a regular model and an irregular model for the analysis of the control results. In the control analysis of the regular model, the CCS and SCS analysis results proved reliable. For the SCS, the errors occur in the control stage and varied according to the control sequence. In the control analysis of the irregular model, the CCS analysis result also proved relatively reliable, and the SCS analysis result with the correction of errors in each stage was found nearly consistent with the target shape after the control. Finally, to investigate an effective control element, the Geiger cable dome was adopted. A set of non-redundant elements was evaluated in the reduced row echelon form of a coefficient matrix of control equations. Important elements for shape control were also evaluated using overlapping elements in the element sets, which were selected based on cable adjustments.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.3 no.2
• /
• pp.48-62
• /
• 1985

A method for the free network adjustment is proposed, based on the application of generalized inverse matrix (g-inverse). If the network adjustment is executed according to the solution with parameters, especially when all coordinates are considered as parameters to keep unity strength, the matrix of normal equation will be singular. This paper discusses the problem of singular matrix and the analysis of accuracy between conventional method and the free network adjustment of trilateration. In case of the adjustment, the RMS errors of adjusted X, Y coordinates are increased to 35.6% in a polygon, central-point figure, and 50.5% in a quadrilateral. In the elements of error ellipse, the RMS errors are decreased by $\pm$24.5% (a) and $\pm$5.0 % (b) in the polygon, $\pm$42.6% (a) and $\pm$49.2% (b) in the quadrilateral. Introduction of free network adjustment, therefore, could be applied to improvement of relative accuracy in the horizontal positioning.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.289-294
• /
• 2002

In this study, performance, flow characteristics and noise of a cross-flow-fan system, used in indoor unit of the split-type air conditioner, were predicted by computational simulation. Triangular elements were used to mesh the calculation domain and quadrilateral elements were attached to the blade surfaces and walls to enhance the simulation quality. The unsteady incompressible Wavier-Stokes equations were solved using a sliding mesh technique on the interface between rotating fan region and the outside. Two stripes of velocity stream inside the cross-flow-fan were shown - the one was due to the eccentric vortex and the other was due to the normal entrance flow. As the flow rate increased, the center of the eccentric vortex moved toward the inner blade tip and rear-guide, and the exiting flow still had velocity variation along the stabilizer, which can increase the noise level. The acoustic pressure was calculated by using Lowson's equation. From the calculated acoustic pressure, it was found that the trailing edge is a dominant of acoustic generation.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2004.10a
• /
• pp.79-82
• /
• 2004

In this study, a new refinement technique for 3-dimensional hexahedral element mesh is proposed, which is aimed at the control of mesh density. With the proposed scheme the mesh is refined adaptively to the elemental error which is estimated by 'a posteriori' error estimator based on the energy norm. A desired accuracy of an analysis i.e. a limit of error defines the new desired mesh density map on the current mesh. To obtain the desired mesh density, the refinement procedure is repeated iteratively until no more elements to be refined exist. In the algorithm, at first the regions of mesh to be refined are defined and, then, the zero-thickness element layers are inserted into the interfaces between the regions. All the meshes in the regions, in which the zero-thickness layers are inserted, are to be regularized in order to improve the shape of the slender elements on the interfaces. This algorithm is tested on a simple shape of 2-d quadrilateral element mesh and 3-d hexahedral element mesh. A numerical example of elastic deformation of a plate with a hole shows the effectiveness of the proposed refinement scheme.

• Computers and Concrete
• /
• v.31 no.3
• /
• pp.223-239
• /
• 2023

Geometric nonlinear performance simulation and analysis of complex modern buildings and industrial products require high-performance shell elements. Balancing multiple aspects of performance in the one geometric nonlinear analysis element remains challenging. We present a new shell element, flat shell DKMGQ-CR (Co-rotational Discrete Kirchhoff-Mindlin Generalized Conforming Quadrilateral), for linear and geometric nonlinear analysis of both thick and thin shells. The DKMGQ-CR shell element was developed by combining the advantages of high-performance membrane and plate elements in a unified coordinate system and introducing the co-rotational formulation to adapt to large deformation analysis. The effectiveness of linear and geometric nonlinear analysis by DKMGQ-CR is verified through the tests of several classical numerical benchmarks. The computational results show that the proposed new element adapts to mesh distortion and effectively alleviates shear and membrane locking problems in linear and geometric nonlinear analysis. Furthermore, the DKMGQ-CR demonstrates high performance in analyzing thick and thin shells. The proposed element DKMGQ-CR is expected to provide an accurate, efficient, and convenient tool for the geometric nonlinear analysis of shells.