分数阶热弹理论下重力场对二维纤维增强介质的影响

基于Sherief等提出的分数阶广义热弹性耦合理论，研究了在热冲击作用下二维纤维增强弹性体的热弹性问题.考虑了重力对二维纤维增强线性热弹性各向同性介质的影响，建立了控制方程.运用正则模态法，经过数值计算，对控制方程进行求解，得到了不同分数阶参数和不同重力场下无量纲温度、位移和应力分量的表达式，以图形的方式展示了变量的分布规律并对结果展开了讨论.研究结果为:重力场和分数阶参数对纤维增强介质的位移及应力有着重要的影响，但对温度的影响有限.     

基于Green-Naghdi模型与三相延迟模型的耦合热弹性波传播

研究了基于G-N模型与三项延迟模型的热弹性波的传播特征.在各种广义热弹性理论中,热位移概念的引入使得G-N模型具有独特性质而引起广泛的关注和应用.基于G-N模型,最近一个三相延迟模型被提出.主要研究了这两种模型下第一和第二声波(耦合热弹性波)的色散和衰减特性,分析了温度场和位移场在两种模型下的振幅比和相位差.结果发现三相延迟模型更为通用和灵活,应得到更多的关注与应用.     

超声共振法测定圆柱体材料弹性常数的影响因素分析

为明确超声共振频谱分析法(RUS)测定圆柱体试样弹性常数的影响因素, 由此针对各向同性圆柱体开展实验, 观察圆柱体试样测量结果的分散性和稳定性, 且主要研究横纵比和有效模态数量对测量结果的影响, 并与拉伸试验的结果进行比较. 以测量结果标准差、变异系数和95%置信度下均值的区间估计3个指标来衡量结果的稳定性和分散性, 且从模态简并角度对实验结果进行分析. 最终结果表明, 当超声共振频谱分析法的有效模态数量在17左右、试样横纵比接近1时, 测量结果的稳定性和分散性较小.     

大展弦比气弹机翼柔性非线性变形的气动力效应分析与风洞试验

从飞行器刚弹耦合动力学模型出发，引入柔性机翼准定常假设，建立大柔性飞行器非线性静气动弹性气动力方程，利用非线性迭代求解思路模拟了柔性飞行器的静气动弹性响应行为，开展了大展弦比飞机静气动弹性风洞试验验证，采用气动力有限基本解与机翼的耦合计算，发现了大柔性飞机大变形状态下载荷及结构变形形式随风速的变化规律.传统基于小变形假设的线性分析方法和刚体分析由于无法考虑气动面随结构变形的曲面气动力因素和结构变形后的非线性刚度特性，均与风洞试验存在一定的误差.对于大展弦比柔性飞机的非线性静气动弹性分析十分必要.      

基于遗传算法的弹性地基加肋板肋梁无网格优化分析

基于遗传算法及一阶剪切理论,提出一种弹性地基上加肋板肋条位置优化的无网格方法.首先,通过一系列点来离散平板及肋条,并用弹簧模拟弹性地基,从而得到加肋板的无网格模型;其次,基于一阶剪切理论及移动最小二乘近似原理导出位移场,求出弹性地基加肋板总势能;再次,根据哈密顿原理导出结构的弯曲控制方程,并通过完全转换法处理边界条件;最后,引入遗传算法和改进遗传算法,以肋条的位置为设计变量、弹性地基板的中心点挠度最小值为目标函数,对肋条位置进行优化达到地基板控制点挠度最小的目的.以不同参数、载荷布置形式的弹性地基加肋板为例,与ABAQUS有限元结果及文献解进行比较.研究表明,采用所提出的无网格模型可有效求解弹性地基上加肋板弯曲问题,结果易收敛,同时基于遗传算法与改进混合遗传算法所提出的无网格优化方法均可有效优化弹性地基加肋板肋条位置,后者计算效率相对较高,只进行了三次迭代便可获得稳定的最优解,此外在优化过程中肋条位置改变时只需要重新计算位移转换矩阵,又避免了网格重构.     

板中热弹波传播: 一种改进的勒让德多项式方法

近年来, 超声导波因其衰减小, 传播距离远和信号覆盖范围广, 成为无损检测领域快速发展的方向之一. 然而, 基于超声导波的高温在线检测和激光超声技术却发展缓慢, 其关键在于热弹耦合波动方程求解难度大、传播与衰减特性研究困难. 作为一种有效的求解方法, 勒让德正交多项式方法已广泛应用于导波传播问题, 但该方法在求解热弹导波传播时存在两个不足, 限制其进一步的发展和应用. 这两个缺陷是: (1)求解过程中大量积分的存在, 致使计算效率低下; (2)仅能处理等热边界条件的热弹导波传播. 针对两项不足之处, 提出一种改进的勒让德正交多项式方法, 以求解分数阶热弹板中的导波传播. 推导求解方法中积分的解析表达式, 以提高计算效率; 引入温度梯度展开式, 发展适合勒让德多项式级数的绝热边界条件处理方法. 与已有文献结果对比表明改进方法的正确性; 与已有方法的计算时间对比说明改进方法的高效性. 最后将改进的方法用于求解分数阶热弹板中的导波传播, 研究分数阶次对频散、衰减曲线和应力、位移、温度分布等的影响.      

高压下β-InSe弹性常数、电子结构的第一性原理研究

基于密度泛函理论的第一性原理计算方法,本文研究了高压对β-InSe弹性常数、机械性能和电子结构的影响.在0～20 GPa范围内,随着压力的增大,β-InSe的晶格常数、晶胞体积逐渐减小,结构参数a/a_0、c/c_0、V/V_0单调减小.在0～12 GPa范围内,弹性模量G、E、B和泊松比v随着压力增大而增大,在16 GPa时大幅减小,G、E、B分别减小了34.9%、53.3%、82.9%.随着压力增大,Se-In和In-In原子之间的电荷密度增大,Se-In原子之间的共价键增强,层间距减小.而且,β-InSe在20 GPa时带隙消失,发生了半导体向半金属的相变.     

一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷，首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点，建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联，建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型，根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处，并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证，结果表明：拟合曲线与试验曲线吻合度较高，所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律，尤其对加速蠕变阶段的非线性特征描述效果很好，验证了模型的合理性．     

形状记忆合金椭圆孔口问题的有限单元法

基于Lagoudas形状记忆合金(SMA)三维本构模型,假设材料为各向同性,推导了SMA平面应力状态的增量型本构方程,继而编写了ABAQUS用户自定义材料(UMAT)子程序,研究了在双向拉伸情况下,外载荷、温度、椭圆孔口长短轴之比对超弹性SMA椭圆孔口板中应力诱发马氏体相变区的影响。数值结果表明:应力诱发马氏体相变首先发生在椭圆孔口长轴端点部位,在外加载荷作用下逐渐扩展到板内,并由内向外形成马氏体相区、相变混合区和奥氏体相区;SMA板内应力诱发马氏体完全相变区面积与施加外载荷成正相关,与温度成负相关;随着椭圆孔口长短轴之比增大,SMA板内应力诱发马氏体完全相变区面积呈现出先减小后增大的趋势;拉应力差值相同时,相较于拉应力沿椭圆孔口长轴方向较大的情况,当拉应力沿椭圆孔口短轴方向较大时,SMA板内完全相变区面积较大,椭圆孔口周边应力集中现象更明显。     

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.