• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4D
• /
• pp.671-678
• /
• 2006

In this study, to describe vehicle-structure dynamic interaction phenomena with 1/4 vehicle model, nonlinear Hertzian contact spring and nonlinear contact damper are adopted. The external loads acting on 1/4 vehicle model are selfweight of vehicle and geometry information of running surface. The constraint equation on contact surface is implemented by the Penalty method with stabilization and the reaction from constraint violation. To describe pitching motion of various vehicles two types of the displacement constraint equations are exerted to connect between car bodies and between bogie frames, i.e., the rigid body connection and the rigid body connection with pin, respectively. For the time integration of dynamic equations of vehicles and structure Newmark time integration scheme is adopted. To reduce the error caused by inadequate time step size, adaptive time-stepping technique is also adopted. Thus, it is expected that more versatile dynamic interaction phenomena can be described by this approach and it can be applied to various railway dynamic problems with low computational cost.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.6 s.40
• /
• pp.73-80
• /
• 2004

A nonlinear time-domain system identification algorithm using measured acceleration data is developed for structural damage assessment. To take account of nonlinear behavior of structural systems, an output error between measured and computed acceleration increments has been defined and a constrained nonlinear optimization problem is solved for optimal structural parameters. The algorithm estimates time-varying properties of stiffness and damping parameters. Nonlinear response of restoring force of a structural system is recovered by using the estimated time-varying structural properties and computed displacement by Newmark-$\beta$ method. In the recovery, no pre-defined model for inelastic behavior has been assumed. In developing the algorithm, noise and incomplete measurement in space and state have been considered. To examine the developed algorithm, numerical simulation and laboratory experimental studies on a three-story shear building have been carried out.

• Computers and Concrete
• /
• v.11 no.1
• /
• pp.51-61
• /
• 2013

This paper describes earthquake response of the Arhavi Highway Tunnel its geometrical properties, 3D finite element model and the linear time history analyses under a huge ground motion considering soil-structure interaction. The Arhavi Highway Tunnel is one of the tallest tunnels constructed in the Black Sea region of Turkey as part of the Coast Road Project. The tunnel has two tubes and each of them is about 1000 m tall. In the study, lineartime history analyses of the tunnel are performed applying north-south, east-west and up accelerations components of 1992 Erzincan, Turkey ground motion. In the time history analyses, Rayleigh damping coefficients are calculated using main natural frequency obtained from modal analysis. Element matrices are computed using the Gauss numerical integration technique. The Newmark method is used in the solution of the equation of motion. Because of needed too much memory for the analyses, the first 10 second of the ground motions, which is the most effective duration, is taken into account in calculations. The results obtained 3D finite element model are presented. In addition, the displacement and stress results are observed to be allowable level of the concrete material during the earthquakes.

• Coupled systems mechanics
• /
• v.2 no.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
• /
• v.75 no.6
• /
• pp.713-722
• /
• 2020

This article aims to illustrate the damped dynamic responses of layered functionally graded (FG) thick 2D beam under dynamic pulse sinusoidal load by using finite element method, for the first time. To investigate the response of thick beam accurately, two-dimensional plane stress problem is assumed to describe the constitutive behavior of thick beam structure. The material is distributed gradually through the thickness of each layer by generalized power law function. The Kelvin-Voigt viscoelastic constitutive model is exploited to include the material internal damping effect. The governing equations are obtained by using Lagrange's equations and solved by using finite element method with twelve -node 2D plane element. The dynamic equation of motion is solved numerically by Newmark implicit time integration procedure. Numerical studies are presented to illustrate stacking sequence and material gradation index on the displacement-time response of cantilever beam structure. It is found that, the number of waves increases by increasing the graduation distribution parameter. The presented mathematical model is useful in analysis and design of nuclear, marine, vehicle and aerospace structures those manufactured from functionally graded materials (FGM).

• Proceedings of the KSR Conference
• /
• 2003.05a
• /
• pp.376-381
• /
• 2003

This study is focused on the dynamic response of curved bridge when the rubber tired AGT vehicles is running with alternative articulations. For the analytic approach, there is necessary for the three dimensional vehicle model with 11 degree of freedom and the three dimensional curved bridge model by means of finite element method. It can be described by conventional Lagrangian formula with respect to the dynamic interactions between vehicles and its met bridge. The formula is implemented by Fortran language on the simulation program designated BADIA II(Bridge-AGT Dynamic Interaction Analysis II). The solutions of the formula are derived by Newmark- ${\beta}$ method. The BADIA II is for the dynamic interactions between vehicle and curved bridge in terms of the roughness of running surface and guide rail. The applicability of the BADIA II is verified in terms of displacement and modal frequency. This study is described that the dynamic interactive behaviors between the rubber tired AGT vehicle and curved bridge in terms of the radius of curvatures of curved bridge, vehicle articulations, vehicle speeds, vehicle weights, flatness of running surface and roughness of guide rail using BADIA II.

• Structural Engineering and Mechanics
• /
• v.73 no.1
• /
• pp.65-73
• /
• 2020

The present work describes a new time marching procedure for structural dynamics analyses. In this novel technique, time integration parameters are automatically evaluated according to the properties of the model. Such parameters are locally defined, allowing the user to input a numerical dissipation property for each element, which defines the amount of numerical dissipation to be introduced. Since the integration parameters are locally defined as a function of the structural element itself, the time marching technique adapts according to the model, providing enhanced accuracy. The new methodology is based on displacement-velocity relations and no computation of accelerations is required. Furthermore, the method is second order accurate, it has guaranteed stability, it is truly self-starting and it allows highly controllable algorithm dissipation in the higher modes. Numerical results are presented and compared to those provided by the Newmark and the Bathe methods, illustrating the good performance of the new time marching procedure.

• Ocean and Polar Research
• /
• v.32 no.4
• /
• pp.361-367
• /
• 2010

The principal objective of this paper was to evaluate the steering characteristics of a tandem tracked vehicle, each side of which features two tandem tracks, when crawling on extremely cohesive soft soil. The tandem tracked vehicle is assumed to be a rigid-body with 6-dof. The dynamic analysis program of the tandem tracked vehicle was developed via Newmark's method and the incremental-iterative method. A terra-mechanics model of extremely cohesive soft soil was implemented according to the relationships of normal pressure to sinkage, of shear resistance to shear displacement, and of dynamic sinkage to shear displacement. In order to simplify the characteristics of contact interaction between track segments and cohesive soft soil, the characteristics of soil are equated to the properties of intact soil. In an effort to evaluate the steering characteristics of a tandem tracked vehicle crawling on extremely cohesive soft soil, a series of dynamic simulations were conducted for a tandem tracked vehicle model with respect to the front and rear steering angle, the steering ratio, and the initial velocity.

• Wind and Structures
• /
• v.10 no.4
• /
• pp.367-382
• /
• 2007

A nonlinear numerical method was developed to assess the stability of suspension bridge catwalks under a wind load. A section model wind tunnel test was used to obtain a catwalk's aerostatic coefficients, from which the displacement-dependent wind loads were subsequently derived. The stability of a suspension bridge catwalk was analyzed on the basis of the geometric nonlinear behavior of the structure. In addition, a full model test was conducted on the catwalk, which spanned 960 m. A comparison of the displacement values between the test and the numerical simulation shows that a numerical method based on a section model test can be used to effectively and accurately evaluate the stability of a catwalk. A case study features the stability of the catwalk of the Runyang Yangtze suspension bridge, the main span of which is 1490 m. Wind can generally attack the structure from any direction. Whenever the wind comes at a yaw angle, there are six wind load components that act on the catwalk. If the yaw angle is equal to zero, the wind is normal to the catwalk (called normal wind) and the six load components are reduced to three components. Three aerostatic coefficients of the catwalk can be obtained through a section model test with traditional test equipment. However, six aerostatic coefficients of the catwalk must be acquired with the aid of special section model test equipment. A nonlinear numerical method was used study the stability of a catwalk under a yaw wind, while taking into account the six components of the displacement-dependent wind load and the geometric nonlinearity of the catwalk. The results show that when wind attacks with a slight yaw angle, the critical velocity that induces static instability of the catwalk may be lower than the critical velocity of normal wind. However, as the yaw angle of the wind becomes larger, the critical velocity increases. In the atmospheric boundary layer, the wind is turbulent and the velocity history is a random time history. The effects of turbulent wind on the stability of a catwalk are also assessed. The wind velocity fields are regarded as stationary Gaussian stochastic processes, which can be simulated by a spectral representation method. A nonlinear finite-element model set forepart and the Newmark integration method was used to calculate the wind-induced buffeting responses. The results confirm that the turbulent character of wind has little influence on the stability of the catwalk.