• Journal of Korean Association for Spatial Structures
• /
• v.3 no.2 s.8
• /
• pp.85-90
• /
• 2003

If we use a third order approximation for the displacement function of beam element in finite element methods, finite element solutions of beams yield nodal displacement values matching to beam theory results to have no connection with the number increasing of elements of beams. It is assumed that, as the member displacement value at beam nodes are correct, the calculation procedure of beam element stiffness matrix have no numerical errors. A the member forces are calculated by the equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$, the member forces at nodes of beams have errors in a moment and a shear magnitudes in the case of smaller number of element. The nodal displacement value of plate subject to the lateral load converge to the exact values according to the increase of the number of the element. So it is assumed that the procedures of plate element stiffness matrix calculations has a error in the fundamental assumptions. The beam methods for the high accuracy ratio solution Is also applied to the plate analysis. The method of reducing a error ratio of member forces and element stiffness matrix in the finite element methods is studied. Results of study were as follows. 1. The matrixes of EI[B] and [K] in the equations of M(x)=EI[B]{q} and M(x) = [K]{q}+{Q} of beams are same. 2. The equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$ for the member forces have a error ratio in a finite element method of uniformly loaded structures, so equilibrium node loads {Q} must be substituted in the equation of member forces as the numerical examples of this paper revealed.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2004.04a
• /
• pp.209-214
• /
• 2004

This paper deals with tipping over of hydraulic excavator's crane work. If excavator lift too heavy weight, excavator tipped up. This is 38% of whole excavator accidents. In this paper, tipping over load which is maximum load of excavator can lift with displacement of excavator links, real load and tipping over rate are calculated with Zero Moment Point. We designed the tipping-over stability criterion algorithm considering the dynamic characteristics to which ZMP theory is applied and discussed the usefulness of the proposed algorithm compared with the moment equilibrium equation through the simulation and the actual test.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.1926-1932
• /
• 2001

In this second part of the paper, the automatic mesh generation and remeshing algorithm using bubble packing method is applied to the nonlinear problem. The remeshing/refinement procedure is necessary in the large deformation process especially because the mesh distortion deteriorates the convergence and accuracy. To perform the nonliear analysis, the transfer of state variables such as displacement and strain is added to the algorithm of Part 1. The equilibrium equation based on total Lagrangian formulation and elasto-viscoplastic model is used. For the numerical experiment, the upsetting process including the contact constraint condition is analyzed by two refinement criteria. And from the result, it is addressed that the present algorithm can generate the refined meshes easily at the largely deformed area with high error.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.11 no.7
• /
• pp.275-286
• /
• 2001

The research presents an analytical theory to calculate the characteristics of the bal bearing with waviness in its rolling elements considering the centrifugal force and gyroscopic moment of bal. The effects of centrifugal force and gyroscopic moment are introduced to the kinematic constraints and force equilibrium equations. and the waviness of rolling elements is modeled by sinusoidal function to calculate the contact force at each ball. The numerical solutions of governing equation of berating due to waviness are calculated by using the Newton-Raphson method. The accuracy of the research is validated by comparing the contact force. contact angle in case of considering the centrifugal force and gyroscopic moment of bal and the contact force and vibration frequencies in cases of considering waviness with the prior researches respectively. It investigates the stiffness, contact force. displacement and vibration frequencies of the ball bearing considering not only the centrifugal force and gyroscopic moment of ball but also the waviness of the rolling elements.

• Computational Structural Engineering
• /
• v.8 no.2
• /
• pp.147-151
• /
• 1995

The time-domain boundary element analysis has been attempted for investigating the displacement and stress in an elastic-viscoelastic compound structure. The subdomain technique has been employed and the whole domain has been divided into two subdomains, which are respectively a homogeneous elastic zone and a homogeneous viscoelastic zone. The boundary element equation has been formulated using the equilibrium and continuity conditions at the common interface. The numerical results of example problem has been presented.

• Structural Engineering and Mechanics
• /
• v.60 no.2
• /
• pp.271-299
• /
• 2016

In this article, based on the higher-order shear deformation plate theory, buckling analysis of a rectangular plate made of functionally graded piezoelectric materials and its effective parameters are investigated. Assuming the transverse distribution of electric potential to be a combination of a parabolic and a linear function of thickness coordinate, the equilibrium equations for the buckling analysis of an FGP rectangular plate are established. In addition to the Maxwell equation, all boundary conditions including the conditions on the top and bottom surfaces of the plate for closed and open circuited are satisfied. Considering double sine solution (Navier solution) for displacement field and electric potential, an analytical solution is obtained for full simply supported boundary conditions. The accurate buckling load of FGP plate is presented for both open and closed circuit conditions. It is found that the critical buckling load for open circuit is more than that of closed circuit in all loading conditions. Furthermore, it is observed that the influence of dielectric constants on the critical buckling load is more than those of others.

• Structural Engineering and Mechanics
• /
• v.53 no.6
• /
• pp.1167-1182
• /
• 2015

In this study, the stresses and electric potential redistributions of a cylinder made from functionally graded piezoelectric material (FGPM) are investigated. All the mechanical, thermal and piezoelectric properties are modeled as power-law distribution of volume fraction. Using the coupled electro-thermo-mechanical relations, strain-displacement relations, Maxwell and equilibrium equations are obtained including the time dependent creep strains. Creep strains are time, temperature and stress dependent, the closed form solution cannot be found for this constitutive differential equation. A semi-analytical method in conjunction with the Mendelson method of successive approximation is therefore proposed for this analysis. Similar to the radial stress histories, electric potentials increase with time, because the latter is induced by the former during creep deformation of the cylinder, justifying industrial application of such a material as efficient actuators and sensors.

• Steel and Composite Structures
• /
• v.12 no.1
• /
• pp.1-11
• /
• 2012

In this paper static analysis of FGM cylinders subjected to internal and external pressure was carried out by a mesh-free method. In this analysis MLS shape functions are used for approximation of displacement field in the weak form of equilibrium equation and essential boundary conditions are imposed by transformation method. Mechanical properties of cylinders were assumed to be variable in the radial direction. Two types of cylinders were analyzed in this work. At first cylinders with infinite length were considered and results obtained for these cylinders were compared with analytical solutions and a very good agreement was seen between them. Then the proposed mesh-free method was used for analysis of cylinders with finite length and two different types of boundary conditions. Results obtained from these analyses were compared with results of finite element analyses and a very good agreement was seen between them.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.3
• /
• pp.237-244
• /
• 2016

This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.