• ETRI Journal
• /
• v.29 no.1
• /
• pp.8-17
• /
• 2007

In this paper, we present a time domain combined field integral equation formulation (TD-CFIE) to analyze the transient electromagnetic response from dielectric objects. The solution method is based on the method of moments which involves separate spatial and temporal testing procedures. A set of the RWG functions is used for spatial expansion of the equivalent electric and magnetic current densities, and a combination of RWG and its orthogonal component is used for spatial testing. The time domain unknowns are approximated by a set of orthonormal basis functions derived from the Laguerre polynomials. These basis functions are also used for temporal testing. Use of this temporal expansion function characterizing the time variable makes it possible to handle the time derivative terms in the integral equation and decouples the space-time continuum in an analytic fashion. Numerical results computed by the proposed formulation are compared with the solutions of the frequency domain combined field integral equation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.415-423
• /
• 2015

In this study, we present a grid refinement model in the lattice Boltzmann method (LBM) for two-dimensional incompressible fluid flow. That is, the model combines the desirable features of the lattice Boltzmann method and stream function-vorticity formulations. In order to obtain an accurate result, very fine grid (or lattice) is required near the solid boundary. Therefore, the grid refinement model is used in the lattice Boltzmann method for stream function-vorticity formulation. This approach is more efficient in that it can obtain the same accurate solution as that in single-block approach even if few lattices are used for computation. In order to validate the grid refinement approach for the stream function-vorticity formulation, the numerical simulations of lid-driven cavity flows were performed and good results were obtained.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.113-140
• /
• 2004

A mixed eight-node hexahedral element formulated via the Hu-Washizu principle as well as the field extrapolation technique is presented. The mixed element with only three translational degrees of freedom at each node can provide extremely accurate and reliable performance for popular benchmark problems such as spacial beams, plates, shells as well as general three-dimensional elasticity problems. Numerical calculations also show that when extremely skewed and coarse meshes and nearly incompressible materials are used, the proposed mixed element can still possess excellent behaviour. The mixed formulation starts with introduction of a parallelepiped domain associated with the given general eight-node hexahedral element. Then, the assumed strain field at the nodal level is constructed via the Hu-Washizu variational principle for that associated parallelepiped domain. Finally, the assumed strain field at the nodal level of the given hexahedral element is established by using the field extrapolation technique, and then by using the trilinear shape functions the assumed strain field of the whole element domain is obtained. All matrices involved in establishing the element stiffness matrix can be evaluated analytically and expressed explicitly; however, a 24 by 24 matrix has to be inverted to construct the displacement extrapolation matrix. The proposed hexahedral element satisfies the patch test as long as the element with a shape of parallelepiped.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.191-213
• /
• 2011

This paper presents an improved 8-node shell element for the analysis of plates and shells. The finite element, based on a refined first-order shear deformation theory, is further improved by the combined use of assumed natural strain method. We analyze the influence of the shell element with the different patterns of sampling points for interpolating different components of strains. Using the assumed natural strain method with proper interpolation functions, the present shell element generates neither membrane nor shear locking behavior even when full integration is used in the formulation. Further, a refined first-order shear deformation theory, which results in parabolic through-thickness distribution of the transverse shear strains from the formulation based on the third-order shear deformation theory, is proposed. This formulation eliminates the need for shear correction factors in the first-order theory. Numerical examples demonstrate that the present element perform better in comparison with other shell elements.

• Steel and Composite Structures
• /
• v.30 no.1
• /
• pp.57-67
• /
• 2019

This paper deals with the static analysis of axially functionally graded rectangular plates. It is assumed that the flexural rigidity of the plate varies exponentially along one of the plate's in-plane dimensions. Both an analytical approach and a numerical method are utilized to solve the problem. The analytical solution is obtained by using the Green's function method. To employ this approach, the adjoint boundary value problem is established. Then, exact solutions for deflection of the plate for different boundary conditions are found. In another way, a finite element formulation for the problem is developed. In order to demonstrate the validity of the Authors' formulation, the results obtained via both mentioned schemes are compared with each other for functionally graded plates and with results of previously published works for homogeneous plates. The effect of plate parameters on the response of the plate is also investigated. To remind the research background, a brief review on the application of Green's function method in plates' analysis and functionally graded plates is also presented.

• Steel and Composite Structures
• /
• v.38 no.4
• /
• pp.463-476
• /
• 2021

The concept of using Steel-concrete (SC) composite walls as retaining walls has recently been introduced by the authors and their effectiveness of resisting out-of-plane loads has also been demonstrated. In this paper, an improved analytical formulation based on partial interaction theory, which has previously been developed by the authors, is presented. The improved formulation considers a new loading condition and also accounts for cracking in concrete to simulate the real conditions. Due to a limited number of test specimens, further finite element (FE)simulations are performed in order to verify the analytical procedure in more detail. It is observed that the results from the improved analytical procedure are in excellent agreement with both experimental and numerical results. Moreover, a detailed parametric study is conducted using the developed FE model to investigate effects of different parameters, such as distance between shear connectors, shear connector length, concrete strength, steel plate thickness, concrete cover thickness, wall's width to thickness ratio, and wall's height to thickness ratio, on the behavior of SC composite walls subjected to out-of-plane loads.

• Ocean Systems Engineering
• /
• v.11 no.1
• /
• pp.59-81
• /
• 2021

The nonlinear formulation using the principle of virtual work-energy for free vibration of a large-sag extensible catenary riser in two dimensions is presented in this paper. A support at one end is hinged and the other is a free-sliding roller in the horizontal direction. The catenary riser has a large-sag configuration in the static equilibrium state and is assumed to displace with large amplitude to the motion state. The total virtual work of the catenary riser system involves the virtual strain energy due to bending, the virtual strain energy due to axial deformation, the virtual work done by the effective weight, and the inertia forces. The nonlinear equations of motion for two-dimensional free vibration in the Cartesian coordinate system is developed based on the difference between the Euler's equations in the static state and the displaced state. The linear and nonlinear stiffness matrices of the catenary riser are obtained and the eigenvalue problem is solved using the Galerkin finite element procedure. The natural frequencies and mode shapes are obtained. The results are validated with regard to the reference research addressing the accuracy and efficiency of the proposed nonlinear formulation. The numerical results for free vibration and the effect of the nonlinear behavior for catenary riser are presented.

• Ocean Systems Engineering
• /
• v.13 no.2
• /
• pp.147-172
• /
• 2023

This study investigates the influence of loading and inflow conditions on tidal turbine performance from a hydrodynamic and hydroacoustic point of view. A boundary element method is utilized for the former to investigate turbine performance at various loading conditions under zero/non-zero yaw inflow. The boundary element method is selected as it has been selected, tested, and validated to be computationally efficient and accurate for marine hydrodynamic problems. Once the hydrodynamic solutions are obtained, such as the time-dependent surface pressures and periodic motion of the turbine blade, they are taken as the known noise sources for the subsequence hydroacoustic analysis based on the Ffowcs Williams-Hawkings formulation given in a form proposed by Farassat. This formulation is coupled with the boundary element method to fully consider the three-dimensional shape of the turbine and the speed of sound in the acoustic analysis. For validations, a model turbine is taken from a reference paper, and the comparison between numerical predictions and experimental data reveals satisfactory agreement in hydrodynamic performance. Importantly, this study shows that the noise patterns and sound pressure levels at both the near- and far-field are affected by different loading conditions and sensitive to the inclination imposed in the incoming flow.

• Structural Engineering and Mechanics
• /
• v.88 no.1
• /
• pp.83-94
• /
• 2023

An accurate and highly efficient inverse element labelled iPCB is developed based on the inverse finite element method (iFEM) for real-time shape estimation of plane-curved structures (such as arch bridges) utilizing onboard strain data. This inverse problem, named shape sensing, is vital for the design of smart structures and structural health monitoring (SHM) procedures. The iPCB formulation is defined based on a least-squares variational principle that employs curved Timoshenko beam theory as its baseline. The accurate strain-displacement relationship considering tension-bending coupling is used to establish theoretical and measured section strains. The displacement fields of the isoparametric element iPCB are interpolated utilizing nonuniform rational B-spline (NURBS) basis functions, enabling exact geometric modelling even with a very coarse mesh density. The present formulation is completely free from membrane and shear locking. Numerical validation examples for different curved structures subjected to different loading conditions have been performed and have demonstrated the excellent prediction capability of iPCBs. The present formulation has also been shown to be practical and robust since relatively accurate predictions can be obtained even omitting the shear deformation contributions and considering polluted strain measures. The current element offers a promising tool for real-time shape estimation of plane-curved structures.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.8
• /
• pp.411-417
• /
• 2008

An alternative formulation of the Helmholtz integral equation derived to express the pressure field explicitly in terms of the velocity vector of a radiating surface is used to solve acoustic radiation and fluid/structure interaction problems. This formulation, derived for arbitrary sources, is similar in form to the Rayleigh's formula for planar sources. Because the surface pressure field is expressed explicitly as a surface integral of the surface velocity, which can be implemented numerically using standard Gaussian quadratures, there is no need to use BEM to solve a set of simultaneous equations for the surface pressure at the discretized nodes. Furthermore the non-uniqueness problem inherent in methods based on Helmholtz integral equation is avoided. Validation of this formulation is demonstrated for some simple geometries.