一个反映土体复杂加卸荷路径的弹黏塑性模型

基于Perzyna 过应力理论,构造了一个土体各向异性弹黏塑性本构模型. 模型引入Wheeler 旋转硬化法则,能够较好的描述土体的初始各向异性和应力诱发各向异性. 借助ABAQUS 软件中UMAT 子程序接口,模型采用隐式积分算法——图形返回算法实现. 通过对一组复杂加卸荷路径的三轴不排水剪切试验(HKMD) 及一组分级加载三轴不排水蠕变试验(Sachville clay) 的模拟,表明本模型能够合理反映土体的率效应、应力松弛、蠕变及土体各向异性现象,验证了模型的有效性.     

泡沫铝率相关本构模型及其在三明治夹芯板冲击吸能特性的应用研究

高孔隙率泡沫铝芯体三明治板具有轻质、高比刚度和减振吸能等优良的力学特性和物理特性,被广泛地应用于碰撞吸能部件上.近年来,高孔隙率泡沫铝在动态压缩下是否具有应变率敏感性成为广大学者的研究焦点.论文建立了横观各向同性率相关本构模型来描述高孔隙率泡沫铝的应变率效应,给出了有限元的计算步骤,基于ABAQUS/Explicit平...     

Nonlinear viscoelastic multi-scale repetitive unit cell model of 3D woven composites with damage evolution

Three-dimensional (3D) textile composites have great potential applications to aircrafts and high speed vehicles because of the high strength/weight ratios and the capabilities of manufacturing complex, net-shape preforms. This paper reports the nonlinear viscoelastic responses and damage mechanisms of one kind of 3D textile composites, named as 3D orthogonal woven composite (3DOWC) under quasi-static tensile loading based on a micro/meso-scale repetitive unit cells (RUCs) model. In the RUCs model, the resin is described with a nonlinear viscoelastic material and the fibers/tows with an elastic material. The damage initiation and propagation in resin are simulated by the post-damage constitutive models with maximum principal theory failure criteria. The fibers/tows impregnated with resin are defined by elastic transverse-isotropic material model with ultimate strengths failure of ‘expanded smeared crack’ both along and perpendicular to fibers/tows axis direction. The engineering parameters and ultimate strengths of homogenized fibers/tows filled with matrix in meso-RUCs model are transferred from the numerical analysis of the micro-RUCs. The results are compared with experimental and theoretical values of RUC deformation and damage initiation and propagation under monotonic axial loading. The methodology of establishing the nonlinear visco-elastic multi-scale model of 3D textile composites without introducing the real fabric architecture in finite element analyses is explained. With the multi-scale RUCs model, the mechanical behaviors of other kinds of 3D textile composites can also be predicted.     

苯乙烯基形状记忆聚合物热力学行为的有限元分析

对苯乙烯基形状记忆聚合物进行了拉伸实验研究,测定了该材料在25℃、30℃、40℃和50℃时的弹性模量和屈服极限.根据实验结果,建立了苯乙烯基形状记忆聚合物的材料参数方程,描述了苯乙烯基形状记忆聚合物在玻璃体转化过程中,材料参数和温度的关系.在假设形状记忆聚合物为各向同性材料的基础上,将Tobushi等建立的热力学本构方程从一维扩展到三维.基于有限元分析软件ABAQUS的二次开发功能,针对上述本构方程和材料参数方程,编写了可供ABAQUS调用的UMAT函数,并对苯乙烯基形状记忆聚合物实现形状记忆效应的高温变形、应力冻结和形状恢复等热力学过程,进行了有限元数值模拟分析.     

连续胶体材料化学-力学耦合有限元理论及其在ABAQUS中UEL程序实现

基于哈密顿原理,得到水凝胶的化学-力学耦合控制方程的等效积分形式和有限元形式。在整体坐标系下推导出用形函数表示的化学-力学耦合应变矩阵和单元刚度矩阵,并且得到在局部坐标系下的离散化形式。结合ABAQUS软件,编制了用户单元子程序UEL,通过数值算例验证了所开发单元的正确性,为在ABAQUS软件中实现各种耦合问题的有限元UEL编程提供了参考依据。     

基于刚度损伤并考虑强度退化的钢筋模型

Clough本构广泛应用于现浇节点的滞回模拟中,针对原始Clough本构中的应力可能会超出材料自身极限强度的问题,提出了一种根据节点刚度损伤来推导钢筋本构的算法。该算法依次建立现浇节点刚度退化与构件中钢筋模量退化间的关系及梁端位移角与钢筋应变的关系,根据以上关系建立钢筋强度退化的表达式,并得到钢筋的本构模型。使用Fortran语言编写自定义纤维单元材料子程序,利用有限元软件ABAQUS进行两个算例验证,得到的模拟滞回曲线具有捏缩现象,形状与试验吻合良好,耗能和延性等方面可与试验接近,在峰值点位移角的预测方面也具备一定的优势,同时,本构在二维和三维模型中均有良好的收敛性,且精度较高,可以应用于实际工程中。     

Dynamic microprogramming and external subroutine calls in a multics-type environment

A scheme is presented for external subroutine calls which makes it possible to have microprogrammed subroutines and to allow subroutines written in different languages. The languages may be compiled to different Directly Executed Languages (DELs) with corresponding microprogrammed interpreters. The scheme makes it possible to intercept each external subroutine call so that such actions as parameter validation, microprogram simulation, and statistics collection can be done. The scheme is intended to be working in a Multics-type environment with late dynamic linking, virtual memory, and ring protection. The scheme is adapted to a model of microprogram development which is also presented.     

蒙特卡罗与离散纵标耦合方法在压水堆堆腔漏束计算中的应用

为克服蒙特卡罗(MC)方法计算时间长和离散纵标(SN)方法复杂几何描述不精确的困难,采用SN-MC耦合计算流程,研究了基于蒙特卡罗方法和离散纵标方法的耦合方法。耦合方法的主要思想是根据离散纵标程序提供的中子角通量,利用接口程序计算出面源的分布概率,然后由修改后的源抽样子程序生成包括上下圆面源和圆柱面源的组合源,提供给蒙特卡罗程序进行计算分析。初步计算结果表明,该耦合方法是正确的,可用于压水堆堆腔漏束的计算分析。     

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

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