考虑颗粒转矩的接触网络诱发各向异性分析

颗粒材料的宏观力学行为与接触网络的组构各向异性密切相关, 根据接触点的滑动与否、转动与否和强弱力情况, 可以将颗粒间的接触系统分为不同的子接触网络. 一般而言, 不同的子接触网络在颗粒体系中的传力机制不同, 对宏观力学响应的贡献也有不同. 采用离散单元法(discrete element method, DEM)模拟了不同抗转动系数$\mu_r$下颗粒材料三轴剪切试验, 分析了剪切过程中不同子接触网络的组构张量的演变规律, 并探究了颗粒抗转动效应对子接触网络各向异性指标演变规律的影响. 研究发现: 剪切过程中转动、非转动接触的组构张量变化不是独立的, 受到颗粒间滑动与否的影响; 非滑动、强接触网络是颗粒间的主要传力结构, 非滑动接触网络的接触法向和法向接触力各向异性均随$\mu_r$的增大而增大, 其对宏观应力的贡献程度随$\mu_r$的增大而减小;强接触网络的接触法向各向异性随$\mu_r$的增大而增大, 但法向接触力各向异性随$\mu_r$的增大无明显变化, 强接触网络对宏观应力的贡献程度在不同$\mu_r$情况下均相同.      

主动脉瓣倾斜角度血流动力学的 PIV 实验研究

瓣叶血栓是主动脉瓣置换术后典型的继发性瓣膜疾病,血流动力学特征异常在其发展过程中至关重要.本文利用粒子图像测速 (particle image velocimetry,PIV) 系统,实验研究了主动脉瓣开口纵向轴线与升主动脉纵向轴线之间倾斜角度 ($\alpha =0^\circ$, $\alpha=5^\circ$,$\alpha =10^\circ$ 和 $\alpha =15^\circ$) 对速度、涡度和黏性剪应力分布等血流动力学特性的影响.研究结果表明:当 $\alpha =0^\circ$ 时,主动脉根部跨瓣血液流动为中心对称流动,而 $\alpha =5^\circ$,$\alpha=10^\circ$ 和 $\alpha =15^\circ$ 时跨瓣血液流动向升主动脉的左冠状动脉一侧倾斜.随着倾斜角度增大,跨瓣血液流动方向倾斜程度增加,血液流动冲击升主动脉壁,损伤内皮细胞导致血栓形成.主动脉瓣倾斜时主动脉窦血液流动速度增大,涡旋也更向主动脉窦底部运动,不利于血液从冠状动脉口流出向心肌供血.同时,主动脉根部的高涡度和高黏性剪应力区域也向升主动脉的左冠状动脉一侧倾斜,主动脉窦的高涡度区域位于主动脉窦底部、高黏性剪应力区域分布于主动脉窦壁面处.主动脉瓣存在倾斜角度时,涡度和黏性剪应力较大,特别是 $\alpha =10^\circ$ 和 $\alpha=15^\circ$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.      

不同壁面取向下超疏水平面直轨道上的气泡滑移

利用特定几何分布的超疏水表面实现气泡定向输运在矿物浮选和生物孵化等领域具有广阔的应用前景, 对平面直线超疏水轨道而言, 其壁面取向是相关工程结构的关键参数, 但超疏水壁面取向对倾斜壁面气泡滑移的影响尚不明确. 本文采用高速阴影成像系统研究了不同壁面取向($-90^\circ\leqslant \beta \leqslant 90^\circ$)及轨道倾角($45^\circ\leqslant \alpha \leqslant 75^\circ$)下, 气泡($D_{eq}=2.4$ mm, $Re=500$ $\sim$ 700, $We=7$ $\sim$ 13)在轨道宽度为2 mm的超疏水直线轨道上的运动特性. 气泡在轨道上的滑移近似为匀速, 形状为具有多脊的半子弹型. 根据气液界面波动程度的不同, 滑移气泡可分为波动型和稳定型, 稳定型气泡只在较小倾角且较大方位角时出现($45^\circ\leqslant \alpha < 70^\circ$, $| \beta | \geqslant 45^\circ$). 根据倾角不同, 滑移速度关于$\beta $有2种变化规律: 当$\alpha \leqslant 65^\circ$, 气泡滑移速度近似为关于$\beta =0^\circ$ 的单峰分布($\beta =0^\circ$时, 气泡滑移速度最大); 当$\alpha \geqslant 70^\circ$, 气泡滑移速度在不同的方位角下基本保持稳定. 气泡的最大滑移速度可达0.66 m/s ($\beta =0^\circ$, $\alpha =70^\circ$), 远大于相同尺度的自由上升气泡($\approx0.25$ m/s), 这主要是壁面浸润性分布和惯性力的耦合效应所致. 轨道取向(方位角$\beta )$及轨道倾角($\alpha )$通过改变气泡沿轨道方向的驱动力和气泡迎风面积影响气泡的滑移速度和气液界面稳定性.      

PIV EXPERIMENTAL STUDY ON THE HEMODYNAMICS OF AORTIC VALVE UNDER VARIED TILTED ANGLES 1)

瓣叶血栓是主动脉瓣置换术后典型的继发性瓣膜疾病,血流动力学特征异常在其发展过程中至关重要.本文利用粒子图像测速 (particle image velocimetry,PIV) 系统,实验研究了主动脉瓣开口纵向轴线与升主动脉纵向轴线之间倾斜角度 ($\alpha =0^\circ$, $\alpha=5^\circ$,$\alpha =10^\circ$ 和 $\alpha =15^\circ$) 对速度、涡度和黏性剪应力分布等血流动力学特性的影响.研究结果表明：当 $\alpha =0^\circ$ 时,主动脉根部跨瓣血液流动为中心对称流动,而 $\alpha =5^\circ$,$\alpha=10^\circ$ 和 $\alpha =15^\circ$ 时跨瓣血液流动向升主动脉的左冠状动脉一侧倾斜.随着倾斜角度增大,跨瓣血液流动方向倾斜程度增加,血液流动冲击升主动脉壁,损伤内皮细胞导致血栓形成.主动脉瓣倾斜时主动脉窦血液流动速度增大,涡旋也更向主动脉窦底部运动,不利于血液从冠状动脉口流出向心肌供血.同时,主动脉根部的高涡度和高黏性剪应力区域也向升主动脉的左冠状动脉一侧倾斜,主动脉窦的高涡度区域位于主动脉窦底部、高黏性剪应力区域分布于主动脉窦壁面处.主动脉瓣存在倾斜角度时,涡度和黏性剪应力较大,特别是 $\alpha =10^\circ$ 和 $\alpha=15^\circ$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.     

双轴压缩下颗粒物质接触力与力链特性研究

为探究双轴压缩过程中颗粒物质接触力与力链特性,采用离散元法研究双轴压缩过程,得到双轴压缩过程中颗粒体系宏观力学特性,即偏应力、体积应变变化规律,并通过对接触力大小分布规律、接触力角度分布规律等细观层次分析,完成接触力大小、力链方向性量化描述,同时探究围压对颗粒物质力学特性的影响。研究表明,双轴压缩过程经历应力强化、应力软化、应力波动三个阶段。在各阶段,强接触与总体接触的各类型接触力大小分布规律可采用高斯函数或指数函数拟合,弱接触的法向与总接触力分布较为均匀;直径较大颗粒上具有接触力集中现象,集中程度随着围压升高而加强;力链网络形态经历"鱼鳞"状-"柱"状-"拱"状-"龟壳"状演变,由初始时基本各向同性,逐渐向最大主应力y向偏转,表现出强烈各向异性,最后又经历小幅度偏转偏离y向。不同围压下力链分量比变化规律证实了该发现,接触力角度分布规律也从力链承载方向验证了此结论。研究结果将为探究双轴压缩中颗粒物质力学特性提供理论基础。     

求解中短半径弯螺杆钻具的纵横弯曲法

介绍了一种求解中短半径($K=1^\circ$/m$ \sim 3^\circ$/m)弯螺杆 钻具的纵横弯曲法, 以指导弯螺杆钻具的设计和使用. 采用变刚度梁柱理论, 并应用等效载荷法对具有初始结构弯角的梁柱进行等效处理, 根据钻具组合支座连续条件及 边界条件建立了中短半径水平井造斜螺杆钻具组合的力学模型, 进而求得钻头侧向力, 依此 预测弯螺杆钻具的造斜能力. 这种分析方法得到了成功应用.     

点载荷作用下密集颗粒物质的传力特性分析

利用颗粒离散元商业软件PFC3D, 模拟了在2m*1m*0.01m容器中直径分别为0.01m, 0.008m和0.006m的颗粒各1*10$^4$个, 受重力作用下的静态密集堆积; 以此为初始条件, 在表层随机选择7个颗粒分别施加5.2*10$^{ - 2}$N(100倍最大颗粒重量)的点载荷, 进行应力传播特点研究. 结果表明: 力的传递在局部范围内呈现很强的各向异性; 应力涨落随着距离的增加呈指数下降; 在大于5倍最大颗粒粒径时, 其分布可以使用弹性力学理论来计算. 探讨了摩擦系数$\mu =0$, 0.2, 1对应力传递的影响, 随着摩擦系数的增加, 各向异性范围减小.     

基于渐进均匀化的平纹编织复合材料低速冲击多尺度方法

针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明:宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.      

高温后高强高性能混凝土双轴压力学性能

利用大型静动真三轴试验机,进行了常温20${^\circ}$C以及200${^\circ}$C$\sim $ 600${^\circ}$C\,6个温度等级高温后高强高性能混凝土在7种应力比双轴压应力状态下的强度与变形试验.测得了双轴主压方向的静态强度、峰值应变与应力应变曲线,剖析了温度和应力比对单、双轴压强度与峰值应变发展趋势的影响规律性以及试件破坏形态. 试验结果表明:随温度的升高,高强高性能混凝土的单轴压减摩强度并不一定降低;双轴压强度相对于单轴压强度的提高倍数取决于应力比、不同温度等级后的高强高性能混凝土``脆硬性'. 提出了带有温度和应力比参数的Kupfer-Gerstle破坏准则公式.      

固体颗粒对半球扰动尾迹影响的实验研究

采用粒子图像测速技术(particle image velocimetry, PIV)研究固体颗粒对放置在平板层流边界层中半球粗糙元尾迹的影响. 实验采集了清水和加入粒径为140 $\mu$m, 220 $\mu$m, 350$\mu$m聚苯乙烯固体颗粒4种工况下二维速度场信息, 基于半球半径的雷诺数为994 ($Re_R=RU/\upsilon$), 固体颗粒的体积浓度为$3.0\times10^{-5}$. 对比清水和两相工况下的平均速度剖面、湍流强度等宏观统计量, 分析固体颗粒对半球尾迹流动宏观特性的影响. 分别利用沿流向不同位置的流向脉动速度的二维空间相关系数和法向脉动速度的功率谱密度函数分析颗粒对尾迹结构演化过程及尾迹结构脱落频率的影响. 结果发现: 与清水相比, 回流区随颗粒粒径增大而逐渐增大; 颗粒使湍流强度增大, 回流区的存在导致在半球后流向位置$2R$前后区域湍流强度呈现不同变化趋势; 颗粒使尾迹结构的流向尺度减小并且随着颗粒粒径的增大先减小后增大; 在尾迹结构运动过程中颗粒的存在促进了尾迹结构的周期性加速和减速运动, 促进作用随着颗粒粒径的增大先增强后减弱; 颗粒的存在促进了尾迹结构的脱落, 脱落频率随颗粒粒径的增大先增大后减小.      

Active control of flow past an elliptic cylinder using an artificial neural network trained by deep reinforcement learning

The active control of flow past an elliptical cylinder using the deep reinforcement learning (DRL) method is conducted. The axis ratio of the elliptical cylinder $\Gamma$ varies from 1.2 to 2.0, and four angles of attack $\alpha=0^\circ, 15^\circ, 30^\circ$, and $45^\circ$ are taken into consideration for a fixed Reynolds number $Re=100$. The mass flow rates of two synthetic jets imposed on different positions of the cylinder $\theta_1$ and $\theta_2$ are trained to control the flow. The optimal jet placement that achieves the highest drag reduction is determined for each case. For a low axis ratio ellipse, i.e., $\Gamma=1.2$, the controlled results at $\alpha=0^\circ$ are similar to those for a circular cylinder with control jets applied at $\theta_1=90^\circ$ and $\theta_2=270^\circ$. It is found that either applying the jets asymmetrically or increasing the angle of attack can achieve a higher drag reduction rate, which, however, is accompanied by increased fluctuation. The control jets elongate the vortex shedding, and reduce the pressure drop. Meanwhile, the flow topology is modified at a high angle of attack. For an ellipse with a relatively higher axis ratio, i.e., $\Gamma\ge1.6$, the drag reduction is achieved for all the angles of attack studied. The larger the angle of attack is, the higher the drag reduction ratio is. The increased fluctuation in the drag coefficient under control is encountered, regardless of the position of the control jets. The control jets modify the flow topology by inducing an external vortex near the wall, causing the drag reduction. The results suggest that the DRL can learn an active control strategy for the present configuration.     

Discrete element modeling of inherently anisotropic granular assemblies with polygonal particles

In the present article, we study the effect of inherent anisotropy, i.e., initial bedding angle of particles and associated voids on macroscopic mechanical behavior of granular materials, by numerical simulation of several biaxial compression tests using the discrete element method (DEM). Particle shape is considered to be irregular convex-polygonal. The effect of inherent anisotropy is investigated by following the evolution of mobilized shear strength and volume change during loading. As experimental tests have already shown, numerical simulations also indicate that initial anisotropic condition has a great influence on the strength and deformational behavior of granular assemblies. Comparison of simulations with tests using oval particles, shows that angularity influences both the mobilized shear strength and the volume change regime, which originates from the interlocking resistance between particles.     

DEM investigation on the effect of sample preparation on the shear behavior of granular soil

The effect of initial fabric anisotropy produced by sample preparation on the shear behavior of granular soil is investigated by performing discrete element method (DEM) simulations of fourteen biaxial tests in drained conditions. Numerical test specimens are prepared by three means: gravitational deposition, multi-layer compression, and isotropic compression, such that different initial inherent soil fabrics are created. The DEM simulation results show that initial fabric anisotropy exerts a considerable effect on the shear behavior of granular soil, and that the peak stress ratio and peak dilatancy increase with an increase in the fabric index a_n that is estimated from the contact orientations. The stress–dilatancy relationship is found to be independent of the initial fabric anisotropy. The anisotropy related to the contact orientation and contact normal force accounts for the main contribution to the mobilized friction angle. Also, the occurrence of contractive shear response in an initial shearing stage is accompanied by the most intense particle rearrangement and microstructural reorganization, regardless of the sample preparation method. Furthermore, the uniqueness of the critical state line in e–log p′ and q–p′ plots is observed, suggesting that the influence of initial fabric anisotropy is erased at large shear strains.     

超音速气流粉碎喷嘴数值模拟

对气流粉碎装置的喷嘴结构和参数进行设计与优化,采用流体动力学软件对所设计喷嘴进行流场模拟,对所喷嘴效果进行检验。分别讨论了锥顶角和内腔造型对超音速喷嘴性能的影响,通过结果比较得出,入口压力3.5MPa、入口直径为6mm的喷嘴为设计的最佳喷嘴。超音速喷嘴锥顶角在8度到12度之间变化时,对喷嘴的性能影响不大,内腔造型为光滑曲面的超音速喷嘴性能最佳。     

椭球颗粒体系剪切过程中自由体积的分布与演化

颗粒材料是一种复杂的多体相互作用体系, 由大量离散的颗粒和其周围的自由体积组成. 虽然颗粒的自由体积与颗粒材料的力学性能和变形特征的相关性已得到证实, 但是由于表征上的困难, 目前对非球形颗粒体系的局部自由体积的认识还不够充分. 本文采用连续离散耦合分析方法进行了不同主轴长度的椭球颗粒试样的三轴剪切数值模拟, 基于Set Voronoi算法对剪切过程中的颗粒试样进行了Voronoi元胞分割, 分析了颗粒试验在剪切过程中自由体积的统计分布特性和演化规律, 研究了颗粒形态对自由体积的影响. 剪切过程中Voronoi元胞的各向异性逐渐增强, 且各项异性增强程度随颗粒非球度的增加而增大, 表明非球颗粒在剪切过程中经历更加强烈的重排列. 具有不同非球度的椭球颗粒体系的局部孔隙比均服从k?Γ分布, 且这个分布仅与颗粒体系的全局孔隙比相关, 不受颗粒形态和剪切状态的影响. 局部孔隙比的波动呈现非对称拉普拉斯分布, 非对称参数刻画了局部自由体积收缩和膨胀的博弈, 其与全局孔隙比呈线性关系.      

横摆振荡翼型动态气动掠效应试验研究

大型风力机设计对获取翼型更加全面、准确的动态载荷提出更高要求, 研究翼型横摆振荡动态气动特性具有重要意义. 借助"电子凸轮"技术和动态数据同步采集手段, 针对翼型动态“掠效应”首次开展了横摆振荡风洞试验研究, 研究表明: 横摆振荡翼型的气动曲线存在明显迟滞效应, 吸力面压力周期性波动是主要诱因, 且随着振荡频率、初始迎角和振幅的增大, 气动迟滞特性均增强; 升力和压差阻力随横摆角变化的迟滞回线呈"W"形, 俯仰力矩迟滞回线呈"M"形, 升力差量迟滞回线呈"$\infty$"形; 负行程下翼型气动力相对于正行程下的更高, 且负行程下翼型气动力随振荡频率的增大而略有增大, 正行程下则明显减小; 升力系数功率谱密度分布在振荡频率倍频处的能量集中的幅值随着振荡频率增大有增大趋势; 吸力面1.2%和40%弦长处压力的滞回特性较强, 是由于翼面剪切层涡和动态分离涡周期性发展、运动、破裂和重建; 振幅为$10^{\circ}$时, 升力迟滞曲线呈"$^{\wedge}$"形, 振幅为$30^{\circ}$ 时, 升力迟滞曲线呈"$^{\wedge\wedge\wedge}$"形.