Search Documents

By Rudrajit Mitra, Tien Dung Le, Chengguo Zhang, Bruce Hebblewhite, Joung Oh

"Geomechanical and geotechnical characteristics of coal seam are important parameters that directly affect the cavability of top coal in the longwall top-coal caving (LTCC) method. The impact of these parameters such as coal strength, coal deformation modulus, coal discontinuities and top-coal thickness on cavability has not been fully understood in previous numerical analyses. The limitation and selection of the numerical modelling codes were two of the main reasons. In this paper, detailed caving simulations are performed using a 2D distinct element method (DEM) code, UDEC, to study the effect of various critical parameters on top-coal cavability. The numerical analysis evidently shows that the coal deformation modulus, coal discontinuities, and top-coal thickness have considerable influences on cavability. Coal strength is also a factor but to a lesser extent. The results presented in this paper contribute to the understanding of top-coal cavability in different mining environments and facilitate the development of an advanced cavability evaluation criterion. INTRODUCTION The longwall top-coal caving (LTCC) method was originally developed in the 1950s in the former Soviet Union and France (Tien, 1998). The method was nearly abandoned in Europe after the mid-1980s; however, it has been significantly improved and widely applied in China since the late 1980s. A conceptual model of the LTCC method is shown in Figure 1 in which a thick coal seam is divided into two sections including the lower or cutting section, and the upper or top coal section. The lower section of coal seam is mined by conventional longwall method. Meanwhile, the upper section is recovered by caving through the face support. The LTCC method is considered to be the most viable option for mining general thick seams. Compared to the multi-slice longwall (MSL) method, LTCC has the major advantages (Zhongming, 2001; cited by Vakili, 2009):"

Jan 1, 2017

By B. Mishra, Q. Shi

Numerical simulations have been used to investigate roof fall in the past decades. Comparing to the continuum based methods, the discrete element method has been well recognized for its excellent performance in simulating large strain problems and fracturing processes in the rock mass. However, the discrete element method has its limitations. When simulating a large-scale problem, a large number of particles are required, which consumes significant computational resources. In the present study, a 3D numerical model coupled with the finite-difference method (FDM) and the discreteelement method (DEM) was created to simulate the laminated roof failure in an eastern coal mine in the US. In the FDM-DEM coupled model, the laminated roof of the entry is represented by an assembly of bonded particles using PFC3D. The laminated roof is simulated by adding parallel discontinuities. The discontinuities are given certain stiffness and strength in the model. The ribs, floor, and far-field surrounding rock are represented by continuum zones using FLAC3D. The strength of the discontinuities is varied sequentially and the fracturing of the laminated roof was analyzed with respect to entry advance. The results showed that the laminated roof fractures along the discontinuities while an intact roof fractures in a dome shape. Increasing the strength of discontinuity could enlarge the tensile zone inside the laminated roof and hence change the crack distribution. In addition, the crack quantity in the laminated roof behind the advancing face decreases with the increment of discontinuity strength while this trend is not occurring ahead of the advancing face. The simulation demonstrated that coupled method using PFC3D/FLAC3D can reproduce the laminated roof failure in underground coal mines.

Mar 2, 2022

By C. O. Karacan, S. J. Schatzel, W. P. Diamond

In coal mining, longwall mining is a preferred method to maximize production by extracting large blocks of coal that have been outlined with a set of development entries. In U.S. mines, longwall panels are typically 305 m (1,000ft) wide (with a continual trend towards even wider panels), and usually over 3,050 m (10,000ft) in length. The increasing size of longwall panels, while helping to increase coal production, may also increase methane emissions due to the exposure of the mining environment to a larger area of fractured, gas-bearing strata. Thus, understanding the impact of increased panel widths on methane emissions, and designing gob gas ventholes and bleeder systems accordingly, can enhance the safety of the underground workforce by reducing their exposure to potentially explosive accumulations of methane/air mixtures. As part of its mine safety research program, NIOSH’s Pittsburgh Research Laboratory has initiated a reservoir modeling effort to better understand the interaction of the various geotechnical factors influencing gas flow within and to the underground longwall mining environment. A focus of this modeling effort has been 1) the prediction of the incremental amount of methane emissions to be expected due to increasing longwall panel widths, and 2) optimizing gob gas venthole completion practices to capture more of the gas in the subsided strata above longwall panels before it can enter the ventilation system of the underground workplace. The history matching phase of the study has been carried out on a 381-m (1,250ft) wide panel in the Pittsburgh Coalbed. Simulations have also been completed for a 442 m (1,450ft) wide panel to estimate the incremental increase in methane emissions and the performance of current gob gas venthole placement configurations. Additional venthole placement configurations were simulated to investigate options for optimizing methane capture from the subsided strata above the wider longwall panel.

By S. J. Schatzel, W. P. Diamond

High methane emissions originating from the active face areas and from the fractured formations overlying and underlying the mined coalbed can adversely affect both safety and productivity in underground coal mines. Since ventilation alone may not be sufficient to control the methane levels in the longwall mining environment, gob gas ventholes have become a standard supplementary methane control option in many mines. As mines progress into deeper and gassier coalbeds, or as longwall panel size increases, ventilation and gob gas ventholes together may not be sufficient to maintain methane levels within statutory limits. To decrease the risk associated with methane emissions under these circumstances, in-seam horizontal methane drainage is often used to reduce the gas content of the coalbed prior to mining. Horizontal methane drainage borehole completion designs, drilling strategies, and degasification lead times may need to be adjusted for site-specific conditions due to mine design, geology, and the gas content of the coalbed. This study investigates different horizontal methane drainage borehole patterns, borehole lengths, and degasification times prior to and during panel extraction to evaluate their effectiveness in reducing methane emissions using a "dynamic" 3D reservoir modeling of a 38 1-m wide longwall panel operating in the Pittsburgh coalbed. Results of this study showed that dual and tri-lateral boreholes are more effective in decreasing emissions and in shielding the entries compared to fewer shorter, cross-panel, horizontal boreholes parallel to the longwall face. Modeling results showed that after 12 months of pre-mining methane drainage, the average longwall face emission rates can be reduced by as much as 10.3 m3/min and 6.8 m3/min using tri- and dual-lateral boreholes, respectively. It was also shown that if pre-mining methane drainage time is short, it is important to continue methane drainage during the panel extraction to maximize reductions in longwall face emissions since additional face emission reductions achieved during this period can be comparable to pre-mining degasification.

By Yves Fautrelle, Kader Zaidat, Hervé Combeau, Miha Zaloznik, Mohammed M'Hamdi, Bachir Saadi, Xiaodong Wang, Knut Omdal Tveito

"Numerical simulations of macrosegregation formation during horizontal solidification of a hypoeutectic Sn-Pb alloy in a rectangular cavity is presented and compared to experimental results. The experiment involves a melting phase, a holding stage at constant temperature with/without application of electromagnetic stirring, and finally solidification. We simulate all stages of the experiment and study the segregation of Pb throughout the experiment. A twophase volume averaged model is used for the numerical simulations, accounting for thermosolutal convection and assuming perfect microscopic mixing (lever rule). The influence of the melting phase and the electromagnetic stirring is studied with different cases to illustrate the consequences on the flow pattern and Pb segregation. The numerical simulations are compared to the experimental data and the impact of chemical inhomogeneities prior to solidification on macrosegregation formation is discussed.IntroductionMany extensive investigations have been performed so far on the solidification of metallic alloys. However, the melting process of metallic alloy has not attracted attention in a same manner. Nevertheless, melting of alloy may put forth peculiar effects as described in the present paper. It has been observed from a benchmark experiment [1, 2] that significant stratification might arise during the melting process of Sn-Pb alloy. That effect was put forth thanks to the analysis of the experimental temperature distribution which indicated that the convective fluid flow was almost zero. The effect was enhanced when the lead composition was higher (above 6wt.%). The present paper is aimed at confirming that the melting of an alloy containing a heavier solute may be accompanied with a stratification of the solute which exerts a stabilizing effect on the convective fluid flow. Results obtained from numerical modeling are compared with experimental data obtained from the benchmark facility."

Jan 1, 2012

By Brunella Recinto

Pile loading tests are carried out with two main objectives: 1) evaluating the ultimate resistance; 2) checking the behaviour of the pile under live load. In the first case and particularly for large diameter bored piles the cost of a loading test may be very high. In the late '80s, a new method to perform static loading tests was proposed. It is based on the use of a sacrificial hydraulic jack placed at the pile toe. The device, internally pressurized, expanding applies an upward force on the pile shaft and an equal but opposite downward force on the soil beneath the pile toe. This device is well known as Osterberg's load cell, shortly O-cell. The obtained results are typically summarized in two separate load-movement curves: one for the shaft side shear resistance and another for the toe resistance. The two curves are used to build what is generally named the "equivalent head load-movement curve" relying upon two hypotheses: - the pile is rigid; - the load movement relationship for the shaft is not affected by the direction of movement. The former assumption can be easily removed, while it is not possible to remove the latter. On the validity of this latter hypothesis, comparisons between conventional head loading tests and toe loading tests could be very helpful. Unfortunately in the literature there are only a few case histories and most of them are also poorly documented. In order to explore limitations and potentialities of the new testing method, a finite element model has been set up. Results of simulations of conventional head loading tests and of O-cell tests are analyzed, discussed and compared.

Jan 1, 2003

By Trevor I. Addenbrooke

Stress reducing piles were used for the first time at the Queen Elizabeth II Conference Centre near the Palace of Westminster in London, UK. The inclusion of these piles directly beneath heavily loaded columns was to reduce the high anticipated bending stress in the foundation raft, and to eliminate local yield of the London Clay beneath the column locations. Load cells were placed in one of the stress reducing piles during construction. The data showed that the pile carried about 30% of the column load. This can be contrasted with conventional piled raft design where the pile carries upwards of 80% of the column load. The build up of load on the pile and pile displacement have been accurately reproduced using coupled finite element analysis. These numerical results demonstrate the value of geometrically simple axially symmetric analyses for modelling a pile-enhanced raft. On site, excavation preceded construction of the raft and subsequent pile loading. Failure to model this construction sequence caused the shear strength to mobilize at the base of the pile first, during loading, whereas following the construction sequence caused the shear strength to mobilize at the top of the pile first.

Jan 1, 2002

By M. Gimpel, R. Galler, T. Antretter, H. Kargl, M. Mikl-Resch, R. Tichy, W. Ecker, G. Pittino

"A computational method based on a continuum-mechanical approach to simulate rock cutting is developed and evaluated. The aim is to obtain the local loading acting at the cutting pick during operation. The explicit FEM is chosen to model the pick-rock interaction. The presented method takes advantage of (i) an advanced material model (ii) a robust contact algorithm and (iii) fast calculations due to parallelization. A damage-plasticity law is chosen and modified with user-subroutines. A user integrated element deletion routine evades critical mesh distortions. Infinite elements at the model boundaries avoid unrealistic deformation wave reflection.INTRODUCTIONExcavation in mining and tunneling is a costly and complex process. Special equipment such as the roadheader (fig. 1) is employed to cut the rock in the excavation process but high material strengths limit their usability (Tockner 1999). Nowadays a challenge is to extend their field of application to harder and more abrasive rock materials. The forces on the pick (fig. 2) are of great importance and determine wear and especially the efficiency of the machine. Figure 2 shows the principle of cutting and the predominant forces. The variation of pick geometry, cutting depth or attack angle will result in an alteration of these forces and determine the efficiency of the process (Tockner 1999). The mechanisms of rock cutting follow the principle of indentation (Kou et al. 1999). Directly under the indenter a hydrostatic stress state evolves causing the material to fail only due to crushing and fragmentation mechanisms. In the close surroundings large strains emerge and tensile cracks are formed but in a further distance the material behaves elastically."

Jan 1, 2015

By M. Gimpel, R. Galler, T. Antretter, H. Kargl, R. Tichy, W. Ecker, M. J. Mikl-Resch, G. Pittino

"A computational method based on a continuum-mechanical approach to simulate rock cutting is developed and evaluated. The aim is to obtain the local loading acting at the cutting pick during operation. The explicit FEM is chosen to model the pick-rock interaction. The presented method takes advantage of (i) an advanced material model (ii) a robust contact algorithm and (iii) fast calculations due to parallelization. A damage-plasticity law is chosen and modified with user-subroutines. A user integrated element deletion routine evades critical mesh distortions. Infinite elements at the model boundaries avoid unrealistic deformation wave reflection.INTRODUCTIONExcavation in mining and tunneling is a costly and complex process. Special equipment such as the roadheader (Figure 1) is employed to cut the rock in the excavation process but high material strengths limit their usability (Tockner 1999). Nowadays a challenge is to extend their field of application to harder and more abrasive rock materials. The forces on the pick (Figure 2) are of great importance and determine wear and especially the efficiency of the machine. Figure 2 shows the principle of cutting and the predominant forces. The variation of pick geometry, cutting depth or attack angle will result in an alteration of these forces and determine the efficiency of the process (Tockner 1999). The mechanisms of rock cutting follow the principle of indentation (Kou et al. 1999). Directly under the indenter a hydrostatic stress state evolves causing the material to fail only due to crushing and fragmentation mechanisms. In the close surroundings large strains emerge and tensile cracks are formed but in a further distance the material behaves elastically."

Jan 1, 2015

By B. Mishra

This paper presents a rock mechanics investigation of the time-dependent behavior of the rock immediately surrounding the excavation in a hypothetical one-entry coal mine based on the geology of the Illinois coalfields Research in the area of time-dependency is almost nonexistent and, therefore, this research was performed to make inroads into this neglected research area. Numerical analysis of a typical one-entry coal mine in an Illinois basin was performed using the commercial finite difference program FLAC3D. The immediate roof of the coal entry was assumed to be a shale formation. Shale in other formations has shown a tendency to creep with time. Therefore, in this study the authors modeled the shale in the roof using time-dependent characteristics. The input properties were taken from laboratory creep tests performed on shale. Model runs were performed at different stages, mimicking the progressive excavation of the entry. The model results showed that stresses in the roof attained the limit of yielding. Using Mohr-Coulomb failure criteria, potential zones of failure were also predicted. The results showed that the current computational capacities are sufficient to explore this previously uncharted territory of ground control.

Mar 1, 2014

By Tian Tang, Sergio D. Felicelli

"The laser cladding is an important surface treatment for enhancing resistance against corrosion, abrasion, and wear. The focus of the present study is to develop a three dimensional finite element model to analyze the temperature history and thermal stress evolution during multi-pass laser cladding process on H13-H13 tool steel sample. A sequentially coupled thermal-stress analysis methodology was adopted in this study. The thermal model was firstly established to investigate the temperature history of the laser cladding process. Then the temperature history was used as input to mechanical model to calculate the thermal stress evolution during the process and the residual stress after the sample completely cooled down to room temperature. The influences of process parameters including the scanning speed and scanning strategies on temperature history, thermal stress evolution, and residual stress were investigated in details.IntroductionLaser cladding is an advanced material process technique, which is often used to improve mechanical properties or increase corrosion resistance, repair worn out parts, and fabricate metal matrix composites. In this process, the injected powder or wire feedstock material and a thin layer of substrate are melted and consolidated by use of a laser in order to coat part of a substrate or fabricate a near-net shape part.Numerous mathematical models have been established for analyzing the effects of process parameters, microstructure evolution, and thermo-mechanical behavior of laser cladding. Zhou et al. [1] established an analytical model to investigate the laser induction hybrid rapid cladding for Ni-based WC composite coatings. Huang et al. [2] also developed an analytical model to theoretically investigate the interaction between the laser beam and the powder stream of laser cladding process. A simple and realistic model of laser cladding was proposed by Picasso et al. [3]. The cladding layer is very sensitive to cracks due to the accumulation of residual stress, caused by non-uniform temperature changes, which have limited the wide application of laser cladding technique in industry. Recently, fmite element method has been increasingly used to simulate the laser cladding process [ 4-6]."

Jan 1, 2014

By R. Chen

Most tunnels built in Sydney encounter similar ground situations; they penetrate either the Hawkesbury Sandstone or the shales that lie beneath the city. Tunnels in these materials are generally sup-ported by rock bolts in combination with a shotcrete lining. A method of analysis is presented to improve our understanding of the long term behaviour of fiber reinforced shotcrete (FRS) materials and the tunnel support systems. A finite element method SAFEA (Semi-analytical FE Analysis) is developed to model the combined creep behaviour of a shotcrete lining and the Hawkesbury Sandstone. To model the time-dependent creep behaviour, the bulk modulus and shear modulus of the materials are represented by logarithmic functions in time. A Laplace transform is then applied to these functions of time and the solution of the problem found from a finite element analysis. Tablot?s method is then used to invert the transformed values and find the solution in real time. The time dependent creep deflection of the numerical tunnel models reflects the interaction of the shotcrete lining and rock which have different creep parameters.

Jan 1, 2008

By Jun-ichiro Yagi

A mathematical model has been developed for simulating the blast furnace operation with carbon composite agglomerates charging based on multi-fluid theory and reaction kinetics. In this model, the behaviours of carbon composite agglomerates are considered on the basis of previously reported conservation equations and experimental research on the reduction rate of carbon composite agglomerates. A series of calculations are performed to examine the effect of charging carbon composite agglomerates. The model calculation gives two-dimensional distributions of process variables and information on the overall operational parameters. It reveals that in-furnace temperature levels significantly decrease and the reduction of carbon composite agglomerates is retarded with carbon composite agglomerates charging ratio. The furnace performance is remarkably improved with the increase in carbon composite agglomerates rate. The productivity tends to increase while coke rate and total reducing agent rate show decreases at the different degree.

Jan 1, 2006

By F. S. Su, P. Z. Pan, X. T. Feng, F. Yan, W. W. Ji

"This paper aims to simulate the spalling failure behavior in hard rock based on a self-developed numerical model, i.e. a rock discontinuous cellular automaton (RDCA). A crack initiation and propagation criterion is proposed based on the plastic strain, which is used to link the continuity and discontinuity during the failure process. The discontinuities of rock, such as cracks and other internal boundaries, are represented by discontinuous displacement functions, i.e. tip enrichment and step enrichment. The partition of unity is used to improve integral precision of the elements intersected internal boundaries, which is traced by a level set method. Arbitrary crack geometries can be conveniently represented and it does not require any remeshing in the modeling of crack initiation, propagation and coalescence. Numerical modeling of excavation in hard rock using RDCA is conducted. The spalling failure process of hard rock induced by excavation is well simulated. The depth and shape of hard rock spalling failure are good agreement with in situ observation.INTRODUCTIONThe excavation in rock mass is more and more popular. At the same time, the disasters, such as rock burst, rock spalling and falling etc, induced by excavation are quite prominent. Numerical modeling is an important way to help us understand the failure mechanism of rocks. Many numerical methods have been used to simulate the excavation of rock mass. For example, the finite element method (FEM) (Karakus and Fowell, 2003), the discrete element method (DEM)(Cundall and Strack, 1979, Labra et al., 2008), Finite difference method (FDM)(Zhu et al., 2009), Element Free Galerkin Method (Belytschko et al., 1994), Discontinuous displacement analysis (DDA) (Shi, 1988) and Boundary element method (BEM) (Katsikadelis, 2002, Chen and Hong, 1992) are commonly found in the stability analysis of excavation in rocks. In the modeling of crack propagation or excavation, most of these methods are tedious since the internal boundaries have to be consistent with the element boundaries. And in many cases, the remesh has to be conducted. In recent years, people are interested in particular types of FEMs, such as the extended finite-element method (XFEM), the particle discretization scheme finite-element method(PDS-FEM(Hori et al., 2005, Oguni et al., 2009)), the combined finite-discrete element method (FEM/DEM) (Munjiza and Wiley, 2004, Cai, 2008) and spectral-element method(Gharti et al., 2012) etc."

Jan 1, 2015

By Go-Gi Lee, Jae-Young Jung, Myung-Duk Seo, Hyun-Na Bae, Seon-Hyo Kim

"Bubbles form at the anode during the electrowinning of magnesium and cause hydrodynamic acceleration and electrical field disturbance. A three-dimensional computational model was developed to investigate the effects of the geometry of the electrochemical chamber and the size distribution of bubbles on the fluid flow in the pilot plant. The results are validated by measurements of the surface characteristics of molten salt and the bubbles formed. Fluid flow at the surface was recorded with a high speed camera, and the distribution of bubble size is determined by inspecting the video images frame by frame. The measured size distribution of the bubbles is entered into the computer simulation of a multi-phase flow. A sloped roof is installed above the electrodes, to enhance the flow of the molten salt, thus resulting in higher current efficiency, lower energy consumption, and higher quality of magnesium metal.IntroductionTitanium's excellent properties of low density, high specific strength, high temperature strength, and exceptional resistance to corrosion are highly beneficial to diverse applications including aerospace, automotive, construction, and medical implants. Although the Kroll process[l] and Hunter process have been used for the extraction of titanium in the 20th century, their complicated engineering stages and high production cost, have led to further research into new processes such as FCC Cambridge, Armstrong and TiRO. Despite these intense research efforts, no process has yet been developed that is more economical than the Kroll process, by which most of the global commercial titanium is still being produced."

Jan 1, 2011

By P. Carlone, A. Astarita, F. Rubino

Cold gas dynamic spray is an innovative technology allowing the deposition of thin metallic layers on a bulk or sheet substrate. Titanium on aluminum deposition is very attractive to the aeronautic industry, due to the enhanced corrosion resistance and wear properties of titanium coating as well as its improved compatibility with CFRP, preserving the reduced cost of aluminum substrate. One main drawback, however, is related to the micro-porosity of the deposed layer, negatively affecting the corrosion barrier performance. In this paper, the effect of selective laser post-treatment on pure grade 2 titanium coatings on AA2020-T3 sheets was experimentally and numerically investigated. Morphological features, microstructure and chemical composition of the treated zone were assessed by means of optical microscopy, scanning electron microscopy and energy dispersive X-ray spectrometry. Microhardness measurements were also carried out to evaluate the mechanical properties of the treated coating. A numerical model of the laser treatment, based on a finite volume scheme, was implemented and solved to simulate the process and discuss the experimental outcomes. Obtained results highlighted the key role played by heat input and dimensional features on the effectiveness of the treatment.

Jan 1, 2015

By S. Bhattacharjee

"Rock failure process is generally recorded or measured using conventional sensors like strain gauges, seismic, micro-seismic, acoustic, electromagnetic and other techniques. The above sensors are either to be pasted on the surface of the rock sample or placed very near to it to capture the response of the sample with the increment of load. This requires precise tools to cut or polish the locations where sensors are to be placed. Moreover, deformation or acoustic or micro-seismic activities are recorded only at the locations where sensors are pasted.Recent development in Digital Image Correlation (DIC) technique has shown promise to measure deformation / strain at every location of the surface of a sample. In DIC, two images (reference and deformed) are taken within a small time interval of the same object are compared to estimate the change in reflectance as a function of displacement, velocity or rate of strain at each pixel while the object is undergoing through a loading process. In this study, extended digital image correlation (X-DIC) technique has been implemented using 4-noded quadrilateral elements with a discontinuity surface in the deformed image. The continuous and enriched displacements are estimated at each node.The X-DIC method is verified using numerically generated images. The study has been further extended to evaluate the performance of X-DIC method on deformation analysis of experimental samples under uni-axial loading condition. The results of this study show that X-DIC method is an effective tool for determination of displacements of the entire surface of a rock sample having a discrete crack or discontinuity in a non-contact manner. This method demonstrates the applicability to investigate rock failure mechanism for the entire sample in a laboratory set up as well as monitoring of pillar failure mechanism in underground mines."

Jan 1, 2015

By F. Guevara

The paper presents results of physical and turbulent three-dimensional numerical simulation of the gas flow in the wind box and reaction shaft of the concentrate burner of a flash furnace. Physical simulation work was carried out in a 1:5 linear scale model, where the velocity field was measured using a Pitot tube. The mathematical model consists of the three-dimensional averaged Navier-Stokes equations, turbulent kinetic energy and dissipation of turbulent kinetic energy using the Renormalization Group (RNG) theory for the k-e turbulence model. Computations were performed using FLUENT, based in a finite- volume code. The geometry of the concentrate burner was separated in two sections: wind box and reaction tower. Results obtained at the exit of the wind box were used as boundary conditions for the reaction shaft. The effect of boundary conditions on the numerical results was also studied. Results of physical and numerical simulations show an asymmetry of the velocity field at the exit of wind box. This asymmetry is possibly due to existence of a separated flow in the back of the exit ring. The agreement of physical and numerical results assess the capacity of the turbulence RNG k-e model to reproduce correctly the behaviour of the velocity field at the exit of the concentrate burner.

Jan 1, 1999

By Mahmoud Shakeri1

Various types of thin-walled structures have been frequently used as energy absorbers to mitigate the adverse effects of impact and hence protect the vehicle structure under consideration when subjected to different loading conditions. Compression tubular structures that buckle in a progressive manner, can be used to advantage, providing inexpensive and versatile, but nevertheless weight and stroke efficient energy absorbers. In this paper quasi-static nonlinear finite element simulations are performed to study the effects of number of corners of the cross-section, wall thickness, plastic hardening rate of material, foam-filler and load angle on deformation mode and mean crushing force of tube. Crush strength increases as the number of corners of the cross-section increases, though it almost saturates for the number of corners beyond 11. Higher plastic hardening rate stabilizes the collapse pattern. Foam-filled tubes seem to have higher critical load angles in oblique loading and display a more stable and gradual decrease in mean load over the transition region. Substantial decrease in the initial peak load and dynamic load at any a given deflection was observed as the load angle increases. The investigation of bitubal arrangements and Multi-cell tubes reveals that these may be preferable to monotubal specimens due to the presence of the inner profiles. In addition, a comparison is conducted between the numerical and experimental results of the mean forces of some special angle elements. Good agreement is obtained.

Oct 1, 2011

By D. D. Tannant

An experimental and numerical study was conducted on the geomechanical behaviour of oil sands during penetration tests in which a steel wedge with an apex angle of either 20 degrees or 30 degrees was forced into compacted oil sands. The boundary conditions simulated in the laboratory tests are similar to shovel penetration into oil sands. The objective of the laboratory tests was to measure the force required to push a steel wedge into compacted oil sands. A wedge penetration test differs from conventional triaxial tests in that much more shear is generated. The laboratory data was needed to improve calibration of micro-mechanical input parameters for a two-dimensional discrete element model of oil sands constructed with Particle Flow Code (PFC). Initial model calibration was performed with triaxial test data and then a PFC model of the wedge penetration test was constructed. The measured force-displacement behaviour of the steel wedge was compared with PFC modelling results. It was found that different combinations of PFC particle and bond parameters could equally well replicate the laboratory triaxial compression data for bitumen-rich oil sands. When the wedge model was run, it was found that the resulting force required to penetrate the wedge into the oil sands was about four to six times higher than that measured. This clearly indicates that further research is required to improve the predictive capability of PFC models of oil sands.

Jan 1, 2002