Optimal control of the rotation of a system of two bodies connected by an elastic rod

Optimal controls of the rotation of a mechanical system consisting of two rigid bodies, joined by an elastic rod, through a specified angle about an axis passing through the centre of mass of one of the bodies are constructed. The problem of the optimal control of the rotation of the system through a given angle with complete suppression of the oscillations of the elastic rod at the minimum of the energy functional of the control moment and the problem of time-optimality for a specified constraint on the energy functional of the control moment are solved.     

Non-linear relations of anisotropic elasticity and a particular postulate of isotropy

A general approach to the construction of six-dimensional images of strain processes is proposed with the introduction of a vector basis which, in special cases, is identical to the well-known bases of A. A. Il’yushin, V. V. Novozhilov and Ye. I. Shemyakin and S. A. Khristianovich. The analysis of the properties of materials is based on the use of the concept of characteristic elastic states, which was introduced in the papers of J. Rychlewski. In the case of an isotropic material and four types of anisotropic materials belonging to the cubic, hexagonal, trigonal and tetragonal systems, characteristic subspaces, corresponding to the multiple eigenvalues of the elasticity tensor are defined in a six-dimensional space. In accordance with Hooke's law, the components of the stress and strain vectors in these subspaces preserve their axial alignment for any of their orthogonal transformations. The particular postulate of isotropy, formulated by Il’yushin, is therefore satisfied by definition within the framework of isotropic characteristic subspaces for linear elastic materials. An extension of the particular postulate to strain processes in non-linear anisotropic materials is proposed, on the basis of which a general form of constitutive relations is obtained containing a minimum number of experimentally determinable material functions.     

Minimum Principle of Complementary Energy for Nonlinear Elastic Cable Networks with Geometrical Nonlinearities

A minimum principle of complementary energy is established for cable networks involving only the stress components as variables with geometrical nonlinearities and nonlinear elastic materials. The minimization problem of total potential energy is reformulated as a variational problem with a convex objective functional and an infinite number of second-order cone constraints; its Fenchel dual problem is shown to coincide with the minimization problem of the complementary energy. It is of interest to note that the obtained complementary energy attains always its minimum value at the equilibrium state irrespective of the stability of the cable networks, contrary to the fact that only stationary principles have been presented for elastic trusses and continua, even in the case of a stable equilibrium state.     

Regularization of optimal design problems for bars and plates,part 1

In this paper, we consider a number of optimal design problems for elastic bars and plates. The material characteristics of rigidity of an elastic nonhomogeneous medium are taken as the control variables. A linear functional of the solutions to the equilibrium boundary-value problem is minimized under additional restrictions upon the control variables.Special variations of the control within a narrow strip provide a necessary condition for a strong local minimum (Weierstrass test). This necessary condition can be used for the detailed analysis of the following problems: bar of extremal torsional rigidity; optimal distribution of isotropic material with variable shear modulus within a plate; and optimal orientation of principal axes of elasticity in an orthotropic plate. For all of these cases, the stationary solutions violate the Weierstrass test and therefore are not optimal. This is because, in the course of the approximation of the optimal solution, the curve dividing zones occupied by materials with different rigidities displays rapid oscillations sweeping over a two-dimensional region. Within this region, the material behaves as a composite medium assembled of materials of the initial class.     

Pure shear axes and elastic strain energy

It is well known that a state of pure shear has distinct setsof basis vectors or coordinate systems: the principal axes,in which the stress is diagonal, and pure shear bases, in whichdiag = 0. The latter is commonly taken as the definition ofpure shear, although a state of pure shear is more generallydefined by tr = 0. New results are presented that characterizeall possible pure shear bases. A pair of vector functions arederived which generate a set of pure shear basis vectors fromany one member of the triad. The vector functions follow froma compatibility condition for the pure shear basis vectors,and are independent of the principal stress values. The complementarytypes of vector basis have implications for the strain energyof linearly elastic solids with cubic material symmetry: fora given state of stress or strain, the strain energy achievesits extreme values when the material cube axes are aligned withprincipal axes of stress or with a pure shear basis. This impliesthat the optimal orientation for a given state of stress iswith one or the other vector basis, depending as the stressis to be minimized or maximized, which involves the sign ofone material parameter.     

Anisotropic creep modeling based on elastic projection operators

The paper presents a unified approach for creep modeling of anisotropic materials, and is specified in more detail to the cases of isotropy, cubic symmetry and transversal isotropy. Thereby an additive decomposition of the elastic and inelastic strain tensors into dilational and isochoric Kelvin modes is assumed. Each of these modes is obtained from fourth order projection operators, resulting from solution of the eigenvalue problem for the fourth order anisotropic elasticity tensor. For simplicity the amount of strain rate for each mode is determined with a Norton type ansatz in terms of an equivalent stress, and the experimental phenomenon of primary creep is taken into account by a back stress tensor of Armstrong‐Frederick type, which is also decompose into Kelvin modes. Two numerical creep simulations investigate the crystal orientation for a compact tension specimen made out of CMSX‐4 superalloy.     

有限位移理论线弹性动力学二类和三类混合变量的最小势作用量原理和驻值余作用量原理及其应用

两个新的概念,即势作用量的概念和余作用量的概念被引入弹性动力学变分原理中.根据势作用量的概念,最小作用量原理（即Hamilton原理）被改称为最小势作用量原理.根据余作用量的概念,首次提出了驻值余作用量原理.考虑边界条件的变化并应用有限位移理论的功的互等定理,导出了以位移和应力为变分变量的二类混合变量的最小势作用量原理及驻值余作用量原理.应用应变势能密度与应力余能密度的关系式于上述二类混合变量作用量原理,导出了以位移、应力和应变为变分变量的三类混合变量的相关作用量原理.最后,应用拉氏乘子法给出了广义势作用量原理及广义余作用量原理,并且应用大挠度梁二类混合变量最小势作用量原理计算了一悬臂梁的受迫振动.     

Applications of anisotropic slipline theory with non-uniform lattice rotation

Anisotropic slipline theory, with non-uniform lattice rotation field, is used to discuss new slipline solutions for the plane strain problems of punch indentation and mode 1 stationary crack in a ductile single crystal with piecewise linear yield locus. The proposed solution allows for both linear dislocation arrays and sectors with bulk dislocation density. Such features provide considerable latitude in the number of allowable stress discontinuities, and their orientation, when compared to the solutions which assume uniform lattice rotation.     

Discussion on relationship between minimal energy and curve shapes

Energy minimization has been widely used for constructing curve and surface in the fields such as computer-aided geometric design, computer graphics. However, our testing examples show that energy minimization does not optimize the shape of the curve sometimes. This paper studies the relationship between minimizing strain energy and curve shapes, the study is carried out by constructing a cubic Hermite curve with satisfactory shape. The cubic Hermite curve interpolates the positions and tangent vectors of two given endpoints. Computer simulation technique has become one of the methods of scientific discovery, the study process is carried out by numerical computation and computer simulation technique. Our result shows that： （1） cubic Hermite curves cannot be constructed by solely minimizing the strain energy; （2） by adoption of a local minimum value of the strain energy, the shapes of cubic Hermite curves could be determined for about 60 percent of all cases, some of which have unsatisfactory shapes, however. Based on strain energy model and analysis, a new model is presented for constructing cubic Hermite curves with satisfactory shapes, which is a modification of strain energy model. The new model uses an explicit formula to compute the magnitudes of the two tangent vectors, and has the properties： （1） it is easy to compute; （2） it makes the cubic Hermite curves have satisfactory shapes while holding the good property of minimizing strain energy for some cases in curve construction. The comparison of the new model with the minimum strain energy model is included.     

立方准晶材料的剪切裂纹问题

通过引入位移函数,使得立方准晶的轴对称弹性问题归结为求解一个高阶偏微分方程.在此基础上,研究了立方准晶含有圆盘状裂纹的剪切问题,得到了此问题弹性场的精确解析解,以及应力强度因子与应变能释放率.     

柔性梁与柔韧板的有限元分析

本文采用有限单元法分析梁、板的大挠度问题,在应变位移关系中考虑了转动引起的中面伸长;在计算应变能时保留了高阶项.应用最小位能原理导出空间柔性梁元、柔韧板元的弹性刚度矩阵、线性和非线性初应力矩阵,算例表明,在不增加存贮量和计算时间的条件下可适当提高解的精度、为排除寄生的刚体位移,应采用拖带坐标系的迭代法.     

On the optimal distribution of elastic moduli of a nonhomogeneous body

The problem of the optimal distribution of elastic moduli is considered for a linearly clastic inhomogeneous body. The cost function is taken to be the work produced by the surface tractions. Necessary conditions for stationary behavior and the Weierstrass condition are obtained. The difference between maximum and minimum problems is underlined, and pecularities connected with different cost functions are discussed.     

Optimization of arches using genetic algorithm

An optimization procedure is presented for the minimum weight and strain energy optimization for arch structures subjected to constraints on stress, displacement and weight responses. Both thickness and shape variables defining the natural line of the arch are considered. The computer program which is developed in this study can be used to optimize thick, thin and variable thickness curved beams/arches. An automated optimization procedure is adopted which integrates finite element analysis, parametric cubic spline geometry definition, automatic mesh generation and genetic algorithm methods. Several examples are presented to illustrate optimal arch structures with smooth shapes and thickness variations. The changes in the relative contributions of the bending, membrane and shear strain energies are monitored during the whole process of optimization.     

The problem of an interface crack with a rigid patch plate on part of its edge

A mixed boundary-value problem is solved for a piecewise-homogeneous elastic body with a rectilinear semi-infinite crack on the line where the materials are joined. A rigid patch plate (a reinforcing plate) of specified shape is attached to the upper edge of the crack on a finite interval adjacent to the crack tip. The edges of the crack are loaded with specified stresses. The body is stretched at infinity by a specified longitudinal stress. External forces with a given principal vector and moment act on the patch plate. The problem reduces to a Riemann-Hilbert boundary-value matrix problem with a piecewise-constant coefficient, the solution of which is explicitly constructed using a Gaussian hypergeometric function. The angle of rotation of the patch plate and the complex potentials describing the stress state of the body are found and the stress state of the body close to the ends of the patch plate, one of which is also simultaneously the crack tip, is investigated. Numerical examples are presented that illustrate the effect of the initial force parameters, the length of the patch plate and other parameters of the body on the angle of rotation of the patch plate and the stress state of the body.     

An exact determination of the effective moduli of elasticity of micro-inhomogeneous medial

An exact solution is obtained of the problem of determining all the macroscopic elastic constants of an anisotropic textured polycrystalline system comprising grains with cubic symmetry of the elastic properties and a discrete distribution of those properties in two equivalent orientations, differing by a rotation about the fourth-order axis through an angle of π/4.     

Optimal material orientation of nonlinear orthotropic materials

The problems of optimal material orientation are studied in the case of orthotropic elastic materials. It is assumed that the stress-strain relation (material behavior) is nonlinear and can be described by a transcendental relation including a logarithmic function. The orientation of the material is established from the condition that the elastic energy density attains its maximal (or minimal) value.Tartu University, Estonia. Translated from Mekhanika Kompozitnykh Materialov, Vol. 35, No. 3, pp. 335–346, May–June, 1999.     

Quasistatic anti-plane motion in the simplest theory of nonlinear viscoelasticity

We consider a very simple model in the framework of differential viscoelastic materials which are isotropic and incompressible. In this model the Cauchy stress tensor is split in an elastic part and a dissipative part. The elastic part is derived from a strain-energy density function only of the first invariant of the Cauchy–Green strain tensor. The dissipative part is like the Navier–Stokes equations: linear in the stretching tensor with a constant viscosity parameter. For this model we provide some time and spatial estimates in the quasistatic approximations for the equations governing anti-plane shear motions. Several explicit examples for specific form of the strain energy are produced. Our results impose analytical restrictions on the mathematical properties of the strain energy to ensure a physical behavior in the creep and recovery experiments. Moreover, we show polynomial decay for the spatial behavior in the class of stress-hardening (or strain-stiffening) materials. For stress-softening materials a Phragmen–Lindelof alternative is proved.     

基于分子动力学模拟的金属构件的弹-塑性分解方法

针对金属多晶材料构件的分子动力学（MD）模拟,本文提出了一种新的弹-塑性分解方法.文章将MD运动轨迹分解为结构变形和热振动,给出了计算结构变形的方法和近似公式;针对金属多晶材料构件的当前构型,给出了基于FCC|BCC晶胞和四原子占位的四面体单元相组合的连续变形函数及计算变形梯度的算法;利用原子-连续关联模型,发展了计算当前构型应力场和弹性张量的算法.分析了当构件承受过大载荷后在材料内部所产生的微观缺陷,并将其分类标定为位错、层错、挛晶界、晶界和空位等;对于层错和挛晶界讨论了在弹性卸载过程中应满足的刚体运动约束方程;利用极小势能原理构造了基于当前构型的弹性卸载算法,进而给出了完整的基于MD模拟的计算弹-塑性应变的算法.最后,针对单晶铜纳米线拉伸的MD模拟,计算了弹-塑性应变场,验证了本文方法的合理性.
本文提出的基于MD模拟的弹-塑性分解方法,为从微观到宏观的多尺度和多模型耦合计算提供了算法支撑.     

Relaxed energy densities for small deformations of membranes

The theory of plane stress for infinitesimal deformations ofelastic membranes commonly yields solutions that are unstablewith respect to out-of-plane perturbations. This can be remediedby replacing the elastic strain energy by a certain relaxedenergy density. This procedure yields tension field theory whenit is appropriate. Minimum energy and minimum complementaryenergy theorems, based on the relaxed energy density, are proved.     

利用近似能量极小构造平面C1三次Hermite插值曲线

研究了利用近似能量极小构造平面$C^1$三次Hermite插值曲线的方法.该方法的主要的目是求出$C^1$三次Hermite插值曲线的最佳切矢.通过将应变能、曲率变化能和组合能的近似函数极小化,得到了求解最佳切矢的线性方程组.通过求解发现,近似曲率变化能极小不存在唯一解, 而近似应变能极小和近似组合能极小由于方程系统的系数矩阵为严格对角占优故都存在唯一解.最后, 通过实例表明了本文方法构造平面$C^1$三次Hermite插值曲线的有效性.