一种非统一黏塑性本构的有限元实现

在与率相关黏塑性本构模型的背应力分量中加入静态回复项通常用来更精确地描述材料应力松弛效应．本文基于包含静态回复项的非统一黏塑性本构模型,提出了对应的隐式应力积分算法．该模型可简化为两个关于偏背应力和等效应力张量的方程组,通过牛顿迭代法进行求解．在应力积分算法中采用了Return-Mapping求解策略．此外还推导出了材料的一致切线刚度矩阵来判断有限元分析收敛性．在此基础上开发出了适用于ABAQUS工程软件的UMAT用户材料子程序．文中选取了高温下镍合金材料在伴随应力松弛的循环加载路径下的试验结果来验证算法的有效性．从比较结果中可以发现：用户子程序即使在大时间增量下也具有良好的收敛性,静态回复项的引入使得模型对应力松弛效应的模拟更加准确．     

基于牛顿-拉夫逊的拉压不同模量问题的数值求解

针对拉压模量不同引起的材料本构非线性,首先,通过引入改进的Heaviside函数将本构方程连续光滑化;然后,基于特征值与特征向量的求导策略,推导有限元求解模型中切线刚度矩阵的列式;最终,提出基于牛顿-拉夫逊迭代格式的拉压不同模量问题的数值求解算法。数值算例验证了本文算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于工程分析中大规模计算问题的求解。     

冲击载荷作用下岩体拉-压损伤破坏的边坡抛掷爆破模拟

为了获得边坡台阶爆破时岩体在冲击载荷作用下破碎和抛掷过程中各物理参量的变化规律,将构建的岩体动态拉-压损伤本构关系嵌入模拟软件进行数值分析。结果表明:振动波三轴合成速率曲线与质点振动位移曲线的变化趋势所对应的时间节点和步长基本吻合,可作为降振减灾的判定指标;岩体最早于0.6 ms左右从坡脚位置产生裂纹,直至12.5 ms左右裂纹裂隙完成扩展,所形成的炮孔粉碎区半径约28 cm;抛掷块状分离现象从炮孔药包的中间部位开始,最大抛掷速度集中于该部位至边坡自由面之间的垂直区域内,边坡自由面抛掷速度小于炮孔周围岩块抛掷速度,导致抛掷过程中形成二次挤压破碎现象;破碎大块主要来源于边坡坡脚、炸药与堵塞物接触面两侧围岩以及台阶顶部自由面处,大块块体直径分布于1.6～2.7 m范围。     

拉压不同弹性模量薄板上圆孔的应力分析

采用弹性理论研究了拉压不同弹性模量薄板上圆孔的孔边应力集中问题．采用广义虎克定律推导出了拉压不同弹性模量薄板上圆孔边的应力平衡方程,并联合利用应力函数及边界条件得到了拉压不同弹性模量薄板上圆孔边的应力表达式．算例分析表明,当薄板材料的拉压弹性模量相差较大时,采用经典弹性理论研究薄板上圆孔的孔边应力是不合适的,当经典弹性理论与拉压不同弹性模量弹性理论的计算结果间的差别超过工程允许误差5%时,应该采用拉压不同弹性模量弹性理论进行计算．     

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

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

不同加载速度下脆性颗粒的破坏特性

对直径为8mm的K9玻璃球进行了加载速度为2×10-7和2×10-6m/s的准静态单轴压缩实验以及加载速度为3.4、7.1和10.6m/s的动态单轴压缩实验,研究了K9玻璃单颗粒破碎强度的Weibull分布特性,结合破碎产物的形貌特征,分析了不同加载速度下脆性材料拉伸破坏机制和剪切破坏机制的转变过程,提出了一种拉剪耦合的时序破坏模型,由此揭示了加载速度与3个破坏区的关系。考虑拉伸和剪切失效准则,应用ABAQUS软件进行数值模拟,并初步验证了该破坏模型的冲击过程。研究结果对于认识脆性颗粒材料的动态破坏具有很好的参考意义。     

卸荷裂隙岩体线弹性阶段的广义变分问题

阐述了卸荷裂隙岩体线弹性阶段本构关系及弹性问题的基本方程,运用弹性问题基本方程中的控制方程和边界条件构造卸荷裂隙岩体线弹性阶段的变分问题.构造卸荷裂隙岩体线弹性阶段的变分问题事先构造了一个适当的最小位能泛函,在最小位能泛函中引入两个待定的拉氏乘子λij和λi,把变分约束条件吸收到泛函中去,从而建立卸荷裂隙岩体线弹性阶段的新泛函.然后将新泛函中的εij、μi、λij、λi作为独立变量,同时考虑新泛函的变分驻值条件识别待定拉氏乘子λij、λz,最后证明卸荷裂隙岩体线弹性阶段双变量的广义变分原理.     

飞秒泵浦不同锥长微结构光纤超连续谱的产生

对飞秒脉冲泵浦下,不同锥长及锥腰直径的微结构光纤的超连续谱产生进行了实验研究。采用“快速低温拉锥方法”,在保持d/Λ不变的情况下,对实验室自制的空气孔间距Λ=6.53 μm,归一化孔径d/Λ=0.79的微结构光纤进行了拉锥,分别得到6,8,10 mm等不同锥长微结构光纤。理论计算表明,随着锥长变长,锥腰直径变小,锥腰处零色散波长向短波移动：未拉锥及6,8和10 mm锥微结构光纤锥腰处零色散波长分别为1 129,885,806和637 nm。利用中心波长为810 nm,重复频率76 MHz,脉宽120 fs的钛蓝宝石飞秒激光器对拉锥后微结构光纤进行了实验研究：锥长为6 mm时,泵浦光中心波长位于整根光纤的正常色散区,锥腰的零色散点附近,内脉冲拉曼散射和级联四波混频是光谱初始展宽的主要因素。泵浦功率达到450 mW时,在可见波段390～461 nm及红外波段1 134～1 512 nm形成-5 dB的平坦宽带连续光谱。泵浦功率达到500 mW时,出现366～2 450 nm覆盖紫外、可见、近红外、中红外的超连续谱,其光谱红蓝移边缘已经接近实验用微结构光纤的传输带宽。锥长为8 mm、泵浦功率为450 mW时,在群速度匹配和群加速度失配的共同影响下,连续谱蓝移边缘达到366 nm,比6 mm锥时蓝移9 nm;锥长为10 mm时,由于锥腰处零色散点移动到可见光区域,可见区光谱仍能满足相位匹配条件。通过级联四波混频效应,在可见区域实现了频率上转换及光谱蓝移。泵浦光功率达到500 mW时,在382～412 nm得到谱宽仅为30 nm,转换效率达到27.7%的频率上转换。     

硅单晶提拉生长中拉速稳定性研究

为解决硅单晶提拉生长过程中拉速不稳定、波动较大的问题,在分析晶体生长界面热量及质量传输的基础上,通过控制温度补偿速率来抑制拉速大幅波动.采取基于遗传算法优化隶属度函数的模糊控制策略,对温度补偿速率进行控制,调节加热功率,使炉内热环境处于适宜晶体生长的范围.实验结果表明,在提高系统控径精度的同时,拉速稳定性也有显著提高,大幅波动明显减少.     

HHV和HVV耦合的有效拉氏量及费曼规则

首先包含Higgs粒子的六维算符用物理场来表示,然后给出包含HHV和HVV相互作用顶点的有效拉氏量,其中V=A,Z,W^±。最后,推导出每个六维算符和有效拉氏量的费曼规则。