• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.281-295
• /
• 2005

An analytical method is presented to solve the elastodynamic problem of finitely long hollow cylinder subjected to torsional impact often occurs in engineering mechanics. The analytical solution is composed of a solution of quasi-static equation satisfied with the non-homogeneous boundary condition and a solution of dynamic equation satisfied with homogeneous boundary condition. The quasi-static solution is obtained directly by solving the quasi-static equation satisfied with the non-homogeneous boundary condition. The solution of the non-homogeneous dynamic equation is obtained by means of finite Hankel transform on the radial variable, r, Laplace transform on time variable, t, and finite Fourier transform on axial variable, z. Thus, the solution for finitely long, hollow cylinder subjected to torsion impact is obtained. In the calculating examples, the response histories and distributions of shear stress in the finitely long hollow cylinder subjected to an exponential decay torsion load are obtained, and the results have been analyzed and discussed. Finally, a dynamic finite element for the same problem is carried out by using ABAQUS finite element analysis. Comparing the analytical solution with the finite element solution, it can be found that two kinds of results obtained by means of two different methods agree well. Therefore, it is further concluded that the analytical method and computing process presented in the paper are effective and accurate.

• 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.

• Structural Engineering and Mechanics
• /
• v.43 no.5
• /
• pp.561-582
• /
• 2012

Despite popularity of FEM in analysis of static and dynamic structural problems and the routine applicability of FE softwares, analytical methods based on simple mathematical relations is still largely sought by many researchers and practicing engineers around the world. Development of such analytical methods for analysis of free vibration of non-prismatic beams is also of primary concern. In this paper a new and simple method is proposed for determination of vibration frequencies of non-prismatic beams under variable axial forces. The governing differential equation is first obtained and, according to a harmonic vibration, is converted into a single variable equation in terms of location. Through repetitive integrations, integral equation for the weak form of governing equation is derived. The integration constants are determined using the boundary conditions applied to the problem. The mode shape functions are approximated by a power series. Substitution of the power series into the integral equation transforms it into a system of linear algebraic equations. Natural frequencies are determined using a non-trivial solution for system of equations. Presented method is formulated for beams having various end conditions and is extended for determination of the buckling load of non-prismatic beams. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained to those obtained using available finite element software.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.280-285
• /
• 2001

As an integral part of the probabilistic fracture mechanics analysis, stress intensity factor calculation scheme for semi-elliptical surface flaws in thin-walled cylinder has been introduced. The approximation solution utilizes the influence coefficients to calculate the stress intensity factor at the crack tip. This method has been compared with other solution methods including 3-D finite element analysis for cooldown boundary condition. The analysis results confirmed that the simplified methods provided sufficiently accurate stress intensity factor values for axial semi-elliptcal flaws on the surface of the reactor pressure vessel.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.

• International Journal of CAD/CAM
• /
• v.3 no.1_2
• /
• pp.19-29
• /
• 2003

Disassembly is a fundamental process needed for component reuse and material recycling in all assembled products. Integral attachments, also known as 'snap' fits, are favored fastening means in design for assembly (DFA) methodologies, but not necessarily a favored choice for design for disassembly. In this paper, design methods of a new class of integral attachments are proposed, where the snapped joints can be disengaged by the application of localized heat sources. The design problem of reversible integral attachments is posed as the design of compliant mechanisms actuated with localized thermal expansion of materials. Topology optimization technique is utilized to obtain conceptual layout of snap-fit mechanisms that realizes a desired deformation of snapped features for joint release. Two design approaches are attempted and design results of each approach are presented, where the geometrical configuration extracted from optimal topologies are simplified to enhance the manufacturability for the conventional injection molding technologies. To maximize the magnitude of deformation, a design scheme has been proposed to include boundary conditions as design variables. Final designs are verified using commercial software for finite element analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.7
• /
• pp.881-889
• /
• 2010

A volume integral equation method (VIEM) is introduced for solving the elastostatic problems related to an unbounded isotropic elastic solid; this solid is subjected to remote uniaxial tension, and it contains multiple interacting isotropic inclusions. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central inclusion is carried out; square and hexagonal packing of the inclusions are considered. The effects of the number of isotropic inclusions and different fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy and efficiency of the method are clarified by comparing the results obtained by analytical and finite element methods. The VIEM is shown to be very accurate and effective for investigating the local stresses in composites containing isotropic fibers.

• Structural Engineering and Mechanics
• /
• v.4 no.5
• /
• pp.469-484
• /
• 1996

Green's functions are obtained in exact closed-forms for the elastic fields in bi-material elastic solids with slipping interface and differing transversely isotropic properties induced by concentrated point and ring force vectors. For the concentrated point force vector, the Green functions are expressed in terms of elementary harmonic functions. For the concentrated ring force vector, the Green functions are expressed in terms of the complete elliptic integral. Numerical results are presented to illustrate the effect of anisotropic bi-material properties on the transmission of normal contact stress and the discontinuity of lateral displacements at the slipping interface. The closed-form Green's functions are systematically presented in matrix forms which can be easily implemented in numerical schemes such as boundary element methods to solve elastic problems in computational mechanics.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.05a
• /
• pp.1038-1043
• /
• 2003

There are many industrial applications including thin-body structures such as fins. For the numerical modeling of radiation of sound from thin bodies, the conventional boundary element method (BEM) using the Helmholtz integral equation fails to yield a reliable solution. Therefore, many researchers have tried to solve the thin-body acoustic problems. In the area of the design sensitivity analysis (DSA) and optimization methods, however, there has been just a few study reported. Especially fur the thin-body acoustics, however, no further study in the DSA and optimization fields has been reported. In this research, the normal derivative integral equation is adopted as an analysis formulation in the thin-body acoustics, and then used for the sizing DSA and optimization. Since the gradient-based method is used for the optimization, it is important to have accurate gradients (design sensitivities) of the objective function and constraints with respect to the design variables. The DSA formulations are derived through chain-ruled derivatives using the finite element method (FEM) and BEM by using the direct differentiation and continuum variation concepts. The proposed approaches are implemented and validated using a numerical example.