• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.39-46
• /
• 1993

The need to develop method which can improve the shape design efficiency using high quality approximation is being brought up. In this study, to perform shape optimal design of three-dimensional continuum structures an efficient approximation method for stress constraints is proposed, based on expanding the nodal forces in Taylor series with respect to shape variables. Numerical examples are performed using the 3-D cantilever beam and fixed-fixed beam and compared with other method to demonstrate the efficiency and convergence rate of the Force Approximation method. It is shown that by taking advantage of this high quality approximation, the total number of finite element analysis required for shape optimization of 3-D continuum structures can be reduced significantly, resulting to the same level of efficiency achieved previously in sizing optimization problems. Also, shape representation by super curve technique applied to obtain optimal shape finds useful method.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.133-140
• /
• 1990

Analysis of cable structures is complex because their force - displacement relationships are highly nonlinear and also because large deformations introduce geometric nonlinearity. Therefore, we must take account their geometric nonlinearity in the analysis and find the equilibrated shape of cable structures. In this paper, to slove these problems, numerical procedures involving geometrical nonlinearity are introduced. They are applicable to general cable net, flexible transmission lines and suspended cable roof. These procedures are divided into two parts; one is to obtain the equilibrated shapes and stresses of the cable structures with uniform load by flexibility iteration method, the other is to analyse the equilibrated structures subjected to nodal external forces by nonlinear finite element method.

• Journal of Power System Engineering
• /
• v.17 no.1
• /
• pp.50-57
• /
• 2013

This describes the formulation for the free vibration of joined conical-cylindrical shells with uniform thickness using the transfer of influence coefficient. This method was developed based on successive transmission of dynamic influence coefficients, which were defined as the relationships between the displacement and the force vectors at arbitrary nodal circles of the system. The two edges of the shell having arbitrary boundary conditions are supported by several elastic springs with meridional/axial, circumferential, radial and rotational stiffness, respectively. The governing equations of vibration of a conical shell, including a cylindrical shell, are written as a coupled set of first order differential equations by using the transfer matrix of the shell. Once the transfer matrix of a single component has been determined, the entire structure matrix is obtained by the product of each component matrix and the joining matrix. The natural frequencies and the modes of vibration were calculated numerically for joined conical-cylindrical shells. The validity of the present method is demonstrated through simple numerical examples, and through comparison with the results of previous researchers.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.6
• /
• pp.1004-1014
• /
• 2001

The fully automatic algorithm from initial finite element mesh generation to remeshing in two dimensional geometry is introduced using bubble packing method (BPM) for finite element analysis. BPM determines the node placement by force-balancing configuration of bubbles and the triangular meshes are made by Delaunay triangulation with advancing front concept. In BPM, we suggest two node-search algorithms and the adaptive/recursive bubble controls to search the optimal nodal position. To use the automatically generated mesh information in FEA, the new enhanced bandwidth minimization scheme with high efficiency in CPU time is developed. In the remeshing stage, the mesh refinement is incorporated by the control of bubble size using two parameters. And Superconvergent Patch Recovery (SPR) technique is used for error estimation. To verify the capability of this algorithm, we consider two elasticity problems, one is the bending problem of short cantilever beam and the tension problem of infinite plate with hole. The numerical results indicate that the algorithm by BPM is able to refine the mesh based on a posteriori error and control the mesh size easily by two parameters.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.9
• /
• pp.609-616
• /
• 2022

The performance of antenna reflectors crucially depends on the faceting error of the surface. The force density method (FDM) has been widely used for the form-finding analysis of the cable nets of reflectors. However, after performing form-finding of some cable nets, the effective reflective area will decrease. In addition, nonlinear deformations of the cable can not be achieved by using the FDM. Thus, an effective form-find methodology is proposed in this research. The whole parts of the cable networks are described by the absolute nodal coordinate formulation. The form-finding analysis of the reflector with standard configuration is performed to validate the proposed methodology. The influence of boundary condition changes on the configuration accuracy of the cable net is investigated.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.28 no.2
• /
• pp.131-136
• /
• 2008

Using ultrasonic standing waves, micro particles submerged or flowing in fluid can be manipulated. Due to acoustic radiation force of ultrasound, particles are forced to move to pressure nodal or antinodal lines. In this work, we propose a method to control the position of micro particle in a flow by adjusting the frequency of the standing wave. To this end, standing wave field generation system including a few millimeter thick micro channel was established using an immersible ultrasonic transducer. The present generation system works valid in a frequency range between 2.0 MHz and 2.5 MHz. We observed the SiC particles in water moved to pressure nodal lines by the standing wave. The effect of the channel thickness and operating frequency was also investigated. Interestingly, it was shown that the operating frequency have a close relation with the location of the pressure nodal line. Consequently, it fan be said that the position of particle movement rail be controlled by adjusting the ultrasound frequency. The maximum range of the controllable position was about 261 micrometers under the given condition. The resulted observations reveal the possibility of various applications of the ultrasonic standing wave to the manipulation of particles submerged in a fluid.

• Structural Engineering and Mechanics
• /
• v.58 no.1
• /
• pp.59-91
• /
• 2016

Laminated composite shells are commonly used in various engineering applications including aerospace and marine structures. In this paper, using semi-analytical finite strip method, the buckling behavior of laminated composite deep as well as thick shells of revolution under follower forces which remain normal to the shell is investigated. The stiffness caused by pressure is calculated for the follower forces subjected to external fibers in thick shells. The shell is divided into several closed strips with alignment of their nodal lines in the circumferential direction. The governing equations are derived based on first-order shear deformation theory which accounts for through thickness-shear flexibility. Displacements and rotations in the middle surface of shell are approximated by combining polynomial functions in the meridional direction as well as truncated Fourier series with an appropriate number of harmonic terms in the circumferential direction. The load stiffness matrix which accounts for variation of loads direction will be derived for each strip of the shell. Assembling of these matrices results in global load stiffness matrix which may be un-symmetric. Upon forming linear elastic stiffness matrix called constitutive stiffness matrix, geometric stiffness matrix and load stiffness matrix, the required elements for the second step analysis which is an eigenvalue problem are provided. In this study, different parameter effects are investigated including shell geometry, material properties, and different boundary conditions. Afterwards, the outcomes are compared with other researches. By considering the results of this article, it can be concluded that the deformation-dependent pressure assumption can entail to decrease the calculated buckling load in shells. This characteristic is studied for different examples.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.24-34
• /
• 1990

In this paper, the formulation of a new simplified finite element is made to analyze the ultimate strength of damaged tubular members subjected to combined axial force and end moment. A damaged tubular member that has the bending deformation and the local dent is modeled by beam elements. Tangent elastic stiffness matrix of a beam element which contains the effect of the geometric nonlinearity is derived by using the updated Lagrangian approach. Here the contribution of the stiffness in the dented area is neglected since its resistance against the external loads is considered to be small. A fully plastic interaction curve of the element under combined loads taking account of the local dent effect is selected as a yielding criterion at each nodal point. Also tangent elasto-plastic stiffness matrix of the element is formulated by plastic node method. Comparison with the present solution and the existing experimental results is made showing that the present method gives quite an accurate solution.

• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.747-769
• /
• 2020

This paper provides a semi-analytical approach to investigate the variations of 3D displacement components, electric potential, stresses, electric displacements and transverse vibration frequencies in laminated piezoelectric composite plates based on the scaled boundary finite element method (SBFEM) and the precise integration algorithm (PIA). The proposed approach can analyze the static and dynamic responses of multilayered piezoelectric plates with any number of laminae, various geometrical shapes, boundary conditions, thickness-to-length ratios and stacking sequences. Only a longitudinal surface of the plate is discretized into 2D elements, which helps to improve the computational efficiency. Comparing with plate theories and other numerical methods, only three displacement components and the electric potential are set as the basic unknown variables and can be represented analytically through the transverse direction. The whole derivation is built upon the three dimensional key equations of elasticity for the piezoelectric materials and no assumptions on the plate kinematics have been taken. By virtue of the equilibrium equations, the constitutive relations and the introduced set of scaled boundary coordinates, three-dimensional governing partial differential equations are converted into the second order ordinary differential matrix equation. Furthermore, aided by the introduced internal nodal force, a first order ordinary differential equation is obtained with its general solution in the form of a matrix exponent. To further improve the accuracy of the matrix exponent in the SBFEM, the PIA is employed to make sure any desired accuracy of the mechanical and electric variables. By virtue of the kinetic energy technique, the global mass matrix of the composite plates constituted by piezoelectric laminae is constructed for the first time based on the SBFEM. Finally, comparisons with the exact solutions and available results are made to confirm the accuracy and effectiveness of the developed methodology. What's more, the effect of boundary conditions, thickness-to-length ratios and stacking sequences of laminae on the distributions of natural frequencies, mechanical and electric fields in laminated piezoelectric composite plates is evaluated.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.