THE PERFORATION PARAMETER OPTIMIZATION FOR MULTI FRACTURE INITIATION IN HYDRAULIC FRACTURING

针对复杂地质近井筒多裂缝起裂的难题, 结合岩体力学与损伤力学相关理论, 建立了近井筒地层流固——损伤力学模型, 采用数值求解获得了射孔壁面的应力及损伤状态, 计算表明:初始地应力最小值越小, 地层的起裂压力越低;射孔方位角越小, 起裂压力越低. 随着射孔方位角增大, 近井筒微裂缝为45° 夹角的拐折裂缝, 起裂压力增大, 基于此, 对于正断层及顺滑断层地应力类型, 取射孔方位角小于30°, 螺旋射孔夹角30°或45°, 对于逆断层地应力类型, 螺旋射孔夹角取45°, 可有效地实现多裂缝起裂.     

含缺陷钢丝多道次拉拔模拟及工艺参数优化

利用有限元软件ABAQUS建立了含缺陷热轧盘条多道次拉拔的分析模型, 研究了道次压缩率、模具顶角和裂纹长度对拉丝安全性的影响,以减少断丝率为目标探讨 道次压缩率和模具顶角的最佳配合方式. 研究结果表明：模具顶角相同时,COD值和$J$积分 值均随着裂纹尺寸的增加而增大;裂纹长度一定时,最佳入线直径随着模具顶角的增加而增 大;最佳模具顶角应在4°- 8°之间. 利用有限元方法和断裂力学理论模拟含缺陷的热轧盘条经过多道次拉丝成型的过程,建立了钢丝经过多道次拉拔模具发生拉拔变形的分析模型,用模拟方法重点研究了道次压缩率、模具顶角以及裂纹长度对拉丝安全性的影响,以减少断丝率为目标探讨道次压缩率和模具顶角的最佳配合方式;研究结果表明：模具顶角相同时,随着裂纹尺寸的增加,J积分值随之增大;裂纹长度一定时,随着模具顶角的增加,入线直径的条件最佳值也在增大;最佳模具顶角的范围应在4°－8°之间。在拉拔前应对线材进行无损探伤从而控制裂纹长度,提高线材使用率和成品质量。研究结果对进一步优化拉丝工艺参数和改进成品PC钢绞线质量提供了有价值的理论指导。     

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

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

主动脉瓣倾斜角度血流动力学的 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$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.      

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$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.     

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.     

一种求解抛物化Navier-Stokes方程的空间推进算法

讨论了抛物化NS方程(parabolized Navier-Stokesequations, PNS)的数学性质,对比分析多种处理流向压力梯度的方法的优缺点. 以此为基础,成功地将LU-SGS隐式时间积分方法推广到PNS方程的流向空间积分上,发展了基于PNS方程的有限体积单次扫描空间推进算法(single-sweep parabolized Navier-Stokesalgorithm, SSPNS). 在该算法中,横向无黏数值通量和黏性通量分别采用混合型迎风格式和中心格式求解. 用SSPNS算法计算了4个典型流场,包括超声速平板流、15$^\circ$楔板压缩高超声速流、带攻角的高超声速锥形流和侧压式高超声速进气道流动. SSPNS计算结果与NASA UPS程序数值结果、文献提供的实验数据及理论分析结果符合得很好.对比研究表明,SSPNS 法与传统时间迭代法相比,二者计算精度相当,而SSPNS计算速度快1~2个量级,存储量至少低1个量级.关键词 抛物化NS方程;空间推进算法;LU-SGS隐式积分方法;超声速/高超声速流动      

一种新的多轴高周疲劳寿命预测模型

分析和讨论3种典型载荷(单轴拉压、纯扭及90^\circ非比例)情况下的5组损伤控制参数,提出了一种以临界面上最大剪切应力幅和最大法向应力的非线性组合作为损伤控制参数的多轴高周疲劳寿命预测模型, 该模型考虑了平均应力对疲劳寿命的影响, 比现有的疲劳预测模型具有更宽的金属材料适用范围. 两种不同类型材料下的多轴非比例试验的预测结果表明,模型的预测结果与试验符合较好.      

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

利用特定几何分布的超疏水表面实现气泡定向输运在矿物浮选和生物孵化等领域具有广阔的应用前景, 对平面直线超疏水轨道而言, 其壁面取向是相关工程结构的关键参数, 但超疏水壁面取向对倾斜壁面气泡滑移的影响尚不明确. 本文采用高速阴影成像系统研究了不同壁面取向($-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 )$通过改变气泡沿轨道方向的驱动力和气泡迎风面积影响气泡的滑移速度和气液界面稳定性.      

锥形涡诱导下建筑物顶面风荷载

通过风洞模拟实验，研究了在锥形涡诱导下建筑物顶面风荷载的特性. 针对一个立方体形模 型，给出了在不同风向角下压力分布的系统结果，进而分析了产生这样结果的原因和流动机 理，以及建筑物顶面的分离流动结构随风向角改变而改变及其演化过程，并指出锥形涡的出 现是建筑物顶面局部出现峰值负压的主要原因. 实验结果还显示，在一定风向角下，由于分 离流动产生的锥形涡结构在屋顶局部可以诱导出时均负压峰值达$-1.0$以上；在风向角 $\beta =30^\circ$时，锥形涡强度达到最大，其诱导的顶面附近局部时均负压峰值可 达$-2.4$, 瞬时风压达$-5.7$以上，是导致建筑物屋顶在强风下损坏的根本原因.     

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.      

基于等温吸附的页岩水分传输特征研究

研究页岩的水分传输特征至关重要,不仅有助于认识页岩的物理化学性质,而且也有助于评价页岩气的吸附扩散和流动能力.本文设计了页岩的水分传输实验装置,采用美国伍德福德和中国南方龙马溪组页岩为研究对象,开展了不同温度、不同湿度下页岩的水分传输实验,研究了页岩的水分传输特征和影响因素.结果表明,页岩的水分吸附属于II型曲线,包含着单分子层吸附、多分子层吸附和毛细凝聚的过程,GAB模型可用于描述页岩的水分吸附过程;水分吸附随着相对压力的增大而增强,有机碳含量和温度对页岩水分吸附起着增强作用,而方解石会抑制页岩的水分吸附;随着相对压力的增大,页岩的水分扩散系数呈现先增大后减小随后增加的趋势,其系数大约在8.73$\times$10$^{ - 9}\sim $5.95$\times $10$^{ - 8 }$m$^{2}$/s之间;伍德福德页岩的等量吸附热均大于龙马溪页岩的等量吸附热,这与其页岩的成熟度有关.研究结果为认识页岩的物理化学性质和力学性能以及评价页岩气的吸附流动能力提供参考依据.      

复合材料水翼水动力与结构强度特性数值研究

针对复合材料水翼存在的流固耦合求解问题,结合其自身特有属性,对复合材料水翼结构变形特性进行了数值仿真计算研究.研究建立了复合材料水翼流固耦合数值计算模型,并将数值计算结果与Zarruk等的实验结果进行对比,验证模型的正确性,得出复合材料水翼尖端扭转角随雷诺数的增加而增加的研究结论.基于数值计算模型,系统地研究了不同铺层角对复合材料水翼水动力特性及强度特性的影响,结果表明:不同铺层角复合材料水翼的尖端扭转角,随铺层角的增大而减小,而其尖端位移量随铺层角的增大先减小后增大.为了削弱工程常数的影响对复合材料水翼变形的影响,研究提出了无量纲扭转角和无量纲位移量,进一步探究复合材料水翼结构的弯扭耦合作用对其变形特性的影响.最后利用蔡$\!$-$\!$-$\!$吴失效准则进行复合材料水翼强度特性的判断和分析,结果表明:不同铺层角复合材料水翼的蔡$\!$-$\!$-$\!$吴系数,随铺层角的增大呈现先减小后增大的趋势,其中0$^\circ$铺层时的复合材料水翼蔡$\!$-$\!$-$\!$吴系数最小,50$^\circ$铺层时的复合材料水翼的蔡$\!$-$\!$-$\!$吴系数最大.      

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

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

基于稀薄效应的微气体径向轴承稳态性能

针对微气体轴承给出参考努森数的定义,根据空气不同温度时的黏度,得到参考努森数的分布范围;考虑气体稀薄效应,给出基于Burgdorfer一阶滑移速度边界的微气体径向轴承润滑Reynolds方程的修正形式; 采用有限差分法求解修正的Reynolds方程,得到不同参考努森数$Kn_0$, 轴承数以及轴颈偏心率情况下轴承的压力分布、无量纲承载能力及偏位角. 数值分析表明:随气体稀薄程度的增强,气体径向轴承的压力明显降低,无量纲承载力降低,而轴承偏位角增大. 当偏心率小于0.6时,轴承偏位角变化平缓,受$Kn_0$数的影响不明显. 当轴承数较小时,气体稀薄程度对轴承的无量纲承载力、偏位角影响较小.      

靶板厚度对卵形弹丸垂直贯穿中等厚度混凝土靶的影响

为研究卵形弹丸贯穿中等厚度混凝土靶体的贯穿规律,开展直径60 mm尖卵形弹丸贯穿不同厚度混凝土靶体的侵彻实验,获得了不同撞击速度的弹丸贯穿不同厚度混凝土靶体的剩余速度规律。结合无网格SPH方法、RHT混凝土本构以及状态方程,对贯穿实验进行数值模拟,对不同工况下的弹丸过载规律以及靶体的损伤过程的分析发现:弹丸贯穿中等厚度混凝土靶体的贯穿过程分为开坑阶段、隧道稳定侵彻阶段以及靶背影响出靶阶段,在相同初始撞击速度下的靶背影响区的厚度随着靶体厚度的增加而增大。实验结果与数值模拟结果对比,表明模型能够有效模拟弹丸贯穿混凝介质问题,研究结果可为贯穿机理的研究提供参考。     

紫藤萝复叶气动特性的风洞实验研究

树叶的形状重构和减阻能力在太阳能帆板、机翼结构、仿生天线设计和新型发电技术等方面具有应用价值.紫藤萝羽状复叶垂直悬挂于风洞中,在风速0~25m/s范围内进行正面和反面迎风测试.发现存在前期稳定、中间过渡和后期稳定3个阶段以及5个临界风速.在前期阶段叶轴随风速弯曲变化剧烈,出现小叶分层飞翼和分层多形状稳定.过渡阶段出现叶轴大幅低频振动和部分小叶小幅高频振动两种不稳定形式.后期出现两层或单一整体稳定,横截面形状分为锥形、楔形和U形.随着风速增大,复叶宽度减小,小叶层数逐步减少,直至出现流线形单一整体.随着雷诺数增大,复叶阻力系数先是快速下降,后又缓慢地趋于常数.复叶Vogel负指数绝对值$\vert \alpha \vert$随小叶数目的增大而增大.反面迎风时$|\alpha|$比正面迎风时大,但随着小叶数目增加两者趋于一致.当复叶旋涡脱落频率与叶轴固有频率接近时,叶轴出现大幅振动.理论分析得到叶轴振动的第二临界风速$V_2/\sqrt{E/\rho}$是$b/l$和$d/l$的函数,其中$E$,$\rho$,$d$和$l$分别为叶轴弹性模量、密度、直径和长度,$b$为变形后的复叶宽度,并由实验数据得到了其变化图.      

聚能射流侵彻页岩储层损伤裂隙形成机制

为研究药型罩对聚能射孔弹侵彻页岩储层的射孔和损伤致裂效果的影响机理，建立了射孔弹-空气-页岩三维模型，设置药型罩的锥角分别为50°、60°、70°和80°，壁厚分别为0.5、1.0和1.5 mm，材料分别为铜、钢、钛和钨。利用ANSYS/LS-DYNA软件进行数值计算，分别从射流速度与形态、页岩射孔效果及页岩孔裂隙形成规律特征等进行系统性分析。研究结果表明：在射孔弹结构中，随着药型罩锥角的减小，射流速度提高、杵体速度降低、侵彻深度增大同时开孔孔径减小。在一定范围内，适当减小药型罩的壁厚，可以提高射流速度、减小杵体质量、增大侵彻深度和开孔倾斜度。药型罩材料对射流速度、杵体结构和页岩射孔效果均有显著影响，其中钨药型罩射孔弹的侵彻深度最大但开孔孔径最小，钛药型罩射孔弹的侵彻深度最小但开孔倾斜度最大，铜比钢药型罩射孔弹的侵彻深度略大但开孔孔径略小。通过研究不同对照组的页岩孔裂隙形成规律特征发现，页岩孔裂隙发育主要发生在杵体对页岩的再扩孔阶段，减小射流初始扩孔孔径、增大杵体直径、提高杵体速度，可以促进页岩孔裂隙发育程度。     

THE EFFECTS OF CAVITY THICKNESS ON SEISMIC BEHAVIOR OF FOAM CORE SANDWICHED WALLS

为研究空腔厚度对现场发泡夹心墙平面内、平面外抗震性能的影响,对20 片现场发泡夹心墙和2 片实心墙进行平面内、平面外受力的数值试验和模型试验,并验证数值试验的正确性和有效性. 对不同空腔厚度的现场发泡夹心墙进行平面内、平面外受力的数值试验研究. 结果表明平面内受力时,对相对位移差的影响最明显,平面外受力时,对承载力的影响最明显. 总体上看,对平面内的影响要大于对平面外的影响. 墙体的协调变形能力随空腔厚度的增加而明显降低,因此建议在高烈度区现场发泡夹心墙的空腔厚度不宜大于100 mm.     

粒料固有各向异性的离散元模拟与细观分析

料在摊铺后形成的颗粒定向排列将导致其力学性质的固有各向异性. 依据粒料的实际不规则形状, 构建了可模拟粒间咬合嵌挤作用的三维离散元复杂形状颗粒; 生成了5 种不同沉积方向的各向异性试件和1种各向同性试件, 对比了各试件在三轴压缩试验中的宏观力学特性差异; 引入组构张量以量化各向异性程度, 利用玫瑰图表达接触法向与接触力的分布特征, 探究了粒料各向异性的细观发展规律. 结果表明: 颗粒长轴愈趋向水平排布, 峰值应力比愈大, 剪缩与剪胀程度愈明显; 相较于各向同性试件, 沉积角$\theta$为料在摊铺后形成的颗粒定向排列将导致其力学性质的固有各向异性. 依据粒料的实际不规则形状, 构建了可模拟粒间咬合嵌挤作用的三维离散元复杂形状颗粒; 生成了5 种不同沉积方向的各向异性试件和1种各向同性试件, 对比了各试件在三轴压缩试验中的宏观力学特性差异; 引入组构张量以量化各向异性程度, 利用玫瑰图表达接触法向与接触力的分布特征, 探究了粒料各向异性的细观发展规律. 结果表明: 颗粒长轴愈趋向水平排布, 峰值应力比愈大, 剪缩与剪胀程度愈明显; 相较于各向同性试件, 沉积角$\theta$为$0^\circ$时, 峰值应力比和最大体积压缩应变分别提高了12.6\%和18.8\%, 其原因在于加载过程中颗粒旋转和滑动百分比更小, 内部调整时间更短、更易被剪密; 固有各向异性对颗粒法向接触力分布的影响不大, 但显著影响接触法向分布特征; 剪切过程中, $\theta$为$90^\circ$时的接触法向各向异性系数先快速减小后逐渐增大, 而$\theta$为$0^\circ$到$60^\circ$时则呈现出增大后稍有回落或趋于稳定的趋势, 且$\theta$ 愈小的试件各向异性系数增大愈快.