Rate‐type Model for electro‐rheological Material Behaviour consistent with Extended Thermodynamics: Application to a steady viscometric flow

Based on extended thermodynamics, a class of rate‐type constitutive equations for electrorheological materials is presented. It is shown unlike many recent works using standard methods of continuum mechanics, that this approach, which consists in postulating expressions which relate the time derivative of the Cauchy stress to the motion, takes well account of the viscoelastic‐plastic effects as well as the coupling between the motion and the electric field. By considering special cases it turns out that the derived nonlinear algebraic model embeds a variety of existing differential models in the literature. The accuracy of the model is illustrated by considering a steady shear flow in which the electric field is applied perpendicularly to the flow direction. Analytical solutions with various parameter variations are obtained and discussed. The model predicts well the behaviour associated with material instabilities observed in viscometric flows for which the existence of the yield stress is required as confirmed by many experiments.     

爆炸作用下钢板层裂的数值分析

用一维运动模型对高能炸药在钢板表面接触爆炸时钢板中应力波的传播及其在钢板自由面上反射后引起的层裂现象进行了数值分析,并和前几年在我国进行的实验结果作了比较.发现当钢板采用流体弹塑性模型,并采用损伤积累层裂准则的情况下,计算所得的主裂片厚度是和实验结果合理地符合的.提出了一个计算主裂片厚度的近似公式,对实验中出现的多层的呈云母状结构的次裂片现象也作出了较满意的解释.     

Analysis of laminated piezoelectric composite plates using an inverse hyperbolic coupled plate theory

This paper presents a non-polynomial coupled plate theory for smart composite structures employing inverse hyperbolic displacement and electric potential functions. The theory is utilized towards analysis of composite piezoelectric plates operating in sensor and actuator modes. Particularly, the following three cases are studied: (i) passive laminated composite structure, (ii) composite piezoelectric plate actuator and (iii) unimorph and bimorph piezoelectric plate sensors. Analytical solutions are obtained for simply supported plates under static electrical and mechanical loads. These results are validated with existing 3D elasticity solutions and compared with other plate theory solutions. Furthermore, parametric studies are performed to determine the effect of loading, span-to-thickness ratio and lamination sequence on the response of the piezoelectric plate. Finally, the theory is applied to a transverse shear sensing device which utilizes transverse shear-electric field coupling in piezoelectric materials. This effect is often ignored in literature.It is observed that the maximum percentage error of the present theory, when compared with 3D results, is less than 3%, which is lower than other higher order plate theories.     

无粘性土的应力矢量本构模型(Ⅰ)——理论

在充分考虑应力的矢量特性基础上,通过将应力矢量的作用效应分解为球应力作用效应与应力比矢量 (偏应力矢量与球应力之比)作用效应的叠加,建立了一个全新的、适用于无粘性土在平面应变和三维条件下的非线性本构模型。该模型可以同时考虑应力的数量和方向变化对变形的影响,既适用于单调静荷作用,也适用于往返动荷作用。     

Characterisation and Modelling of Piezoceramic Actuator Hystereses

Piezoelectric ceramics are often used as actuators in smart structures technology. In the vast majority of papers dealing with this topic only linear constitutive relations are used. However, the electric field-strain relations of such actuators show hysteretic behaviour, which means that the piezoelectric coupling coefficient is not constant. In this study the hysteresis of a mechanically unconstrained actuator is obtained using the Michelson interferometry. The hysteretic behaviour is modelled by a Preisach model. Using these experimental data, for the modelling of an active structure with embedded piezoelectric actuators the actual coupling coefficient can then be determined with the help of the Preisach model. With this procedure the actuation strain of an embedded actuator, including the physical nonlinearities, can be calculated using the material characteristics obtained for an unconstrained actuator. For an experimental validation of the method outlined above, a Lead Zirconate Titanate (PZT) actuator is characterised experimentally and then glued to a cantilever beam. Then, the tip displacement of the actuated beam is determined experimentally and simulated numerically using the above method. The experimental and numerical results agree reasonably well if the shear lag due to the bonding layer between the actuator and the structure is taken into consideration. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Piezoelectric Finite Shell Element for the Geometrically Nonlinear Analysis of Sensors

A geometrically nonlinear finite element formulation to analyze piezoelectric shell structures is presented. The formulation is based on the mixed field variational functional of Hu–Washizu. Within this variational principle the independent fields are displacements, electric potential, strains, electric field, stresses and dielectric displacements. The mixed formulation allows an interpolation of the strains and the electric field through the shell thickness, which is an essential advantage when using a three dimensional material law. It is remarked that no simplification regarding the constitutive relation is assumed. The normal zero stress condition and the normal zero dielectric displacement condition are enforced by the independent resultant stress and resultant dielectric displacement fields. The shell structure is modeled by a reference surface with a four node element. Each node possesses six mechanical degrees of freedom, three displacements and three rotations, and one electrical degree of freedom, which is the difference of the electric potential through the shell thickness. The developed mixed hybrid shell element fulfills the in–plane, bending and shear patch tests, which have been adopted for coupled field problems. A numerical investigation of a smart antenna demonstrates the applicability of the piezoelectric shell element under the consideration of geometrical nonlinearity. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A solid shell finite element formulation to simulate the behaviour of thin dielectric elastomer structures

To analyse the behaviour of thin structures of dielectric elastomer (DE) material a solid shell finite element is presented. The main characteristics of DEs are a non-linear hyper elastic behaviour, the quasi-incompressibility, and the ability to transform electric energy into mechanical work. Applying a voltage to thin DE structures may produce large elongation strains of 120-380%. These large strains, the efficient electro-mechanical coupling, and the light weight make DEs very attractive for the usage in actuators. Thus, there is a need for detailed research. With respect to the electro-mechanical coupling a constitutive model is presented. An electric stress tensor and a total stress tensor are introduced by considering the electrical body force and couple in the balance of linear momentum and angular momentum, respectively. The governing equations are derived and embedded in the solid shell formulation. The element formulation is based on a Hu-Washizu mixed variational principle using six independent fields: displacements, electric potential, strains, electric field, mechanical stresses, and dielectric displacements. It allows large deformations and accounts for physical nonlinearities to capture two of the main characteristics of DEs. The shell element could be applied for the modelling of arbitrary curved thin structures. The ability of the present element formulation is demonstrated in several examples. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Bending of sinusoidal functionally graded piezoelectric plate under an in-plane magnetic field

This work investigates the bending of a simply supported functionally graded piezoelectric plate under an in-plane magnetic field. The extended sinusoidal plate theory for piezoelectric plate is adopted. The governing equations are derived by the principle of the virtual work considering the Lorentz magnetic force obtained from the Maxwell's relation. The effect of magnetic field, electric loading and gradient index on the displacement, electric potential, stress and electric displacement are numerically presented and discussed in detail. These conclusions will be of particular interest to the future analysis of piezoelectric plate in magnetic field.     

CFRP修复缺陷钢板应力解析模型

在使用碳纤维复合材料(carbon fiber reinforced polymer, CFRP)修复钢结构腐蚀缺陷的修复方式中,CFRP应力及胶层应力是确定碳纤维修复结构承载能力的关键。基于平截面假设,得到弯矩作用下应力与应变分布;基于胶层剪切模型,得到胶层剪应力与CFRP和钢板位移间的关系;基于力的平衡,得到CFRP和钢板的应力关系。结合得到的各种材料之间关系,推导出轴力和弯矩联合作用状态下CFRP双面修复钢板的CFRP与胶层应力分布解析解。采用数值分析对CFRP双侧粘贴修复缺陷钢板进行分析,分析结果与解析结果具有一致性,同时获得了CFRP双侧粘贴修复缺陷钢板的应力分布特点,以及构件可能发生破坏的位置,为计算构件极限承载力提供了基础。     

Micro-mechanical analyses of the effect of stress triaxiality and the Lode parameter on ductile damage and failure

The paper deals with the effect of different stress states on plastic deformations, damage and fracture of ductile materials. To be able to model these effects a continuum damage model has been introduced taking into account the dependence of stress-state on the constitutive equations. The model is based on the introduction of damaged and fictitious undamaged configurations. All parameters appearing in the constitutive equations are stress-state-dependent which can be characterized by the stress intensity, the stress triaxiality and the Lode parameter. Only experiments are not adequate enough to determine all constitutive parameters. Thus, additional series of three-dimensional micro-mechanical simulations of representative volume elements have been performed to get more insight in the complex damage mechanisms. These simulations cover a wide range of stress triaxialities and Lode parameters in tension, shear and compression domains. After all, the results from the micro-mechanical simulations are used to suggest the damage equations and to identify corresponding parameters. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Constitutive equations of non-isothermal polymer melt

Constitutive equations of non-isothermal polymer melt are presented by the analysis of entropic free energy contribution of the macromolecular chains, which are treated as elastic dumbbell models. With describing non-isothermal dumbbell spring, as the function of temperature, the non-linear elastic coefficient expression causes the appearance of temperature gradient in stress constitutive equations. Following the constitutive equation of Hookean dumbbell model, non-isothermal stress constitutive equations of FENE and FENE-P models are derived. In deriving process of constitutive equations, the second moment approximation is used to closure FENE model. Using the non-isothermal constitutive equations, numerical simulations of polymer flow through shear cavity and planar contraction cavity are presented. And the distributions of correlative stress functions and the effects of different temperatures on stress functions are discussed. The present results are shown to explore the non-isothermal constitutive equations of elastic dumbbell models, and to search more accurately describing way of non-isothermal polymer melt.     

A 1D constitutive law for hysteresis effects in piezoceramics

This paper is concerned with a macroscopic constitutive law for domain switching effects, which occur in piezoelectric ceramics. The thermodynamical framework of the law is based on two scalar valued functions: the Helmholtz free energy and a switching surface. In common usage, the remanent polarization and the remanent strain are employed as internal variables. The novel aspect of the present work is to introduce an irreversible electric field, which serves besides the irreversible strain as internal variable. The irreversible electric field has only theoretical meaning, but leads to advantages within the finite element implementation, where displacement and the electric potential are the nodal degrees of freedoms. A common assumption is a one-to-one relation between the irreversible strain and the polarization. This simplification is not employed in the present paper. To accomplish enough space for the polarization, resulting from an applied electric field, the irreversible strains are additively split and a special hardening function is introduced. This balances the amount of space and the domain switching due to polarization. The constitutive model reproduces the ferroelastic and the ferroelectric hysteresis as well as the butterfly hysteresis for piezoelectric ceramics and it accounts for the mechanical depolarization effect. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

关于损伤张量的阶次

本文首先讨论了较为广泛的连续介质材料的应力变形本构关系,得到了通常以泛函表示的应力变形本构关系的张量表达式.以此为基础,研究了各向异性材料各向异性损伤时,无论从连续介质力学模型出发还是从缺陷模型出发,描述损伤的张量都存在最高阶次的限制;指出了在什么条件下,损伤变量可用低于最高阶次的张量来描述.     

Peridynamics model for flexoelectricity and damage

A flexoelectric peridynamic (PD) theory is proposed. In the PD framework, the formulation introduces a nanoscale flexoelectric coupling that entails non-uniform strain in centrosymmetric dielectrics. This potentially enables PD modeling of a large class of phenomena in solid dielectrics involving cracks, discontinuities etc. wherein large strain gradients are present and the classical electromechanical theory based on partial differential equations do not directly apply. PD electromechanical equations, derived from Hamilton's principle, satisfy the global balance laws. Linear PD constitutive equations reflect the electromechanical coupling effect, with the mechanical force state affected by the polarization state and the electrical force state in turn by the displacement state. An analytical solution to the PD electromechanical equations is presented for the static case when a point mechanical force and a point electric force act in an infinite 3D solid dielectric. A parametric study on how different length scales influence the response is undertaken. In addition, the model is extended to incorporate damage using phase field – an order parameter, supplemented with a PD bond breaking criterion to study flexoelectric effects in damage and fracture problems. To demonstrate the performance of our proposal, we first simulate, considering small flexoelectricity effect and no damage, an externally pressured 2D flexoelectric disk subjected to a potential difference between the inner and outer surfaces and compare the results with existing solutions in the literature. Next, we simulate a plate with a central pre-crack under tension considering damage and flexoelectricity effects, and study the effect of various constitutive parameters on the damage evolution. We also furnish a classical derivation of phase field based flexoelectricity in Appendix I.     

集成结构振动主动控制和抑制

采用一种新的压电板单元,建立了含有分布压电传感元件和执行元件的集成结构的有限元动力模型。研究了这种集成结构在常增益负速度反馈控制规律作用下,振动的主动控制与抑制的问题,并提出了集成结构振动主动控制和抑制的一般方法。最后,提供了数值示例,说明本文提出方法的有效性。     

An efficient iteration approach for nonlinear boundary value problems in 2D piezoelectric semiconductors

Based on initial nonlinear constitutive equations, we establish the extended displacement and traction boundary integral equations for a piezoelectric medium with a volume electric charge, along with electron and electric current density boundary integral equations for a conductor with a volume electric current. Then, an iterative approach is proposed for investigation of boundary value problems in two-dimensional piezoelectric semiconductors (PSCs). Compared with extended displacements obtained by finite element analysis, this approach is validated via a rectangular PSC under extended external loads. Furthermore, as a numerical example, extended displacements across an elliptical hole in a rectangular PSC are investigated. It is shown that there is a stress concentration near the elliptical hole, which is closely dependent on its shape.     

The role of magnetic fields in the strong stabilization of a hybrid magneto‐elastic structure

In this paper, we are concerned with a model for the magneto–elastic interactions of a three‐dimensional elastic body and a two‐dimensional flexible plate, which is attached to the flat flexible part of the surface of the body. Both the solid body and the plate are permeated by magnetic fields. The mathematical model is analyzed from the point of view of existence and uniqueness and stabilization.It turns out that, in the presence of the magnetic fields in the solid and the plate, strong stabilization can be achieved under viscous damping in the plate in one direction that is determined by the nature of the primary magnetic fields in the body and the plate. Copyright © 2011 John Wiley & Sons, Ltd.     

3D nonlinear variable strain-rate-dependent-order fractional thermoviscoelastic dynamic stress investigation and vibration of thick transversely graded rotating annular plates/discs

In the present article, the idea of using the variable-order fractional-derivative thermoviscoelastic constitutive laws in dynamic stress and vibration analysis of the engineering structures, the required implementation backgrounds, and the relevant numerical solution procedures are investigated for the first time. In this regard, dynamic 3D stress and displacement fields and radial/transverse vibrations of transversely graded viscoelastic spinning thick plates/discs exposed to sudden thermoelastic loads are investigated. Instead of using the approximate plate theories, the exact thermoviscoelasticity theory is employed in the development of the governing equations. Since the variable fractional order is dependent on the localized deformation rates, the resulting thermoviscoelastic integro-differential equations are nonlinear. These equations are solved by utilizing a combination of the second-order backward/central/forward finite difference discretization of the spatial and time domains, numerical evaluation and updating of the Caputo-type fractional derivatives, updating the growing number of terms of the governing equations, and Picard's iterations. Various edge conditions are considered. Finally, comprehensive sensitivity analyses and various 3D plots are presented and discussed regarding the effects of the variable fractional order of the constitutive law, time variations of the nonuniformly distributed transverse loads, and edge conditions on the distributions and damping of the resulting displacement and stress components.     

Flow-induced anisotropic viscosity in short fiber reinforced polymers

The commonly used flow models for fiber reinforced polymers often neglect the flow induced mechanical anisotropy of the suspension. With an increasing fiber volume fraction, this plays, however, an important role. There are some models which count on this effect, they are, however, phenomenological and require a fitted model parameter. In this paper, a micromechanically based constitutive law is proposed which considers the flow induced anisotropic viscosity of the fiber suspension. The introduced viscosity tensor can handle arbitrary anisotropy of the fluid-fiber mixture depending on the actual fiber orientation distribution. A homogenization method for unidirectional structures in contribution with orientation averaging is used to determine the effective viscosity tensor. The motion of rigid ellipsoidal fibers induced by the flow of the matrix material is described by Jeffery's equation. A numerical implementation of the introduced model is applied to representative flow modes. The calculated stress values are analyzed in transient and stationary flow cases. They show a less pronounced anisotropic viscous behaviour in every investigated case compared to the results obtained by the use of the Dinh-Armstrong constitutive law. The reason for the qualitative difference is that the presented model depends on the complete orientation information, while the other one is linear in the fourth-order fiber orientation tensor. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)