• Earthquakes and Structures
• /
• 제26권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
• /
• 제23권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
• /
• 제56권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
• /
• 제8권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
• /
• 제63권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
• /
• 제64권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.

• 대한기계학회논문집A
• /
• 제39권2호
• /
• pp.143-156
• /
• 2015

함수구배재료의 모드 III 균열이 물성치 구배방향과 다른 방향으로 비정상적으로 전파할 때 전파균열선단부근의 응력 및 변위장에 대하여 연구하였다. 함수구배재료는 밀도가 일정한 상태에서 전단탄성계수가 선형적으로 변화하는 경우와 밀도와 전단탄성계수가 지수형적으로 변화하는 경우로 가정했다. 조화함수의 해를 얻기 위하여 일반적인 편미분방정식의 동적평형방정식을 라플라스 방정식으로 변환하였다. 라플라스 방정식으로부터 균열속도 변화률, 응력확대계수의 변화률 등에 의존되는 응력장과 변위장을 근접해법으로 얻었다. 본 연구에서 얻어진 응력장과 변위장을 사용하여 재료의 비 균질성, 균열속도의 변화률, 응력확대계수의 변화률 등을 고려한 상태에서 균열이 임의의 방향으로 전파할 때 균열선단부근의 응력 및 변위 그리고 응력확대계수에 대하여 연구하였다.

• 대한기계학회논문집A
• /
• 제22권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
• /
• 제48권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
• /
• 제27권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.