一类铁电模型的参变量变分原理及其应用

将参变量变分原理引入铁电问题。对一类借用了经典弹塑性理论中的概念和方法的多轴铁电模型建立基于Helmholtz自由能的参变量变分原理,可以有效处理传统变分原理中由非关联流动法则或屈服面不考虑材料系数变化所引起的切线模量非对称困难。相应于参变量变分原理,引入参数二次规划算法,可获得具有可靠数值稳定性的一套铁电算法。将该算法应用于一个具体的铁电模型,数值计算结果表明本文方法的有效性。     

解Drucker-Prager塑性问题的二阶锥互补法

基于经典弹塑性理论中多数屈服准则具有凸锥数学结构的事实,将在大规模计算中更具潜力的锥规划法引入弹塑性分析。考虑到弹塑性流动理论有关联与非关联之分,本文提出利用锥型互补法求解弹塑性问题。具体以Drucker-Prager弹塑性模型为例,首先利用最大塑性功耗散原理和圆锥对偶理论等工具,建立了弹塑性本构方程的等价二阶锥互补模型;然后,基于参变量变分原理和有限元技术,建立了弹塑性增量分析的二阶锥线性互补模型;最后,利用一类半光滑Newton算法求解。数值算例表明了本文方法的有效性。     

参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

Cosserat连续体模型中本构参数对应变局部化模拟结果影响的数值分析

基于所发展的压力相关弹塑性Cosserat连续体模型及相应的数值方法,以一维剪切层及二维平板压缩问题为例,数值分析了Cosserat连续体模型中的本构参数Cosserat剪模、软化模量及内部长度参数对应变局部化数值模拟结果的影响.结果表明在一定取值范围内,Cosserat剪模对数值模拟结果几乎没有影响,并给出了具体数值计算时的取值范围;软化模量绝对值越大,后破坏段的荷载-位移曲线越陡,计算得到的剪切带宽度越窄;内部长度参数越大,后破坏段的荷载-位移曲线越平缓,计算得到的剪切带越宽.     

基于Cosserat理论的小孔应力集中问题的有限元分析

近来由于细观力学的发展和对材料尺度效应的研究使得Cosserat理论受到了关注.本文给出了基于Cosserat理论的有限元八节点等参元格式,数值计算了单向均匀拉伸小孔应力集中问题,研究了Cosserat理论中有关参数对应力集中因子和尺度效应的影响.计算结果与理论解十分吻合.这表明,基于Cosserat理论的平面八节点等参元适用于求解基于Cosserat理论的平面问题.     

参数二次规划法在计算力学中的应用(一)

绪言 一系列论文发表了关于参数变分原理在弹性接触问题和弹塑性分析中的理论和应用。所提出的参数最小势能原理和参数最小余能原理的特点是构造的能量泛函中含有两类变量:一类是状态变量,它们和经典变分原理中的一样;另一类是控制变量,它们不参加变分,但控制着变分过程,使问题的非线性本构特性(或状态方程)得以满足。参数变分原理     

改进的Prandtl—Reuss理论与有限元杂交／混合变刚度法

1 引言常规弹塑性有限元法基于传统的流动法则与Prandtl-Reuss 理论.由位能原理导出的某些迭代算法在计算效率和收敛速度方面是偏于保守的.再者,由于Mises 材料的不可压缩性,它的应力不能由变形唯一地确定。导致常规算法数值求解的极大困难.采用某些特殊的修正技术或引用广义变分原理的协调模型及杂交模型是克服这种数值奇性的有效途径.在本文中,作者改进了传统的Prandtl-Reuss 理论,提出等向强化Mises 材料的拟流动准则及应变增量假设的新想法.并引入文[7]修正余能原理的杂交/混合模型.由此建立     

基于参变量变分原理的非均质多边形微结构材料弹塑性分析

为了更好地模拟复合材料及含夹杂非均质材料等的宏观弹塑性力学性能,简化有限元建模时间和减少有限元模拟计算量。本文基于参变量变分原理,提出了一种采用任意多边形弹塑性单元进行结构非线性分析的参数二次规划算法,给出了参变量最小势能原理以及最终的二次规划模型,并在有限元分析与优化设计软件系统JIFEX上进行了程序实现。数值算例证明了本文方法的正确与可行性。     

基于界带模型的碳纳米管声子谱的辛分析

针对碳纳米管声子谱的数值计算方法研究,基于对偶体系和辛几何算法提出了一套全新的计算方法和相应的界带结构模型,通过将碳纳米管模拟成不同的结构力学模型,利用分析结构力学中的振动理论来计算碳纳米管的色散关系.理论框架包括:周期结构的变分原理、周期结构中波的传播分析、子结构方法、界带理论和声子色散关系的基本算法.数值算例验证了理论和算法的有效性,而且也指出了针对碳纳米管的声子谱的计算,界带模型相对于其它传统模型存在着一定的优势.     

一种新的弹塑性杂交应力元法

本文基于一个改进的弹塑性的Hellinger/Reis■ner 混合变分原理构造了一种用于解弹塑性问题的四节点等参杂交应力元.新的模型中,在单元内增加了等效应力增量、塑性等效应变增量及不协调位移变量,从而使单元内的屈服准则及流动法则平均得到满足,不协调位移改进了单元应力精度.计算表明,新的模型可以提高弹塑性杂交法的精度和计算效率.     

PARAMETRIC VARIATIONAL PRINCIPLE BASED ELASTIC-PLASTIC ANALYSIS OF COSSERAT CONTINUUM

A new algorithm is developed based on the parametric variational principle for elastic-plastic analysis of Cosserat continuum. The governing equations of the classic elastic-plastic problem are regularized by adding rotational degrees of freedom to the conventional translational degrees of freedom in conventional continuum mechanics. The parametric potential energy principle of the Cosserat theory is developed, from which the finite element formulation of the Cosserat theory and the corresponding parametric quadratic programming model are constructed. Strain localization problems are computed and the mesh independent results are obtained.     

平面Cosserat模型有限元分析的4和8节点单元与分片试验研究

经典连续体理论不包括物质内部尺度,当考虑应变软化问题时,有限元结果对网格具有很强的依赖性。与经典连续介质力学理论不同,Cosserat连续体模型在传统平动自由度的基础上添加了一独立的旋转自由度,在本构模型中引入了内尺度参数。本文研究了基于Cosserat理论的平面4和8节点等参元以及8（4）节点线、角位移混合插值等参单元,给出Cosserat单元分片试验的实施过程。最后将单元运用到小孔应力集中问题的分析当中,通过计算结果与理论解的比较,表明了4和8节点以及8（4）节点等参元的适用性,为问题的非线性分析打下基础。     

三维Cosserat连续体模型与微结构尺寸相关效应的有限元分析

分析了三维Cosserat连续体理论中的应力应变特征,推导了三维Cosserat连续体的有限元方程,基于ABAQUS计算软件提供的用户单元子程序(UEL)接口编写了弹性Cosserat连续体三维20节点有限元程序,并分析了微悬臂梁自由端的挠度问题和微杆扭转问题。通过与基于经典连续体理论的解析解及有限元数值计算结果进行比较,表明所发展的三维Cosserat连续体有限元能有效地模拟微结构尺寸相关效应问题,即随着微结构尺寸与材料内部长度参数的接近,基于Cosserat连续体有限元分析得到的微梁的挠度以及微杆的转角与经典连续体的解析解及有限元解相比越来越小;反之,Cosserat连续体有限元的计算结果与经典连续体的解析解及有限元数值解较为一致。     

Computational multiscale methods for granular materials

The fine-scale heterogeneity of granular material is characterized by its polydisperse microstructure with randomness and no periodicity. To predict the mechanical response of the material as the microstructure evolves, it is demonstrated to develop computational multiscale methods using discrete particle assembly-Cosserat continuum modeling in micro- and macro- scales, respectively. The computational homogenization method and the bridge scale method along the concurrent scale linking approach are briefly introduced. Based on the weak form of the Hu-Washizu variational principle, the mixed finite element procedure of gradient Cosserat continuum in the frame of the second-order homogenization scheme is developed. The meso-mechanically informed anisotropic damage of effective Cosserat continuum is characterized and identified and the microscopic mechanisms of macroscopic damage phenomenon are revealed.     

梯度增强Cosserat连续体的广义Hill定理

基于经典Cauchy连续体的Hill定理,在平均场理论的框架下导出了梯度增强Cosserat连续体细、宏观均匀化方法的广义Hill定理。在梯度增强Cosserat连续体中,不仅宏观样条点上的应变和应力张量,而且它们的梯度均作用于与该样条点相关联的细观表征元（RVE）。依据此广义Hill定理,对梯度增强Cosserat连...     

Transversely isotropic material: nonlinear Cosserat versus classical approach

We consider a specific case of unidirectional reinforced material under applied tensile load. The reinforcement of the material is inclined with 45° to the direction of the tensile resultant. Different approaches are discussed: one experiment and three computational models. Two models use the classical Cauchy continuum theory whereas the third computational model is based on a Cosserat continuum. It is well known that test specimen with inclination between unidirectional reinforcement and tensile direction show, besides Poissons effect, additional deformation perpendicular to the load direction. The classical transversely isotropic continuum theory predicts this deformation as typical S-shape. In the Cosserat continuum the orientation of the inner structure is incorporated. Thus, structural parameters influence the deformation. With the proposed geometrically non-linear Cosserat model classical and non-classical behaviour can be modelled. In the non-classical case, the transverse deformation is not described by one S-shape but by multiple S-shaped modes. The additional rotational parameters in the Cosserat continuum are responsible for the non-classical behaviour which is due to non-symmetric strain.     

Finite element modeling of micropolar-based phononic crystals

The performance of a Cosserat/micropolar solid as a numerical vehicle to represent dispersive media is explored. The study is conducted using the finite element method with emphasis on Hermiticity, positive definiteness, principle of virtual work and Bloch–Floquet boundary conditions. The periodic boundary conditions are given for both translational and rotational degrees of freedom and for the associated force- and couple-traction vectors. Results in terms of band structures for different material cells and mechanical parameters are provided.     

Asymptotic expansion homogenization of discrete fine-scale models with rotational degrees of freedom for the simulation of quasi-brittle materials

Discrete fine-scale models, in the form of either particle or lattice models, have been formulated successfully to simulate the behavior of quasi-brittle materials whose mechanical behavior is inherently connected to fracture processes occurring in the internal heterogeneous structure. These models tend to be intensive from the computational point of view as they adopt an “a priori” discretization anchored to the major material heterogeneities (e.g. grains in particulate materials and aggregate pieces in cementitious composites) and this hampers their use in the numerical simulations of large systems. In this work, this problem is addressed by formulating a general multiple scale computational framework based on classical asymptotic analysis and that (1) is applicable to any discrete model with rotational degrees of freedom; and (2) gives rise to an equivalent Cosserat continuum. The developed theory is applied to the upscaling of the Lattice Discrete Particle Model (LDPM), a recently formulated discrete model for concrete and other quasi-brittle materials, and the properties of the homogenized model are analyzed thoroughly in both the elastic and the inelastic regime. The analysis shows that the homogenized micropolar elastic properties are size-dependent, and they are functions of the RVE size and the size of the material heterogeneity. Furthermore, the analysis of the homogenized inelastic behavior highlights issues associated with the homogenization of fine-scale models featuring strain-softening and the related damage localization. Finally, nonlinear simulations of the RVE behavior subject to curvature components causing bending and torsional effects demonstrate, contrarily to typical Cosserat formulations, a significant coupling between the homogenized stress–strain and couple-curvature constitutive equations.