Force chains, friction, and flow: Behavior of granular media across length scales

We study the behavior of granular materials at three length scales. At the smallest length scale, the grain-scale, we study inter-particle forces and "force chains". Inter-particle forces are the natural building blocks of constitutive laws for granular materials. Force chains are a key signature of the heterogeneity of granular systems. Despite their fundamental importance for calibrating grain-scale numerical models and elucidating constitutive laws, inter-particle forces have not been fully quantified in natural granular materials. We present a numerical force inference technique for determining inter-particle forces from experimental data and apply the technique to two-dimensional and three-dimensional systems under quasi-static and dynamic load. These experiments validate the technique and provide insight into the quasi-static and dynamic behavior of granular materials.

At a larger length scale, the mesoscale, we study the emergent frictional behavior of a collection of grains. Properties of granular materials at this intermediate scale are crucial inputs for macro-scale continuum models. We derive friction laws for granular materials at the mesoscale by applying averaging techniques to grain-scale quantities. These laws portray the nature of steady-state frictional strength as a competition between steady-state dilation and grain-scale dissipation rates. The laws also directly link the rate of dilation to the non-steady-state frictional strength.

At the macro-scale, we investigate continuum modeling techniques capable of simulating the distinct solid-like, liquid-like, and gas-like behaviors exhibited by granular materials in a single computational domain. We propose a Smoothed Particle Hydrodynamics (SPH) approach for granular materials with a viscoplastic constitutive law. The constitutive law uses a rate-dependent and dilation-dependent friction law. We provide a theoretical basis for a dilation-dependent friction law using similar analysis to that performed at the mesoscale. We provide several qualitative and quantitative validations of the technique and discuss ongoing work aiming to couple the granular flow with gas and fluid flows.

Regional variations in upper mantle structure beneath North America

Several types of seismological data, including surface wave group and phase velocities, travel times from large explosions, and teleseismic travel time anomalies, have indicated that there are significant regional variations in the upper few hundred kilometers of the mantle beneath continental areas. Body wave travel times and amplitudes from large chemical and nuclear explosions are used in this study to delineate the details of these variations beneath North America.

As a preliminary step in this study, theoretical P wave travel times, apparent velocities, and amplitudes have been calculated for a number of proposed upper mantle models, those of Gutenberg, Jeffreys, Lehman, and Lukk and Nersesov. These quantities have been calculated for both P and S waves for model CIT11GB, which is derived from surface wave dispersion data. First arrival times for all the models except that of Lukk and Nersesov are in close agreement, but the travel time curves for later arrivals are both qualitatively and quantitatively very different. For model CIT11GB, there are two large, overlapping regions of triplication of the travel time curve, produced by regions of rapid velocity increase near depths of 400 and 600 km. Throughout the distance range from 10 to 40 degrees, the later arrivals produced by these discontinuities have larger amplitudes than the first arrivals. The amplitudes of body waves, in fact, are extremely sensitive to small variations in the velocity structure, and provide a powerful tool for studying structural details.

Most of eastern North America, including the Canadian Shield has a Pn velocity of about 8.1 km/sec, with a nearly abrupt increase in compressional velocity by ~ 0.3 km/sec near at a depth varying regionally between 60 and 90 km. Variations in the structure of this part of the mantle are significant even within the Canadian Shield. The low-velocity zone is a minor feature in eastern North America and is subject to pronounced regional variations. It is 30 to 50 km thick, and occurs somewhere in the depth range from 80 to 160 km. The velocity decrease is less than 0.2 km/sec.

Consideration of the absolute amplitudes indicates that the attenuation due to anelasticity is negligible for 2 hz waves in the upper 200 km along the southeastern and southwestern margins of the Canadian Shield. For compressional waves the average Q for this region is > 3000. The amplitudes also indicate that the velocity gradient is at least 2 x 10^-3 both above and below the low-velocity zone, implying that the temperature gradient is < 4.8°C/km if the regions are chemically homogeneous.

In western North America, the low-velocity zone is a pronounced feature, extending to the base of the crust and having minimum velocities of 7.7 to 7.8 km/sec. Beneath the Colorado Plateau and Southern Rocky Mountains provinces, there is a rapid velocity increase of about 0.3 km/sec, similar to that observed in eastern North America, but near a depth of 100 km.

Complicated travel time curves observed on profiles with stations in both eastern and western North America can be explained in detail by a model taking into account the lateral variations in the structure of the low-velocity zone. These variations involve primarily the velocity within the zone and the depth to the top of the zone; the depth to the bottom is, for both regions, between 140 and 160 km.

The depth to the transition zone near 400 km also varies regionally, by about 30-40 km. These differences imply variations of 250 °C in the temperature or 6 % in the iron content of the mantle, if the phase transformation of olivine to the spinel structure is assumed responsible. The structural variations at this depth are not correlated with those at shallower depths, and follow no obvious simple pattern.

The computer programs used in this study are described in the Appendices. The program TTINV (Appendix IV) fits spherically symmetric earth models to observed travel time data. The method, described in Appendix III, resembles conventional least-square fitting, using partial derivatives of the travel time with respect to the model parameters to perturb an initial model. The usual ill-conditioned nature of least-squares techniques is avoided by a technique which minimizes both the travel time residuals and the model perturbations.

Spherically symmetric earth models, however, have been found inadequate to explain most of the observed travel times in this study. TVT4, a computer program that performs ray theory calculations for a laterally inhomogeneous earth model, is described in Appendix II. Appendix I gives a derivation of seismic ray theory for an arbitrarily inhomogeneous earth model.

Imaging the earth with ambient noise and earthquakes

In this thesis, I develop the velocity and structure models for the Los Angeles Basin and Southern Peru. The ultimate goal is to better understand the geological processes involved in the basin and subduction zone dynamics. The results are obtained from seismic interferometry using ambient noise and receiver functions using earthquake- generated waves. Some unusual signals specific to the local structures are also studied. The main findings are summarized as follows:

(1) Los Angeles Basin

The shear wave velocities range from 0.5 to 3.0 km/s in the sediments, with lateral gradients at the Newport-Inglewood, Compton-Los Alamitos, and Whittier Faults. The basin is a maximum of 8 km deep along the profile, and the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the center decreasing to 1.3 at the edges, and is in approximate isostatic equilibrium. This "high-density" (~1 km spacing) "short-duration" (~1.5 month) experiment may serve as a prototype experiment that will allow basins to be covered by this type of low-cost survey.

(2) Peruvian subduction zone

Two prominent mid-crust structures are revealed in the 70 km thick crust under the Central Andes: a low-velocity zone interpreted as partially molten rocks beneath the Western Cordillera – Altiplano Plateau, and the underthrusting Brazilian Shield beneath the Eastern Cordillera. The low-velocity zone is oblique to the present trench, and possibly indicates the location of the volcanic arcs formed during the steepening of the Oligocene flat slab beneath the Altiplano Plateau.

The Nazca slab changes from normal dipping (~25 degrees) subduction in the southeast to flat subduction in the northwest of the study area. In the flat subduction regime, the slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes a normal dipping geometry. The flat part closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the western half of the flat slab, which indicates the lack of melting and thus explains the cessation of the volcanism above. The velocity turns to normal values before the slab steepens again, indicating possible resumption of dehydration and ecologitization.

(3) Some unusual signals

Strong higher-mode Rayleigh waves due to the basin structure are observed in the periods less than 5 s. The particle motions provide a good test for distinguishing between the fundamental and higher mode. The precursor and coda waves relative to the interstation Rayleigh waves are observed, and modeled with a strong scatterer located in the active volcanic area in Southern Peru. In contrast with the usual receiver function analysis, multiples are extensively involved in this thesis. In the LA Basin, a good image is only from PpPs multiples, while in Peru, PpPp multiples contribute significantly to the final results.

Interpretación litoestratigráfica y paleogeográfica del Jurásico marino del este Cantabria-oeste de Vizcaya en base a diagrafías de sondeo y su correlación con series tipo

Este trabajo se centra en el estudio del Jurásico marino (Lías y Dogger) de un área de la parte central de La Cuenca Vasco-Cantábrica situada en el este de Cantabria y oeste de Vizcaya donde sus afloramientos son muy escasos pero se dispone de datos de sondeos de exploración petrolífera (Fig 1A y B), poniendo especial interés en la localización y distribución lateral de los niveles de rocas madre de hidrocarburos (black shales) y de potenciales almacenes (unidades carbonatadas fracturadas o dolomitizadas). La sucesión del Jurásico de la Cuenca Vasco-Cantábrica (CVC) está formada por dos unidades diferenciadas por edad y ambiente sedimentario. Por una parte tenemos los materiales del “Jurásico marino” (Robles et al., 1989) que representan la mayor parte de la sucesión (Lías y Dogger) y por otro lado, tenemos los materiales del “Jurásico continental” pertenecientes exclusivamente al Tithoniense Superior y que se engloban en las facies Purbeck que abarcan hasta el Berriasiense (Rat 1962).

Estudo do conforto humano em pisos mistos (aço-concreto) submetidos a ações humanas rítmicas.

Este trabalho de pesquisa objetiva o estudo do comportamento dinâmico de pisos mistos (aço-concreto), em edificações de andares múltiplos, sob o ponto de vista de conforto humano, quando essas estruturas encontram-se submetidas às atividades rítmicas provenientes dos seres humanos. A definição das ações dinâmicas atuantes sobre os modelos estruturais foi feita com base em resultados experimentais, com os indivíduos praticando atividades rítmicas e não rítmicas associadas à ginástica aeróbica e saltos à vontade. Os modelos estruturais investigados baseiam-se em edificações mistas de andares múltiplos. O sistema estrutural é do tipo misto (aço-concreto), composto por vigas de aço em seção do tipo I e laje de concreto armado. A análise fundamenta-se na modelagem computacional dos sistemas estruturais, através do Método dos Elementos Finitos (MEF). São empregadas técnicas usuais de discretização, por meio do emprego do programa ANSYS. Uma análise paramétrica foi desenvolvida sobre três modelos estruturais, com dois, três e quatro pavimentos. Os valores das acelerações máximas encontradas na análise são confrontados e comparados com os limites propostos por recomendações internacionais. Os resultados obtidos mostram que os limites recomendados em diversas normas de projeto foram ultrapassados. Esses resultados demonstram que atividades rítmicas oriundas dos seres humanos podem gerar acelerações de pico elevadas, violando critérios de projeto, no que concerne ao conforto humano. Foi observado também que estas ações dinâmicas podem comprometer o conforto humano em pisos adjacentes, próximos do local onde a carga dinâmica está sendo efetivamente aplicada.

Marine Fisheries Stock Assessment Improvement Plan: report of the National Marine Fisheries Service National Task Force for Improving Fish Stock Assessments

This report argues for greatly increased resources in terms of data collection facilities and staff to collect, process, and analyze the data, and to communicate the results, in order for NMFS to fulfill its mandate to conserve and manage marine resources. In fact, the authors of this report had great difficulty defining the "ideal" situation to which fisheries stock assessments and management should aspire. One of the primary objectives of fisheries management is to develop sustainable harvest policies that minimize the risks of overfishing both target species and associated species. This can be achieved in a wide spectrum of ways, ranging between the following two extremes. The first is to implement only simple management measures with correspondingly simple assessment demands, which will usually mean setting fishing mortality targets at relatively low levels in order to reduce the risk of unknowingly overfishing or driving ecosystems towards undesirable system states. The second is to expand existing data collection and analysis programs to provide an adequate knowledge base that can support higher fishing mortality targets while still ensuring low risk to target and associated species and ecosystems. However, defining "adequate" is difficult, especially when scientists have not even identified all marine species, and information on catches, abundances, and life histories of many target species, and most associated species, is sparse. Increasing calls from the public, stakeholders, and the scientific community to implement ecosystem-based stock assessment and management make it even more difficult to define "adequate," especially when "ecosystem-based management" is itself not well-defined. In attempting to describe the data collection and assessment needs for the latter, the authors took a pragmatic approach, rather than trying to estimate the resources required to develop a knowledge base about the fine-scale detailed distributions, abundances, and associations of all marine species. Thus, the specified resource requirements will not meet the expectations of some stakeholders. In addition, the Stock Assessment Improvement Plan is designed to be complementary to other related plans, and therefore does not duplicate the resource requirements detailed in those plans, except as otherwise noted.

Seismic performance of geomembrane placed on a slope in a MSW landfill cell -A centrifuge study

Geomembranes are one of the most commonly used geosynthetics in landfill liner systems. They retain the leachate produced by the waste and prevent leakage. Geomembranes may experience harsh environmental conditions such as extreme temperatures or earthquake loading. Earthquake loading can be an extreme loading case for landfills located in seismic regions. This study, based on dynamic centrifuge testing, investigates the effects of simulated earthquake loading on the tension experienced bythe geomembrane on a landfill slope. The landfill modeled in the dynamic centrifuge test was a municipal solid waste (MSW) landfill cell with a single geomembrane-clay liner system (45° side slope and 10 m slope length). The paper shows that moderate earthquake loading (base acceleration between 0.1g to 0.2g) can result in transient increases of around 20% in geomembrane tension, with permanent tension increases of around 5%.

The influence of peel angle on the mechanics of peeling flexible adherends with arbitrary load-extension characteristics

The peel test is commonly used to determine the strength of adhesive joints. In its simplest form, a thin flexible strip which has been bonded to a rigid surface is peeled from the substrate at a constant rate and the peeling force which is applied to the debonding surfaces by the tension in the tape is measured. Peeling can be carried out with the peel angle, i.e. the angle made by the peel force with the substrate surface, from any value above about 10° although peeling tests at 90 and 180° are most common. If the tape is sufficiently thin for its bending resistance to be negligibly small then as well as the debonding or decohesion energy associated with the adhesive in and around the point of separation, the relation between the peeling force and the peeling angle is influenced both by the mechanical properties of the tape and any pre-strain locked into the tape during its application to the substrate. The analytic solution for a tape material which can be idealised as elastic perfectly-plastic is well established. Here, we present a more general form of analysis, applicable in principle to any constitutive relation between tape load and tape extension. Non-linearity between load and extension is of increasing significance as the peel angle is decreased: the model presented is consistent with existing equations describing the failure of a lap joint between non-linear materials. The analysis also allows for energy losses within the adhesive layer which themselves may be influenced by both peel rate and peel angle. We have experimentally examined the application of this new analysis to several specific peeling cases including tapes of cellophane, poly-vinyl chloride and PTFE. © 2005 Elsevier Ltd. All rights reserved.

A normal force distance regulation scheme for near-field optical microscopy

A near-field optical microscope (NFOM) has been developed that combines the features of a near-field optical microscope and an atomic force microscope. Improved control over tip-sample separation has led to improved optical imaging and independent surface topography information. The tip oscillation is normal to the sample plane thereby reducing lateral forces - important for nonperturbative imaging of soft samples. Both topographic images and reflection near-field optical images are presented which demonstrate the capability of the system. © 1996 American Institute of Physics.

Modelling effects of axial load on behaviour of pile groups in laterally-spreading soil

Previous research into the behaviour of piled foundations in laterally-spreading soil deposits has concentrated on pile groups that carry small or negligible axial loads. This paper presents dynamic centrifuge test results for 2 x 2 pile groups with bending and geometric properties similar to real 0.5 m diameter tubular steel and solid circular reinforced-concrete field piles. Axial loads applied represented upper-bounds on typical working loads. The simultaneous scaling of the relevant properties controlling both lateral and axial behaviour allows comparisons to be drawn regarding the particular mechanisms of failure that would dominate for each type of pile. Flexible reinforced-concrete piles which tend to carry lower loads were found to be dominated by lateral effects, while steel piles, which are much stiffer and usually carry greater loads are dominated by settlement considerations. © 2006 Taylor & Francis Group.

Seismic behaviour of anchored quay walls with dry backfill

The seismic behaviour of anchored sheet pile walls is a complex soil-structure interaction problem. Damaged sheet pile walls are very expensive to repair and their seismic behaviour needs to be investigated in order to understand their possible mechanisms of failure. The research described in this paper involves both centrifuge testing and Finite Element (FE) analyses aimed at investigating the seismic behaviour of an anchored sheet pile wall in dry sand. The model wall is tied to the backfill with two tie rods connected to an anchor beam. The accelerations of the sheet pile wall, the anchor beam and the soil around the wall were measured using miniature piezoelectric accelerometers. The displacement at the tip of the wall was also measured. Stain gauges at five different locations on the wall were used to measure the bending moments induced in the the wall. The anchor forces in the tie rods were also measured using load cells. The results from the centrifuge tests were compared with 2-D, plane strain FE analyses conducted using DIANA-SWANDYNE II and the observed seismic behaviour was explained in the light of these findings. © 2011 Taylor & Francis.

Preservation of Sardinella longiceps in iced and chilled seawater. Part 1 - Changes during storage with particular reference to bacterial load and nitrogenous compounds

Oil sardines in prime condition were subjected to onboard chilling. Two lots were chilled in CSW (samples C and CI), a third lot was chilled in crushed ice (sample I) and a fourth lot left not iced on deck (Sample AI). Upon landing sample AI was iced and sample CI was removed from the CSW and iced. All the four samples were kept in a chilled room for storage studies. The fish chilled and stored in CSW recorded the least, and the fish subjected to delayed icing, the highest values for all the indices of spoilage namely, free amino nitrogen, trimethylamine (TMA) and total volatile base nitrogen (TVBN). The total psychrophilic bacterial number also showed a similar trend. The organoleptic assessment of the cooked samples revealed C I, CI, AI to be the order of preference throughout the storage. This assessment was found to hold good for the rest of the parameters as well. The CSW held fishes were found to be distinctly superior to the iced ones for the first five days of storage. Such a marked prevalence in quality for five days would suffice for the fish to fetch a premium in the market over other landings of the same fish whether chilled or not chilled. Chilling on board in CSW and icing the same after landings, did not show encouraging results.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.

3D FE Analyses of Buried Pipeline with Elbows Subjected to Lateral Loading

This study investigates the interaction between soil and pipeline in sand subjected to lateral ground displacements with emphasis on the peak force exerted to a bended elbow-pipe. A series of three-dimensional (3D) finite-element (FE) analyses were performed in both opening and closing modes of the elbow section for different initial pipe bending angles. To model the mechanical behavior of sands, two soil models were adopted: Mohr-Coulomb and Nor-Sand soil model. Investigations also included the effects of pipe embedment depth and soil density. Results show that the opening mode exhibits higher ultimate forces and greater localized deformations than the closing mode. Nondimensional charts that account for pipeline location, bending angle, and soil density are developed. Soil-spring pipeline analyses of an elbow-pipe were performed using modified F-δ soil-spring models based on the 3D FE results and were compared to the findings of conventional spring model analyses using the standard two-dimensional soil-spring model. Results show that the pipe strain does not change in the closing mode case. However, in the opening mode case, the pipe strain computed by the modified analysis is larger than that by the conventional analysis and the difference is more pronounced when the pipe stiffness is stiffer. © 2011 American Society of Civil Engineers.

Modelling of liquefaction-induced instability in pile

Centrifuge testing has been undertaken to investigate instability failure of pile groups during seismic liquefaction, with specific reference to the 'top-down' propagation of liquefaction during the earthquake and to account for initial imperfections in pile geometry. The results of these tests were used to validate numerical models within the finite element program ABAQUS, based on the popular p-y analysis method. Pseudostatic classical and post-buckling analyses were conducted to examine the collapse behaviour of the pile groups and were found to give reasonable predictions of collapse load and conservative predictions of the associated deflection conditions. This numerical model was compared to currently published methods which were found to over-predict collapse loads. The resulting insights into the collapse of axially loaded pile groups revealed that the failure load is strongly dependent on both the depth of liquefaction propagation and initial imperfections, which reduce the collapse load.