力学学报

姓名
邮箱
手机号码
标题
留言内容
验证码

2017, 49(4).

PDF

马天雪, 苏晓星, 董浩文, 汪越胜, 张传增

声光子晶体是一种同时具有光子和声子带隙的周期性结构.声光子晶体为同时控制电磁波和弹性波的传播提供了一个系统的平台，并在光学、声学及声光多功能器件、腔光力学等领域展现了十分广阔的应用前景.论文首先介绍了声光子晶体的基本概念，包括声光子晶体的材料、空间周期分类、能带结构的计算方法等；阐述了不同体系下声光子晶体双重带隙的特性，以及拓扑优化方法在声光子晶体带隙优化方面的应用；然后简要介绍了腔光力学，以及计算声光耦合作用的准静态方法和光力耦合系数方法，并针对当前各种声光子晶体结构中声光耦合作用的研究进行了阐述；还进一步介绍了声光子晶体波导和传感器的相关研究；最后，基于当前声光子晶体的研究进展对未来的研究方向进行了展望，其中涉及到增强声光子晶体谐振腔的声光相互作用、三维声光子晶体的研究、声光超构材料的设计、声光子晶体器件设计与应用等.

2017, 49(4): 743-757. doi: 10.6052/0459-1879-17-130

HTML PDF

超弹性薄膜与可压缩基底双层结构表面失稳分析

黄春阳, 唐山, 彭向和

当上层超弹性硬质薄膜和下层可膨胀基底构成的双层结构受压时，薄膜的自由表面可通过形成褶皱降低系统能量.研究表明，上下两层的模量比不同时，上层弹性硬质薄膜将表现出不同的表面失稳模式.本文提出了一种新颖的方法可有效抑制双层软材料的表面失稳，即改变基底材料的泊松比，这种方法同时适用于不具有应变硬化的软材料.首先基于Neo-Hookean模型发展了小变形条件下双层结构表面失稳的理论模型，通过半解析的方法得到了表面失稳的临界应变；然后通过有限元计算与模拟，进一步验证了负泊松比基底可延缓表面失稳.结果表明：（1）当双层结构基底泊松比为正且趋于0.5（不可压缩）时，双层结构在较小的压缩应变下出现表面失稳；（2）当基底的泊松比为负且趋于-1时，可被压缩至46%而不出现表面失稳，即可膨胀基底能有效抑制薄膜的表面失稳.本文发展的方法及主要结果可为延展性电子器件的设计提供指导.

2017, 49(4): 758-762. doi: 10.6052/0459-1879-17-161

HTML PDF

基于微结构动态演化机制的单晶镍基高温合金晶体塑性本构及其有限元模拟

熊骏, 李振环, 朱亚新, 黄敏生

因其优异的高温力学性能，镍基单晶高温合金在航空航天和能源等领域得到了广泛的应用.镍基单晶高温合金优异的高温性能来源于其特有的两相微结构.基于代表体胞模型及分块均匀化方法，以位错密度为主要内变量，发展了一个包含两相微结构和位错演化信息的单晶镍基高温合金塑性行为的本构模型.该本构模型充分考虑了镍基单晶合金中位错在基体相和沉淀增强相中的多种演化机制，例如，基体位错八面体滑移、立方滑移、位错攀移、交滑移、位错弓出、位错切过沉淀增强相以及位错Kear-Wilsdolf（K-W）锁形成与解锁等.在商用有限元软件ABAQUS的框架下，编制了UMAT用户材料子程序.利用该用户子程序，对单晶和多晶镍基高温合金在不同温度、不同加载方向下的单调塑性、循环塑性、蠕变等典型行为进行了计算模拟.结果表明：该晶体塑性本构模型能“统一地”刻画镍基高温合金在不同温度、不同方向下的多种变形行为，并与实验结果具有良好的一致性.

2017, 49(4): 763-781. doi: 10.6052/0459-1879-17-183

HTML PDF

航行体水下发射流固耦合效应分析

杜特专, 王一伟, 黄晨光, 廖丽涓

对于水下发射过程来说，掌握水动力载荷形成机理与结构响应特征是一个亟待解决的问题.研究该问题需要考虑含相变的复杂多相流动，变约束的结构运动以及这二者之间的耦合效应.本文采用松耦合的方法，以流体求解器为主体，将自编的固体结构程序接入流体求解器中，在每个时间步长内分别对流体动力学方程和固体结构动力学方程进行求解，通过流固界面之间的数据交换实现耦合计算.其中，流体求解器基于雷诺平均纳维斯托克斯方程，采用单流体模型处理多相流问题，引入空化模型描述空化相变，采用修正的湍流模型模拟混合物的湍流效应，并采用动网格技术处理移动边界问题.航行体的刚体运动和结构振动分开求解.结构求解器采用等效梁模型描述结构的振动，通过坐标变换给出了随体坐标系下的结构振动方程，求解方法采用时域积分法.所建立的流固耦合方法不仅能够捕捉到自然空化的演化情况，还可获得航行体所受水动力、结构振动响应以及截面的弯矩，获得了实验的验证.基于该方法研究了结构刚度、发射速度对空泡溃灭与结构振动耦合效应的影响规律.结果表明，同步溃灭是影响结构载荷的主要因素，包括溃灭压力幅值，溃灭压力作用位置，以及溃灭压力与结构振动的相位关系.

2017, 49(4): 782-792. doi: 10.6052/0459-1879-16-401

HTML PDF

改进巴西试验:从平台巴西圆盘到切口巴西圆盘

王启智, 李炼, 吴礼舟, 黄润秋

巴西试验是测试混凝土、陶瓷、岩石等脆性或准脆性材料拉伸强度的标准方法，在材料科学和土木工程中有广泛的应用，研究该方法的改进在科学和工程界都受到日益增长的关注.我们对巴西试验的原始试样即巴西圆盘（Brazilian disc，BD）提出一种新的构型：切口巴西圆盘（grooved Brazilian disc，GBD），利用GBD加载直径两端的窄而浅的切口，消除完整巴西圆盘对径压缩加载时可能非中心起裂的缺点，这一点与我们过去提出的平台巴西圆盘（flattened Brazilian disc，FBD）相同.GBD的构型更普适，因为它涵盖了巴西圆盘和平台巴西圆盘.对于求解GBD切口顶点及其正前方近邻的应力，提出两个近似解析模型，分别求出对应的压应力和拉应力，用叠加原理最后推出的公式能够定性地预测应力的变化趋势，以及试样几何参数的影响.近似解析公式，有限元数值计算，对比实验的结果都证实，对径压缩圆盘的几何与表面边界条件的一个微小的改变，对巴西试验试样的几何稳定性和断裂过程产生显著和有利的影响.利用推出的切口应力公式进行试样几何参数的优化选择，可以使GBD比BD和FBD更具优越性.

2017, 49(4): 793-801. doi: 10.6052/0459-1879-17-031

HTML PDF

转捩位置对全动舵面热气动弹性的影响

刘成, 叶正寅, 叶坤

高超声速附面层的转捩预测一直是流体力学研究中的难点，转捩前后物面的摩擦系数和传热系数会发生改变，转捩位置的不同会影响到飞行器表面热环境，进而使得飞行器的气动弹性特性发生显著变化.鉴于高超声速附面层转捩预测的不确定性，本文分析了转捩位置对高超声速全动舵面热气动弹性的影响.首先分别用层流模型和湍流模型求解N-S方程，得到气动热环境，并对气动热进行参数化；然后在不同转捩位置情况下构造出不同转捩位置的热分布模型，基于此种温度分布，结合热应力和材料属性的影响分析结构的热模态，将结构模态插值到气动网格上，采用基于CFD的当地流活塞理论进行气动弹性分析.以M=6，H=15 km的某舵面为对象进行计算，结果表明：（1）随着转捩位置向后缘移动，结构频率上升，结构颤振速度呈增大趋势，转捩位置的变化能够带来颤振临界速度最大6%的变化量；（2）当转捩位置位于舵轴附近时，结构的颤振特性变化剧烈.通过刚度特性的分解和分析发现，导致颤振特性变化的主要因素在于舵轴的刚度特性变化，舵轴的影响量占整个结构刚度特性变化量的80%以上.

2017, 49(4): 802-810. doi: 10.6052/0459-1879-17-016

HTML PDF

疏液壁面上亲液杂质对滑移特性的影响

鲍路瑶, 罗凯, 文俊, 胡海豹

纳米流动系统具有高效、经济等优势，在众多领域具有广泛的应用前景.因该类系统具有极高的表面积体积比，致使界面滑移效应对流动具有显著影响.论文采用非平衡分子动力学模拟方法，研究了疏液壁面表层混入少量亲液原子时纳米通道内液体的滑移特性，并基于分子动理论解释了其影响机制.研究发现，亲液杂质（均布或集中）对液体法向密度振荡程度影响较弱，但会显著改变壁面附近类固体层的分布和液体滑移规律；随亲液杂质占比增加，液体类固体现象更趋明显，壁面处液体接触密度也线性增大，但通道内液体平均速度逐渐降低，滑移长度也迅速减小；相比于集中的亲液杂质，均匀分布亲液杂质对滑移的弱化效应更强.如当亲液杂质占比为28%时，其滑移长度比单纯疏液表面的降低率从50%扩大至56%.基于分子动理论的分析发现，亲液杂质会导致杂质原子附近第一液体层内的原子发生跃迁的能垒加大，即弱化了液体原子的流向跃迁行为，从而降低了滑移量；相比于集中杂质，均匀分布的杂质还会降低固液间的非公度性，致使滑移特性破坏更严重.

2017, 49(4): 811-817. doi: 10.6052/0459-1879-16-368

HTML PDF

柔性圆柱涡激振动流体力系数识别及其特性

徐万海, 马烨璇, 罗浩, 栾英森

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前，海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据，存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用，建立了柔性圆柱涡激振动流体力模型，运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数，设计并开展了拖曳水池模型实验，实验用柔性圆柱模型的质量比为1.82，长径比为195.5.通过与刚性圆柱流体力系数对比，深入分析了柔性圆柱流体力系数的特性.结果表明：柔性圆柱在一阶模态控制区，流体力系数随约化速度变化趋势与刚性圆柱大致相似；二阶模态控制区，升力系数和脉动阻力系数显著增大；附加质量系数在响应频率较低时与振动位移的相关性增强；当响应频率较低时，振动位移较大区域为能量耗散区，当响应频率较高时，振动位移较大区域为能量输入区.

2017, 49(4): 818-827. doi: 10.6052/0459-1879-16-263

HTML PDF

低渗透煤层气藏中气-水两相不稳定渗流动态分析

刘文超, 刘曰武

针对低渗透煤层渗流问题，考虑启动压力梯度及其引起的动边界和动边界内吸附气解吸作用的渗流模型研究目前仅限于单相流，而更符合实际的气-水两相渗流动边界模型未见报道.本文综合考虑了煤层吸附气的解吸作用、气-水两相渗流、非达西渗流、地层应力敏感等影响因素，进行了低渗透煤层的气-水两相渗流模型研究.采用了试井技术中的“分相处理”方法，修正了两相渗流的综合压缩系数和流度，并基于含气饱和度呈线性递减分布的假设，建立了煤层气藏的气-水两相渗流耦合模型.该数学模型不仅可以描述由于低渗透煤层中渗流存在启动压力梯度而产生的可表征煤层有效动用范围随时间变化的移动边界，还可以描述煤层有效动用范围内吸附气的解吸现象以及吸附气解吸作用所引起的煤层含气饱和度的上升；为了提高模型精度，控制方程还保留了二次压力梯度项.采用了稳定的全隐式有限差分方法进行了模型的数值求解，并验证了数值计算方法的正确性，获得了模型关于瞬时井底压力与压力导数响应的双对数特征曲线，由此分析了各渗流参数的敏感性影响.本文研究结果可为低渗透煤层气藏开发的气-水两相流试井技术提供渗流力学的理论基础.

2017, 49(4): 828-835. doi: 10.6052/0459-1879-17-034

HTML PDF

不同拉压特性的厚壁圆筒极限内压统一解

赵均海, 姜志琳, 张常光, 曹雪叶

厚壁圆筒在实际工程领域中应用广泛，若能精确计算出极限内压，对预防事故发生，降低风险有重要意义.工程中存在许多材料，其拉压强度和拉压模量均存在差异，这些差异对极限内压的大小有显著影响.以往研究表明，仅考虑拉压强度与拉压模量的一个方面，计算结果与实际情况存在一定的误差.本文基于双剪统一强度理论，综合考虑中间主应力效应及材料拉压强度和拉压模量的不同，推导了内压作用下厚壁圆筒的弹、塑性状态的应力分布及弹性极限内压、塑性极限内压与安定极限内压的统一解，通过与其他文献对比分析验证了本文计算结果的正确性，分析了半径比、统一强度理论参数、拉压强度比与拉压模量系数对弹性极限内压、塑性极限内压及安定极限内压的影响.结果表明：统一解均随半径比和统一强度理论参数的增大而增大，随拉压强度比的增大而减小，弹性极限内压随材料拉压模量系数的增大而减小，当壁厚增加到一定值后，安定极限内压随材料拉压模量系数的增大而减小；材料的拉压模量不同、拉压强度差异对厚壁圆筒的安定性影响显著，考虑中间主应力效应可使材料的潜能得到更充分发挥，极限内压随半径比的变化规律可为选择合理壁厚提供参考，该结论可为厚壁圆筒的工程应用提供理论依据.

2017, 49(4): 836-847. doi: 10.6052/0459-1879-17-006

HTML PDF

半无限板边缘裂纹的权函数解法与评价

童第华, 吴学仁, 胡本润, 陈勃

权函数法是求解裂纹体在任意受载条件下的应力强度因子和裂纹面位移等断裂力学参量的高效、高精度方法，与有限元等数值方法相比，在求解效率和可靠性方面均具有明显优势.针对半无限板边缘裂纹，系统分析了在国际断裂力学界较有代表性的Wu-Carlsson、Glinka-Shen和Fett-Munz三种解析形式的权函数法，进而以在远端均匀加载下的半无限板边缘裂纹面位移Wigglesworth解析解导得的权函数及其对应的格林函数解（即裂纹面受一对单位集中力作用下的应力强度因子）为基准，沿整个裂纹长度对3种权函数的精度逐点进行比较，并与文献中基于其他方法求得的权函数做了广泛对比，包括Bueckner，Hartranft-Sih以及Wigglesworth利用不同解析方法推导出的高精度的权函数.研究了3种参考载荷（均布/正反向线性分布应力、集中力）及其不同组合，以及裂纹嘴位移的几何条件对权函数精度的影响.结果表明，基于一种参考载荷下的裂纹面张开位移比基于两种参考载荷下的应力强度因子所得到的权函数具有更高的精度，而且后一种方法的精度明显受到所选参考载荷组合的影响；裂纹面位移在裂纹嘴处三阶导数等于零的条件对基于一个参考解的权函数精度的改进效果较小.最后给出了利用各种权函数方法计算得到的4种载荷条件下的应力强度因子，并对结果进行了比较.

2017, 49(4): 848-857. doi: 10.6052/0459-1879-17-024

HTML PDF

弹塑性微凸体侧向接触相互作用能耗

高志强, 傅卫平, 王雯, 康维超, 吴洁蓓, 刘雁鹏

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设，忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用，这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面，从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力，获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理，从而获取加、卸载法向合力与总变形之间的关系，由于法向分力产生的塑性变形及切向分力产生的摩擦，导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分，获得一个周期的微凸体接触能耗，包括应变能耗及摩擦能耗.仿真分析表明：微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大，能耗越大，且在弹性阶段最为明显.在弹性阶段，仅存在侧向的摩擦能耗，故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段，由于微凸体接触能耗为应变能耗，且接触角对其能耗影响甚微，故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型，本文提出的模型与Bartier的实验结果更吻合.

2017, 49(4): 858-869. doi: 10.6052/0459-1879-17-103

HTML PDF

纯铜后继屈服面的测试与晶体塑性模型模拟

张克实, 黄世鸿, 刘贵龙, 陆大敏

通过纯铜薄壁圆管试样的实测和晶体塑性模拟，用单试样法和多试样法对分别经历拉伸、扭转和组合拉扭变形的试样后继屈服面进行研究.考虑预变形方式、测点数目、测试顺序和指定平移应变等不同条件，对后继屈服面测定结果差异及屈服面内凹现象进行探讨.在此基础上，比较了单试样和多试样两种方法的合理性与有效性.数值模拟采用能反映Bauschinger效应的晶体塑性模型，试样有限元模型每个单元的晶体取向均随机生成，能反映多晶材料变形特征.模拟试验加载过程与真实试验一致.研究表明：（1）采用本文方法可再现真实试验过程，模拟后继屈服面测试展示的现象与实测相近，证实了方法的有效性和合理性；（2）模拟测试与实测均发现，薄壁圆管组合拉扭加载测得的后继屈服面可能出现内凹，单试样法测得屈服面的内凹现象尤为显著；（3）若试验材料的材质比较一致，用多试样法测试后继屈服比用单试样法更合理.

2017, 49(4): 870-879. doi: 10.6052/0459-1879-17-074

HTML PDF

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

武守信, 魏吉瑞, 杨舒蔚

以非局部塑性理论为基础，应用状态空间理论，通过局部和非局部两个状态空间的塑性能量耗散率等效原理，提出了一种求解应变局部化问题的新方法，以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程，并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法，对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明，塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛，应变局部化区域的尺寸只与材料内尺度有关，而对有限元网格的大小不敏感.对于一维问题，当有限元网格尺寸减小时，数值解收敛于解析解.对于二维剪切带局部化问题，数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时，剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明，所提方法可以克服经典连续介质力学模型导致的网格相关性问题，从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C⁰连续性，因而容易在现有的有限元程序中实现而无需对程序作大的修改.

2017, 49(4): 880-893. doi: 10.6052/0459-1879-16-330

HTML PDF

透射边界条件在波动谱元模拟中的实现:二维波动

邢浩洁, 李鸿晶

将邢浩洁和李鸿晶提出的多次透射公式（multi-transmitting formula，MTF）的谱元格式应用于均匀介质中线弹性SH波动问题的谱元模拟.假定紧邻人工边界的一层谱单元为具有直线边界的四边形单元，以保证每个人工边界节点都唯一对应一条指向内域的离散网格线.人工边界节点在某时刻的位移由该离散网格线上的节点在前若干时刻的位移确定，按照MTF谱元格式进行计算.通过平面波以一定角度传播的外源问题算例和点源脉冲自由扩散的内源问题算例，验证了方法的可行性以及对实际复杂波动问题的适用性.通过不同类型初值问题算例，在时域内分析了插值多项式阶次、人工波速和透射阶次三个参数对反射误差的影响.结果表明：插值多项式阶次较高的格式会表现出更好的精度，但总体上对反射误差的影响较小；人工波速对反射误差具有显著影响，当人工波速小于介质物理波速时反射误差较大，而当人工波速等于或稍大于介质物理波速时反射误差处于较低水平；透射阶次对反射误差具有决定性影响，表现在不失稳的情形下提高透射阶次能够迅速降低反射误差，但内源问题从三阶MTF开始出现飘移失稳，外源问题从二阶MTF开始出现轻微的飘移失稳.

2017, 49(4): 894-906. doi: 10.6052/0459-1879-16-393

HTML PDF

具有曲率和惯性非线性以及材料内阻的旋转复合材料轴的主共振

任勇生, 姚东辉

旋转复合材料轴作为一类典型的转子动力学系统，在先进直升机和汽车动力驱动系统中有着广阔的应用前景.研究旋转复合材料轴的非线性振动特性具有重要的理论与实用价值.然而，目前有关旋转轴的非线性振动研究仅限于各向同性金属材料轴，很少考虑材料内阻的影响.本文研究具有材料内阻的旋转非线性复合材料轴的主共振.非线性来源于不可伸长复合材料轴的大变形引起的非线性曲率和非线性惯性，材料内阻来源于复合材料的黏弹性.动力学建模计入转动惯量和陀螺效应.基于扩展的Hamilton原理，导出具有偏心激励的旋转复合材料轴的弯-弯耦合非线性振动偏微分方程组.采用Galerkin法将偏微分方程离散化为常微分方程，采用多尺度法对常微分方程进行摄动分析，导出主共振响应的解析表达式.对内阻、外阻、铺层角、长径比、铺层方式和偏心距进行数值分析，研究上述参数对旋转非线性复合材料轴的稳态受迫振动响应行为的影响.研究发现，角铺设复合材料轴的内阻系数随着铺层角的增大而增大；内阻对主共振响应特性的影响主要体现在对抑制振幅和改变频率响应的稳定性方面；发生在正进动固有频率附近的主共振响应具有典型的硬弹簧非线性特性.本文提出的模型能够用于描述旋转复合材料轴的主共振特性，是对不可伸长旋转金属轴非线性动力学模型的重要推广.

2017, 49(4): 907-919. doi: 10.6052/0459-1879-17-002

HTML PDF

两自由度舵-轴系统振动三维效应修正模型

王人凤, 尤云祥, 陈科, 段金龙

考虑到小展弦比舵所存在的三维效应，利用附加质量系数ε和环量系数δ对经典Theodorsen两自由度运动方程进行修正，并与经典颤振实验结果进行比较，验证了修正后两自由度运动方程的适用性.质量比μ的不同会引起两自由度舵-轴系统振动V-g曲线形态的差异，故根据V-g曲线形状的不同将系统的振动分为第一类振动和第二类振动，其对应情况下可能发生的颤振为第一类颤振和第二类颤振.利用修正后的两自由度颤振理论模型分析了支撑刚度k_h、扭转刚度k_α、舵弦向重心位置x_α和初始攻角AOA对舵-轴系统颤振特性的影响规律，并通过开展相关实验对理论计算值进行验证，实验结果与计算值吻合良好.计算结果表明，k_h，k_α，x_α和AOA对颤振速度V_F存在显著影响，它们可以分别在一定的取值范围内导致系统发生第二类颤振.并且，V_F随k_h的增大单调增大，随k_α和x_α的增大先增大再减小，随AOA的增大则逐渐减小.其中，令V_F存在非零值的x_α取值范围狭小，反映了系统振动形态对x_α的敏感性.因此，在设计阶段避免将x_α设置在这个狭小的范围内可以降低颤振的发生几率.另一方面，由于V_F对k_h和k_α的反应缓慢，一旦颤振发生就可以通过将刚性轴锁紧来消除颤振效应.

2017, 49(4): 920-928. doi: 10.6052/0459-1879-17-042

HTML PDF

曲线加筋Kirchhoff-Mindlin板自由振动分析

刘璟泽, 姜东, 韩晓林, 费庆国

相比传统加筋板，曲线加筋板能够更充分地发挥材料力学性能.在加筋板力学分析中，厚板通常采用Reissner-Mindlin理论，然而当板厚较薄时易出现剪切自锁，离散的Kirchhoff-Mindlin理论采用假设剪切应变场可避免该问题.针对曲线加筋Kirchhoff-Mindlin板自由振动分析，采用离散的Kirchhoff-Mindlin三角形单元和Timoshenko曲梁单元分别模拟板和加强筋，根据板的位移插值函数及筋板交界面的位移协调条件，建立基于板单元位移自由度的有限元方程.为了验证方法的有效性和准确性，采用直线加筋薄板、曲线加筋薄板和厚板3种模型进行算例研究，通过收敛性和精度分析来选择合理的有限元网格密度.直线加筋薄板前20阶固有频率均与文献结果吻合良好；曲线加筋板算例中，本文方法满足收敛条件的板单元数目为2469，Nastran模型板单元数目为6243；本文所得曲线加筋板固有频率与Nastran计算结果最大误差为3.4%.研究结果表明，本文方法无需筋板单元共节点，可使用较少的有限元网格数量，并能够保证计算精度；在离散Kirchhoff-Mindlin三角形板单元基础上构造Timoshenko梁单元可同时适用于曲线加筋薄板与厚板自由振动分析.

2017, 49(4): 929-939. doi: 10.6052/0459-1879-17-041

HTML PDF

静动组合载荷下混凝土率效应机理及强度准则

路德春, 李萌, 王国盛, 杜修力

混凝土结构在受到动载荷作用之前，通常已承受着初始静载荷的作用.大量关于混凝土应变率效应的研究均没有考虑初始静载荷对动强度的影响，会导致过高地估计混凝土的动强度，使混凝土结构设计偏于危险.本文通过分析混凝土材料在静动组合载荷下的率效应机理，给出了初始静载荷的定义.在此基础上，推导了混凝土材料参数与初始静载荷和应变率的表达式，提出了建立静动组合强度准则的一般方法.通过材料参数反映初始静载荷与应变率的联合影响，给出了由初始有效静载荷、动态黏聚强度和摩擦强度共同组成的混凝土动态强度，将广义非线性强度准则发展为静动组合多轴强度准则.建立的强度面在相同初始静载荷下随应变率的增大向外扩张，在相同应变率下随初始静载荷的增大向里收缩，即混凝土的强度在相同初始静载荷下随应变率的增大而增大，在相同应变率下随初始静载荷的增大而减小.此外，当初始静载荷和应变率不变时，加载路径对混凝土材料的应变率效应无影响，但会影响混凝土材料的静水压力效应，即当初始静载荷和应变率固定不变时，静动组合强度面的位置和大小即可确定，不同加载路径下强度的不同是由于静水压力效应导致的.最后利用多组混凝土材料静动组合强度试验对建立的静动组合强度准则进行了验证.

2017, 49(4): 940-952. doi: 10.6052/0459-1879-17-013

HTML PDF

基于细观混凝土模型的时间逆转损伤成像方法

赵国旗, 仇亚萍, 骆英, 冯侃

提出了一种针对混凝土结构损伤检测的时间逆转损伤成像方法.以检测混凝土结构中与骨料尺寸相近的微小损伤为目的，引入细观混凝土随机骨料模型，该模型将混凝土结构视为由水泥浆基底、骨料及粘接层组成的三相复合材料，基于Monte Carlo随机样本原理并结合真实试件的骨料级配曲线建立.在数值模拟分析中，将生成含损伤的细观模型导入有限元分析软件进行超声波场模拟，同时采用自适应性强的时间逆转模型（timereversed model，TRM）进行损伤定位.TRM分为正向检测和逆时成像两个部分：正向检测过程得到包含损伤的一系列散射回波信号，从数值角度进行时间反演并作为逆时过程的输入信号；逆时成像过程选用等效弹性参数模型，几何尺寸与随机骨料模型相同，时反信号在相应几何位置同时加载形成时反波场，时反波场在损伤位置会发生干涉叠加从而导致能量峰值的出现，通过确定干涉峰值时刻，并获取该时刻对应原始波场以及小波变换能量场完成成像.与原始数据波场图相比，小波变换处理成像结果消除了杂波干扰，成像结果更加清晰.进一步对等效弹性参数的取值进行讨论，并且在骨料尺寸范围内调整损伤大小，结果显示成像结果匹配度高，对于非均质混凝土结构的损伤检测能很好满足损伤定位需求.由此证明，时间逆转成像方法对于具有复杂结构的混凝土材料的损伤检测具有较好的适用性.

2017, 49(4): 953-960. doi: 10.6052/0459-1879-17-007

HTML PDF

基于分子动力学-格林函数法的微凸体接触数值分析

黄仕平, 吴杰, 胡俊亮, 郑恒斌, 王卫锋

表面接触是摩擦的先决条件，其真实接触面积、压应力大小、空间分布等一直是接触力学关注的核心问题.采用分子动力学-格林函数法（GFMD）模拟粗糙面的接触过程，验证了其在大规模接触分析中的高效及准确性，同时探讨了由微球体组成的粗糙面的接触力学特性，并分析了分子尺度下的结果和传统力学模型计算结果的差异.结果表明，单个微凸体接触结果和分子动力学-格林函数法模拟所得非常接近，误差在5%以内.数值模拟发现，在微凸体高度符合高斯分布的情况下，接触面积和接触力成线性关系；在相同接触面积下，微凸体模型得出的接触力偏高，是上限值.微凸体模型没有考虑微凸体间的相互影响，实际是高估了弹性体的刚度；实际接触过程中微凸体相互影响，微凸体对临域形变影响尤其大，使接触区域更加离散.GFMD模型可以准确计算数十亿量级别分子、原子接触过程中真实接触面积及分布，为后续摩擦、滑移等分析提供可靠的参考.

2017, 49(4): 961-967. doi: 10.6052/0459-1879-17-084

HTML PDF

力学处2016年度重点项目结题审查情况简要介绍

詹世革, 张攀峰, 白坤朝, 王建山

对2016年度国家自然科学基金委员会数理科学部力学科学处结题的12项重点项目结题审查情况进行了简要介绍.给出了12项2016年度结题重点项目清单以及经评审专家组讨论形成的结题审查评价意见.

2017, 49(4): 968-972. doi: 10.6052/0459-1879-17-250

HTML PDF

留言板

2017年第49卷第4期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

留言板

2017年 第49卷 第4期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

2017年第49卷第4期