磁电复合材料中拓扑磁结构的力学调控

磁性斯格明子是在一些铁磁材料中存在的一种重要拓扑磁结构,由于其具有独特的磁-电-力-热多场耦合特性,在未来新型自旋电子器件中有着广泛的应用前景。然而,磁性斯格明子一般需要在外加磁场下才能稳定存在,极大地限制了其在自旋电子器件中的实际应用。本文基于实空间下磁电材料的相场模拟,发现铁电和铁磁复合薄膜中铁电斯格明子可以通过界面变形来稳定铁磁斯格明子。由于力电耦合效应,铁电层中铁电斯格明子的非均匀分布极化在界面产生周期性的非均匀界面变形。界面变形通过力磁耦合效应,使铁磁层中的磁性斯格明子在没有外加磁场的条件下能够稳定存在。本文的研究结果表明,基于磁电复合材料中的力-电-磁耦合效应,通过优化设计复合材料中不同组元的结构,可以实现拓扑磁结构的力学调控,从而为设计基于拓扑磁结构的新型自旋电子器件提供了新思路。     

磁电复合材料中拓扑磁结构的力学调控

磁性斯格明子是在一些铁磁材料中存在的一种重要拓扑磁结构,由于其具有独特的磁-电-力-热多场耦合特性,在未来新型自旋电子器件中有着广泛的应用前景。然而,磁性斯格明子一般需要在外加磁场下才能稳定存在,极大地限制了其在自旋电子器件中的实际应用。本文基于实空间下磁电材料的相场模拟,发现铁电和铁磁复合薄膜中铁电斯格明子可以通过界面变形来稳定铁磁斯格明子。由于力电耦合效应,铁电层中铁电斯格明子的非均匀分布极化在界面产生周期性的非均匀界面变形。界面变形通过力磁耦合效应,使铁磁层中的磁性斯格明子在没有外加磁场的条件下能够稳定存在。本文的研究结果表明,基于磁电复合材料中的力-电-磁耦合效应,通过优化设计复合材料中不同组元的结构,可以实现拓扑磁结构的力学调控,从而为设计基于拓扑磁结构的新型自旋电子器件提供了新思路。     

纵横载荷作用下功能梯度梁的弯曲和过屈曲

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(FGM)梁非线性弯曲和过屈曲问题的控制方程,并获得了该问题的精确解;据此解研究了梯度材料性质、外载荷、横向剪切变形以及边界条件等因素对功能梯度材料梁非线性力学行为的影响,分析中假设功能梯度材料性质只沿梁厚度方向,并按成分含量的幂指数函数形式变化。结果表明:纵横载荷共同作用下,功能梯度梁的弯曲构形将有无限多个;随着梯度指数的增大,梁的变形减小,临界载荷升高;随着长高比的增大,横向剪切变形的影响减小。     

铁电材料的疲劳失效行为

在航空航天、核能发电等重大装备技术领域, 作为高温传感/驱动/能量收集器件的敏感材料——铋层状结构铁电(BLSF)陶瓷在复杂载荷环境下的疲劳失效问题严重限制着器件寿命和可靠性的提高. 本文以BLSF陶瓷的应用需求为背景, 围绕铁电材料的疲劳裂纹扩展与电畴极化翻转及其相互作用机制等关键问题, 综述了铁电材料在热、力、电三种载荷及其耦合作用下疲劳失效行为的研究现状, 并根据当前铁电材料的一些新发展、新应用对其未来研究方向进行了展望, 旨在为高性能、长寿命铁电/压电器件设计提供参考.      

大变形下初始斜交异性本构方程

采用材料主轴法,建立了初始斜交异性材料在变形构形（Euler描述）下的斜交异性本构方程,以及在初始构形（Lagrange描述）下的形式。具体给出了斜交异性线弹性材料方程的显式,它在Lagrange描述下形式简洁,可方便地用于有限元计算。文中指出,在变形构形下是线弹性的材料,在Lagrange描述下其本构方程一般已成为非线性,我们称之为本构转换非线性。这种非线性在实际的有限元计算中还未引起重视。为理论简明,本构方程是对二维给出的。     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

力电耦合载荷作用下铁电陶瓷破坏行为的云纹干涉方法研究

本文介绍了如何用云纹干涉法实时地观察铁电陶瓷在力载荷和电载荷共同作用下裂尖的破坏行为．测量了三点弯试验中由电场和应力集中导致的裂尖的变形场．对变形的云纹图分析表明：当极化方向与裂纹扩展方向一致，且都与电场方向垂直，裂尖附近的正应变随电场的增加而增加，应变集中现象比较突出，电场促进和加速了裂纹的扩展．     

含随机分布缺陷的多晶铁电材料的有效电弹性能

针对铁电材料含有随机分布缺陷以及微结构在外场作用下发生变化的特点,建立起一个细观统计力学模型,考虑到缺陷和酶极化转动之间的相互影响,利用细观力学方法－Eshelby等效夹杂法和Mori-Tanaka的平均理论,具体分析了所含缺陷以及外场对铁电材料有效电弹性能和模量的影响。针对BaTiO3铁电陶瓷的有效电弹性能与常数的预报结果表明了缺陷的存在将增强材料的压电性能。     

铁电陶瓷的单电畴力学模型

物理实验证实铁电陶瓷的畴变是一个复杂的物理渐变过程,因此论文中引入了电畴翻转体积分数来表征这一过程,以电畴翻转体积分数为内变量,建立了考虑电畴翻转体积分数的铁电材料本构关系.接着提出了单电畴力学模型,把极化后的铁电试样看作横观各向同性单电畴的均匀材料.利用考虑电畴翻转体积分数的铁电材料本构关系,根据部分实验数据,拟合出电畴翻转驱动力与电畴翻转体积分数之间的强化关系曲线.然后根据标定出的强化关系曲线计算铁电试样在力电载荷下的响应.理论预测与实验结果符合良好.     

乘波构形和乘波飞行器研究综述

乘波构形的特点是高升阻比,下表面上的流动是均匀的,因此是推进系统/机身一体化设计的理想候选构形.乘波飞行器是源于乘波构形的高超音速飞行器,利用了乘波构形的高升阻比,并可为吸气发动机提供已知的均匀流场.本文比较全面地总结了乘波构形的生成方法和乘波飞行器的设计方法,介绍了乘波构形的优化方法及影响因素,给出了优化的乘波构形, 并介绍了乘波飞行器的研究进展,提出了今后的研究重点.      

A research note on certain loading directions in rate independent plasticity

Starting with the strain space formulation of rate independent elastic-plastic materials, the notion of “work-instability”, i.e. the possibility of an elastic plastic material undergoing spontaneous plastic deformation in a workless manner, is discussed. Some simple criteria for a material to be work stable are presented. An example of a material that satisfies Prandtl-Reuss type equations with strain softening is presented and it is demonstrated that for plane homogeneous deformations, the material shows the possibility of spontaneous plastic flow even if the in-plane strains are held fixed, if it softens beyond a critical amount.     

On some connections between equivalent single material and mixture theory models for fiber reinforced hyperelastic materials

There are two approaches that can be used to model the large strain mechanical response of material systems in which elastic fibers are embedded in an elastic matrix. In the first approach, a fiber reinforced material undergoing large deformation is homogenized in the sense that it is assumed to act as an equivalent single material that is transversely isotropic and hyperelastic. Both constituents then share a common reference configuration, which is typically assumed to be a natural or stress-free configuration for the equivalent single material. The stress depends on a single deformation gradient defined with respect to the natural configuration.In the second approach, the fiber/matrix system is treated as a mixture, with the matrix and the fibrous constituents having their own reference configurations and material symmetries. The total stress depends on the deformation gradients and material symmetries for both constituents, defined with respect to their reference configurations.Under appropriate assumptions, the constitutive theory developed using mixture theory can coincide with the constitutive theory assuming an equivalent single material that is transversely isotropic and hyperelastic. This paper explores the connection between the two approaches by considering the various reference configurations and material symmetries.     

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

A FE formulation for elasto-plastic materials with planar anisotropic yield functions and plastic spin

In this article a stress integration algorithm for shell problems with planar anisotropic yield functions is derived. The evolution of the anisotropy directions is determined on the basis of the plastic and material spin. It is assumed that the strains inducing the anisotropy of the pre-existing preferred orientation are much larger than subsequent strains due to further deformations. The change of the locally preferred orientations to each other during further deformations is considered to be neglectable. Sheet forming processes are typical applications for such material assumptions. Thus the shape of the yield function remains unchanged. The size of the yield locus and its orientation is described with isotropic hardening and plastic and material spin.The numerical treatment is derived from the multiplicative decomposition of the deformation gradient and thermodynamic considerations in the intermediate configuration. A common formulation of the plastic spin completes the governing equations in the intermediate configuration. These equations are then pushed forward into the current configuration and the elastic deformation is restricted to small strains to obtain a simple set of constitutive equations. Based on these equations the algorithmic treatment is derived for planar anisotropic shell formulations incorporating large rotations and finite strains. The numerical approach is completed by generalizing the Return Mapping algorithm to problems with plastic spin applying Hill’s anisotropic yield function. Results of numerical simulations are presented to assess the proposed approach and the significance of the plastic spin in the deformation process.     

A class of exact solutions of the ideal plasticity equations

A class of exact solutions is formulated for the ideal plasticity equations in the case of plane deformation. The solutions describe the plastic state of various wedges, notched half-planes, and domains in the form of funnels produced by detonation, the plastic state of domains having a horn configuration, etc. Many of the solutions have natural boundaries comprising envelopes of slip lines. The equations for the boundaries, slip lines, and stresses are presented in explicit form.     

Generalization of strain-gradient theory to finite elastic deformation for isotropic materials

This paper concerns finite deformation in the strain-gradient continuum. In order to take account of the geometric nonlinearity, the original strain-gradient theory which is based on the infinitesimal strain tensor is rewritten given the Green–Lagrange strain tensor. Following introducing the generalized isotropic Saint Venant–Kirchhoff material model for the strain-gradient elasticity, the boundary value problem is investigated in not only the material configuration but also the spatial configuration building upon the principle of virtual work for a three-dimensional solid. By presenting one example, the convergence of the strain-gradient and classical theories is studied.     

Measurements of Inelastic Strain Evolution of Single Solder Grain Subject to Nominal Shear Loading

Experimental Mechanics - A novel E-shaped symmetric configuration is proposed and implemented to document deformation evolutions of single-grain solder joints with different orientations. The...     

Numerical simulation of the transient shape of the red blood cell in microcapillary flow

The transient shape of a red blood cell (RBC) in a microcapillary flow is simulated under different initial conditions, including various axis orientations and centroid locations, using the LBM-DLM/FD method, which is derived from the lattice Boltzmann method and the distributed Lagrange multiplier/fictitious domain method. Although the terminal velocity is not sensitive to the initial configuration, the evolution of the velocity and the shape are determined by the initial conditions. The parachute and the slipper shape are the most probable shapes for a deformed RBC in the flow. An RBC with an initial axis orientation of 90 degrees exhibits a more complicated deformation. RBCs have a tendency to move to the centerline of a tube if an offset between the RBC centroid and the centerline exists. Our numerical results are validated by experiments, and some details beyond the experiment are provided.     

An optical method to assess electromechanical coupling in ferroelectric ceramics

A new technique is described, which allows the assessment of elastic and inelastic regions around a macroscopic defect in ferroelectric-ferroelastic ceramics. The accuracy and robustness of the method are demonstrated on a PZT plate with a centered hole subjected to uni-axial compressive stresses. From the electrical potential distribution on the sample surface, the mechanical response of the material is obtained at different load levels.     

Torsion of an Elastomeric Cylinder Undergoing Microstructural Changes

A constitutive theory which accounts for scission and cross linking processes in polymers during deformation is used to analyze the torsion of a circular bar. In each increment of deformation at a material element of the torsion bar, some volume fraction of material undergoes scission and then re-cross links to form a new network with a new reference state. The scission process reduces the ability of the material to transmit stress. The newly formed networks restore the ability of the material to transmit stress. The total stress is assumed to be the superposition of the stress in the remainder of the original material, determined by its deformation from its original configuration, and the stress in each newly formed network, determined by the deformation in that network from the configuration at which it formed.The interaction of this material response with the inhomogeneous deformation during torsion is studied. The analysis shows the evolution of regions of original and modified material, the softening effects associated with the process of scission and re-cross linking and the occurrence of residual stress and deformation on removal of load.