• The Journal of the Acoustical Society of Korea
• /
• v.22 no.1E
• /
• pp.11-18
• /
• 2003

Non-singular boundary element method (BEM) codes are developed in acoustics application. The BEM code is then used to calculate unknown boundary surface normal displacements and surface pressures from known exterior near field pressures. And then the calculated surface normal displacements and surface pressures are again applied to the BEM in forward in order to calculate reconstructed field pressures. The initial exterior near field pressures are very well agreed with the later reconstructed field pressures. Only the same number of boundary surface nodes (1178) are used for the initial exterior pressures which are at first calculated by Finite Element Method (FEM) and BEM. Pseudo-inverse technique is, used for the calculation of the unknown boundary surface normal displacements. The structural object is a tuning fork with 128.4 ㎐ resonant. The boundary element is a quadratic hexahedral element (eight nodes per element).

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

Non-singular boundary element method (BEM) codes are developed in acoustics application. The BEM code is then used to calculate unknown boundary surface normal displacements and surface pressures from known exterior near Held pressures. And then the calculated surface normal displacements and surface pressures are again applied to the BEM in forward in order to calculate reconstructed field pressures. The initial exterior near field pressures are very well agreed with the later reconstructed field pressures. Only the same number of boundary surface nodes (1178) are used far the initial exterior pressures which are initially calculated by Finite Element Method (FEM) and BEM. Pseudo-inverse technique is used for the calculation of the unknown boundary surface normal displacements. The structural object is a tuning fork with 128.4 Hz resonant. The boundary element is a quadratic hexahedral element (eight nodes per element).

• Journal of Ocean Engineering and Technology
• /
• v.23 no.5
• /
• pp.1-8
• /
• 2009

This paper is a continuation of the recent development on Hermite-based divergence free element method and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square enclosure with temperature difference across the two sides. The basis functions for the velocity field consist of the Hermite function and its curl while the basis functions for the temperature field consists of the Hermite function and its gradients. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x and y derivatives. This paper presents numerical results for Ra = 105, and compares with those from a stabilized finite element method developed by Illinca et al. (2000). The comparison has been done on 32 by 32 uniform elements and the degree of approximation of elements used for the stabilized finite element are linear (Deg. 1) and quadratic (Deg. 2). The numerical results from both methods show well agreements with those of De vahl Davi (1983).

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.8
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• Journal of Magnetics
• /
• v.16 no.2
• /
• pp.120-124
• /
• 2011

The design of an air compressor-driving quadratic linear actuator in a fuel cell BOP system is studied using orthogonal techniques. The approach utilizes an orthogonal array for design of 'experiments', i.e. the scheme for numerical simulations using a finite element method. Eco-friendly energy is increasingly important due to the depletion of fossil fuels and environmental pollution. Among the new energy sources, fuel cell is spotlighted as renewable energy because it produces few dusts. The air compressor performance is directly related to the efficiency of the fuel cell BOP system has high power consumption. In this paper, an optimized technique using an orthogonal matrix is applied to the design problem to improve the performance of quadratic linear actuator.

• Structural Engineering and Mechanics
• /
• v.69 no.3
• /
• pp.269-281
• /
• 2019

This work addresses a three-dimensional nonlinear structural analysis of the constructive phases of a cable-stayed segmental concrete bridge using The Finite Element Method through ANSYS, version 14.5. New subroutines have been added to ANSYS via its UPF customization tool to implement viscoelastoplastic constitutive equations with cracking capability to model concrete's structural behavior. This numerical implementation allowed the use of three-dimensional twenty-node quadratic elements (SOLID186) with the Element-Embedded Rebar model option (REINF264), conducting to a fast and efficient solution. These advantages are of fundamental importance when large structures, such as bridges, are modeled, since an increasing number of finite elements is demanded. After validating the subroutines, the bridge located in Rio de Janeiro, Brazil, and known as "Ponte do Saber" (Bridge of Knowledge, in Portuguese), has been numerically modeled, simulating each of the constructive phases of the bridge. Additionally, the data obtained numerically is compared with the field data collected from monitoring conducted during the construction of the bridge, showing good agreement.

• Computational Structural Engineering
• /
• v.5 no.2
• /
• pp.105-118
• /
• 1992

Static performance was compared for the triangular plate elements through some numerical experiments. Four Kirchhoff elements and six Mindlin elements were selected for the comparison. Numerical tests were executed for the problems of rectangular plates with regular and distorted meshes, rhombic plates, circular plates and cantilever plates. Among the Kirchhoff 9 DOF elements, the discrete Kirchhoff theory element was the best. Element distortion and the aspect ratio were shown to have negligible effects on the displacement behaviour. The Specht's element resulted in better results than the Bergan's but it was sensitive to the aspect ratio. The element based on the hybrid stress method also resulted in good results but it assumed to be less reliable. Among the linear Mindlin elements, the discrete shear triangle was the best in view of reliability, accuracy and convergence. Since the thin plate behaviour of it was as good as the DKT element, it can be used effectively in the finite element code regardless of the thickness. As a quadratic Mindlin element, the MITC7 element resulted in best results in almost all cases considered. The results were at least as good as those of doubly refined meshes of linear elements.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.291-300
• /
• 1998

This paper demonstrates how ambient vibration measurements at a limited number of locations can be effectively utilized to estimate parameters of a finite element model of a large-scale structural system involving a large number of elements. System identification using ambient vibration measurements presents a challenge requiring the use of special identification techniques, which ran deal with very small magnitudes of ambient vibration contaminated by noise without the knowledge of input farces. In the present study, the modal parameters such as natural frequencies, damping ratios, and mode shapes of the structural system were estimated by means of appropriate system identification techniques including the random decrement method. Moreover, estimation of parameters such as the stiffness matrix of the finite element model from the system response measured by a limited number of sensors is another challenge. In this study, the system stiffness matrix was estimated by using the quadratic optimization involving the computed and measured modal strain energy of the system, with the aid of a sensitivity relationship between each element stiffness and the modal parameters established by the second order inverse modal perturbation theory. The finite element models thus identified represent the actual structural system very well, as their calculated dynamic characteristics satisfactorily matched the observed ones from the ambient vibration test performed on a large-scale structural system subjected primarily to ambient wind excitations. The dynamic models identified by this study will be used for design of an active mass damper system to be installed on this structure fer suppressing its wind vibration.

• International Journal of Fuzzy Logic and Intelligent Systems
• /
• v.5 no.2
• /
• pp.119-123
• /
• 2005

This paper describes a fuzzy-based system for analyzing the stress intensity factors (SIFs) of three-dimensional (3D) cracks. 3D finite element method(FEM) was used to obtain the SIF for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy theory. Nodes are generated by the bucketing method, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. The singular elements such that the mid-point nodes near crack front are shifted at the quarter-points, and these are automatically placed along the 3D crack front. The complete FE model is generated, and a stress analysis is performed. The SIFs are calculated using the displacement extrapolation method. The results were compared with those surface cracks in homogeneous materials. Also, this system is applied to analyze cladding effect of surface cracks in inhomogeneous materials.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.587-614
• /
• 2014

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.