An Approach to the Modeling of Physical Ageing in Rubbery Polymers

The oxidation reactions responsible for physical ageing are thermally activated processes and yield cross-linking similar to the vulcanization process and chain-scission reactions. Chain scission, is responsible for the relaxation behaviour under constant stretch, whereas cross-linking is responsible for hardening of the material. The decrease in the stretch at break upon cross-linking motivates the so called network alteration proposed in this contribution. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Approach to Multiphase Flow Modeling

We review the phase field (otherwise called diffuse interface) model for fluid flows, where all quantities, such as density and composition, are assumed to vary continuously in space. This approach is the natural extension of van der Waals?? theory of critical phenomena both for one-component, two-phase fluids and for partially miscible liquid mixtures. The equations of motion are derived, assuming a simple expression for the pairwise interaction potential. In particular, we see that a non-equilibrium, reversible body force appears in the Navier-Stokes equation, that is proportional to the gradient of the generalized chemical potential. This, so called Korteweg, force is responsible for the convection that is observed in otherwise quiescent systems during phase change. In addition, in binary mixtures, the diffusive flux is modeled using a Cahn-Hilliard constitutive law with a composition-dependent diffusivity, showing that it reduces to Fick??s law in the dilute limit case. Finally, the results of several numerical simulations are described, modeling, in particular, a) mixing, b) spinodal decomposition, c) nucleation, d) enhanced heat transport, e) liquid-vapor phase separation.     

A Regularized Sharp Interface Model for Martensitic Phase Transformations at Large Strains

The macroscopic mechanical behavior of many functional materials crucially depends on the formation and evolution of their microstructure. When considering martensitic shape memory alloys, this microstructure typically consists of laminates with coherent twin boundaries. We suggest a variational-based phase field model for the dissipative evolution of microstructure with coherence-dependent interface energy and construct a suitable gradient-extended incremental variational framework for the proposed dissipative material. We use our model to predict laminate microstructure in martensitic CuAlNi. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Modeling of Brittle and Ductile Fracture

The phase field modeling of brittle fracture was a topic of intense research in the last few years and is now well-established. We refer to the work [1-3], where a thermodynamically consistent framework was developed. The main advantage is that the phase-field-type diffusive crack approach is a smooth continuum formulation which avoids the modeling of discontinuities and can be implemented in a straightforward manner by multi-field finite element methods. Therefore complex crack patterns including branching can be resolved easily. In this paper, we extend the recently outlined phase field model of brittle crack propagation [1-3] towards the analysis of ductile fracture in elastic-plastic solids. In particular, we propose a formulation that is able to predict the brittle-to-ductile failure mode transition under dynamic loading that was first observed in experiments by Kalthoff and Winkler [4]. To this end, we outline a new thermodynamically consistent framework for phase field models of crack propagation in ductile elastic-plastic solids under dynamic loading, develop an incremental variational principle and consider its robust numerical implementation by a multi-field finite element method. The performance of the proposed phase field formulation of fracture is demonstrated by means of the numerical simulation of the classical Kalthoff-Winkler experiment that shows the dynamic failure mode transition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formulation and Numerical Implementation of Ferroelectrics at Large Strains

In recent years, an increasing interest has been shown in functional materials such as ferroelectric polymers. For such materials, viscous effects and electric polarizations cause hysteresis phenomena accompanied with possibly large remanent strains and rotations. Ferroelectric polymers have many attractive characteristics. They are light, inexpensive, fracture tolerant, and pliable. Furthermore, they can be manufactured into almost any conceivable shape and their properties can be tailored to suit a broad range of requirements. In this work, continuous and discrete variational formulations are exploited for the treatment of the non-linear dissipative response of ferroelectric polymers under electrical loading. The point of departure is a general internal variable formulation that determines the hysteretic response of this class of materials in terms of an energy storage and a rate-dependent dissipation function. The ferroelectric constitutive assumptions, which account for specific problems arising in the geometric nonlinear setting, are discussed. With regard to the choice of the internal variables, a critical factor is the kinematic assumptions. Here, we propose the multiplicative decomposition of the local deformation gradient into reversible and remanent parts, where the latter is characterized by a metric tensor. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase Field Modeling of Fracture in Plates and Shells

The numerical modeling of failure mechanisms in plates and shells due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by a diffusive crack modeling, which is based on the introduction of a crack phase field. In this paper, we extend ideas recently outlined in [1, 2] towards the phase field modeling of fracture in dimension-reduced continua with application to Kirchhoff plates and shells. The introduction of history fields, containing the maximum reference energy obtained in history, provides a very transparent representation of the coupled balance equations and allows the construction of an extremely robust operator split technique. The performance of the proposed models is demonstrated by means of representative numerical examples. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

纯弯曲矩形截面梁Ⅰ型单边裂纹端部的应力应变场及裂纹失稳扩展临界应力的计算

本文应用文[1]的分析方法,研究了纯弯曲矩形载面梁Ⅰ型单边裂纹端部的应力应变场,给出了裂纹尖端的应力应变分量和计算裂纹端部弹性变形区和变形强化区宽度的公式以及计算裂纹失稳扩展临界应力的方程组。最后用计算实例对裂纹失稳扩展临界应力方程组进行了验证,最大误差不超过0.18%.     

线性硬化材料中稳恒扩展裂纹尖端场的粘塑性解

采用弹粘塑性力学模型,对线性硬化材料中平面应变扩展裂纹尖端场进行了渐近分析．假设人工粘性系数与等效塑性应变率的幂次成反比,通过量级匹配表明应力和应变均具有幂奇异性,奇异性指数由粘性系数中等效塑性应变率的幂指数唯一确定．通过数值计算讨论了Ⅱ型动态扩展裂纹尖端场的分区构造随各材料参数的变化规律．结果表明裂尖场构造由硬化系数所控制而与粘性系数基本无关．弱硬化材料的二次塑性区可以忽略,而较强硬化材料的二次塑性区和二次弹性区对裂尖场均有重要影响．当裂纹扩展速度趋于零时,动态解趋于相应的准静态解;当硬化系数为零时便退化为HR(Hui-Riedel)解．     

Towards Phase Field Modeling of Ductile Fracture in Gradient-Extended Elastic-Plastic Solids

The modeling of failure in ductile metals must account for complex phenomena at a micro-scale as well as the final rupture at the macro-scale. Within a top-down viewpoint, this can be achieved by the combination of a micro-structure-informed elastic-plastic model with a concept for the modeling of macroscopic crack discontinuities. In this context, it is important to account for material length scales and thermo-mechanical coupling effects due to dissipative heating. This can be achieved by the construction of non-standard, gradient-enhanced models of plasticity with a full embedding into continuum thermodynamics [1,2]. The modeling of macroscopic cracks can be achieved in a convenient way by recently developed continuum phase field approaches to fracture based on regularized crack discontinuities. This avoids the use of complex discretization methods for crack discontinuities, and can account for complex crack patterns within a pure continuum formulation. Moreover, the phase field modeling of fracture is related to gradient theories of continuum damage mechanics, and fits nicely the structure of constitutive models for gradient plasticity. The main focus of this work is the extensions to gradient thermoplasticity and phase field formulation of ductile fracture, conceptually in line with the work [3]. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On The Formulatio and Numerical Implementation of Dissipative Electro-Mechanics at Large Strains

In the last year an increasing interest in functional materials such as ferroelectric polymers and ceramics has been shown. For those materials viscous effects or electric polarizations cause hysteresis phenomena accompanied with possibly large remanent strains. This paper outlines aspects of the formulation and numerical implementation of dissipative electro-mechanics at large strains. In the first part, we focus on the geometric nature of dissipative electro-mechanics. In a second part, we discuss constitutive assumptions which account for specific problems arising in the geometric nonlinear setting. This concerns the definition of objective energy storage and dissipation functions with suitable symmetries and (weak) convexity properties. With regard to the choice of the internal variables entering these functions, a critical point are kinematic assumptions. Here, we investigate the multiplicative decomposition of the local deformation gradient into reversible and remanent parts as well as the introduction of a remanent metric. In a third part, we summarize details of the constitutive updates as well as the finite element formulations, both on the basis of compact incremental variational principles. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

磁场对液体-微伸缩固体介面波传播的影响

在一个传播理想的非粘性液体半空间,和一个传播理想的微伸缩弹性固体半空间之间,研究介面处纵波的反射和透射.在两个半空间中,满足介面处必需的边界条件下,得到控制方程的适当解,是一组以不同反射和透射波振幅比表示的5个非齐次方程.以水和铝-环氧树脂合成材料介面为实际例子,用Gauss消除法的Fortran程序求解方程组.考虑存在和不存在外加横向磁场两种情况,在某些入射角范围内,计算振幅比的数值解.最后用图形给出横向磁场对不同反射和透射波振幅比的影响.     

Variational Phase Field Modeling of Fracture Caused by Fluid Transport in Poroelastic Solids

The numerical modeling of failure mechanisms due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by recently developed diffusive crack modeling concepts, which are based on the introduction of a crack phase field. Such an approach is conceptually in line with gradient-extended continuum damage models which include internal length scales. In this paper, we extend our recently outlined mechanical framework [1–3] towards the phase field modeling of fracture in the coupled problem of fluid transport in deforming porous media. Here, extremely complex crack patterns may occur due to drying or hydraulic induced fracture, the so called fracking. We develop new variational potentials for Biot-type fluid transport in porous media at finite deformations coupled with phase field fracture. It is shown, that this complex coupled multi-field problem is related to an intrinsic mixed variational principle for the evolution problem. This principle determines the rates of deformation, fracture phase field and fluid content along with the fluid potential. We develop a robust computational implementation of the coupled problem based on the potentials mentioned above and demonstrate its performance by the numerical simulation of complex fracture patterns. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Variational Phase Field Modeling of Laminate Deformation Microstructure in Finite Gradient Crystal Plasticity

The macroscopic mechanical behavior of many materials crucially depends on the formation and evolution of their microstructure. In this work, we consider the formation and evolution of laminate deformation microstructure in plasticity. Inspired by work on the variational modeling of phase transformation [5] and building on related work on multislip gradient crystal plasticity [9], we present a new finite strain model for the formation and evolution of laminate deformation microstructure in double slip gradient crystal plasticity. Basic ingredients of our model are a nonconvex hardening potential and two gradient terms accounting for geometrically necessary dislocations (GNDs) by use of the dislocation density tensor and regularizing the sharp interfaces between different kinematically coherent plastic slip states. The plastic evolution is described by means of a nonsmooth dissipation potential for which we propose a new regularization. We formulate a continuous gradient-extended rate-variational framework and discretize it in time to obtain an incremental-variational formulation. Discretization in space yields a finite element formulation which is used to demonstrate the capability of our model to predict the formation and evolution of laminate deformation microstructure in f.c.c. Copper with two active slip systems in the same slip plane. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

守恒相场模型弱解的Blow—up

本文讨论一类守恒相场模型弱解的性态,证明当ａ２ｐ－１＜０及初始数据充分大时解在有限时刻Ｂｌｏｗ－ｕｐ     

泊松比对静止平面应变裂纹尖端的理想塑性应力场的影响

在裂纹尖端的理想塑性应力分量都只是θ的函数的条件下,利用平衡方程和含有泊松比的Mises屈服条件,本文导出了静止平面应变裂纹尖端的理想塑性应力场的一般解析表达式.将这些一般解析表达式用于具体裂纹,我们就可以得到静止平面应变Ⅰ型、Ⅱ型及Ⅰ-Ⅱ复合型裂纹尖端的理想塑性应力场的解析表达式,这些表达式含有泊松比.     

长程选举模型的平均场极限

本文研究长程选举模型的平均场极限，利用对偶关系和特征函数方法证得长程选举模型的平均场极限满足下列微分方程：     

含埋藏椭圆形裂纹的金属构件脉冲放电瞬间温度场分析

对带有椭圆埋藏裂纹金属构件在脉冲放电瞬间的温度场进行了理论分析．在带有裂纹问题温度场的理论求解中,根据相似性原理,将电流通过带有椭圆面裂纹导体类比于流体流过障碍物的情形,从而引入了裂纹的边界条件,得到域内电流密度的分布,进而得到了放电瞬间埋藏椭圆裂纹尖端的温度场表达式,为空间裂纹电磁热止裂技术的实际应用提供了理论基础．     

Error Analysis of SAV Finite Element Method to Phase Field Crystal Model

In this paper, we construct and analyze an energy stable scheme by combining the latest developed scalar auxiliary variable (SAV) approach and linear finite element method (FEM) for phase field crystal (PFC) model, and show rigorously that the scheme is first-order in time and second-order in space for the $L^2$ and $H^{-1}$ gradient flow equations. To reduce efficiently computational cost and capture accurately the phase interface, we give a simple adaptive strategy, equipped with a posteriori gradient estimator, i.e., $L^2$ norm of the recovered gradient. Extensive numerical experiments are presented to verify our theoretical results and to demonstrate the effectiveness and accuracy of our proposed method.