Search

The main objective of this research is to formulate and couple technologies for modeling discrete and continuous media using real particle morphologies. To that end, two coupled formulations based on virtual modeling technologies of single real particles with another one called real particle packing technique are presented. The first formulation employs Fourier descriptors’ theory to virtually achieve the morphology and construct a repository of real particle geometries. The second formulation is a particle packing method, supported by advancing front techniques combined with dynamic methods. This method presents a stochastic formulation and allows the packing of particle systems following continuous, discrete and empirical statistical distributions. The coupling of both techniques is a very efficient tool to achieve discrete or continuous media geometries to solve engineering problems. Three different examples are developed to illustrate the usefulness of the formulations. The first one is a discrete angle-of-repose problem involving clusters of spheres (real particle morphologies are described with groups of spheres); in the second example the same angle-of-repose problem is resolved with real particles. In the third case, which involves continuous medium mechanics, a small-scale road engineering problem is modeled, specifically, the testing of an asphalt concrete.

Borehole breakout initiation, progression, and stabilization are modeled using a semi-analytical method based on Melentiev’s graphical conformal mapping procedure, and Kolosov–Muskhelishvili complex stress potentials. The only input data required are the elastic moduli of the rock, the Coulomb strength parameters (cohesion and angle of internal friction), and the far-field stresses. The stresses around the borehole wall are computed, the region in which the rock has failed is then “removed”, creating a new borehole shape. This process is iterated until a shape is obtained for which the breakout will progress no further, and a stable state has been reached. This modeling shows that stresses around the flank of the breakout evolve so as to reduce the propensity for shear failure, which helps to explain why the breakout width remains relatively constant throughout the process, even as the breakout region deepens radially. The failed area around the borehole becomes smaller and more localized, as the breakout tip sharpens and deepens. Using the Mogi–Coulomb failure criterion, a good match is obtained between the modeled breakout geometry and the geometry observed by Herrick and Haimson in laboratory experiments on an Alabama limestone. The new method leads to a correlation between breakout geometry, rock strength properties, and in situ stress. The paper ends with a critical discussion of the possibility of inferring the in situ stress state from observed breakout geometries.

This paper presents an experimental investigation into the structural and material response of ambient-dry and wet clay-brick/lime-mortar masonry elements. In addition to cyclic tests on four large-scale masonry walls subjected to lateral in-plane displacement and co-existing compressive gravity load, the study also includes complementary tests on square masonry panels under diagonal compression and cylindrical masonry cores in compression. After describing the specimen details, wetting method and testing arrangements, the main results and observations are provided and discussed. The results obtained from full-field digital image correlation measurements enable a detailed assessment of the material shear-compression strength envelope, and permit a direct comparison with the strength characteristics of structural walls. The full load-deformation behaviour of the large-scale walls is also evaluated, including their ductility and failure modes, and compared with the predictions of available assessment models. It is shown that moisture has a notable effect on the main material properties, including the shear and compression strengths, brick–mortar interaction parameters, and the elastic and shear moduli. The extent of the moisture effects is a function of the governing behaviour and material characteristics as well as the interaction between shear and precompression stresses, and can lead to a loss of more than a third of the stiffness and strength. For the large scale wall specimens subjected to lateral loading and co-existing compression, the wet-to-dry reduction was found to be up to 20% and 11% in terms of stiffness and lateral strength, respectively, whilst the ductility ratio diminished by up to 12%. Overall, provided that the key moisture-dependent material properties are appropriately evaluated, it is shown that analytical assessment methods can be reliably adapted for predicting the response, in terms of the lateral stiffness, strength and overall load-deformation, for both dry and wet masonry walls.

Knowledge regarding microseismic characteristics associated with longwall coal mining is crucial in evaluating the potential for underground mining hazards. Although microseismicity is induced by mining activities, it still remains uncertain as to what extent mining activities influence the spatial, temporal, and magnitude characteristics of microseismicity. To establish a thorough understanding of the relationship between microseismic characteristics and mining activities, a 27-month long microseismic monitoring campaign was conducted around a highly stressed coal zone and eight producing longwall panels at Coal Mine Velenje in Slovenia. Each microseismic event was classified to be associated with the producing longwall panel that triggered it, and the microseismic response to multi-panel longwall top coal caving face advance was analysed. Monitoring data have shown that locations of microseismic events coincided with stress concentrated regions. It was established that both seismic count and energy-intensive regions associated with coal mining in different panels are spatially connected, but they do not fully overlap with mined-out or stress concentrated areas. In addition, microseismic event counts frequency was found to be well correlated with mining intensity, while seismic energy magnitude and spatial distribution are poorly correlated with the same. Therefore, microseismic characteristics could not be explained solely by the mining-induced stress transfer and mining intensity, but are believed to be dominated by pre-existing natural fractures throughout the coal seam. Analyses of these observations helped the development of a conceptual seismic-generation model, which provides new insights into the causes of microseismicity in coal mining.

The final publication is available at Springer via http://dx.doi.org/10.1007/s10237-020-01291-x This work explores the use of an embedded computational fluid dynamics method to study the type B aortic dissection. The use of the proposed technique makes it possible to easily test different intimal flap configurations without any need of remeshing. To validate the presented methodology, we take as reference test case an in vitro experiment present in the literature. This experiment, which considers several intimal flap tear configurations (number, size and location), mimics the blood flow in a real type B aortic dissection. We prove the correctness and suitability of the presented approach by comparing the pressure values and waveform. The obtained results exhibit a remarkable similarity with the experimental reference data. Complementary, we present a feasible surgical application of the presented computer method. The aim is to help the clinicians in the decision making before the type B aortic dissection surgical fenestration. The capabilities of the proposed technique are exploited to efficiently create artificial reentry tear configurations. We highlight that only the radius and center of the reentry tear need to be specified by the clinicians, without any need to modify neither the model geometry nor the mesh. The obtained computational surgical fenestration results are in line with the medical observations in similar clinical studies.

In granular soils grain crushing reduces dilatancy and stress obliquity enhances crushability. These are well-supported specimen-scale experimental observations. In principle, those observations should reflect some peculiar micromechanism associated with crushing, but which is it? To answer that question the nature of crushing-induced particle-scale interactions is here investigated using an efficient DEM model of crushable soil. Microstructural measures such as the mechanical coordination number and fabric are examined while performing systematic stress probing on the triaxial plane. Numerical techniques such as parallel and the newly introduced sequential probing enable clear separation of the micromechanical mechanisms associated with crushing. Particle crushing is shown to reduce fabric anisotropy during incremental loading and to slow fabric change during continuous shearing. On the other hand, increased fabric anisotropy does take more particles closer to breakage. Shear-enhanced breakage appears then to be a natural consequence of shear-enhanced fabric anisotropy. The particle crushing model employed here makes crushing dependent only on particle and contact properties, without any pre-established influence of particle connectivity. That influence does not emerge, and it is shown how particle connectivity, per se, is not a good indicator of crushing likelihood.

This paper deals with the estimation of fragility functions for acceleration-sensitive components of buildings subjected to earthquake action. It considers ideally coherent pulses as well as real non-pulselike ground-motion records applied to continuous building models formed by a flexural beam and a shear beam in tandem. The study advances the idea of acceleration-based dimensionless fragility functions and describes the process of their formulation. It demonstrates that the mean period of the Fourier Spectrum, $$T_m$$ , is associated with the least dispersion in the predicted dimensionless mean demand. Likewise, peak ground acceleration, PGA-, and peak ground velocity, PGV-based length scales are found to be almost equally appropriate for obtaining efficient ‘universal’ descriptions of maximum floor accelerations. Finally, this work also shows that fragility functions formulated in terms of dimensionless $$\varPi$$ -terms have a superior performance in comparison with those based on conventional non-dimensionless terms (like peak or spectral acceleration values). This improved efficiency is more evident for buildings dominated by global flexural type lateral deformation over the whole intensity range and for large peak floor acceleration levels in structures with shear-governed behaviour. The suggested dimensionless fragility functions can offer a ‘universal’ description of the fragility of acceleration-sensitive components and constitute an efficient tool for a rapid seismic assessment of building contents in structures behaving at, or close to, yielding which form the biggest share in large (regional) building stock evaluations.

© 2018, Springer Nature B.V. The Groningen gas field is one of the largest in the world and has produced over 2000 billion m3 of natural gas since the start of production in 1963. The first earthquakes linked to gas production in the Groningen field occurred in 1991, with the largest event to date being a local magnitude (ML) 3.6. As a result, the field operator is leading an effort to quantify the seismic hazard and risk resulting from the gas production operations, including the assessment of liquefaction hazard. However, due to the unique characteristics of both the seismic hazard and the geological subsurface, particularly the unconsolidated sediments, direct application of existing liquefaction evaluation procedures is deemed inappropriate in Groningen. Specifically, the depth-stress reduction factor (rd) and the magnitude scaling factor relationships inherent to existing variants of the simplified liquefaction evaluation procedure are considered unsuitable for use. Accordingly, efforts have first focused on developing a framework for evaluating the liquefaction potential of the region for moment magnitudes (M) ranging from 3.5 to 7.0. The limitations of existing liquefaction procedures for use in Groningen and the path being followed to overcome these shortcomings are presented in detail herein.

The main constrain for effective gas drainage in coal mines is the low permeability nature of coal reservoirs. As coal mining activities are extending to deeper subsurface, the ever-increasing in situ stress conditions is anticipated to result in much lower permeability and more challenges for gas emission control in coal mines. In recent years, hydraulic slotting using high-pressure waterjet along underground gas drainage boreholes, as a general solution to stimulate low permeability coal seams, has become increasingly favourable. This paper presents a systematic investigation into the sensitivity of borehole slotting performance to a number of field and operational parameters. A wide range of geomechanical properties, in situ stress conditions, slot geometry and spacing of multiple slots were considered in a series of numerical simulations. The relations between these key parameters and the failure zone size/volume induced by slotting were quantified. The effect of different parameters in improving slotting performance has also been ranked, which provides theoretical base for mine operators to optimise slotting operations.

The growth of fractures around an undercut of a block cave is simulated. A finite element based approach is used, in which fractures are represented as non-planar 3D surfaces that grow in response to boundary stresses and interaction. A new mesh is recreated at each step to compute the displacement field. Stress intensity factors are computed around fracture tips using a technique that computes the interaction integral over a virtual disk. Fracture geometry is updated using Paris and Schöllmann propagation laws, and a geometric fracture pattern ensues from the simulation. The growth of fractures is examined in the lower 20 m of a mine at 812 m depth. The growth of 30, 60 and 90 fractures is examined. Realistic extraction schedules for over 100 draw points control the rate of mass extraction. The effect of rock bridges as overburden stress shields is investigated. Bridges are modelled by constraining the vertical displacement of the top boundary. This case is compared to a Neumann-type overburden stress boundary condition in which the overburden is felt throughout the top of the cave. In both cases, fractures grow to form a dome shape above and around the cave during extraction. For the case of a fixed top boundary, fracture growth is observed away from the cave, while in the direct overburden stress case, fractures tend to grow close to the cave. Over-arching fractures concentric to the undercut continue to grow as the cave progresses.

In Japan and Turkey, predominantly traditional design concept of pile group is still adopted even though most pile foundations are piled raft foundations in reality. This situation may be attributed to the still limited knowledge about behaviour of piled rafts subjected to the combination of static and dynamic loads. Hence, in this study, a series of shaking table tests on 3-pile piled raft and 3-pile pile group models were carried out in saturated sand at 1-g field under the combinations of static vertical and dynamic horizontal loading, where the input motion was sinusoidal wave having a frequency of 20 Hz and 1.50 and 6.00 m/s(2) amplitudes. The results showed that, raft-model ground contact in piled raft is very effective to reduce settlement compared to pile group even in liquefied soil. Moreover, even in piled raft, the horizontal resistance reduces if liquefaction of the ground surface occurs, due to the loss of the raft base resistance.

An empirical joint constitutive model (JCM) that captures the rough wall interaction behaviour of individual fractures associated with roughness characteristics observed in laboratory experiments is combined with the solid mechanical model of the finite-discrete element method (FEMDEM). The combined JCM-FEMDEM formulation gives realistic fracture behaviour with respect to shear strength, normal closure, and shear dilatancy and includes the recognition of fracture length influence as seen in experiments. The validity of the numerical model is demonstrated by a comparison with the experimentally established empirical solutions. A 2D plane strain geomechanical simulation is conducted using an outcrop-based naturally fractured rock model with far-field stresses loaded in two consecutive phases, i.e. take-up of isotropic stresses and imposition of two deviatoric stress conditions. The modelled behaviour of natural fractures in response to various stress conditions illustrates a range of realistic behaviour including closure, opening, shearing, dilatancy, and new crack propagation. With the increase in stress ratio, significant deformation enhancement occurs in the vicinity of fracture tips, intersections, and bends, where large apertures can be generated. The JCM-FEMDEM model is also compared with conventional approaches that neglect the scale dependency of joint properties or the roughness-induced additional frictional resistance. The results of this paper have important implications for understanding the geomechanical behaviour of fractured rocks in various engineering activities., Rock Mechanics and Rock Engineering, 49 (12), ISSN:1434-453X, ISSN:0723-2632

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $$\text{ V }_\mathrm{S30}$$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $$\text{ M }_\mathrm{w}$$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $$\text{ R }_\mathrm{epi}$$ , and hypocentral distance, $$\text{ R }_\mathrm{hyp}$$ ) and finite-fault (distance to the surface projection of the rupture, $$\text{ R }_\mathrm{JB}$$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.