• Earthquakes and Structures
• /
• v.26 no.2
• /
• pp.163-174
• /
• 2024

Dynamic equilibrium equations for finite element analysis were derived for the free field one-dimensional shear wave propagation through the horizontally layered soil deposits with the elastic half-space. We expressed Rayleigh's viscous damping consisting of mass and stiffness proportional terms. We considered two cases where damping matrices are defined in the total and relative displacement fields. Two forms of equilibrium equations are presented; one in terms of total motions and the other in terms of relative motions. To evaluate the performance of new equilibrium equations, we conducted two sets of site response analyses and directly compared them with the exact closed-form frequency domain solution. Results show that the base shear force as earthquake load represents the simpler form of equilibrium equation to be used for the finite element method. Conventional finite element procedure using base acceleration as earthquake load predicts exact solution reasonably well even in soil deposits with unrealistically high damping.

• Geomechanics and Engineering
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2020

Slope displacement and factor of safety (FOS) of a slope are two aspects that reflect the stability of a slope. However, the traditional limit equilibrium (LE) methods only give the result of the slope FOS and cannot be used to solve for the slope displacement. Therefore, developing a LE method to obtain the results of the slope FOS and slope displacement has significance for engineering applications. Based on a force-displacement coupled mode, this work expands the LE Morgenstern-Price (M-P) method. Except for the mechanical equilibrium conditions of a sliding body adopted in the traditional M-P method, the present method introduces a nonlinear model of the shear stress and shear displacement. Moreover, the energy equation satisfied by a sliding body under a small slope displacement is also applied. Therefore, the double solutions of the slope FOS and horizontal slope displacement are established. Furthermore, the flow chart for the expanded LE M-P method is given. By comparisons and analyses of slope examples, the present method has close results with previous research and numerical simulation methods, thus verifying the feasibility of the present method. Thereafter, from the parametric analysis, the following conclusions are obtained: (1) the shear displacement parameters of the soil affect the horizontal slope displacement but have little effect on the slope FOS; and (2) the curves of the horizontal slope displacement vs. the minimum slope FOS could be fitted by a hyperbolic model, which would be beneficial to obtain the horizontal slope displacement for the slope in the critical state.

• Structural Engineering and Mechanics
• /
• v.56 no.1
• /
• pp.137-156
• /
• 2015

In this paper, an analytical solution of displacement, strain and stress field for rotating thick-walled cylinder made of functionally graded material subjected to the uniform external magnetic field and thermal field in plane strain state has been studied. Stress, strain and displacement field as a function of radial coordinates considering magneto-thermo-elasticity are derived analytically. According to the Maxwell electro-dynamic equations, Lorentz force in term of displacement is obtained in cylindrical coordinates. Also, symmetric temperature distribution along the thickness of hollow cylinder is obtained by solving Fourier heat transfer equation in cylindrical coordinates. Using equation of equilibrium and thermo-mechanical constitutive equations associated with Lorentz force, a second-order inhomogeneous differential equation in term of displacement is obtained and will be solved analytically. Except Poisson's ratio, other mechanical properties such as elasticity modulus, density, magnetic permeability coefficient, heat conduction coefficient and thermal expansion coefficient are assumed to vary through the thickness according to a power law. In results analysis, non-homogeneity parameter has been chosen arbitrary and inner and outer surface of cylinder are assumed to be rich metal and rich ceramic, respectively. The effect of rotation, thermal, magnetic field and non-homogeneity parameter of functionally graded material which indicates percentages of cylinder's constituents are studied on displacement, Von Mises equivalent stress and Von Mises equivalent strain fields.

• Journal of Ship and Ocean Technology
• /
• v.8 no.3
• /
• pp.50-60
• /
• 2004

In this paper, a continuum-based design sensitivity analysis (DSA) method is developed for the weakly coupled thermo-elasticity problems. The temperature and displacement fields are described in a common domain. Boundary value problems such as an equilibrium equation and a heat conduction equation in steady state are considered. The direct differentiation method of continuum-based DSA is employed to enhance the efficiency and accuracy of sensitivity computation. We derive design sensitivity expressions with respect to thermal conductivity in heat conduction problem and Young's modulus in equilibrium equation. The sensitivities are evaluated using the finite element method. The obtained analytical sensitivities are compared with the finite differencing to yield very accurate results. Extensive developments of this method are useful and applicable for the optimal design problems incorporating welding and thermal deformation problems.

• Structural Engineering and Mechanics
• /
• v.63 no.2
• /
• pp.251-257
• /
• 2017

A local damage identification method using measured structural static displacement is proposed in this study. Based on the residual force vector deduced from the static equilibrium equation, residual strain energy (RSE) is introduced, which can localize the damage in the element level. In the case of all the nodal displacements are used, the RSE can localize the true location of damage, while incomplete displacement measurements are used, some suspicious damaged elements can be found. A model updating method based on static displacement response sensitivity analysis is further utilized for accurate identification of damage location and extent. The proposed method is verified by two numerical examples. The results indicate that the proposed method is efficient for damage identification. The advantage of the proposed method is that only limited static displacement measurements are needed in the identification, thus it is easy for engineering application.

• Structural Engineering and Mechanics
• /
• v.64 no.2
• /
• pp.205-212
• /
• 2017

Based on the coupled elastic bending deformation features and relationships between the internal force and deformation of pre-twisted Euler beam, the generalized strain, the equivalent constitutive equation and the equilibrium equation of pre-twisted Euler beam are developed. Based on the properties of the dual-antisymmetric matrix, the general solution of pre-twisted Euler beam is obtained. By comparison with ANSYS solution by using straight Beam-188 element based on infinite approach strategy, the results show that the developed method is available for pre-twisted Euler beam and also provide an accuracy displacement interpolation function for the subsequent finite element analysis. The effect of pre-twisted angle on the mechanical property has been investigated.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.2
• /
• pp.143-156
• /
• 2015

The stress and displacement fields at the crack tip were studied during the unsteady propagation of a mode III crack in a direction that was different from the property graduation direction in functionally graded materials (FGMs). The property graduation in FGMs was assumed based on the linearly varying shear modulus under a constant density and the exponentially varying shear modulus and density. To obtain the solution of the harmonic function, the general partial differential equation of the dynamic equilibrium equation was transformed into a Laplace equation. Based on the Laplace equation, the stress and displacement fields, which depended on the time rates of change in the crack tip speed and stress intensity factor, were obtained through an asymptotic analysis. Using the stress and displacement fields, the effects of the angled property variation on the stresses, displacements, and stress intensity factors are discussed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.2
• /
• pp.278-288
• /
• 1998

An efficient and useful method to improve the local stress field in a displacement-formulated finite element solution has been proposed using the theory of conjugate approximations for a stress field and the Loubignac's iterative method for a displacement field. Validity of the proposed method has been tested through three test examples, to improve the stress field and displacement field in the whole domain and the local regions. As a result of analysis on the test examples, it is found that the stress field in the local regions are approximated to those in the whole domain within a few iterations which have satisfied the original finite element equilibrium equation. In addition, it is found that the local stress field are by far better approximated to the exact stress field than the displacement-based stress field with the reduction of the finite-element mesh-size.

• Structural Engineering and Mechanics
• /
• v.48 no.6
• /
• pp.849-878
• /
• 2013

This paper consists of two parts, which broadly examines solution techniques abilities for the structures with geometrical nonlinear behavior. In part I of the article, formulations of several well-known approaches will be presented. These solution strategies include different groups, such as: residual load minimization, normal plane, updated normal plane, cylindrical arc length, work control, residual displacement minimization, generalized displacement control, modified normal flow, and three-parameter ellipsoidal, hyperbolic, and polynomial schemes. For better understanding and easier application of the solution techniques, a consistent mathematical notation is employed in all formulations for correction and predictor steps. Moreover, other features of these approaches and their algorithms will be investigated. Common methods of determining the amount and sign of load factor increment in the predictor step and choosing the correct root in predictor and corrector step will be reviewed. The way that these features are determined is very important for tracing of the structural equilibrium path. In the second part of article, robustness and efficiency of the solution schemes will be comprehensively evaluated by performing numerical analyses.

• Steel and Composite Structures
• /
• v.27 no.1
• /
• pp.13-25
• /
• 2018

The main objective of this study is to develop a dual approach for geometrically nonlinear finite element analysis of plane truss structures. The geometric nonlinearity is considered using the Total Lagrangian formulation. The nonlinear solution is obtained by introducing and minimizing an objective function subjected to displacement-type constraints. The proposed method can fully trace the whole equilibrium path of geometrically nonlinear plane truss structures not only before the limit point but also after it. No stiffness matrix is used in the main approach and the solution is acquired only based on the direct classical stress-strain formulations. As a result, produced errors caused by linearization and approximation of the main equilibrium equation will be eliminated. The suggested algorithm can predict both pre- and post-buckling behavior of the steel plane truss structures as well as any arbitrary point of equilibrium path. In addition, an equilibrium path with multiple limit points and snap-back phenomenon can be followed in this approach. To demonstrate the accuracy, efficiency and robustness of the proposed procedure, numerical results of the suggested approach are compared with theoretical solution, modified arc-length method, and those of reported in the literature.