构造杂交应力单元的柔度矩阵H对角化方法

证明了杂交元柔度矩阵 H非奇异的充分必要条件是假设应力模式线性无关 ;以及等价应力模式形成相同的杂交元。在此基础上建立了假设应力模式的 Hilbert子空间 ,从而可以利用 Schmidt方法简单地得到等价的正交应力模式 ,实现了柔度矩阵 H对角化 ,使得杂交元形成过程中完全避免了繁杂的矩阵求逆运算 ,提高了杂交元分析的计算效率 ,特别在柔度矩阵不容易显式求逆的材料非线性分析中更具有实际意义     

Response of an elastic Bingham fluid to oscillatory shear

The response of an elastic Bingham fluid to oscillatory strain has been modeled and compared with experiments on an oil-in-water emulsion. The newly developed model includes elastic solid deformation below the yield stress (or strain), and Newtonian flow above the yield stress. In sinusoidal oscillatory deformations at low strain amplitudes the stress response is sinusoidal and in phase with the strain. At large strain amplitudes, above the yield stress, the stress response is non-linear and is out of phase with strain because of the storage and release of elastic recoverable strain. In oscillatory deformation between parallel disks the non-uniform strain in the radial direction causes the location of the yield surface to move in-and-out during each oscillation. The radial location of the yield surface is calculated and the resulting torque on the stationary disk is determined. Torque waveforms are calculated for various strains and frequencies and compared to experiments on a model oil-in-water emulsion. Model parameters are evaluated independently: the elastic modulus of the emulsion is determined from data at low strains, the yield strain is determined from the phase shift between torque and strain, and the Bingham viscosity is determined from the frequency dependence of the torque at high strains. Using these parameters the torque waveforms are predicted quantitatively for all strains and frequencies. In accord with the model predictions the phase shift is found to depend on strain but to be independent of frequency.Notation A plate strain amplitude (parallel plates) - A _R plate strain amplitude at disk edge (parallel disks) - G elastic modulus - m torque (parallel disks) - M normalized torque (parallel disks) = 2m/R ³₀ - N ratio of viscous to elastic stresses (parallel plates) =µ A/ ₀ ratio of viscous to elastic stresses (parallel disks) =µ A _R/₀ - r normalized radial position (parallel disks) =r/R - r radial position (parallel disks) - R disk radius (parallel disks) - t normalized time = t — /2 - t time - _E elastic strain - _P plate strain (displacement of top plate or disk divided by distance between plates or disks) - _PR plate strain at disk edge (parallel disks) - ₀ yield strain - _E normalized elastic strain = _E/₀ - _P normalized plate strain = _P/₀ - _PR normalized plate strain at disk edge (parallel disks) = _PR/₀ - ₀ normalized plate strain amplitude (parallel plates) =A/ ₀ — normalized plate strain amplitude at disk edge (parallel disks) =A _R/₀ - phase shift between _P andT (parallel plates) — phase shift between _PR andM (parallel disks) - µ Bingham viscosity - stress - ₀ yield stress - T normalized stress =/ ₀ - frequency     

Dynamical analysis of cavitation for a transversely isotropic incompressible hyper-elastic medium: Periodic motion of a pre-existing micro-void

In this paper, the dynamical cavitation behavior is analyzed for a sphere composed of a class of transversely isotropic incompressible hyper-elastic materials, where there is a pre-existing micro-void in the interior of the sphere. A second-order non-linear ordinary differential equation that governs the motion of the initial micro-void is obtained by using the boundary conditions. On analyzing the qualitative properties of the solutions of the differential equation, some interesting conclusions are proposed. It is proved that the number of equilibrium points of the differential equation depends on the values of the material parameters, and that the phase diagrams of the equation are closed, smooth and convex trajectories. For any prescribed surface tensile dead-loads, the motion of the initial micro-void undergoes a non-linear periodic oscillation. The dependence of the periodic motion of the initial micro-void on material parameters and the radius of the initial micro-void is examined, and numerical results are also provided. It is worth pointing out that the conclusions in this paper can be used to describe approximately the physical implications of the dynamical formation of a cavity in the sphere.     

Dynamic dip testing as a method to assess the condensate drainage behavior from the air-side surface of compact heat exchangers

A new method to assess the condensate drainage behavior of the air-side surface of compact heat exchangers—dynamic dip testing—is introduced. The new method is shown to provide highly repeatable data for real-time drainage. Results from experiments with more than 20 flat-tube and round-tube-and-fin heat exchangers are presented, and the data clearly show geometrical effects such as the impact of the tube type on condensate drainage. By comparing the results from dip testing to wind-tunnel experiments for the same heat exchangers, we find dip testing can serve as a powerful tool for assessing the condensate retention behavior. The coils retaining the most and the least condensate in a steady-state wind-tunnel test, likewise held the most and the least in a dip test. However, different amounts of water are retained on the air-side surface during dip tests and wind-tunnel tests. A model based on gravity, surface tension and drag effects is developed to help understand and predict the drainage behavior of heat exchangers. The new model and experimental approach are useful in screening heat exchangers for condensate retention and for assessing off-cycle drainage behavior.     

The “decaying springs” model for irreversibility in polymer melts

In this work, entanglements in a polymer melt are modeled as a system of parallel springs which form and decay spontaneously. The springs are assumed to be nonlinear, and a certain fraction of them is torn apart by a certain strain.Based on these assumptions, a model of behavior in simple shear is developed. This model is shown to predict a behavior comprising that of a Wagner fluid, and is generalized to a tensorial model of single integral type. The integrand depends on a product of a material function, modeling reversible behavior, and a material functional which takes irreversible processes into account.Irreversibility of network disentanglement, which may occur when deformation changes or reverses direction, can be modeled in this way. It is shown that the two well-known Wagner constitutive equations with and without irreversibility assumptions are special cases of the model developed. In case of a deformation which does not change directions, the new material function and the material functional are multiplied to yield Wagner's damping function.When the rate of spring formation is a function of temperature, the developed model is shown to predict thermorheologically simple behavior. A constitutive equation for non-isothermal flow of polymers is developed with this assumption.     

有初内力的板的内共振

本文利用非线性时空有限方法和样条有限元技术对具有初内力的板的非线性频响特性进行了分析，计算了在不同初始内力下方板的大振幅自由振动、有阻尼强迫振动和矩形板的内共振。     

Low-frequency dilatational wave propagation through unsaturated porous media containing two immiscible fluids

An analytical theory is presented for the low-frequency behavior of dilatational waves propagating through a homogeneous elastic porous medium containing two immiscible fluids. The theory is based on the Berryman–Thigpen–Chin (BTC) model, in which capillary pressure effects are neglected. We show that the BTC model equations in the frequency domain can be transformed, at sufficiently low frequencies, into a dissipative wave equation (telegraph equation) and a propagating wave equation in the time domain. These partial differential equations describe two independent modes of dilatational wave motion that are analogous to the Biot fast and slow compressional waves in a single-fluid system. The equations can be solved analytically under a variety of initial and boundary conditions. The stipulation of “low frequency” underlying the derivation of our equations in the time domain is shown to require that the excitation frequency of wave motions be much smaller than a critical frequency. This frequency is shown to be the inverse of an intrinsic time scale that depends on an effective kinematic shear viscosity of the interstitial fluids and the intrinsic permeability of the porous medium. Numerical calculations indicate that the critical frequency in both unconsolidated and consolidated materials containing water and a nonaqueous phase liquid ranges typically from kHz to MHz. Thus engineering problems involving the dynamic response of an unsaturated porous medium to low excitation frequencies (e.g., seismic wave stimulation) should be accurately modeled by our equations after suitable initial and boundary conditions are imposed.     

Insight into the physics of foam densification via numerical simulation

Foamed materials are increasingly finding application in engineering systems on account of their unique properties. The basic mechanics which gives rise to these properties is well established, they are the result of collapsing the foam microstructure. Despite a basic understanding, the relationship between the details of foam microstructure and foam bulk response is generally unknown. With continued advances in computational power, many researchers have turned to numerical simulation to gain insight into the relationship between foam microstructure and bulk properties. However, numerical simulation of foam microscale deformation is a very challenging computational task and, to date, simulations over the full range of bulk deformations in which these materials operate have not been reported. Here a particle technique is demonstrated to be well-suited for this computational challenge, permitting simulation of the compression of foam microstructures to full densification. Computations on idealized foam microstructures are in agreement with engineering guidelines and various experimental results. Dependencies on degree of microstructure regularity and material properties are demonstrated. A surprising amount of porosity is found in fully-densified foams. The presence of residual porosity can strongly influence dynamic material response and hence needs to be accounted for in bulk (average) constitutive models of these materials.     

Constitutive equations for martensitic reorientation and detwinning in shape-memory alloys

A crystal-inelasticity-based constitutive model for martensitic reorientation and detwinning in shape-memory alloys (SMAs) has been developed from basic thermodynamics principles. The model has been implemented in a finite-element program by writing a user-material subroutine. We perform two sets of finite-element simulations to model the behavior of polycrystalline SMAs: (1) The full finite-element model where each finite element represents a collection of martensitic microstructures which originated from within an austenite single crystal, chosen from a set of crystal orientations that approximates the initial austentic crystallographic texture. The macroscopic stress-strain responses are calculated as volume averages over the entire aggregate: (2) The Taylor model (J. Inst. Metals 62 (1938) 32) where an integration point in a finite element represents a material point which consist of sets of martensitic microstructures which originated from within respective austenite single-crystals. Here the macroscopic stress-strain responses are calculated through a homogenization scheme.Experiments in tension and compression were conducted on textured polycrystalline Ti-Ni rod initially in the martensitic phase by Xie et al (Acta Mater. 46 (1998) 1989). The material parameters for the constitutive model were calibrated by fitting the tensile stress-strain response from a full finite-element calculation of a polycrystalline aggregate to the simple tension experiment. With the material parameters calibrated the predicted stress-strain curve for simple compression is in very good accord with the corresponding experiment. By comparing the simulated stress-strain response in simple tension and simple compression it is shown that the constitutive model is able to predict the observed tension-compression asymmetry exhibited by polycrystalline Ti-Ni to good accuracy. Furthermore, our calculations also show that the macroscopic stress-strain response depends strongly on the initial martensitic microstructure and crystallographic texture of the material.We also show that the Taylor model predicts the macroscopic stress-strain curves in simple tension and simple compression reasonably well. Therefore, it may be used as a relatively inexpensive computational tool for the design of components made from shape-memory materials.     

化学镀耐磨自润滑Ni-P复合镀层的摩擦磨损性能

利用化学复合镀技术制备了Ni-P-碳纳米管(Ni-P-CNTs)和Ni-P-无机类富勒烯WS2[Ni-P-(IF-WS2)]复合镀层,考察了复合镀层的减摩抗磨性能.结果表明,Ni-P-CNTs和Ni-P-(IF-WS2)复合镀层的减摩抗磨性能优于化学镀Ni-P和Ni-P-石墨镀层,其原因在于Ni-P-CNTs化学复合镀层中的碳纳米管具有优异力学性能和同轴石墨纳米管结构,而Ni-P-(IF-WS2)化学复合镀层中的IF-WS2具有封闭层状类富勒烯球形结构,二者均具有优异的自润滑性能.