• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.5 s.149
• /
• pp.568-577
• /
• 2006

With the increase of ship size and speed, the loading on the propeller is increasing, which in turn increases the rotational speed in the propeller slipstream. The rudder placed in the propeller slip stream is therefore subject to severe cavitation with the increased angle of attack due to the increased rotational induction speed of the propeller. In the present paper the surface panel method, which has been proved useful in predicting the sheet cavitation on the propeller blade, is applied to solve the cavity boundary value problem on the rudder. The problem is then solved numerically by discretizing the rudder and cavity surface elements of the quadrilateral panels with constant strengths of sources and dipoles. The strengths of the singularities are determined satisfying the boundary conditions on the rudder and cavity surfaces. The extent of the cavity, which is unknown a priori, is determined by iterative procedure. Series of numerical experiments are performed increasing the degree of complexity of the rudder geometry and oncoming flows from the simple hydrofoil case to the real rudder in the circumferentially averaged propeller slipstream. Numerical results are presented with experimental results.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.151-159
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory) concept such as the effect of compression softening in cracked concrete, and macroscopic and rotating crack models etc. was presented for the nonlinear behaviour of structural concrete. In this category, tangential or secant material stiffnesses for cracked concrete were also defined and discussed in view of the iterative solution schemes for nonlinear equations. Considering the computational efficiency and the ability of modelling the post-ultimate behaviour as major concerns, the incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Nonlinear Analysis of structrual Concrete by FEM : Monotonic Loading) developed baed on the CFT constitutive relationships and the incremetal solution strategy described enables the predictions of strength and deformation capacities in a full range. crack patterns and their corresponding widths, and yield extents of reinforcement. As the verfication purpose of NASCOM, the prediction of Cervenka's panel test results including the load resistance and the deformation history was made. A limited number of predictions indicate a good correlation in a general sense.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.3
• /
• pp.187-195
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory), considering the effect of compression softening in cracked concrete, and macro-scopic and rotating crack models etc., was presented for the nonlinear behaviour of structural concrete. Considering the computational efficency and the ability of modelling the post-ultimate behaviour as major concerns, the Incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Non1inear Analysis of Structural Concrete by FEM : Monotonic Loading) developed in this way enables the predictions of strength and deformation capacities in a full range, crack patterns and their corresponding widths, and yield extents of reinforcement. As the verification purpose of NASCOM, the predictions were made for Bhide's Panel(PB21) and Leonhardt's deep beam tests. The predicted results shows somewhat stiff behaviour for the panel test, and vice versa for deep beam tests. More refining process would be necessary hereafter in terms of more accurately simulating the effects of tension-stiffening and compression softening in concrete.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.90-100
• /
• 1993

A finite element formulation of vibration control of a laminated plate with piezoelectric sensor/ actuators is presented. Classical lamination theory with the induced strain actuation and Hamilton's principle are used to formulate the equations of motion of the system. The total charge developed on the sensor layer is calculated from the direct piezoelectric equation. The equations of motion and the total charge are discretized with 4 node, 12 degrees of freedom quadrilateral plate bending elements with one electrical degree of freedom. The mass and stiffness of the piezoelectric layer are introduced by treating them as another layer in laminated plate. Piezoelectric sensor/actuators are distributed, but discrete due to the geometry of electrodes. By defining an i.d. number of electrode for each element, modelling of electrodes with variable geometry can be achieved. The static response of a piezoelectric bimorph beam to electrical loading and sensor voltage to given displacement are calculated. For a laminated plate under the negative velocity feedback control, the direct time response by the Newmark-.betha. method and damped frequencies and modal damping ratios by modal state space analysis are derived.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.6 no.1
• /
• pp.3-12
• /
• 1988

The objective of this paper is to study the characteristics of combination method to correct both angle and distance errors simultaneously based on the least square adjustment methods. Changing the standard errors of distance and angle, the simulation errors of triangultion, trilateration, and combination result in some 39.8%, 33.9%, and 26.3% respectively. As the above, combination method shows more consistent accuracy than other methods. When considering the weight factor about error elements with independence, the diminishing rate of simulated average standard error represents a various change in each method. But considering them simultaneously, it shows a remarkable rate of diminishing 75.5%, 74.1%, and 69.2% in each method. And also, by growing the weight factor, accuracy of triangulation method is growing, whereas that of trilateration is diminishing. Therefore, determining the reasonable weight factors of distance and angle errors simultaneously in the analytical combination method, this method is expected to be one of more accurate and more effective methods for determining horizontal positions on the earth.

• Coupled systems mechanics
• /
• v.11 no.5
• /
• pp.439-458
• /
• 2022

Over-coating is one of the most popular engineering practices to strengthen Reinforced Concrete (RC) structures, due to the relative quickness and ease of construction. It consists of an external coat bonded to the outer surface of the structural RC element, either by the use of chemical adhesives, mechanical anchor bolts or simply mortar injection. In contrast to these constructive advantages, the numerical estimation of the bearing capacity of the strengthened reinforced concrete element is still complicated, not only for the complexity of modelling a flexible membrane or plate attached to a quasi-rigid solid, but also for the difficulties that raise of simulating any potential delamination between both materials. For these reasons, the standard engineering calculations used in the practice remain very approximated and clumsy. In this work, we propose the formulation of a new 2D solid-layer finite element capable to link a solid body with a flexible thin layer, as it were the "skin" of the body, allowing the potential delamination between both materials. In numerical terms, this "skin" element is intended to work as a transitional region between a solid body (modelled with a classical formulation of a standard quadrilateral four-nodes element) and a flexible coat layer (modelled with cubic beam element), dealing with the incompatibility of Degrees-Of-Freedom between them (two DOF for the solid and three DOF for the beam). The aim of the solid-layer element is to simplify the mesh construction of the strengthened RC element being aware of two aspects: a) to prevent the inappropriate use of very small solid elements to simulate the coat; b) to improve the numerical estimation of the real bearing capacity of the strengthened element when the coat is attached or detached from the solid body.

• Coupled systems mechanics
• /
• v.12 no.6
• /
• pp.481-501
• /
• 2023

Over-coating is one of the most popular engineering practices to strengthen Reinforced Concrete (RC) structures, due to the relative quickness and ease of construction. It consists of an external coat bonded to the outer surface of the structural RC element, either by the use of chemical adhesives, mechanical anchor bolts or simply mortar injection. In contrast to these constructive advantages, the numerical estimation of the bearing capacity of the strengthened reinforced concrete element is still complicated, not only for the complexity of modelling a flexible membrane or plate attached to a quasi-rigid solid, but also for the difficulties that raise of simulating any potential delamination between both materials. For these reasons, the standard engineering calculations used in the practice remain very approximated and clumsy. In this work, we propose the formulation of a new 2D solid-layer finite element capable to link a solid body with a flexible thin layer, as it were the "skin" of the body, allowing the potential delamination between both materials. In numerical terms, this "skin" element is intended to work as a transitional region between a solid body (modelled with a classical formulation of a standard quadrilateral four-nodes element) and a flexible coat layer (modelled with cubic beam element), dealing with the incompatibility of Degrees-OfFreedom between them (two DOF for the solid and three DOF for the beam). The aim of the solid-layer element is to simplify the mesh construction of the strengthened RC element being aware of two aspects: a) to prevent the inappropriate use of very small solid elements to simulate the coat; b) to improve the numerical estimation of the real bearing capacity of the strengthened element when the coat is attached or detached from the solid body.

• Journal of Korean Society of Steel Construction
• /
• v.19 no.2
• /
• pp.233-245
• /
• 2007

Previous finite element studies have shown that AASHTO Standard load distribution factor (LDF) equations appear to be conservative for longer spans and larger girder spacing, but too permissible for short spans and girder spacing. AASHTO LRFD specification defines the distribution factor equation for girder spacing, span length, slab thickness, and longitudinal stiffness. However, this equation requires an iterative procedure to correctly determine the LDF value due to an initially unknown longitudinal stiffness parameter. This study presents a simplified LDF equation for interior and exterior girders of two-span continuous I-girder bridges that does not require an iterative design procedure. The finite element method was used to investigate the effect of girder spacing, span length, slab thickness, slab width, and spacing and size of bracing. The computer program, GTSTRUDL, was used to idealize the bridge superstructures as the eccentric beam model, the concrete slab by quadrilateral shell elements, steel girders by space frame members, and the composite action between these elements by rigid links. The distribution factors obtained from these analyses were compared with those from the AASHTO Standard and LRFD methods. It was observed through the parametric studies that girder spacing, span length, and slab thickness were the dominant parameters compared with others. The LRFD distribution factor for the interior girder was found to be conservative in most cases, whereas the factor for the exterior girder to be unconservative in longer spans. Furthermore, a regression analysis was performed to develop simplified LDF formulas. The formulas developed in this study produced LDF values that are always conservative to those from the finite element method and are generally smaller than the LDF values obtained from the AASHTO LRFD specification. The proposed simplified equation will assist bridge engineers in predicting the actual LDF in two-span continuous I-girder bridges.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.