Variational formulation of a rod under torsional vibration for crack identification

A variational formulation is developed for the torsional vibration of a cylindrical shaft with a circumferential crack. The work is compared with existing methods. The Hu–Washizu–Barr variational formulation was used to develop the differential equation and the boundary conditions of the cracked rod. The general variational principle and the independent assumptions about the displacement, the momentum, the strain and the stress fields of the cracked rod, and the equations of motion for a uniform rod in torsional vibration, are derived. The crack was modelled as a continuous flexibility using the displacement field in the vicinity of the crack, found with fracture mechanics methods. Rayleigh quotient was used to approximate the natural frequencies of the cracked rod. Independent evaluations of crack identification methods in rotating shafts are reported and compared with methods using the continuous crack flexibility theory.     

大理岩的高应变率动态劈裂实验

首次把平台巴西圆盘试样引入动态劈裂试验。利用直径100mm的分离式Hopkinson杆对大理岩巴西圆盘和平台巴西圆盘试样进行了动态劈裂实验。结合有限元分析,得到了大理岩的动态劈裂破坏的拉伸强度。分析了巴西圆盘和平台巴西圆盘的典型破坏方式。结果表明,大理岩的动态拉伸强度随着应变率的提高而增加。利用圆盘中心粘贴的应变片来测大理岩等脆性材料的动态拉伸强度,是一种简便高效的试验方法。和巴西圆盘相比,平台巴西圆盘具有更大的优越性和更好的测量效果。     

Deduction of Spinal Loading from Vertebral Body Surface Strain Measurements

The lumbar intervertebral disc, the apparent nexus of low back pain, undergoes biomechanical changes during its degeneration which are as yet poorly understood. In an effort to ultimately examine in vivo daily activity loads across intervertebral discs, we engaged in the following methodological study. The aim of this research was to correlate vertebral body surface strains with the loads across a lumbar spine segment. Rosette strain gages were affixed anterolaterally on L4 and L5 in a macaque monkey model. These tissues were loaded axially and with sagittal plane moments and the principal strains were compared with the applied loads. Predictable axial and sagittal plane loading profiles were found for similar strain measurements and the system was found to be robust through freezing and thawing. These results support future research aimed at quantifying the in vivo disc mechanics of healthy and degenerate tissues in an attempt to develop prevention or intervention strategies to ease those afflicted with low back pain.     

Bifurcation and chaotic response of a cracked rotor system with viscoelastic supports

Cracks appearing in the shaft of a rotary system are one of the main causes of accidents for large rotary machine systems. This research focuses on investigating the bifurcation and chaotic behavior of a rotating system with considerations of various crack depth and rotating speed of the system’s shaft. An equivalent linear-spring model is utilized to describe the cracks on the shaft. The breathing of the cracks due to the rotation of the shaft is represented with a series truncated time-varying cosine series. The geometric nonlinearity of the shaft, the masses of the shaft and a disc mounted on the shaft, and the viscoelasticity of the supports are taken into account in modeling the nonlinear dynamic rotor system. Numerical simulations are performed to study the bifurcation and chaos of the system. Effects of the shaft’s rotational speed, various crack depths and viscosity coefficients on the nonlinear dynamic properties of the system are investigated in detail. The system shows the existence of rich bifurcation and chaos characteristics with various system parameters. The results of this research may provide guidance for rotary machine design, machining on rotary machines, and monitoring or diagnosing of rotor system cracks.     

Rock deformation equations and application to the study on slantingly installed disc cutter

At present the mechanical model of the interac- tion between a disc cutter and rock mainly concerns indentation experiment, linear cutting experiment and tunnel boring machine （TBM） on-site data. This is not in line with the actual rock-breaking movement of the disc cutter and impedes to some extent the research on the rock-breaking mechanism, wear mechanism and design theory. Therefore, our study focuses on the interaction between the slantingly installed disc cutter and rock, developing a model in accordance with the actual rock-breaking movement. Displacement equations are established through an analysis of the velocity vector at the rock-breaking point of the disc cutter blade; the func- tional relationship between the displacement parameters at the rock-breaking point and its rectangular coordinates is established through an analysis of micro-displacement vectors at the rock-breaking point, thus leading to the geometric equations of rock deformation caused by the slantingly installed disc cutter. Considering the basically linear relationship between the cutting force of disc cutters and the rock deformation before and after the leap break of rock, we express the constitutive relations of rock deformation as generalized Hooke＇s law and analyze the effect of the slanting installa- tion angle of disc cutters on the rock-breaking force. This will, as we hope, make groundbreaking contributions to the development of the design theory and installation practice of TBM.     

PBX炸药的拉伸断裂损伤分析

采用巴西实验作为间接拉伸加载手段,研究了某PBX(polymer-bonded explosive)炸药试样拉伸作用下的断裂损伤特性.发展了PBX炸药的光学制备技术,获得了试样在光学显微镜下的细观表面形貌和断裂损伤形貌,结合高速摄影和数字相关分析技术获得了试样的形变和破坏过程.实验结果表明该PBX炸药的静态拉伸强度低于...     

岩石断裂韧度试样CCNBD临界应力强度因子的全新数值标定

国际岩石力学学会（ISRM）提出的用人字型切槽巴西圆盘（Cracked Chevron Notched Brazilian Disc—CCNBD）测试岩石I型断裂韧度所必需的量纲为一的临界应力强度因子Y＊min的精度需要进一步改进。本文对CCNBD试样的Y＊min进行了系统的重新标定,CCNBD的几何参数有效范围为0.44≤αB≤1.04、0≤α0≤0.69、0.4≤α1≤0.8。采用有限单元法对CCNBD做三维数值分析,得到了435种CCNBD试样的Y＊min值标定结果;在此基础上全面修正了CCNBD试样Y＊min计算公式中u和v的取值表;给出了不用查表直接确定CCNBD试样Y＊min值的近似表达式,该表达式计算结果与标定结果的相对误差绝对值在1.87%以内。对于ISRM建议的CCNBD标准试样的Y＊min值：ISRM标定值0.84比本文结果0.957小12.2%;分片合成法标定值0.947比本文结果小1.0%;子模型法标定值0.943比本文结果小1.5%。本文特别强调了任何CCNBD试样Y＊min的取值必须在它对应的上限和下限所限定的范围内,这一要求对判断Y＊min标定值是否合理是很重要的。     

Stress intensity factors for the Brazilian disc with a short central crack: Opening versus closing cracks

Closed form expressions are obtained for the stress intensity factors (SIFs) in case of a Brazilian disc with a short central crack, the length of which does not exceed one fifth of the disc radius. The disc is loaded by uniform radial pressure along two finite symmetric arcs of its periphery. The solution is achieved using the method of complex potentials introduced by Kolosov and Muskhelishvili. The advantage of the expressions obtained is that they are valid both for cracks under opening mode as well as for closing cracks. For the first case (opening cracks) the results of the present study are compared with existing approximate solutions and it is concluded that the agreement is excellent as long as the length of the crack remains relatively small compared to the radius of the disc. Regarding the case of a closing mode crack the procedure proposed here (based on a recent alternative approach of the cracked Brazilian disc) leads to a physically acceptable deformed crack shape instead to an unnatural crack with overlapped lips. At the same moment the dependence of the SIFs on the properties of the material is eliminated.     

Three dimensional analytical solution for finite circular cylinders subjected to indirect tensile test

This paper derives a new three-dimensional (3-D) analytical solution for the indirect tensile tests standardized by ISRM (International Society for Rock Mechanics) for testing rocks, and by ASTM (American Society for Testing and Materials) for testing concretes. The present solution for solid circular cylinders of finite length can be considered as a 3-D counterpart of the classical two dimensional (2-D) solutions by Hertz in 1883 and by Hondros in 1959. The contacts between the two steel diametral loading platens and the curved surfaces of a cylindrical specimen of length H and diameter D are modeled as circular-to-circular Hertz contact and straight-to-circular Hertz contact for ISRM and ASTM standards respectively. The equilibrium equations of the linear elastic circular cylinder of finite length are first uncoupled by using displacement functions, which are then expressed in infinite series of some combinations of Bessel functions, hyperbolic functions, and trigonometric functions. The applied tractions are expanded in Fourier–Bessel series and boundary conditions are used to yield a system of simultaneous equations. For typical rock cylinders of 54 mm diameter subjected to ISRM indirect tensile tests, the contact width is in the order of 2 mm (or a contact angle of 4°) whereas for typical asphalt cylinders of 101.6 mm diameter subjected to ASTM indirect tensile tests the contact width is about 10 mm (or a contact angle of 12°). For such contact conditions, 50 terms in both Fourier and Fourier–Bessel series expansions are found sufficient in yielding converged solutions. The maximum hoop stress is always observed within the central portion on a circular section close to the flat end surfaces. The difference in the maximum hoop stress between the 2-D Hondros solution and the present 3-D solution increases with the aspect ratio H/D as well as Poisson’s ratio ν. When contact friction is neglected, the effect of loading platen stiffness on tensile stress in cylinders is found negligible. For the aspect ratio of H/D = 0.5 recommended by ISRM and ASTM, the error in tensile strength may be up to 15% for both typical rocks and asphalts, whereas for longer cylinders with H/D up to 2 the error ranges from 15% for highly compressible materials, and to 60% for nearly incompressible materials. The difference in compressive radial stress between the 2-D Hertz solution or 2-D Hondros solution and the present 3-D solution also increases with Poisson’s ratio and aspect ratio H/D. In summary, the 2-D solution, in general, underestimates the maximum tensile stress and cannot predict the location of the maximum hoop stress which typically locates close to the end surfaces of the cylinder.     

On numerical modelling of embedded subsonic flow

An analytical model has been developed for computing embedded subsonic flow in rocket plumes from underexpanded axisymmetric supersonic nozzles. Numerical procedures based on the analysis have been incorporated in a simplified, non-reacting exhaust structure program and calculations for representative plume conditions performed. The technique is numerically stable and has provided satisfactory predictions of Mach–disc associated embedded subsonic flow.