Correlation between the rock mass properties and maximum horizontal stress:A case study of overcoring stress measurements

Understanding the mechanical properties of the lithologies is crucial to accurately determine the horizontal stress magnitude.To investigate the correlation between the rock mass properties and maximum horizontal stress,the three-dimensional(3D)stress tensors at 89 measuring points determined using an improved overcoring technique in nine mines in China were adopted,a newly defined characteristic parameter C_(ERP)was proposed as an indicator for evaluating the structural properties of rock masses,and a fuzzy relation matrix was established using the information distribution method.The results indicate that both the vertical stress and horizontal stress exhibit a good linear growth relationship with depth.There is no remarkable correlation between the elastic modulus,Poisson's ratio and depth,and the distribution of data points is scattered and messy.Moreover,there is no obvious relationship between the rock quality designation(RQD)and depth.The maximum horizontal stress σ_(H) is a function of rock properties,showing a certain linear relationship with the C_(ERP)at the same depth.In addition,the overall change trend of σ_(H) determined by the established fuzzy identification method is to increase with the increase of C_(ERP).The fuzzy identification method also demonstrates a relatively detailed local relationship betweenσ_H and C_(ERP),and the predicted curve rises in a fluctuating way,which is in accord well with the measured stress data.

An effective stress-based DSC model for predicting hydromechanical shear behavior of unsaturated collapsible soils subjected to initial shear stress

Evaluation of hydromechanical shear behavior of unsaturated soils is still a challenging issue. The time and cost needed for conducting precise experimental investigation on shear behavior of unsaturated soils have encouraged several investigators to develop analytical, empirical, or semi-empirical models for predicting the shear behavior of unsaturated soils. However, most of the previously proposed models are for specimens subjected to the isotropic state of stress, without considering the effect of initial shear stress. In this study, a hydromechanical constitutive model is proposed for unsaturated collapsible soils during shearing, with consideration of the effect of the initial shear stress. The model implements an effective stress-based disturbed state concept (DSC) to predict the stress-strain behavior of the soil. Accordingly, material/state variables were defined for both the start of the shearing stage and the critical state of the soil. A series of laboratory tests was performed using a fully automated unsaturated triaxial device to verify the proposed model. The experimental program included 23 suction-controlled unsaturated triaxial shear tests on reconstituted specimens of Gorgan clayey loess wetted to different levels of suctions under both isotropic and anisotropic stress states. The results show excellent agreement between the prediction by the proposed model and the experimental results.

Rheological properties and concentration evolution of thickened tailings under the coupling effect of compression and shear

Cemented paste backfill(CPB)is a key technology for green mining in metal mines,in which tailings thickening comprises the primary link of CPB technology.However,difficult flocculation and substandard concentrations of thickened tailings often occur.The rheological properties and concentration evolution in the thickened tailings remain unclear.Moreover,traditional indoor thickening experiments have yet to quantitatively characterize their rheological properties.An experiment of flocculation condition optimization based on the Box-Behnken design(BBD)was performed in the study,and the two response values were investigated:concentration and the mean weighted chord length(MWCL)of flocs.Thus,optimal flocculation conditions were obtained.In addition,the rheological properties and concentration evolution of different flocculant dosages and ultrafine tailing contents under shear,compression,and compression-shear coupling experimental conditions were tested and compared.The results show that the shear yield stress under compression and compression-shear coupling increases with the growth of compressive yield stress,while the shear yield stress increases slightly under shear.The order of shear yield stress from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Under compression and compression-shear coupling,the concentration first rapidly increases with the growth of compressive yield stress and then slowly increases,while concentration increases slightly under shear.The order of concentration from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Finally,the evolution mechanism of the flocs and drainage channels during the thickening of the thickened tailings under different experimental conditions was revealed.

Degradation and resynthesis of chlorophyll during increased oxidative stress and prolonged darkness differ between annual and perennial flax(Linum L.)

Among plants,there is considerable variation in lifespan:annuals live less than one year,whereas perennials live for several years,with the longest-living perennial having survived 43,600 years.As proposed by the Disposable Soma Theory,this lifespan variation among plants likely reflects differential investment of limited energy and nutrient resources,with perennials investing more energy and nutrients into biomolecular maintenance compared to annuals in order to ensure persistence over multiple seasons.Such differential investment may be particularly important during periods of exogenous stress,which are known to accelerate biomolecular damage.The present study evaluated this hypothesis using annual and perennial flax(Linum L.)subjected to two exogenous stressors—increased oxidative stress(i.e.,foliar H2O2spraying)and complete prolonged darkness.As chlorophyll has been shown to exhibit degradation in response to changes in environmental conditions,we utilized changes in chlorophyll levels during and after periods of exogenous stress to evaluate our hypotheses.We predicted that i)perennials would exhibit a slower rate of chlorophyll degradation during exposure to exogenous stressors compared to annuals,and ii)perennials would exhibit a faster rate of chlorophyll resynthesis following such exposure compared to annuals.Chlorophyll levels before,during,and after exposure to both exogenous stressors were measured in two separate trails,once using image colour analysis and once using spectrophotometry.While chlorophyll degradation rates in response to oxidative stress did not differ between annuals and perennials,contrary to our predictions,chlorophyll resynthesis rates following such exposure were significantly higher in perennials,as predicted.When plants were subjected to complete prolonged darkness,chlorophyll degradation rates were significantly lower in perennials than annuals,as predicted;however,when plants were subsequently reintroduced to natural photoperiod,chlorophyll resynthesis rates did not consistently differ between annuals and perennials,though they tended to be higher in the latter,as predicted.Overall,our study illuminates that evolutionary transitions between life history strategies in plants have been accompanied by physiological modifications to chlorophyll dynamics that permit perennial species to better maintain chlorophyll levels—and thus photosynthetic energy acquisition-in the face of exogenous stressors,which likely underlies their capacity to survive for multiple growing seasons.Future studies should explore whether other key biomolecules(e.g.,proteins,DNA)are also better maintained in perennial plants,especially in the face of exogenous stress.

Evaluation of the coefficient of lateral stress at rest of granular materials under repetitive loading conditions

Although the internal stress state of soils can be affected by repetitive loading,there are few studies evaluating the lateral stress(or K_(0))of soils under repetitive loading.This study investigates the changes in K_(0) and directional shear wave velocity(V_(s))in samples of two granular materials with different particle shapes during repetitive loading.A modified oedometer cell equipped with bender elements and a diaphragm transducer was developed to measure the variations in the lateral stress and the shear wave velocity,under repetitive loading on the loading and unloading paths.The study produced the following results:(1)Repetitive loading on the loading path resulted in an increase in the K_(0) of test samples as a function of cyclic loading number(i),and(2)Repetitive loading on the unloading path resulted in a decrease in K_(0) according to i.The shear wave velocity ratio(i.e.V_(s)(HH)/V_(s)(VH),where the first and second letters in parentheses corresponds to the directions of wave propagation and particle motion,respectively,and V and H corresponds to the vertical and horizontal directions,respectively)according to i supports the experimental observations of this study.However,when the tested material was in lightly over-consolidated state,there was an increase in K_(0) during repetitive loading,indicating that it was the initial K_(0),rather than the loading path,which is responsible for the change in K_(0).The power model can capture the variation in the K_(0) of samples according to i.Notably,the K_(0)=1 line acts as the boundary between the increase and decrease in K_(0) under repetitive loading.

A method for extracting the preseismic gravity anomalies over the Tibetan Plateau based on the maximum shear strain using GRACE data

The occurrence of earthquakes is closely related to the crustal geotectonic movement and the migration of mass,which consequently cause changes in gravity.The Gravity Recovery And Climate Experiment(GRACE)satellite data can be used to detect gravity changes associated with large earthquakes.However,previous GRACE satellite-based seismic gravity-change studies have focused more on coseismic gravity changes than on preseismic gravity changes.Moreover,the noise of the north–south stripe in GRACE data is difficult to eliminate,thereby resulting in the loss of some gravity information related to tectonic activities.To explore the preseismic gravity anomalies in a more refined way,we first propose a method of characterizing gravity variation based on the maximum shear strain of gravity,inspired by the concept of crustal strain.The offset index method is then adopted to describe the gravity anomalies,and the spatial and temporal characteristics of gravity anomalies before earthquakes are analyzed at the scales of the fault zone and plate,respectively.In this work,experiments are carried out on the Tibetan Plateau and its surrounding areas,and the following findings are obtained:First,from the observation scale of the fault zone,we detect the occurrence of large-area gravity anomalies near the epicenter,oftentimes about half a year before an earthquake,and these anomalies were distributed along the fault zone.Second,from the observation scale of the plate,we find that when an earthquake occurred on the Tibetan Plateau,a large number of gravity anomalies also occurred at the boundary of the Tibetan Plateau and the Indian Plate.Moreover,the aforementioned experiments confirm that the proposed method can successfully capture the preseismic gravity anomalies of large earthquakes with a magnitude of less than 8,which suggests a new idea for the application of gravity satellite data to earthquake research.

Size-dependent thermomechanical vibration characteristics of rotating pre-twisted functionally graded shear deformable microbeams

A three-dimensional(3D)thermomechanical vibration model is developed for rotating pre-twisted functionally graded(FG)microbeams according to the refined shear deformation theory(RSDT)and the modified couple stress theory(MCST).The material properties are assumed to follow a power-law distribution along the chordwise direction.The model introduces one axial stretching variable and four transverse deflection variables including two pure bending components and two pure shear ones.The complex modal analysis and assumed mode methods are used to solve the governing equations of motion under different boundary conditions(BCs).Several examples are presented to verify the effectiveness of the developed model.By coupling the slenderness ratio,gradient index,rotation speed,and size effect with the pre-twisted angle,the effects of these factors on the thermomechanical vibration of the microbeam with different BCs are investigated.It is found that with the increase in the pre-twisted angle,the critical slenderness ratio and gradient index corresponding to the thermal instability of the microbeam increase,while the critical material length scale parameter(MLSP)and rotation speed decrease.The sensitivity of the fundamental frequency to temperature increases with the increasing slenderness ratio and gradient index,and decreases with the other increasing parameters.Moreover,the size effect can suppress the dynamic stiffening effect and enhance the Coriolis effect.Finally,the mode transition is quantitatively demonstrated by a modal assurance criterion(MAC).

The Effects of the Longitudinal Axis of Loading upon Bending, Shear and Torsion of a Thin-Walled Cantilever Channel Beam

Three aluminium channel sections of US standard extruded dimension are mounted as cantilevers with x-axis symmetry. The flexural bending and shear that arise with applied axial torsion are each considered theoretically and numerically in terms of two longitudinal axes of loading not coincident with the shear centre. In particular, the warping displacements, stiffness and stress distributions are calculated for torsion applied to longitudinal axes passing through the section’s centroid and its web centre. The stress conversions derived from each action are superimposed to reveal a net sectional stress distribution. Therein, the influence of the axis position upon the net axial and shear stress distributions is established compared to previous results for each beam when loading is referred to a flexural axis through the shear centre. Within the net stress analysis is, it is shown how the constraint to free warping presented by the end fixing modifies the axial stress. The latter can be identified with the action of a ‘bimoment’ upon each thin-walled section.

Wave-current bottom shear stresses and sediment re-suspension in the mouth bar of the Modaomen Estuary during the dry season

On the basis of the measurement data pertaining to waves, current, and sediment in February 2012 in the mouth bar of the Modaomen Estuary, the Soulsby formulae with an iterative method are applied to calculating bottom shear stresses （BSS） and their effect on a sediment resuspension. Swell induced BSS have been found to be the most important part of the BSS. In this study, the correlation coefficient between a wavecurrent shear stress and SSC is 0.86, and that between current shear stresses and SSC is only 0.40. The peaks of the SSC are consistent with the height and the BSS of the swell. The swell is the main mechanism for the sediment re-suspension, and the tidal current effect on sediment re-suspension is small. The peaks of the SSC are centered on the high tidal level, and the flood tide enhances the wave shear stresses and the SSC near the bottom. The critical shear stress for sediment re-suspension at the observation station is between 0.20 and 0.30 N/m2. Tidal currents are too weak to stir up the bottom sediment into the flow, but a WCI （wave-current interaction） is strong enough to re-suspend the coarse sediment.

Mathematical Model and Experiment Validation of Fluid Torque by Shear Stress under Influence of Fluid Temperature in Hydro-viscous Clutch

The current design of hydro-viscous clutch（HVC） in tracked vehicle fan transmission mainly focuses on high-speed and high power. However, the fluid torque under the influence of fluid temperature can not be predicted accurately by conventional mathematical model or experimental research. In order to validate the fluid torque of HVC by taking the viscosity-temperature characteristic of fluid into account, the test rig is designed. The outlet oil temperature is measured and fitted with different rotation speed, oil film thickness, oil flow rate, and inlet oil temperature. Meanwhile, the film torque can be obtained. Based on Navier-Stokes equations and the continuity equation, the mathematical model of fluid torque is proposed in cylindrical coordinate. Iterative method is employed to solve the equations. The radial and tangential speed distribution, radial pressure distribution and theoretical flow rate are determined and analyzed. The models of equivalent radius and fluid torque of friction pairs are introduced. The experimental and theoretical results indicate that tangential speed distribution is mainly determined by the relative rotating speed between the friction plate and the separator disc. However, the radial speed distribution and pressure distribution are dominated by pressure difference at the lower rotating speed. The oil film fills the clearance and the film torque increases with increasing rotating speed. However, when the speed reaches a certain value, the centrifugal force will play an important role on the fluid distribution. The pressure is negative at the outer radius when inlet flow rate is less than theoretical flow, so the film starts to shrink which decreases the film torque sharply. The theoretical fluid torque has good agreement with the experimental data. This research proposes a new fluid torque mathematical model which may predict the film torque under the influence of temperature more accurately.

Effects of bottom shear stresses on the wave-induced dynamic response in a porous seabed:PORO-WSSI (shear) model

When ocean waves propagate over the sea floor,dynamic wave pressures and bottom shear stresses exert on the surface of seabed.The bottom shear stresses provide a horizontal loading in the wave-seabed interaction system,while dynamic wave pressures provide a vertical loading in the system.However,the bottom shear stresses have been ignored in most previous studies in the past.In this study,the effects of the bottom shear stresses on the dynamic response in a seabed of finite thickness under wave loading will be examined,based on Biot＇s dynamic poro-elastic theory.In the model,an ＂u-p＂ approximation will be adopted instead of quasi-static model that have been used in most previous studies.Numerical results indicate that the bottom shear stresses has certain influences on the wave-induced seabed dynamic response.Furthermore,wave and soil characteristics have considerable influences on the relative difference of seabed response between the previous model（without shear stresses） and the present model（with shear stresses）.As shown in the parametric study,the relative differences between two models could up to 10% of p0,depending on the amplitude of bottom shear stresses.

A new analytical solution for calculation the displacement and shear stress of fully grouted rock bolts and numerical verifications

In presence of difficult conditions in coal mining roadways, an adequate stabilization of the excavation boundary is required to ensure a safe progress of the construction. The stabilization of the roadways can be improved by fully grouted rock bolt, offering properties optimal to the purpose and versatility in use. Investigations of load transfer between the bolt and grout indicate that the bolt profile shape and spacing play an important role in improving the shear strength between the bolt and the surrounding strata. This study proposes a new analytical solution for calculation displacement and shear stress in a fully encapsulated rock bolt in jointed rocks. The main characteristics of the analytical solution consider the bolt profile and jump plane under pull test conditions. The performance of the proposed analytical solution, for three types of different bolt profile configurations, is validated by ANSYS software. The results show there is a good agreement between analytical and numerical methods. Studies indicate that the rate of displacement and shear stress from the bolt to the rock exponentially decayed. This exponential reduction in displacement and shear stress are dependent on the bolt characteristics such as: rib height, rib spacing, rib width and grout thickness, material and joint properties.

In search of potential earthquake source regions in the Chinese mainland in the light of ambient shear stress field

Earth media are incomplete media.There exist many cracks in it. The achievements of fracture mechanics showthat the strength of the incomplete materials will be much lower than that of the complete materials. We consider that earthquake occurrence is the result of unstable propagation of a crack in crust media in proper conditionand the earthquake rupture is the phenomenon of a failure by fast fracture under applied low shear stress. It hasalready been explained by fracture mechanics.The occurrence of failure by fast fracture is necessarily associated with the presence of high level concentration of local stress and strain. The elastic/plastic stress analysis in cracked pieces by Dugdale indicates that thestate of stress at the tip of a crack takes a very important role to crack propagation. A plastic zone has necessarilyformed in the tip of a crack due to stress concentration. Therefore, the dislocations st the tip of a crack are naturally a plastic displacement, rather than elastic one. The plastic displacement, where τ0 is appliedshear stress which is equivalent to initial or tectonic shear stress when the quake occurs, a is the half length of acrack, It is the rigidity,τy is the yield stresses in shear. The main seismic dislocations take place exactly at theends of the crack where the plastic zone had been formed. SO, a critical assumption is adopted, i. e. we assumethe dislocation D(1,,t) as formula (5) in text. The maximum earthquake dislocation, whereL is the fault length. If p is taken the value in the upper crust, μ=33 GPa; and τy is taken the average valuegiven from laboratories,τy= 30 MPa. Thus, according to observation values of Dmax and L, using the formula,one can estimate the initial shear stresses for large earthquakes. Computations show that the initial shear stressesfor large earthquakes all over the world are about 5-20 MPa which have some differences between regions.We further research the characteristics of source spectra and have derived the dependent relation of bodywave magnitude mb on the shear stress τ0 and seismic moment M, as formula (11)in text. Thus, the formulaprovides a POssibility of computation of large amount of tectonic shear stress values from seismic data. We consider that the tectonic shear stress field is a main factor which controls the earthquake occurrence. The regions withhigh tectonic shear stress values are considered to be prone to occur great earthquakes (Ms>6) and called earthquake hazard regions. Based on this criterion, τ0 values for all earthquakes with mb≥3. 8 all over China since1987 have been computed, and the great earthquake hazard regions with magnitude ranges have been zoned inthe Chinese mainland.During April 1992 -January 31, 1994, there were 9 Ms≥6 earthquakes which occurred in the Chinesemainland, 8 earthquakes of the 9 had fallen into the regions delineated by us prior to the earthquake occurrence,with only one failure. This new approach as a method for medium--term prediction of strong earthquakes hasbeen proved by practice to be an efficient one.It has good physical bases and bright prospect and worth furtherresearch. Received February 7,1994 1 Accepted February 10, 1995.Contribution No. 95A0061, Institute of Geophysics,SSB, China.

Analysis on Interface Shear Stress of Thermally InsulatedOcean Pipelines Under Installation

It has been proved that the thermally insulated ocean pipeline has advantages over the conventional pipe-in-pipe pipeline. The risk of using the thermally insulated pipeline is that the exterior layers covering the steel pipe may be. pulled off if the shear stress on the interface induced by the pullout fore from the tensioner is greater than the binding fore between two neighboring layers during installation. This paper develops a procedure to calculate the shear stress on the interface. The binding force between two neighboring layers can be determined with full scale model tests. The safety of the thermally insulated pipe under installation can then be checked by comparison of the interface shear stress with the binding force.

Influence of maximum principal stress direction on the failure process and characteristics of D-shaped tunnels

To investigate the failure process and characteristics of D-shaped tunnels under different maximum principal stress directions θ, true-triaxial tests were conducted on cubic sandstone samples with a through D-shaped hole. The test results show that the failure process can be divided into 4 periods:calm, buckling deformation, gradual buckling and exfoliation of rock fragment, and formation of a Vshaped notch. With an increase in θ from 0° to 90°, the size of the rock fragments first decreases and then increases, whereas the fractal dimension of the rock fragments first increases and then decreases. Meanwhile, the failure position at the left side shifts from the sidewall to the corner and finally to the floor, whereas the failure position at the right side moves from the sidewall to the spandrel and finally to the roof, which is consistent with the failure position in underground engineering. In addition, the initial vertical failure stress first decreases and then increases. By comparing the results,the failure severities at different maximum principal stress directions can be ranked from high to low in the following order: 90°>60°>30°>45°>0°.

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.

Dynamic modeling and control of extracellular ATP concentration on vascular endothelial cells via shear stress modulation

A new dynamic model for cell-deformation-induced adenosine triphosphate （ATP） release from vascular endothelial cells （VECs） is proposed in this paper to quantify the relationship between the ATP concentration at the surface of VECs and blood flow-induced shear stress. The simulation results demonstrate that ATP concentration at the surface of VECs predicted by the proposed new dynamic model is more consistent with the experimental observations than those by the existing static and dynamic models. Furthermore, it is the first time that a proportional-integral-derivative （PID） feedback controller is applied to modulate extracellular ATP concentration. Three types of desired ATP concentration profiles including constant, square wave and sinusoid are obtained by regulating the wall shear stress under this PID control. The systematic methodology utilized in this paper to model and control ATP release from VECs via adjusting external stimulus opens up a new scenario where quantitative investigations into the underlying mechanisms for many biochemical phenomena can be carded out for the sake of controlling specific cellular events.

Estimation of Bed Shear Stresses in the Pearl River Estuary

Mean and fluctuating velocities were measured by use of a pulse coherent acoustic Doppler profiler （PC-ADP） and an acoustic Doppler velocimeter in the tidal bottom boundary layer of the Pearl River Estuary. The bed shear stresses were estimated by four different methods： log profile （LP）, eddy correlation （EC）, turbulent kinetic energy （TKE）, and inertial dissipation （ID）. The results show that （a） all four methods for estimating bed stresses have advantages and disadvantages, and they should be applied simultaneously to obtain reliable frictional velocity and to identify potential sources of errors; （b） the LP method was found to be the most suitable to estimate the bed stresses in non-stratified, quasi-steady, and homogeneous flows; and （c） in the estuary where the semi-diurnal tidal current is dominant, bed shear stresses exhibit a strong quarter-diurnal variation.

Stationary Flow of Blood in a Rigid Vessel in the Presence of an External Magnetic Field: Considerations about the Forces and Wall Shear Stresses

The magnetohydrodynamics laws govern the motion of a conducting fluid, such as blood, in an externally applied static magnetic field B0. When an artery is exposed to a magnetic field, the blood charged particles are deviated by the Lorentz force thus inducing electrical currents and voltages along the vessel walls and in the neighboring tissues. Such a situation may occur in several biomedical applications: magnetic resonance imaging (MRI), magnetic drug transport and targeting, tissue engineering… In this paper, we consider the steady unidirectional blood flow in a straight circular rigid vessel with non-conducting walls, in the presence of an exterior static magnetic field. The exact solution of Gold (1962) (with the induced fields not neglected) is revisited. It is shown that the integration over a cross section of the vessel of the longitudinal projection of the Lorentz force is zero, and that this result is related to the existence of current return paths, whose contributions compensate each other over the section. It is also demonstrated that the classical definition of the shear stresses cannot apply in this situation of magnetohydrodynamic flow, because, due to the existence of the Lorentz force, the axisymmetry is broken.

The Analysis of Wall Shear Stress Modulated by Acute Exercise in the Human Common Carotid Artery with an Elastic Tube Model

Assessment of the magnitude and pattern of wall shear stress(WSS)in vivo is the prerequisite for studying the quantitative relationship between exercise-induced WSS and arterial endothelial function.In the previous studies,the calculation of the WSS modulated by exercise training was primarily based upon the rigid tube model,which did not take non-linear effects of vessel elastic deformation into consideration.In this study,with an elastic tube model,we estimated the effect of a bout of 30-minute acute cycling exercise on the WSS and the flow rate in the common carotid artery according to the measured inner diameter,center-line blood flow velocity,heart rates and the brachial blood pressures before and after exercise training.Furthermore,the roles of exerciseinduced arterial diameter and blood flow rate in the change of WSS were also determined.The numerical results demonstrate that acute exercise significantly increases the magnitudes of blood flow rate and WSS.Moreover,the vessel elastic deformation is a non-negligible factor in the calculation of the WSS induced by exercise,which generates greater effects on the minimum WSS than the maximum WSS.Additionally,the contributions of exercise-induced variations in blood flow rate and diameter are almost identical in the change of the mean WSS.