• Journal of Theoretical and Applied Mechanics
• /
• 제3권1호
• /
• pp.45-65
• /
• 2002

A constitutive model on oorthotropic thermo-elasto-viscoplasticity for fiber-reinforced composite materials Is illustrated, and their thermomechanical responses are predicted with the fully-coupled finite element formulation. The unmixing-mixing scheme can be adopted with the multipartite matrix method as the constitutive model. Basic assumptions based upon the composite micromechanics are postulated, and the strain components of thermal expansion due to temperature change are included In the formulation. Also. more than two sets of mechanical variables, which represent the deformation states of multipartite matrix can be introduced arbitrarily. In particular, the unmixing-mixing scheme can be used with any well-known isotropic viscoplastic theory of the matrix material. The scheme unnecessitates the complex processes for developing an orthotropic viscoplastic theory. The governing equations based on fully-coupled thermomechanics are derived with constitutive arrangement by the unmixing-mixing concept. By considering some auxiliary conditions, the Initial-boundary value problem Is completely set up. As a tool of numerical analyses, the finite element method Is used with isoparametric Interpolation fer the displacement and the temperature fields. The equation of mutton and the energy conservation equation are spatially discretized, and then the time marching techniques such as the Newmark method and the Crank-Nicolson technique are applied. To solve the ultimate nonlinear simultaneous equations, a successive iteration algorithm is constructed with subincrementing technique. As a numerical study, a series of analyses are performed with the main focus on the thermomechanical coupling effect in composite materials. The progress of viscoplastic deformation, the stress-strain relation, and the temperature History are careful1y examined when composite laminates are subjected to repeated cyclic loading.

• Geomechanics and Engineering
• /
• 제5권2호
• /
• pp.143-164
• /
• 2013

Calculation of seismic displacements in reinforced slopes plays a crucial role in appropriate design of these structures however current analytical methods result indifferent values for permanent displacements of the slope. In this paper, based on limit equilibrium and using the horizontal slices method, a new formulation has been proposed for estimating the seismic displacements of a reinforced slope under earthquake records. In this method, failure wedge is divided into a number of horizontal slices. Assuming linear variations for tensile forces of reinforcements along the height of the slope, the coefficient of yield acceleration has been estimated. The simplicity of calculations and taking into account the frequency content of input triggers are among the advantages of the present formulation. Comparison of the results shows that the yield acceleration calculated by the suggested method is very close to the values resulted from other techniques. On the other hand, while there is a significant difference between permanent displacements, the values obtained from the suggested method place somehow between those calculated by the other techniques.

• Advances in aircraft and spacecraft science
• /
• 제1권3호
• /
• pp.291-310
• /
• 2014

An advanced model for the linear flutter analysis is introduced in this paper. Higher-order beam structural models are developed by using the Carrera Unified Formulation, which allows for the straightforward implementation of arbitrarily rich displacement fields without the need of a-priori kinematic assumptions. The strong form of the principle of virtual displacements is used to obtain the equations of motion and the natural boundary conditions for beams in free vibration. An exact dynamic stiffness matrix is then developed by relating the amplitudes of harmonically varying loads to those of the responses. The resulting dynamic stiffness matrix is used with particular reference to the Wittrick-Williams algorithm to carry out free vibration analyses. According to the doublet lattice method, the natural mode shapes are subsequently used as generalized motions for the generation of the unsteady aerodynamic generalized forces. Finally, the g-method is used to conduct flutter analyses of both isotropic and laminated composite lifting surfaces. The obtained results perfectly match those from 1D and 2D finite elements and those from experimental analyses. It can be stated that refined beam models are compulsory to deal with the flutter analysis of wing models whereas classical and lower-order models (up to the second-order) are not able to detect those flutter conditions that are characterized by bending-torsion couplings.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.177-189
• /
• 2015

This paper presents an all-direct domain decomposition approach for large-scale structural analysis. The proposed approach achieves computational robustness and efficiency by enforcing the compatibility of the displacement field across the sub-domain boundaries via local Lagrange multipliers and augmented Lagrangian formulation (ALF). The proposed domain decomposition approach was compared to the existing FETI approach in terms of the computational time and memory usage. The parallel implementation of the proposed algorithm was described in detail. Finally, a preliminary validation was attempted for the proposed approach, and the numerical results of two- and three-dimensional problems were compared to those obtained through a dual-primal FETI approach. The results indicate an improvement in the performance as a result of the implementing the proposed approach.

• Steel and Composite Structures
• /
• 제29권3호
• /
• pp.363-377
• /
• 2018

An efficient and free of shear locking finite element model is developed and employed to study free vibration of tapered bidirectional functionally graded material (BFGM) beams. The beam material is assumed to be formed from four distinct constituent materials whose volume fraction continuously varies along the longitudinal and thickness directions by power-law functions. The finite element formulation based on the first-order shear deformation theory is derived by using hierarchical functions to interpolate the displacement field. In order to improve efficiency and accuracy of the formulation, the shear strain is constrained to constant and the exact variation of the cross-sectional profile is employed to compute the element stiffness and mass matrices. A comprehensive parametric study is carried out to highlight the influence of the material distribution, the taper and aspect ratios as well as the boundary conditions on the vibration characteristics. Numerical investigation reveals that the proposed model is efficient, and it is capable to evaluate the natural frequencies of BFGM beams by using a small number of the elements. It is also shown that the effect of the taper ratio on the fundamental frequency of the BFGM beams is significantly influenced by the boundary conditions. The present results are of benefit to optimum design of tapered FGM beam structures.

• 한국지반환경공학회 논문집
• /
• 제2권3호
• /
• pp.89-99
• /
• 2001

본 연구에서는 가장 널리 사용되고 있는 지반모델인 Modified Cam-Clay 모델을 이용하여 피에조콘 관입 및 소산시험의 수치해석을 수행하였다. Modified Cam-Clay 모델 및 관련 유한요소 식들을 피에조콘 관입의 대변형 현상을 고려하기 위하여 Updated Lagrangian frame에서 formulation 하였다. 유한요소 해석 결과 얻어진 콘 관입저항치, 간극수압 및 소산곡선을 양산지역 현장 시험결과와 비교 분석하였다. 수치해석 결과는 시험결과와 비교적 잘 부합하는 것으로 고찰되었으나 보다 현실에 근접한 simulation을 위하여 연속적인 깊은 관입의 적절한 수치해석적 modeling이 요구되는 것으로 고찰되었다.

• 대한기계학회논문집A
• /
• 제30권1호
• /
• pp.18-25
• /
• 2006

The purpose of this research work is to demonstrate a successful application of hybrid-mixed formulation and nodeless degrees of freedom in developing a very accurate in-plane curved beam element for free vibration analysis. To resolve the numerical difficulties due to the spurious constraints, the present element, based on the Hellinger-Reissner variational principle and considering the effect of shear deformation, employed consistent stress parameters corresponding to cubic displacement polynomials with additional nodeless degrees. The stress parameters were eliminated by the stationary condition, and the nodeless degrees were condensed by Guyan Reduction. Several numerical examples indicated that the property of the mass matrix as well as that of the stiffness matrix have a great effect on the numerical performance. The element with consistent mass matrix produced best results on convergence and accuracy in the numerical analysis of Eigenvalue problems. Also, the higher-order mixed curved beam element showed a superior numerical behavior for the free vibration analyses.

• Coupled systems mechanics
• /
• 제2권4호
• /
• pp.375-387
• /
• 2013

Soils consist of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled with water and air. Hence, soil behaviour must be analyzed by incorporating the effects of the transient flow of the pore-fluid through the voids, and therefore requires a two-phase continuum formulation for saturated porous media. The present paper presents briefly the Biot's basic theory of dynamics of saturated porous media with u-P formulation to determine the responses of pore fluid and soil skeleton during cyclic loading. Kelvin elements are attached to transmitting boundary. The Pastor-Zienkiewicz-Chan model has been used to describe the inelastic behavior of soils under isotropic cyclic loadings. Newmark-Beta method is employed to discretize the time domain. The response of fluid-saturated porous media which are subjected to time dependent loads has been simulated numerically to predict the liquefaction potential of a semi-infinite saturated sandy layer using finite-infinite elements. A settlement of 17.1 cm is observed at top surface. It is also noticed that liquefaction occurs at shallow depth. The mathematical advantage of the coupled finite element analysis is that the excess pore pressure and displacement can be evaluated simultaneously without using any empirical relationship.

• Structural Engineering and Mechanics
• /
• 제21권5호
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• Structural Engineering and Mechanics
• /
• 제57권4호
• /
• pp.733-762
• /
• 2016

Based on a reduced displacement field, a layer-wise (LW) formulation is developed for analysis of thick shell panels which is subjected to axial tension. Employing the principle of minimum total potential energy, the local governing equations of thick panel which is subjected to axial extension are obtained. An analytical method is developed for solution of the governing equations for various edge conditions. The governing equations are solved for free and simply supported edge conditions. The interlaminar stresses in the panel are investigated by means of Hooke's law and also by means of integration of the equilibrium equations of elasticity. Dependency of the result upon the number of numerical layers in the layerwise theory (LWT) is studied. The accuracy of the numerical results is validated by comparison with the results of the finite element method and with other available results in the open literature and good agreement is seen between the results. Numerical results are then presented for the distribution of interlaminar normal and shear stresses within the symmetric and un-symmetric cross-ply thick panels with free and simply supported boundaries. The effects of the geometrical parameters such as radius to thickness and width to thickness ratio are investigated on the distribution of the interlaminar stresses in thick panels.