滑轨导向式静/动态双轴拉伸实验技术

基于液压伺服高速加载系统,发展了一种材料双轴拉伸力学性能测试技术。利用锥面接触导向驱动方法,把加载锤竖直方向的驱动力转化为水平方向的双轴驱动力,从而实现对十字形试样平面双轴加载。借助有限元数值模拟手段优化了锥面接触角和十字形试样尺寸。当接触锥角为45°时,既有较好的水平驱动转化效率,同时又保持较小的接触力,确保水平驱动加载各组件在弹性变形范围内,可多次重复使用。确定了加载臂狭缝个数、狭缝与减薄区边缘长度和标距段厚度等试样设计关键参数,在十字形试样测试标距段内实现了均匀变形。设计了测力夹持一体化导杆和非接触光学全场应变测试系统,准确获得了试样的应力和应变。利用此平面双轴拉伸加载装置,开展2024-T351铝合金板单轴拉伸实验和激光探测同步性验证实验,验证装置设计的可行性;开展铝合金板材在不同加载速率下的双轴拉伸实验,得到在双轴加载下铝合金板材应力应变曲线,并与单轴加载下实验结果进行了对比分析。     

爆轰加载下银气溶胶源项实验研究

为研究钚气溶胶生成规律及源项分布,提出了一种密封环境下的金属气溶胶源项实验方法,利用爆炸容器及气溶胶采样器等实验设备,开展了钚替代材料银的气溶胶源项生成实验,分析了空气动力学直径小于10 μm的银气溶胶源项分布特性,并与外场扩散试验中钚气溶胶数据进行了对比。进一步采用不同炸药及加载装置,对相同银片样品进行加载实验,研究了加载峰值压力对银气溶胶源项分布的影响。研究结果表明,金属银可作为替代材料研究爆轰加载条件下钚气溶胶的源项参数,在一定实验条件下,二者归一化积累质量分布具有较高的一致性;爆轰加载产生的银气溶胶中,较小粒径气溶胶质量含量较高。气溶胶总量与加载峰值压力呈二次函数分布,当峰值压力超过某个临界值时,气溶胶总量将达到极值;此外,通过对同发次实验的3次不同采样数据进行分析,发现了较明显的气溶胶凝并和沉降现象。     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象, 该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径. 为探究颗粒毛细效应的影响因素, 采用离散元方法, 模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性, 考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律. 结果表明, 在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下, 管径与粒径比$D/d = 3.33$时, 管内颗粒堵塞严重, 使得颗粒上升缓慢,并造成颗粒柱中断; $D/d = 8.33$时, 起初毛细上升高度增加迅速, 随后毛细上升高度的增大逐渐减缓, 管内颗粒在管径方向几乎不存在速度梯度; $D/d =15$时, 随着颗粒毛细上升高度的增大, 管内颗粒柱分离为速度截然不同的两层, 上层颗粒在管径方向几乎不存在速度梯度, 而下层颗粒存在明显的速度梯度.研究还发现, 在毛细效应能够发生的管径范围内, 存在一个对应于颗粒最终毛细上升高度最大值的临界管径, 当管径小于临界管径时, 颗粒最终毛细上升高度随管径的增大而增大, 当管径大于临界管径时, 颗粒最终毛细上升高度随管径的增大而趋于减小; 增大容器宽度,临界管径有所增大; 增大振幅、适当提高频率能够有效促进临界管径的增大.      

爆炸加载下金属柱壳破片软回收技术研究

针对爆炸加载下金属柱壳膨胀断裂破片软回收的研究需求,本文通过理论分析和初步的数值模拟设计了由低密度聚氨酯泡沫与水介质为主体的回收装置。与传统单一材料为主的回收装置相比,该回收装置既能在破片高速阶段将低阻抗聚氨酯泡沫对破片的冲击压力减小到约为水对破片冲击压力的1/3,又使破片速度全程持续地较大幅度衰减,还能在破片低速阶段又能充分利用水介质密度大的优势,减小以聚氨酯泡沫单一材料为主的回收装置尺寸。依托该装置开展了炸药加载下304不锈钢柱壳膨胀断裂回收实验。通过测量回收池外壁速度、检查实验后的回收池外观,发现回收池池壁和底部完好,可以重复使用;通过对回收破片称重统计,破片回收率超过85%,破片内外界面辨识度高,破片表面车刀纹清晰可见,内部可见多条未贯穿的裂纹。表明该回收装置对破片的冲击损伤显著降低。根据破片断口和表面信息,推测了破片在金属柱壳的大致位置。本文最后初步给出了回收破片的平均厚度及质量分布等相关信息的统计结果。     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象,该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径.为探究颗粒毛细效应的影响因素,采用离散元方法,模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性,考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律.结果表明,在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下,管径与粒径比D/d=3.33时,管内颗粒堵塞严重,使得颗粒上升缓慢,并造成颗粒柱中断; D/d=8.33时,起初毛细上升高度增加迅速,随后毛细上升高度的增大逐渐减缓,管内颗粒在管径方向几乎不存在速度梯度; D/d=15时,随着颗粒毛细上升高度的增大,管内颗粒柱分离为速度截然不同的两层,上层颗粒在管径方向几乎不存在速度梯度,而下层颗粒存在明显的速度梯度.研究还发现,在毛细效应能够发生的管径范围内,存在一个对应于颗粒最终毛细上升高度最大值的临界管径,当管径小于临界管径时,颗粒最终毛细上升高度随管径的增大而增大,当管径大于临界管径时,颗粒最终毛细上升高度随管径的增大而趋于减小;增大容器宽度,临界管径有所增大;增大振幅、适当提高频率能够有效促进临界管径的增大.     

基于近场动力学的单晶冰弹性各向异性数值模拟方法

天然冰材料在变形与破坏行为上的各向异性特征是冰与结构相互作用中产生复杂载荷过程的关键诱因, 而天然冰各向异性的根源则在于单晶冰的各向异性. 目前, 学术界针对单晶冰各向异性的数值模拟方法研究仍较为缺乏. 为了准确再现天然冰材料的特殊力学性质, 本文基于近场动力学理论, 提出了一种单晶冰弹性各向异性的数值模拟方法. 该方法的核心思想是将单晶冰杨氏模量沿不同加载方向的变化规律引入到近场动力学力密度向量的影响函数中. 以前人实验测试得到的杨氏模量值为参考, 通过开展与C轴呈0°, 45°和90°三个加载方向的单晶冰单轴压缩数值模拟实验, 提出了针对该影响函数的修正和辅助参数标定方法, 最终在15°, 30°, 60°和75°等其他四个加载方向进行了验证. 结果表明: 本文提出的针对影响函数的修正与参数标定方法, 能够较为便捷地找到数值模型杨氏模量与参考杨氏模量相一致的影响函数最优解, 即本文提出的基于影响函数的近场动力学数值模拟方法, 能够合理、准确地模拟单晶冰的弹性各向异性行为. 本文研究成果可为后续多晶冰各向异性数值模拟方法的建立提供基础性参考.      

各向异性介质波动问题的时域边界元法及实验验证

对前已建立的各向异性介质波动问题的时域边界元计算模型进行了实验验证．用单向纤维增强光弹性复合材料模拟正交异性介质，用冲击加载，加载方向与纤维方向分别成0°，90°及45°角度，进行了正交异性动态光弹性实验及动态应变测量，并同时对该模型进行了时城边界元计算．将时域边界元方法计算出的应力分量代入正交异性动态光弹性的动态应力-光性定律，得到双折射条纹级数随时间的变化曲线，将其与动态光弹性实验的结果进行比较；此外，由动态电测获得的应变响应曲线推算出应力时程；与时城边界元计算出的应为响应曲线也进行了比较．两种情况下，时域边界元的计算成果均与实验成果吻合较好，从而证明该各向异性介质波动问题的时域边界元计算模型具有较好的精度及稳定性，可用于各向异性介质的动态问题特别是波传播问题的分析研究．     

砂岩球体法向恢复系数实验研究

法向恢复系数是岩崩块石运动分析的关键参数,其取值直接决定了块石的运动轨迹。本文中采用自行设计的碰撞实验装置和声频采样技术,测定了砂岩球体碰撞的法向恢复系数,研究了粒径、碰撞速度、含水状态和板的弹性特性4个因素对恢复系数的影响。结果表明:砂岩球体法向恢复系数存在复杂的尺寸效应,恢复系数随粒径的增大先增大后减小;碰撞过程中存在的黏弹性耗能机理和弹塑性损伤耗能机理共同作用产生了复杂的尺寸效应;受砂岩非均质特性的作用,粒径较小时,恢复系数的速度效应较明显(随速度增大而增大),粒径较大时速度对恢复系数的影响消失;砂岩饱和使黏弹性耗能和弹塑性损伤耗能增加,使恢复系数比风干时低;等效弹性模量对恢复系数的影响较大,等效弹性模量越大,法向恢复系数越小。     

双轴应力状态下正交异性动态光弹性应力——光性定律研究

基于静态下Ｈｙｅｒ和Ｌｉｕ表述的正交异性应力－光性定律，在前文中，提出了正交异性光弹性复合材料的动态应力－光性定律并证明了其在单轴应力状态下的正确性。本文旨在进一步考察在双轴应力状态下正交异性动态应力－光性定律的适用性，采用的方法是对纤维增强光弹性复合材料制作的平板模型施加冲击荷载，加载方向与材料纤维方向分别成０°、９０°及４５°角，同时进行正交异性动态光弹性实验和动态应变测量，另外，对该模型进行相应的各向异性介质时域边界元计算。把动态应变测量推算出的应力分量以及时域边界元计算出的应力分量分别代入正交异性动态应力－光性定律，得到随时间变化的双折射条纹级数历程，将其与正交异性动态光弹性实验的结果进行比较。实验及计算结果表明，在三个加载方向下，由这三种方法得到的双折射条纹级数历程均吻合良好，从而证明了前文提出的正交异性动态应力－光性定律在双轴应力状态下的正确性。     

探头在准二维颗粒介质中慢速运行的受阻实验研究

测量了圆锥探头缓慢压入准二维颗粒介质过程中所受阻力随深度的变化,发现阻力曲线在不同区域呈现不同的增长规律,曲线全过程存在两次增长规律的转变.本实验条件下观测到的现象表明:当填充颗粒总高度较高时,阻力曲线存在两个拐点,前部拐点是由颗粒回填所致,拐点深度约65mm(与探头尺寸有关),且不随颗粒填充高度变化,后部拐点是由于底部边界影响所致;后部拐点出现的深度Z2随着颗粒填充高度Zmax的增加线性增加,容器底部边界开始影响探头受力时探头距离底部的距离H也会随着颗粒填充高度Zmax的增加而增加,Z2-Zmax和H-Zmax关系拟合直线的斜率均约为0.5.     

竖直振动管中颗粒毛细效应的离散元模拟

将一根细管插入填充有颗粒的静止容器中并对管施加竖直振动,颗粒将在管内发生上升运动,并最终稳定在一定高度,这一现象与液体毛细效应类似,被称为颗粒毛细效应.为探究颗粒毛细效应过程中伴随的颗粒尺度动力学行为及机理,基于离散元方法建立颗粒运动模型,对颗粒毛细效应动力学过程和特性开展数值模拟研究.模拟再现了文献中实验得到的颗粒毛细效应全过程,给出了管内颗粒柱高度随时间的演变规律,结果表明,受到颗粒系统参数的影响,本模拟条件下颗粒毛细效应过程呈现单周期上升、倍周期上升和倍周期稳定三个阶段,在倍周期上升阶段颗粒柱上升速度逐渐减小,平缓过渡到稳定阶段.在此基础上,分析了管内颗粒速度场和填充率分布随时间的演变特性,揭示了颗粒毛细效应过程中由容器传输到管内的颗粒的占比分布.研究发现,管内不同高度位置颗粒的运动并不同步,随着管的振动,管内出现速度波,速度波的传播引起管内颗粒出现膨胀和压缩交替的情况,从而管内颗粒填充率随时间发生周期性波动;在上升阶段,越接近管壁由容器传输到管内的颗粒占比越大,在稳定阶段,管内上层颗粒的对流引起容器传输到管内的颗粒占比发生反转.      

Representative elementary volume analysis of polydisperse granular packings using discrete element method

The representative elementary volume (REV) for three-dimensional polydisperse granular packings was determined using discrete element method simulations. Granular mixtures of various sizes and particle size distributions were poured into a cuboid chamber and subjected to uniaxial compression. Findings showed that the minimum REV for porosity was larger compared with the REV for parameters such as coordination number, effective elastic modulus, and pressure ratio. The minimum REV for porosity and other parameters was found to equal 15, 10, and 5 times the average grain diameter, respectively. A study of the influence of sample size on energy dissipation in random packing of spheres has also confirmed that the REV size is about 15 times the average grain diameter. The heterogeneity of systems was found to have no effect on the REV for the parameters of interest for the narrow range of coefficient of uniformity analyzed in this paper. As the REV approach is commonly applied in both experimental and numerical studies, determining minimum REV size for polydisperse granular packings remains a crucial issue.     

The contribution of frictional contacts to the shear strength of coarse glass bead powders and slurries

This paper investigates the use of the shear vane technique as a means of determining frictional and cohesive interparticle force contributions to the shear strength of coarse glass bead powders and slurries. To this end, the shear strength of 203-μm glass beads in air and slurried in water and kaolinite suspensions was determined as a function of vane immersion depth, vane geometry, and container size. Both vane immersion depth and container diameter are found to significantly impact the shear strength measured using the vane technique. An equation describing interparticle frictional and cohesive contributions to shear vane measurements was derived in an effort to describe the experimental results. A Janssen stress distribution model for granular materials forms the basis for this equation and appears to explain the behavior of shear strength measurements at varying immersion depths. The presence of the Janssen stress distribution can affect the interpretation of shear vane results. Rather than shear strength being a material property, as is the case with flocculated colloid slurries and polymer solutions, shear strength becomes a process property where vane depth, container size, and container material can result in significant measurement variations. Such parameters should be considered before using the shear vane results on applications involving granular material components.US Department of Energy Environmental Management Science Program project no. 90162.     

Segregation of particulate solids: Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method (DEM) can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable.An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investigated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behavior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

Segregation of particulate solids＂ Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method （DEM） can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable. An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investi- gated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behav- ior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

Experimental study on the influence of bed height on convection cell formation

In this study, we have investigated experimentally the convection phenomenon of granular matter in a vertical shaker device. The motions of the granular materials were recorded by a high-speed camera. Image processing technology and a particle tracking method were employed to measure the convection strength, convection size, and granular temperature. The results show that the initial bed height, dimensionless vibrating acceleration, and vibrating frequency have significant influence on the convection mechanism. An interesting finding that occurs with increasing the bed height is the existence of two peaks and one low convection strength value which occurs due to the growth of the solid-like region. We also find that the two peak and one low values become greater when higher dimensionless vibrating acceleration is applied to the granular bed.     

High-amplitude elastic solitary wave propagation in 1-D granular chains with preconditioned beads: Experiments and theoretical analysis

Elastic solitary waves resulting from Hertzian contact in one-dimensional (1-D) granular chains have demonstrated promising properties for wave tailoring such as amplitude-dependent wave speed and acoustic band gap zones. However, as load increases, plasticity or other material nonlinearities significantly affect the contact behavior between particles and hence alter the elastic solitary wave formation. This restricts the possible exploitation of solitary wave properties to relatively low load levels (up to a few hundred Newtons). In this work, a method, which we term preconditioning, based on contact pre-yielding is implemented to increase the contact force elastic limit of metallic beads in contact and consequently enhance the ability of 1-D granular chains to sustain high-amplitude elastic solitary waves. Theoretical analyses of single particle deformation and of wave propagation in a 1-D chain under different preconditioning levels are presented, while a complementary experimental setup was developed to demonstrate such behavior in practice. The experimental results show that 1-D granular chains with preconditioned beads can sustain high amplitude (up to several kN peak force) solitary waves. The solitary wave speed is affected by both the wave amplitude and the preconditioning level, while the wave spatial wavelength is still close to 5 times the preconditioned bead size. Comparison between the theoretical and experimental results shows that the current theory can capture the effect of preconditioning level on the solitary wave speed.     

Nonlinear vibration analysis of laminated composite Mindlin micro/nano-plates resting on orthotropic Pasternak medium using DQM

The nonlocal nonlinear vibration analysis of embedded laminated microplates resting on an elastic matrix as an orthotropic Pasternak medium is investigated. The small size effects of micro/nano-plate are considered based on the Eringen nonlocal theory. Based on the orthotropic Mindlin plate theory along with the von Kármán geometric nonlinearity and Hamilton's principle, the governing equations are derived. The differential quadrature method (DQM) is applied for obtaining the nonlinear frequency of system. The effects of different parameters such as nonlocal parameters, elastic media, aspect ratios, and boundary conditions are considered on the nonlinear vibration of the micro-plate. Results show that considering elastic medium increases the nonlinear frequency of system. Furthermore, the effect of boundary conditions becomes lower at higher nonlocal parameters.     

Load-cell-design developments by numerical and experimental methods

Design data for a family of ring-shaped elastic elements are derived using analytical, numerical and experimental methods. The development of a family of elastic elements for load cells operating at medium to low load level is considered in detail, using several approaches. Snapes analyzed were circular ring with integral bosses and ‘square’ ring. The effects of thickness, boss size and fillet radius (reduced to nondimensional form) on strain/load and defiection/load relationships are evaluated within a given sample space, considering a two-dimensional problem. Thin-ring and curved-beam theory support the analytical approach, with correction terms covering some departures from theoretical shape. Numerical analysis with the finite-element method was mainly used to evaluate the effects of boss size and shape. Accordingly, the required grid patterns were developed using a mesh-generating method capable of covering the range of parameters considered. A set of configurations was selected, enabling single and combined effects to be analyzed with little entanglement between estimates. For the square ring, frame theory and a limit solution for the case of the plate with central hole are resorted to. Photoelastic tests were performed, covering a range of shapes by sequential remachining of the model, and equal contigurations were analyzed with the finite-element method. Strain-gage tests on actual elastic elements were made, considering also the results of a series of tests carried out at the Metrology Institute some time ago. Results obtained yield comprehensive information on the effect of several parameters on the load-output transfer function of ring-shaped elastic elements.