• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Transactions of Materials Processing
• /
• v.17 no.3
• /
• pp.161-169
• /
• 2008

In the present work, the two-back stress model is proposed and continuum damage mechanics (CDM) is incorporated into the plastic constitutive relation in order to describe the plastic deformation localization and the damage evolution in a deforming continuum body. Coupling between damage mechanics and isothermal rate independent plasticity is performed using the kinematic hardening rule, which in turn is formulated by combining the nonlinear Armstrong-Frederick rule and the Phillips rule. The numerical analyses are carried out within h deformation theory. It is noted that the damage evolution within a work piece accelerates the plastic deformation localization such that the material with lower hardening exponent results in a rapid shear band formation. Moreover, the results from the numerical analysis reflected closely with the micro-structures around the fractured regime. The effects of the various hardening parameters on deformation localization are also investigated. As the nonlinear strain rate description in the back stress evolution becomes dominant, the strain localization becomes intensified as well as the damage evolution.

• Structural Engineering and Mechanics
• /
• v.5 no.1
• /
• pp.1-19
• /
• 1997

A numerical procedure for dynamic analysis of structures including lateral-torsional coupling, axial force effect and soil-structure interaction is presented in this study. A simple soil-structure system model has been designed for microcomputer applications capable of reflecting both kinematic and inertial soil-foundation interaction as well as the effect of this interaction on the superstructure response. A parametric study focusing on inertial soil-structure interaction is carried out through a simplified nine-degree of freedom building model with different foundation conditions. The inertial soil-structure interaction and axial force effects on a 20-storey building excited by an Australian earthquake is analysed through its top floor displacement time history and envelope values of structural maximum displacement and shear force.

• Coupled systems mechanics
• /
• v.10 no.6
• /
• pp.521-544
• /
• 2021

In this work, we present the development of a 3D lattice-type model at microscale based upon the Voronoi-cell representation of material microstructure. This model can capture the coupling between mechanic and electric fields with non-linear constitutive behavior for both. More precisely, for electric part we consider the ferroelectric constitutive behavior with the possibility of domain switching polarization, which can be handled in the same fashion as deformation theory of plasticity. For mechanics part, we introduce the constitutive model of plasticity with the Armstrong-Frederick kinematic hardening. This model is used to simulate a complete coupling of the chosen electric and mechanics behavior with a multiscale approach implemented within the same computational architecture.

• International Journal of Aeronautical and Space Sciences
• /
• v.9 no.1
• /
• pp.66-75
• /
• 2008

In this work, we explore the feasibility of roll-pitch-yaw integrated autopilots for high angle-of-attack missiles. An investigation of the aerodynamic characteristics of a surface-to-air missile is presented, which reveals the strong effects of cross coupling between the longitudinal and lateral dynamics. Robust control techniques based on $H_{\infty}$ synthesis are employed to design roll-pitch-yaw integrated autopilots. The performance of the proposed roll-pitch-yaw integrated controller is tested in high-fidelity nonlinear five-degree-of-freedom simulations accounting for kinematic cross-coupling effects between the lateral and longitudinal channels. Against nonlinearity and cross-coupling effects of the missile dynamics, the integrated controller demonstrates superior performance when compared with the controller designed in a decoupled manner.

• Coupled systems mechanics
• /
• v.8 no.1
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• Korean Journal of Applied Biomechanics
• /
• v.31 no.4
• /
• pp.259-269
• /
• 2021

Objective: The purpose of this systematic review and meta-analysis was to investigate age-related bimanual coordination functions in older adults. Method: Thirteen studies that compared bimanual coordination functions in older adults with those in healthy young adults qualified for this meta-analysis. We additionally categorized 21 total comparisons from the 13 qualified studies into two types of task-related moderator variables: (1) kinematic versus kinetic movements and (2) symmetry versus asymmetry movements. Results: Random effects model meta-analysis found that older adults revealed significant bimanual coordination impairments as compared with young adults (Hedges's g = -0.771; p < .0001; I² = 74.437%). We additionally confirmed specific bimanual coordination deficits using two moderator variables: 1) kinematic (Hedges's g = -0.884; p < .0001; I² = 0.000%) and kinetic (Hedges's g = -0.666; p = .023; I² = 86.170%). 2) symmetry (Hedges's g = -0.712; p = .001; I² = 74.291%) and asymmetry (Hedges's g = -0.817; p < .0001; I² = 76.322%). The moderator variable analysis indicated older adults indicated bimanual coordination deficits in the upper extremities than healthy young adults while performing kinematic bimanual coordination tasks and asymmetry coordination tasks. Conclusion: These findings suggest that developing motor rehabilitation programs based on asymmetric bimanual movement task for enhancing interlimb coordination functions of older adults may be crucial for increasing their independence in everyday activities. Given that elderly revealed the deficits in lower extremities coordination when older adults perform gait, posture, and balance, future studies should estimate lower limb coordination functions in elderly people.

• Structural Engineering and Mechanics
• /
• v.63 no.1
• /
• pp.37-46
• /
• 2017

This paper is concerned with the numerical study on the resonance response of a rigid spar-type floating platform in coupled heave and pitch motion. Spar-type floating platforms, widely used for supporting the offshore structures, offer an economic advantage but those exhibit the dynamically high sensitivity to external excitations due to their shape at the same time. Hence, the investigation of their dynamic responses, particularly at resonance, is prerequisite for the design of spar-type floating platforms which secure the dynamic stability. Spar-type floating platform in 2-D surface wave is assumed to be a rigid body having 2-DOFs, and its coupled dynamic equations are analytically derived using the geometric and kinematic relations. The motion-variance of the metacentric height and the moment of inertia of floating platform are taken into consideration, and the hydrodynamic interaction between the wave and platform motions is reflected into the hydrodynamic force and moment and the frequency-dependent added masses. The coupled nonlinear equations governing the heave and pitch motions are solved by the RK4 method, and the frequency responses are obtained by the digital Fourier transform. Through the numerical experiments to the wave frequency, the resonance responses and the coupling in resonance between heave and pitch motions are investigated in time and frequency domains.

• Wind and Structures
• /
• v.32 no.3
• /
• pp.267-281
• /
• 2021

Bending-twisting coupling induced in big composite wind turbine blades is one of the passive control mechanisms which is exploited to mitigate loads incurred due to deformation of the blades. In the present study, flutter characteristics of bend-twist coupled blades, designed for load alleviation in wind turbine systems, are investigated by time-domain analysis. For this purpose, a baseline full GFRP blade, a bend-twist coupled full GFRP blade, and a hybrid GFRP and CFRP bend-twist coupled blade is designed for load reduction purpose for a 5 MW wind turbine model that is set up in the wind turbine multi-body dynamic code PHATAS. For the study of flutter characteristics of the blades, an over-speed analysis of the wind turbine system is performed without using any blade control and applying slowly increasing wind velocity. A detailed procedure of obtaining the flutter wind and rotational speeds from the time responses of the rotational speed of the rotor, flapwise and torsional deformation of the blade tip, and angle of attack and lift coefficient of the tip section of the blade is explained. Results show that flutter wind and rotational speeds of bend-twist coupled blades are lower than the flutter wind and rotational speeds of the baseline blade mainly due to the kinematic coupling between the bending and torsional deformation in bend-twist coupled blades.

• Korean Journal of Applied Biomechanics
• /
• v.16 no.3
• /
• pp.149-163
• /
• 2006

The purpose of this study was to compare the lower extremity's joint and segment coupling patterns between forward and backward running in subjects who were twelve healthy males. Three-dimensional kinematic data were collected with Qualisys system while subjects ran to forward and backward. The thigh internal/external rotation and tibia internal/external rotation, thigh flexion/extension and tibia flexion/extension, tibia internal/external rotation and foot inversion/eversion, knee internal/external rotation and ankle inversion/eversion, knee flexion/extension and ankle inversion/eversion, knee flexion/extension and ankle flexion/extension, and knee flexion/extension and tibia internal/external rotation coupling patterns were determined using a vector coding technique. The comparison for each coupling between forward and backward running were conducted using a dependent, two-tailed t-test at a significant level of .05 for the mean of each of five stride regions, midstance(1l-30%), toe-off(31-50%), swing acceleration(51-70%), swing deceleration(71-90), and heel-strike(91-10%), respectively. 1. The knee flexion/extension and ankle flexion/extension coupling pattern of both foreward and backward running over the stride was converged on a complete coordination. However, the ankle flexion/extension to knee flexion/extension was relatively greater at heel-strike in backward running compared with forward running. At the swing deceleration, backward running was dominantly led by the ankle flexion/extension, but forward running done by the knee flexion/extension. 2. The knee flexion/extension and ankle inversion/eversion coupling pattern for both running was also converged on a complete coordination. At the mid-stance. the ankle movement in the frontal plane was large during forward running, but the knee movement in the sagital plane was large during backward running and vice versa at the swing deceleration. 3. The knee flexion/extension and tibia internal/external rotation coupling while forward and backward run was also centered on the angle of 45 degrees, which indicate a complete coordination. However, tibia internal/external rotation dominated the knee flexion/extension at heel strike phase in forward running and vice versa in backward running. It was diametrically opposed to the swing deceleration for each running. 4. Both running was governed by the ankle movement in the frontal plane across the stride cycle within the knee internal/external rotation and tibia internal/external rotation. The knee internal/external rotation of backward running was greater than that of forward running at the swing deceleration. 5. The tibia internal/external rotation in coupling between the tibia internal/external rotation and foot inversion/eversion was relatively great compared with the foot inversion/eversion over a stride for both running. At heel strike, the tibia internal/external rotation of backward running was shown greater than that of forward(p<.05). 6. The thigh internal/external rotation took the lead for both running in the thigh internal/external rotation and tibia internal/external rotation coupling. In comparison of phase, the thigh internal/external rotation movement at the swing acceleration phase in backward running worked greater in comparison with forward running(p<.05). However, it was greater at the swing deceleration in forward running(p<.05). 7. With the exception of the swing deceleration phase in forward running, the tibia flexion/extension surpassed the thigh flexion/extension across the stride cycle in both running. Analysis of the specific stride phases revealed the forward running had greater tibia flexion/extension movement at the heel strike than backward running(p<.05). In addition, the thigh flexion/extension and tibia flexion/extension coupling displayed almost coordination at the heel strike phase in backward running. On the other hand the thigh flexion/extension of forward running at the swing deceleration phase was greater than the tibia flexion/extension, but it was opposite from backward running. In summary, coupling which were the knee flexion/extension and ankle flexion/extension, the knee flexion/extension and ankle inversion/eversion, the knee internal/external rotation and ankle inversion/eversion, the tibia internal/external rotation and foot inversion/eversion, the thigh internal/external rotation and tibia internal/external rotation, and the thigh flexion/extension and tibia flexion/extension patterns were most similar across the strike cycle in both running, but it showed that coupling patterns in the specific stride phases were different from average point of view between two running types.