一种考虑孔洞形状变化的损伤细观模型及其在孔洞闭合过程中的应用

本文在一种特殊的坐标系下,建立了非线性的基体材料,有限大的椭球体中含椭球形孔洞的损伤细观模型,考虑了孔洞形状的影响.得出的粘性约束方程(或称屈服面方程)除应力∑_ij,孔隙度f,幂硬化指数m外,还与孔洞的形状有关.通过曲线拟合的方法,对Gurson方程进行了修正,使之适合于非线性的基体材料、变形状孔洞的情形.最后将此模型用于分析非线性材料内部孔洞的闭合过程.     

Shape Optimization of Elasto-Plastic Bodies

Existence of an optimal shape of a deformable body made from a physically nonlinear material obeying a specific nonlinear generalized Hooke's law (in fact, the so called deformation theory of plasticity is invoked in this case) is proved. Approximation of the problem by finite elements is also discussed.     

Shape design sensitivity analysis of plates and plane elastic solids under unilateral constraints

The local dependence of static response and eigenvalues on the shape of plates and plane elastic solids is characterized. The so-called material derivative method is used. The shape sensitivity analysis includes, besides linear problems, nonlinear problems with unilateral conditions, e.g., the frictionless contact problem for an elastic body on a rigid foundation. The results on shape sensitivity analysis can be used to obtain expressions for variations of integral functionals that arise in structural optimization problems.The authors are indebted to Professor N. Olhoff and Dr. M. P. Bendsøe for stimulating discussions and valuable comments on design sensitivity analysis.     

基于曲率插值的大变形梁单元

线性梁单元的形函数在单元大转动时会引起虚假应变,不适用于几何非线性分析．传统的几何非线性梁单元由于位移插值和转角插值的相干性,常常引起剪切闭锁等问题．该文 提出了一种平面大变形梁单元,通过单元域内的曲率插值以及曲率与节点位移之间的函数关系,将单元节点力和节点位移表示为节点曲率的函数．由于曲率插值本质上是对梁的应变进行插值,保证了单元任意刚体运动不会产生虚假的节点力;且将梁的截面形心位移表示为曲率的函数,避免了传统单元中的剪切闭锁问题．因而所提方法特别适用于梁的几何非线性分析．数值算例说明了所提方法的正确性和有效性.     

Finite element approximation of a shape optimization problem for von karman system

An optimal shape design problem of an elastic body described by a system of two nonlinear elliptic equations is considered. The problem is to find the boundary of the domain occupied by the body in such a way that the stiffnes of the system in the equilibrium state is minimized.

It is assumed that the volume of the body is constant. Moreover, the function describing the boundary of the domain and its gradient are bounded.     

Variational approach to contact problems in nonlinear elasticity

We use variational methods to study problems in nonlinear 3-dimensional elasticity where the deformation of the elastic body is restricted by a rigid obstacle. For an assigned variational problem we first verify the existence of constrained minimizers whereby we extend previous results. Then we rigorously derive the Euler-Lagrange equation as necessary condition for minimizers, which was possible before only under strong smoothness assumptions on the solution. The Lagrange multiplier corresponding to the obstacle constraint provides structural information about the nature of frictionless contact. In the case of contact with, e.g., a corner of the obstacle, we derive a qualitatively new contact condition taking into account the deformed shape of the elastic body. By our analysis it is shown here for the first time rigorously that energy minimizers really solve the mechanical contact problem. Received: 20 October 2000 / Accepted: 7 June 2001 / Published online: 5 September 2002     

A Geometric Theory of Nonlinear Morphoelastic Shells

Many thin three-dimensional elastic bodies can be reduced to elastic shells: two-dimensional elastic bodies whose reference shape is not necessarily flat. More generally, morphoelastic shells are elastic shells that can remodel and grow in time. These idealized objects are suitable models for many physical, engineering, and biological systems. Here, we formulate a general geometric theory of nonlinear morphoelastic shells that describes both the evolution of the body shape, viewed as an orientable surface, as well as its intrinsic material properties such as its reference curvatures. In this geometric theory, bulk growth is modeled using an evolving referential configuration for the shell, the so-called material manifold. Geometric quantities attached to the surface, such as the first and second fundamental forms, are obtained from the metric of the three-dimensional body and its evolution. The governing dynamical equations for the body are obtained from variational consideration by assuming that both fundamental forms on the material manifold are dynamical variables in a Lagrangian field theory. In the case where growth can be modeled by a Rayleigh potential, we also obtain the governing equations for growth in the form of kinetic equations coupling the evolution of the first and the second fundamental forms with the state of stress of the shell. We apply these ideas to obtain stress-free growth fields of a planar sheet, the time evolution of a morphoelastic circular cylindrical shell subject to time-dependent internal pressure, and the residual stress of a morphoelastic planar circular shell.     

Shape and topology optimization of contact problems by level set methods

This paper deals with the numerical solution of a topology and shape optimization problems of an elastic body in unilateral contact with a rigid foundation. The contact problem with the prescribed friction is considered. The structural optimization problem consists in finding such shape of the boundary of the domain occupied by the body that the normal contact stress along the contact boundary of the body is minimized. In the paper shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization. Level set based numerical algorithm for the solution of the shape optimization problem is proposed. Level set method is used to describe the position of the boundary of the body and its evolution on a fixed mesh. This evolution is governed by Hamilton – Jacobi equation. The speed vector field driving the propagation of the boundary of the body is given by the shape derivative of a cost functional with respect to the free boundary. Numerical examples are provided. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of roll grinding system dynamics with rotor equations and speed control

A nonlinear dynamical model for paper machine roll grinding process is investigated through a group of delay differential equations with one constant time-delay. In this model, the time-delay effect is originated from shape error traces on the surface of the roll. The contact interaction of the roll and the grindstone is based on the wear theory, and the lateral deformations of the roll as a simply supported continuous beam element inside a rotational coordinate frame and the rotational rigid body vibration system are considered. The PD-controllers of the roll and the grindstone drives are also included. The numerical simulations for time-history responses provide a view of the stability of this grinding process for the design, analysis and verification of industrial roll grinding measurements in future.     

Surface tension effects on interaction between two fluids near a wall

Interaction between two fluids near a fixed solid surface ismodelled, with surface tension acting as an important influenceon the assumed planar motion. The two fluids are immiscible,incompressible and have small density and viscosity ratios;the heavier more viscous body of fluid is approaching the solidsurface and the other fluid is lying as a thin layer in between.In the so-called supercritical range where, for both fluids,inviscid forces dominate over viscous ones, a pair of pressure–shaperelations is found which leads to a nonlinear integro-differentialequation for the unknown interface shape. Analysis, computationand comparisons are applied to the equation. Travelling-statesolutions are found of periodic and non-periodic form, includinginteresting cases which exhibit parabolic growth of the layerthickness in the far field.     

Nonlinear Alfvén/Beltrami waves—An integrable structure built around the Casimir

We formulate a nonlinear wave equations that describe amplitude and pitch modulations of one-dimensional Alfvén waves propagating on a dispersive nonlinear plasma. The well-known fact that the ideal Alfvén wave can propagate on a homogeneous ambient magnetic field with conserving an arbitrary wave shape of any amplitude is explained by invoking the Casimirs stemming from a “topological defect” (or, a kernel) in the Poisson bracket operator of the ideal magnetohydrodynamic (MHD) system. Including the Hall term, however, the Alfvén waves are affected by the dispersive effect, and the aforementioned simplicity of the ideal Alfvén waves is greatly lost; an arbitrary wave can no longer propagate with a constant shape. Yet, we observe an integrable structure in the nonlinear modulation (induced by a compressible motion) of the Alfvén waves, which is described as nonlinear deformation of “Beltrami vortex” pertaining to the Casimirs.     

两个同心旋转圆柱之间的两种流体的交界面几何形状问题

研究两个同心旋转圆柱之间的两种流体的交界面几何形状问题．利用张量分析工具,给出了忽略耗散能量影响下交界面几何形状是一种能量泛函的临界点,其对应的Euler-Lagrange方程是1个非线性椭圆边值问题．对于粘性引起的耗散能量不能忽略的情况下,同样给出了1个带有耗散能量的能量泛函,其临界点是交界面几何形状,相应的Euler-Lagrange方程也是1个二阶的非线性椭圆边值问题．这样,交界面几何形状问题转化为求解非线性椭圆边值问题．     

Predicting the elastoplastic response of fiber-reinforced metal matrix composites

This paper aims to investigate the effect of microstructure parameters (such as the cross-sectional shape of fibers and fiber volume fraction) on the stress–strain behavior of unidirectional composites subjected to off-axis loadings. A micromechanical model with a periodic microstructure is used to analyze a representative volume element. The fiber is linearly elastic, but the matrix is nonlinear. The Bodner–Partom model is used to characterize the nonlinear response of the fiber-reinforced composites. The analytical results obtained show that the flow stress of composites with square fibers is higher than with circular or elliptic ones. The difference in the elastoplastic response, which is affected by the fiber shape, can be disregarded if the fiber volume fraction is smaller than 0.15. Furthermore, the effect of fiber shape on the stress–strain behavior of the composite can be ignored if the off-axis loading angle is smaller than 30°.     

Determination of the coefficient of thermal conductivity and heat capacity per unit volume in a multilayer medium

We consider the inverse problem of simultaneously determining two time-dependent thermophysical characteristics—the coefficient of thermal conductivity and the heat capacity per unit volume—for a body having the shape of a layer situated between two other layers with known thermophysical characteristics. The necessary measurements are carried out on their outside boundaries. The problem is reduced to a system of nonlinear equations for which the existence of a solution is established by using Schauder's fixed-point theorem. We find conditions that guarantee that the solution of the inverse problem is unique. Translated fromMatematichni Metodi ta Fiziko-Mekhanichni Polya, Vol. 40, No. 2, 1997, pp. 153–159.     

具有强结构阻尼的非线性梁方程的整体动力学和控制

在本文中，我们考虑了高维具有强结构阻尼和全指数Ｂａｌａｋｒｉｓｈｎａｎ－Ｔａｙｌｏｒ阻尼的非线性固定边界可伸展的弹性梁方程，得到它的吸收集和平坦惯性流形的存在性．基于无控制方程的惯性流形的存在性，得到了相应的溢出问题的有限维反馈镇定控制．进而，此结果关于结构参数的不确定性是鲁棒的．     

Shape-Aware Matching of Implicit Surfaces Based on Thin Shell Energies

A shape sensitive, variational approach for the matching of surfaces considered as thin elastic shells is investigated. The elasticity functional to be minimized takes into account two different types of nonlinear energies: a membrane energy measuring the rate of tangential distortion when deforming the reference shell into the template shell, and a bending energy measuring the bending under the deformation in terms of the change of the shape operators from the undeformed into the deformed configuration. The variational method applies to surfaces described as level sets. It is mathematically well-posed, and an existence proof of an optimal matching deformation is given. The variational model is implemented using a finite element discretization combined with a narrow band approach on an efficient hierarchical grid structure. For the optimization, a regularized nonlinear conjugate gradient scheme and a cascadic multilevel strategy are used. The features of the proposed approach are studied for synthetic test cases and a collection of geometry processing examples.     

刚性目标形状反演的一种非线性最优化方法

发展了从声散射场的远场分布的信息来再现声刚性目标形状反问题的一种非线性最优化方法,它是通过独立地求解一个不适定的线性系统和一个适定的非线性最小化问题来实现的。对反问题的非线性和不适定性的这种分离式数值处理,使所建立方法的数值实现是非常容易和快速的,因为在确定声刚性障碍物形状的非线性最优化步中,只需求解一个只有一个未知函数的小规模的最小平方问题。该方法的另一个特别的性质是,只需要远场分布的一个Fourier系数,即可对未知的刚性目标作物形设别。进而提出了数值实现该方法的一种两步调整迭代算法。对具有各种形状的二维刚性障碍物的数值试验保证了本算法是有效和实用的。     

强非线性广义 Boussinesq 方程孤波解的波形分析及求解

上海理工大学理学院\quad 上海 200093该文建立了强非线性广义 Boussinesq 方程的耗散项、波速、渐进值与波形函数的导数之间的关系.利用适当变换和待定假设方法,作者求出了上述广义 Boussinesq 方程的扭状或钟状孤波解,还求出了以前文献中未曾提到过的余弦函数的周期波解.进一步给出了波速对波形影响的结论,即:``好'广义 Boussinesq 方程的行波当波速由小变大时,波形由钟状孤波变成余弦函数周期波解;``坏'广义 Boussinesq 方程的行波当波速由小变大时,波形由余弦函数周期波解变成钟状孤波.     

Prediction of the paper path in a laser printer by finite element analysis

A conventional problem of paper handling in a page-printing printer is to design a reliable paper path. Sometimes, this also affects the printing quality. The nonlinear theory of elastica has often been used to model paper deflection shape and to predict the paper path between some guiding rollers. The present paper proposes another method that can solve this problem more efficiently. Considering a piece of paper as a nonlinear beam, finite element analysis with geometric nonlinearity and manipulating the contact gap element between paper and guiding surfaces can be used to simulate the large-deflection behavior of paper. In this paper, the paper path in the vicinity of an OPC (organic photoconductive) drum is investigated with a general purpose finite element package. The effect of the electrostatic force between the OPC drum and the paper is considered, and the deformed shape of the paper for different forces lengths is obtained. According to those results an appropriate paper guide is added for an exact paper path.     

Global dynamics and control of a nonlinear body equation with strong structural damping

A nonlinear hinged extensible elastic body equation with strong structural damping and Balakrishnan-Taylor damping of full exponent is studied as a general model for large space structures of higher dimensions. In this paper, the absorbing sets and flat inertial manifold are obtained for this nonlinear body equation. The control spillover problem associated with the stabilization of this equation is resolved by constructing a linear finite dimensional feedback, control based on the existence of inertial manifolds of the uncontrolled equation. Moreover, the results obtained are robust with respect to the uncertainty in structural parameters. Supported by the National Natural Science Foundation of China (No. 19701023)