Lateral impact response and parametric studies of axially loaded square concrete filled steel tube columns

This paper presents a numerical study on the response of axially loaded slender square concrete filled steel tube (CFST) columns under low velocity lateral impact loading. A finite element analysis (FEA) model was developed using the explicit dynamic nonlinear finite element code LS -DYNA in which the strain rate effects of both steel and concrete, contact between steel tube and concrete and confinement effect provided by the steel tube for the concrete were considered. The model also benefited from a relatively recent feature of LS-DYNA for applying a pre-loading in the explicit solver. The developed numerical model was verified for its accuracy and adequacy by comparing the results with experimental results available in the literature. The verified model was then employed to conduct a parametric study to investigate the influence of axial load level, impact location, support conditions, and slenderness ratio on the response of the CFST columns. A good agreement between the numerical and experimental results was achieved. The model could reasonably predict the impact load-deflection history and deformed shape of the column at the end of the impact event. The results of the parametric study showed that whilst impact location, axial load level and slenderness ratio can have a significant effect on the peak impact force, residual lateral deflection and maximum lateral deflection, the influence of support fixity is minimal. With an increase of axial load to up to a certain level, the peak force increases; however, a further increase in the axial load causes a decrease in the peak force. Both residual lateral deflection and maximum lateral deflection increase as axial load level increases. Shifting the impact location towards the supports increases the peak force and reduces both residual and maximum lateral deflections. A rise in slenderness ratio decreases the peak force and increases the residual and maximum lateral deflections.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

R2P ten years after the World Summit: Explaining ongoing contestation over Pillar III

This article examines how and why contrasting interpretations of the international community’s role in preventing and responding to mass atrocity crimes continue to exist a decade after the Responsibility to Protect (R2P) was unanimously endorsed at the 2005 World Summit. Building on recent critical constructivist insights into the fluid, dynamic nature of norms, it advances two main arguments. The first is that continuing contestation over R2P’s third pillar is a product of a combination of internal and external sources of norm dynamism. R2P’s inherently complex normative structure, coupled with several external factors, including the broader normative environment, norm implementation experiences and a shift in global power towards the BRICS, have contributed to a period of renewed contestation and triggered attempts to re-formulate R2P thorough Brazil’s ‘Responsibility while Protecting’ (RwP) proposal and China’s semi-official ‘Responsible Protection’ concept. The second central argument is that such contestation is affecting R2P’s distinct normative prescriptions in different ways. While resistance to the implementation of coercive pillar III measures is currently impeding the normative progress of that component of the norm, this contestation has not prevented consensual pillar II assistance from becoming more deeply embedded in international practice and discourse.

Particle field imaging using digital in-line holography

Digital holography is the direct recording of holograms using a CCD camera and is an alternative to the use of a film or a plate. In this communication in-line digital holographic microscopy has been explored for its application in particle imaging in 3D. Holograms of particles of about 10 mu m size have been digitally reconstructed. Digital focusing was done to image the particles in different planes along the depth of focus. Digital holographic particle imaging results were compared with conventional optical microscope imaging. A methodology for dynamic analysis of microparticles in 3D using in-line digital holography has been proposed.

The hemodynamic effects of in-tandem carotid artery stenosis: Implications for carotid endarterectomy

Objectives: It remains controversial whether patients with severe disease of the internal carotid artery and a coexisting stenotic lesion downstream would benefit from a carotid endarterectomy (CEA) of the proximal lesion. The aim of this study was to simulate the hemodynamic and wall shear effects of in-tandem internal carotid artery stenosis using a computational fluid dynamic (CFD) idealized model to give insight into the possible consequences of CEA on these lesions. Methods: A CFD model of steady viscous flow in a rigid tube with two asymmetric stenoses was introduced to simulate blood flow in arteries with multiple constrictions. The effect of varying the distance between the two stenoses, and the severity of the upstream stenosis on the pressure and wall shear stress (WSS) distributions on the second plaque, was investigated. The influence of the relative positions of the two stenoses was also assessed. Results: The distance between the plaques was found to have minimal influence on the overall hemodynamic effect except for the presence of a zone of low WSS (range -20 to 30 dyne/cm2) adjacent to both lesions when the two stenoses were sufficiently close (<4 times the arterial diameter). The upstream stenosis was protective if it was larger than the downstream stenosis. The relative positions of the stenoses were found to influence the WSS but not the pressure distribution. Conclusions: The geometry and positions of the lesions need to be considered when considering the hemodynamic effects of an in-tandem stenosis. Low WSS is thought to cause endothelial dysfunction and initiate atheroma formation. The fact that there was a flow recirculation zone with low WSS in between the two stenoses may demonstrate how two closely positioned plaques may merge into one larger lesion. Decision making for CEA may need to take into account the hemodynamic situation when an in-tandem stenosis is found. CFD may aid in the risk stratification of patients with this problem.

Dynamics of a class of social interaction systems

A mathematical model of social interaction in the form of two coupler! first-order non-linear differential equations, forms the topic of this study. This non-conservative model io representative of such varied social interaction problems as coexisting sub-populations of two different species, arms race between two rival countries and the like. Differential transformation techniques developed elsewhere in the literature are seen to be effective tools of dynamic analysis of this non-linear non-conservative mode! of social interaction process.

Corneal nerves in refractive surgery and dry eye

This study aimed to investigate the morphology and function of corneal sensory nerves in 1) patients after corneal refractive surgery and 2) patients with dry eye due to Sjögren's syndrome. A third aim was to explore the possible correlation between cytokines detected in tears and development of post-PRK subepithelial haze. The main methods used were tear fluid ELISA analysis, corneal in vivo confocal microscopy, and noncontact esthesiometry. The results revealed that after PRK a positive correlation exists between the regeneration of subbasal nerves and the thickness of regenerated epithelium. Pre- or postoperative levels of the tear fluid cytokines TGF-β1, TNF-α, or PDGF-BB did not correlate with the development of corneal haze objectively estimated by in vivo confocal microscopy 3 months after PRK. After high myopic LASIK, a discrepancy between subjective dry eye symptoms and objective signs of dry eye was observed. The majority of patients reported ongoing dry eye symptoms even 5 years after LASIK, although no objective clinical signs of dry eye were apparent. In addition, no difference in corneal sensitivity was observed between these patients and controls. Primary Sjögren's syndrome patients presented with corneal hypersensitivity, although their corneal subbasal nerve density was normal. However, alterations in corneal nerve morphology (nerve sprouting and thickened stromal nerves) and an increased number of antigen-presenting cells among subbasal nerves were observed, implicating the presence of an ongoing inflammation. Based on these results, the relationship between nerve regeneration and epithelial thickness 3 months after PRK appears to reflect the trophic effect of corneal nerves on epithelium. In addition, measurement of tear fluid cytokines may not be suitable for screening patients for risk of scar (haze) formation after PRK. Presumably, at least part of the symptoms of "LASIK-associated dry eye" are derived from aberrantly regenerated and abnormally functioning corneal nerves. Thus, they may represent a form of corneal neuropathy or "phantom pain" rather than conventional dry eye. Corneal nerve alterations and inflammatory findings in Sjögren's syndrome offer an explanation for the corneal hypersensitivity or even chronic pain or hyperalgesia often observed in these patients. In severe cases of disabling chronic pain in patients with dry eye or after LASIK, when conventional therapeutic possibilities fail to offer relief, consultation of a physician specialized in pain treatment is recommended.

Ultrasonic wave characteristics of nanorods via nonlocal strain gradient models

This paper studies an ultrasonic wave dispersion characteristics of a nanorod. Nonlocal strain gradient models (both second and fourth order) are introduced to analyze the ultrasonic wave behavior in nanorod. Explicit expressions are derived for wave numbers and the wave speeds of the nanorod. The analysis shows that the fourth order strain gradient model gives approximate results over the second order strain gradient model for dynamic analysis. The second order strain gradient model gives a critical wave number at certain wave frequency, where the wave speeds are zero. A relation among the number of waves along the nanorod, the nonlocal scaling parameter (e(0)a), and the length of the nanorod is obtained from the nonlocal second order strain gradient model. The ultrasonic wave characteristics of the nanorod obtained from the nonlocal strain gradient models are compared with the classical continuum model. The dynamic response behavior of nanorods is explained from both the strain gradient models. The effect of e(0)a on the ultrasonic wave behavior of the nanorods is also observed. (C) 2010 American Institute of Physics.

Nonlinear aeroelasticity of rotating and flapping wings - a review

The interaction between large deflections, rotation effects and unsteady aerodynamics makes the dynamic analysis of rotating and flapping wing a nonlinear aeroelastic problem. This problem is governed by nonlinear periodic partial differential equations whose solution is needed to calculate the response and loads acting on vehicles using rotary or flapping wings for lift generation. We look at three important problems in this paper. The first problem shows the effect of nonlinear phenomenon coming from piezoelectric actuators used for helicopter vibration control. The second problem looks at the propagation on material uncertainty on the nonlinear response, vibration and aeroelastic stability of a composite helicopter rotor. The third problem considers the use of piezoelectric actuators for generating large motions in a dragonfly inspired flapping wing. These problems provide interesting insights into nonlinear aeroelasticity and show the likelihood of surprising phenomenon which needs to be considered during the design of rotary and flapping wing vehicle

A study of symbolic processing and computational aspects in helicopter dynamics

Even research models of helicopter dynamics often lead to a large number of equations of motion with periodic coefficients; and Floquet theory is a widely used mathematical tool for dynamic analysis. Presently, three approaches are used in generating the equations of motion. These are (1) general-purpose symbolic processors such as REDUCE and MACSYMA, (2) a special-purpose symbolic processor, DEHIM (Dynamic Equations for Helicopter Interpretive Models), and (3) completely numerical approaches. In this paper, comparative aspects of the first two purely algebraic approaches are studied by applying REDUCE and DEHIM to the same set of problems. These problems range from a linear model with one degree of freedom to a mildly non-linear multi-bladed rotor model with several degrees of freedom. Further, computational issues in applying Floquet theory are also studied, which refer to (1) the equilibrium solution for periodic forced response together with the transition matrix for perturbations about that response and (2) a small number of eigenvalues and eigenvectors of the unsymmetric transition matrix. The study showed the following: (1) compared to REDUCE, DEHIM is far more portable and economical, but it is also less user-friendly, particularly during learning phases; (2) the problems of finding the periodic response and eigenvalues are well conditioned.

Schlieren visualization of shock wave phenomena over a missile-shaped body at hypersonic Mach numbers

The flow over a missile-shaped configuration is investigated by means of Schlieren visualization in short-duration facility producing free stream Mach numbers of 5.75 and 8. This visualization technique is demonstrated with a 41 degrees full apex angle blunt cone missile-shaped body mounted with and without cavity. Experiments are carried out with air as the test gas to visualize the flow field. The experimental results show a strong intensity variation in the deflection of light in a flow field, due to the flow compressibility. Shock stand-off distance measured with the Schlieren method is in good agreement with theory and computational fluid dynamic study for both the configurations. Magnitude of the shock oscillation for a cavity model may be greater than the case of a model without cavity. The picture of visualization shows that there is an outgoing and incoming flow closer to the cavity. Cavity flow oscillation was found to subside to steady flow with a decrease in the free stream Mach number.

Hidden order in serrated flow of metallic glasses

We report results of statistical and dynamic analysis of the serrated stress-time curves obtained from compressive constant strain-rate tests on two metallic glass samples with different ductility levels in an effort to extract hidden information in the seemingly irregular serrations. Two distinct types of dynamics are detected in these two alloy samples. The stress-strain curve corresponding to the less ductile Zr65Cu15Ni10Al10 alloy is shown to exhibit a finite correlation dimension and a positive Lyapunov exponent, suggesting that the underlying dynamics is chaotic. In contrast, for the more ductile Cu47.5Zr47.5Al5 alloy, the distributions of stress drop magnitudes and their time durations obey a power-law scaling reminiscent of a self-organized critical state. The exponents also satisfy the scaling relation compatible with self-organized criticality. Possible physical mechanisms contributing to the two distinct dynamic regimes are discussed by drawing on the analogy with the serrated yielding of crystalline samples. The analysis, together with some physical reasoning, suggests that plasticity in the less ductile sample can be attributed to stick-slip of a single shear band, while that of the more ductile sample could be attributed to the simultaneous nucleation of a large number of shear bands and their mutual interactions. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Modeling and optimization of a natural convection-based agricultural dryer

This study presents development of a computational fluid dynamic (CFD) model to predict unsteady, two-dimensional temperature, moisture and velocity distributions inside a novel, biomass-fired, natural convection-type agricultural dryer. Results show that in initial stages of drying, when material surface is wet and moisture is easily available, moisture removal rate from surface depends upon the condition of drying air. Subsequently, material surface becomes dry and moisture removal rate is driven by diffusion of moisture from inside to the material surface. An optimum 9-tray configuration is found to be more efficient than for the same mass of material and volume of dryer. A new configuration of dryer, mainly to explore its potential to increasing uniformity in drying across all trays, is also analyzed. This configuration involves diverting a portion of hot air before it enters over the first tray and is supplied directly at an intermediate location in the dryer. Uniformity in drying across trays has increased for the kind of material simulated.

An efficient hybrid approach for design of automotive wheel bearings

Wheel bearings play a crucial role in the mobility of a vehicle by minimizing motive power loss and providing stability in cornering maneuvers. Detailed engineering analysis of a wheel bearing subsystem under dynamic conditions poses enormous challenges due to the nonlinearity of the problem caused by multiple factional contacts between rotating and stationary parts and difficulties in prediction of dynamic loads that wheels are subject to. Commonly used design methodologies are based on equivalent static analysis of ball or roller bearings in which the latter elements may even be represented with springs. In the present study, an advanced hybrid approach is suggested for realistic dynamic analysis of wheel bearings by combining lumped parameter and finite element modeling techniques. A validated lumped parameter representation serves as an efficient tool for the prediction of radial wheel load due to ground reaction which is then used in detailed finite element analysis that automatically accounts for contact forces in an explicit formulation.