应力约束下薄板结构的拓扑优化

研究了应力约束下薄板结构的拓扑优化问题，分析了极限应力的影响，建立了拓扑优化的数学模型，讨论了若干优化过程中的技术问题，最后，进行了实例计算。     

高压干气密封扭转变形结构因素影响分析

基于圆环理论以及气体润滑理论,以复杂截面机械密封环为研究对象,建立密封扭转变形无量纲分析模型;提出了干气密封端面膜压引起端面扭转变形的计算方法,并基于该法,采用多因素优化法研究密封环结构和辅助密封圈安放位置等因素对其扭转变形的影响规律,以期使结构影响规律化.结果表明:采用解析法与集中法计算膜压扭转变形是可行的;各结构参数对扭转变形的影响规律及程度不同;扭转角随台阶半径比π2、长度比π3、O型圈位置的长度比π5的增大而增大,随宽度比π4、内径ri的增大而减少;扭转角在π4≤3.5时随长宽比π1的增大而减小,当π43.5时,随π1的增大先增大后减小,在π1=2~3间形成最大值;π1值相比π2值,对扭转角影响不大;相比π1、π2、π3和π4的影响,ri和π5影响较小,相比密封结构,O型圈位置影响较小.     

简单柔性机械臂的构造性逆动力学方法

本文提出了适于简单柔性机械臂的一种基于谱分析的构造性逆动力学方法。数值模拟和实验结果表明，这种方法提供的开环输入力矩，不仅可以实现机械臂特定的刚体转动，而且可以有效地减弱柔性对机械臂精确定位的影响。     

A theory of large-strain isotropic thermoplasticity based on metric transformation tensors

Summary A formulation of isotropic thermoplasticity for arbitrary large elastic and plastic strains is presented. The underlying concept is the introduction of a metric transformation tensor which maps a locally defined six-dimensional plastic metric onto the metric of the current configuration. This mixed-variant tensor field provides a basis for the definition of a local isotropic hyperelastic stress response in the thermoplastic solid. Following this fundamental assumption, we derive a consistent internal variable formulation of thermoplasticity in a Lagrangian as well as a Eulerian geometric setting. On the numerical side, we discuss in detail an objective integration algorithm for the mixed-variant plastic flow rule. The special feature here is a new representation of the stress return and the algorithmic elastoplastic moduli in the eigenvalue space of the Eulerian plastic metric for plane problems. Furthermore, an algorithm for the solution of the coupled problem is formulated based on an operator split of the global field equations of thermoplasticity. The paper concludes with two representative numerical simulations of thermoplastic deformation processes.     

A stream tube model for miscible flow

A simple theoretical model is described for deriving a 1-dimensional equation for the spreading of a tracer in a steady flow at the field scale. The originality of the model is to use a stochastic appoach not in the 3-dimensional space but in the 1-D space of the stream tubes. The simplicity of calculation comes from the local relationship between permeability and velocity in a 1-D flow. The spreading of a tracer front is due to local variations in the cross-sectional area of the stream tubes, which induces randomness in travel time. The derived transport equation is averaged in the main flow direction. It differs from the standard dispersion equation. The roles of time and space variables are exchanged. This result can be explained by using the statistical theory of Continuous Time Random Walk instead of a standard Random Walk. However, the two equations are very close, since their solutions have the same first and second moments. Dispersivity is found to be equal to the product of the correlation length by the variance of the logarithm of permeability, a result similar to Gelhar's macrodispersion.Nomenclature A total cross-section area of the sample - C (resident) concentration of tracer - D,D ^* dispersion coefficient - F flux of tracer - G probability distribution function for permeability in the stream-tube segments - I tracer intensity (mass crossing a surface per unit time) - K permeability - L length of the medium - M number of stream tubes in the medium - N number of segments along a stream tube - P pressure - Q total flow rate in the sample - a length of an elementary stream-tube segment - g probability distribution function for permeability in the space - i, j indices, tube numbers - q flow rate in each stream tube - s variable cross-section area of a stream tube - t, t time - u front velocity - x space variable in the flow direction - small local variation in time - , _t longitudinal, transverse dispersivity - porosity of the porous medium - correlation length in the permeability field - viscosity of the fluid - time for filling an elementary stream tube segment - standard deviation of a stochastic variable - probability distribution of arrival times (Gaussian)     

Stochastic modeling for starting-time of phase evolution of random seismic ground motions

In response to the challenge inherent in classical high-dimensional models of random ground motions, a family of simulation methods for non-stationary seismic ground motions was developed previously through employing a wave-group propagation formulation with phase spectrum model built up on the frequency components' starting-time of phase evolution. The present paper aims at extending the formulation to the simulation of non-stationary random seismic ground motions. The ground motion records associated with N—S component of Northridge Earthquake at the type-II site are investigated. The frequency components' starting-time of phase evolution of is identified from the ground motion records, and is proved to admit the Gamma distribution through data fitting. Numerical results indicate that the simulated random ground motion features zero-mean, non-stationary, and non-Gaussian behaviors, and the phase spectrum model with only a few starting-times of phase evolution could come up with a sound contribution to the simulation.     

Parameter identification for cracks in elastic plate structures based on remote strain fields

A method for the detection of cracks in plate structures is presented. In contrast to most of the common monitoring concepts taking advantage of the reflection of elastic waves at crack faces, the presented approach is based on the strain measured at different locations on the surface of the structure. This allows both the identification of crack position parameters, such as length, location and angles with respect to a reference coordinate system and the calculation of stress intensity factors (SIF). The solution of the direct problem is performed on the basis of the BFM (body force method). The inverse problem is solved applying the particle swarm optimization (PSO) algorithm. The BFM is based on the principle of linear superposition which allows the calculation of the strain field in a cracked body. The strain at an arbitrary point in the structure is replaced by the strain provided by body force doublets in the uncracked structure. The doublets as well as external loads are parameters which have to be determined solving the inverse problem by minimizing a fitness function, which is defined by a square sum of residuals between measured strain distributions and computed ones for an assumed crack. The PSO algorithm applied to the fitness function operates on the basis of a swarm of candidate solutions. Once knowing loading and crack parameters, the SIF can be determined.     

A model for the deformations and thermodynamics of liquids involving a dislocation kinetics

A model for the deformation and thermodynamics of liquids is developed that depends on dislocation kinetics. The approach uses concepts from statistical mechanics to model a stochastic evolution equation for a scalar dislocation density function. The dislocation density is used in an idealized model for the discrete discontinuous deformation due to dislocation motion and dislocation creation kinetics. The total deformation functional for a liquid is modelled as a continuum deformation of an idealized lattice structure plus the discontinuous deformation due to dislocation kinetics. This results in a thermodynamic model that has an elastic response from the continuum lattice structure and a fluid response from the dislocation kinetics.In the thermodynamics, a generalized internal energy functional is assumed to exist and to have a dependence on the functions of entropy, continuum lattice strain, scalar dislocation density, velocity, and mass density. The continuum lattice strain is termed the recoverable strain and its conjugate variable is the thermodynamic stress. The conjugate variable to the scalar dislocation density is the thermodynamic chemical potential for a dislocation configuration, somewhat analogous to Gibbs' treatment of chemical potential for various mass species.This model implies that a liquid and a crystalline solid have analogous deformation and thermodynamic responses. Their differences appear in the dislocation densities and in the dislocation chemical potentials. To illustrate the deformation response analogy, some solutions are developed for simple laminar shear flows. Also, using some concepts primarily from Kuhlmann-Wilsdorf's melting model, a definition for a specific dislocation creation heat equivalent is given. This thermodynamic formalism suggests that the melting process can be modelled as the consequence of a continuous change in the dislocation density function.Work performed under the auspices of the U.S. Department of Energy by the Lawrence Livermore National Laboratory under contract No. W-7405-ENG-48.     

Solution of the plane stress problems of strain-hardening materials described by power-law using the complex pseudo-stress function

In the present paper,the compatibility equation for the plane stress problems of power-law materials is transformed into a biharmonic equation by introducing the so-calledcomplex pseudo-stress function,which makes it possible to solve the elastic-plastic planestress problems of strain hardening materials described by power-law using the complexvariable function method like that in the linear elasticity theory.By using this generalmethod,the close-formed analytical solutions for the stress,strain and displacementcomponents of the plane stress problems’of power-law materials is deduced in the paper,which can also be used to solve the elasto-plastic plane stress problems of strain-hardeningmaterials other than that described by power-law.As an example,the problem of a power-law material infinite plate containing a circular hole under uniaxial tension is solved byusing this method,the results of which are compared with those of a known asymptoticanalytical solution obtained by the perturbation method.