Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inf...Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inferior by design but complete by function. To understand the loading response, this paper numerically investigated structural remodeling of mature sheep femur using evolutionary structural optimization method (ESO). Femur images from Computed Tomography scanner were used to determine the elastic modulus variation and subsequently construct finite element model of the femur with stiffest elasticity measured. Major muscle forces on dominant phases of healthy sheep gait were imposed on the femur under static mode. ESO was applied to progressively alter the remodeling of numerically simulated femur from its initial to final design by iteratively removing elements with low strain energy density (SED). The computations were repeated with two different mesh sizes to test the convergence. The elements within the medullary canal had low SEDs and therefore were removed during the optimization. The SEDs in the remaining elements varied with angle around the circumference of the shaft. Those elements with low SED were inefficient in supporting the load and thus fundamentally explained how bone remodels itself with less stiff inferior tissue to meet load demand. This was in line with the Wolff’s law of transformation of bone. Tissue growth and remodeling process was found to shape the sheep femur to a mechanically optimized structure and this was initiated by SED in macro-scale according to traditional principle of Wolff’s law.展开更多
The most commonly used method for assessing the hydraulic erodibility of rock is Annandale's method.This method is based on a correlation between the erosive force of flowing water and the capacity of rock resista...The most commonly used method for assessing the hydraulic erodibility of rock is Annandale's method.This method is based on a correlation between the erosive force of flowing water and the capacity of rock resistance. This capacity is evaluated using Kirsten's index, which was initially developed to evaluate the excavatability of earth materials. For rocky material, this index is determined according to certain geomechanical factors related to intact rock and rock mass, such as compressive strength of intact rock, rock block size, discontinuity shear strength and relative block structure. To quantify the relative block structure, Kirsten(1982) developed a mathematical expression that accounts for the shape and orientation of the blocks relative to the direction of flow. Kirsten's initial concept for assessing the relative block structure considers that the geological formation is mainly fractured by two joint sets forming an orthogonally fractured system. An adjusted concept is proposed to determine the relative block structure when the fractured system is non-orthogonal where the angle between the planes of the two joint sets is greater or less than 90°. An analysis of the proposed relative block structure rating shows that considering a non-orthogonally fractured system has a significant effect on Kirsten's index and, as a consequence, on the assessment of the hydraulic erodibility of rock.展开更多
Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in ...Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.展开更多
Shale samples of Silurian Longmaxi Formation in the Changning area of the Sichuan Basin, SW China, were selected to carry out scanning electron microscopy, CT imaging, high-pressure mercury injection, low-temperature ...Shale samples of Silurian Longmaxi Formation in the Changning area of the Sichuan Basin, SW China, were selected to carry out scanning electron microscopy, CT imaging, high-pressure mercury injection, low-temperature nitrogen adsorption and imbibition experiments to compare the hydration characteristics of montmorillonite and illite, analyze the main factors affecting the water block removal of shale, and reveal the mechanisms of pore structure evolution during shale hydration. The hydration characteristics of shale are closely related to the composition of clay minerals, the shale with high illite content is not susceptible to hydration and thus has limited room for pore structure improvement;the shale with high montmorillonite is susceptible to hydration expansion and thus has higher potential of pore structure improvement by stimulation;the shale with high illite content has stronger imbibition in the initial stage, but insufficient diffusion ability, and thus is likely to have water block;the shale with high montmorillonite content has weaker imbibition in the initial stage but better water diffusion, so water blocking in this kind of shale can be removed to some degree;the shale reservoir has an optimal hydration time, when it is best in physical properties, but hydration time too long would cause damage to the reservoir, and the shale with high illite content has a shorter optimal hydration time;inorganic cations can inhibit the hydration of clay minerals and have stronger inhibition to illite expansion, especially K^(+);for the reservoir with high content of montmorillonite, the cation content of fracturing fluid can be lowered to promote the shale hydration;fracturing fluid with high K^(+) content can be injected into reservoirs with high illite content to suppress hydration.展开更多
In this paper, we propose and analyze a mathematical model to study the dynamics of a fishery resource system with stage structure in an aquatic environment that consists of two zones namely unreserved zone (fishing p...In this paper, we propose and analyze a mathematical model to study the dynamics of a fishery resource system with stage structure in an aquatic environment that consists of two zones namely unreserved zone (fishing permitted) and reserved zone (fishing is strictly prohibited). In this model we introduce a stage structure in which predators are split into two kinds as immature predators and mature predators. It is assumed that immature predators cannot catch the prey and their foods are given by their parents (mature predators). It is also assumed that the fishing of immature predators prohibited in the unreserved zone and predator species are not allowed to enter inside the reserved zone. The local and global stability analysis has been specified. Biological and Bionomical equilibriums of the system are derived. Mathematical formulation of the optimal harvesting policy is given and its solution is derived in the equilibrium case by using Pontryagin’s maximum principle.展开更多
文摘Live bone inherently responds to applied mechanical stimulus by altering its internal tissue composition and ultimately biomechanical properties, structure and function. The final formation may structurally appear inferior by design but complete by function. To understand the loading response, this paper numerically investigated structural remodeling of mature sheep femur using evolutionary structural optimization method (ESO). Femur images from Computed Tomography scanner were used to determine the elastic modulus variation and subsequently construct finite element model of the femur with stiffest elasticity measured. Major muscle forces on dominant phases of healthy sheep gait were imposed on the femur under static mode. ESO was applied to progressively alter the remodeling of numerically simulated femur from its initial to final design by iteratively removing elements with low strain energy density (SED). The computations were repeated with two different mesh sizes to test the convergence. The elements within the medullary canal had low SEDs and therefore were removed during the optimization. The SEDs in the remaining elements varied with angle around the circumference of the shaft. Those elements with low SED were inefficient in supporting the load and thus fundamentally explained how bone remodels itself with less stiff inferior tissue to meet load demand. This was in line with the Wolff’s law of transformation of bone. Tissue growth and remodeling process was found to shape the sheep femur to a mechanically optimized structure and this was initiated by SED in macro-scale according to traditional principle of Wolff’s law.
基金the organizations that have funded this project:the Natural Sciences and Engineering Research Council of Canada(Grant No.498020-16)Hydro-Quebec(NC-525700)Mitacs Accelerate Program(Grant Ref.IT10008)
文摘The most commonly used method for assessing the hydraulic erodibility of rock is Annandale's method.This method is based on a correlation between the erosive force of flowing water and the capacity of rock resistance. This capacity is evaluated using Kirsten's index, which was initially developed to evaluate the excavatability of earth materials. For rocky material, this index is determined according to certain geomechanical factors related to intact rock and rock mass, such as compressive strength of intact rock, rock block size, discontinuity shear strength and relative block structure. To quantify the relative block structure, Kirsten(1982) developed a mathematical expression that accounts for the shape and orientation of the blocks relative to the direction of flow. Kirsten's initial concept for assessing the relative block structure considers that the geological formation is mainly fractured by two joint sets forming an orthogonally fractured system. An adjusted concept is proposed to determine the relative block structure when the fractured system is non-orthogonal where the angle between the planes of the two joint sets is greater or less than 90°. An analysis of the proposed relative block structure rating shows that considering a non-orthogonally fractured system has a significant effect on Kirsten's index and, as a consequence, on the assessment of the hydraulic erodibility of rock.
文摘Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.
基金Supported by the Science and Technology Planning Project of Sichuan Province,China(2020YJ0135)National Natural Science Foundation of China(51874250)。
文摘Shale samples of Silurian Longmaxi Formation in the Changning area of the Sichuan Basin, SW China, were selected to carry out scanning electron microscopy, CT imaging, high-pressure mercury injection, low-temperature nitrogen adsorption and imbibition experiments to compare the hydration characteristics of montmorillonite and illite, analyze the main factors affecting the water block removal of shale, and reveal the mechanisms of pore structure evolution during shale hydration. The hydration characteristics of shale are closely related to the composition of clay minerals, the shale with high illite content is not susceptible to hydration and thus has limited room for pore structure improvement;the shale with high montmorillonite is susceptible to hydration expansion and thus has higher potential of pore structure improvement by stimulation;the shale with high illite content has stronger imbibition in the initial stage, but insufficient diffusion ability, and thus is likely to have water block;the shale with high montmorillonite content has weaker imbibition in the initial stage but better water diffusion, so water blocking in this kind of shale can be removed to some degree;the shale reservoir has an optimal hydration time, when it is best in physical properties, but hydration time too long would cause damage to the reservoir, and the shale with high illite content has a shorter optimal hydration time;inorganic cations can inhibit the hydration of clay minerals and have stronger inhibition to illite expansion, especially K^(+);for the reservoir with high content of montmorillonite, the cation content of fracturing fluid can be lowered to promote the shale hydration;fracturing fluid with high K^(+) content can be injected into reservoirs with high illite content to suppress hydration.
文摘In this paper, we propose and analyze a mathematical model to study the dynamics of a fishery resource system with stage structure in an aquatic environment that consists of two zones namely unreserved zone (fishing permitted) and reserved zone (fishing is strictly prohibited). In this model we introduce a stage structure in which predators are split into two kinds as immature predators and mature predators. It is assumed that immature predators cannot catch the prey and their foods are given by their parents (mature predators). It is also assumed that the fishing of immature predators prohibited in the unreserved zone and predator species are not allowed to enter inside the reserved zone. The local and global stability analysis has been specified. Biological and Bionomical equilibriums of the system are derived. Mathematical formulation of the optimal harvesting policy is given and its solution is derived in the equilibrium case by using Pontryagin’s maximum principle.
