水和水蒸汽热力学物性参数及导数的计算与应用

本文针对现代换热设备设计、两相流动与传热研究的需要,介绍了可供科研使用的水和水蒸汽物性计算方法,完善了计算程序。根据超临界锅炉设计中的超临界水的物性计算,两相流研究和换热设备动态特性研究中对物性导数计算的需要,提供了两个算例,并进行了定性分析。     

Chemically and mechanically driven creep due to generation and annihilation of vacancies with non-ideal sources and sinks

The classical concept of Nabarro creep is extended for a general dislocation microstructure. The specific mechanism of the creep consists in generation and annihilation of vacancies at dislocation jogs acting as non-ideal sources and sinks for vacancies. This mechanism causes the climb of dislocations, allowing for local volume and shape change. The final kinetic equations, relating the dislocation microstructure and the local stress state to the creep rate, are derived by means of the thermodynamic extremal principle. Closed-form equations for the creep rate are derived for isotropic polycrystals. Based on the model the creep rate in the ferritic P-91 type steel at very low applied stress is evaluated and compared with experiment.     

Saturated Compressible and Incompressible Porous Solids: Macro- and Micromechanical Approaches

In this paper two complementary approaches are used to describe the mechanical behavior of saturated compressible and incompressible porous solids. The macroscopic investigation is based on the mixture theory, restricted by the volume fraction concept. In the micromechanical approach, a hierarchy of conditionally ensemble averaged fluid and solid phase momentum balance equations are derived for a simple model of quasi-static liquid saturated porous media. The ensemble averaged equations for both the phases agree remarkably well with the macroscopic results. A micromechanical basis for Terzhagi's effective stress concept is presented. In addition, an expression for additional partial solid stress modifying the effective stress principle, to account for deformability of solid materials, is also derived.     

Reordering and incomplete preconditioning in serial and parallel adaptive mesh refinement and coarsening flow solutions

The effects of reordering the unknowns on the convergence of incomplete factorization preconditioned Krylov subspace methods are investigated. Of particular interest is the resulting preconditioned iterative solver behavior when adaptive mesh refinement and coarsening (AMR/C) are utilized for serial or distributed parallel simulations. As representative schemes, we consider the familiar reverse Cuthill–McKee and quotient minimum degree algorithms applied with incomplete factorization preconditioners to CG and GMRES solvers. In the parallel distributed case, reordering is applied to local subdomains for block ILU preconditioning, and subdomains are repartitioned dynamically as mesh adaptation proceeds. Numerical studies for representative applications are conducted using the object‐oriented AMR/C software system libMesh linked to the PETSc solver library. Serial tests demonstrate that global unknown reordering and incomplete factorization preconditioning can reduce the number of iterations and improve serial CPU time in AMR/C computations. Parallel experiments indicate that local reordering for subdomain block preconditioning associated with dynamic repartitioning because of AMR/C leads to an overall reduction in processing time. Copyright © 2011 John Wiley & Sons, Ltd.     

Fracture and stability of materials and structural members with cracks: approaches and results

A new alternative approach to fracture problems for materials and structural elements with cracks is set out. It is based on the mechanism of local instability near defects. The approach is used to study the fracture of materials compressed along interacting cracks and the fracture of thin structural members with cracks under tension with allowance for local buckling.__________Translated from Prikladnaya Mekhanika, Vol. 40, No. 12, pp. 18–64, December 2004.     

Displacement fields and self-energies of circular and polygonal dislocation loops in homogeneous and layered anisotropic solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order −1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.     

多幅动态等和线及等差线的分离和图象处理

提出用夹层光载波法分离多幅动态等和线与等差线,将VAX机及Model75用于图象的自动处理,求出用每一载面分离了的主应力。     

岩石细观损伤破坏的观测研究

利用扫描电镜（ Ｓ Ｅ Ｍ ）,对香港白岗岩等岩样在单轴压缩状态下,进行细观结构（μｍ 尺度）观察,研究其微损伤的萌生、扩展、连接直至破坏的行为,分析了白岗岩等岩石的细观破坏机制及其与宏观力学行为的关系,得到了一些重要的感性认识和若干观察记录．为认识岩爆发生的机理及复杂应力状态下岩石损伤破坏的理论分析提供了实验依据     

Spinning rods, elliptical disks and solid ellipsoidal bodies: Elastic and plastic stresses and limit spins

The analysis of the stresses in one-, two- and three-dimensional spinning bodies is discussed in a systematic and comprehensive way. First elastic solutions are derived for rods, for elliptical-shaped flat disks and for ellipsoidal solid bodies spinning about their sideways axes. Then the spins for first plastic yield are found in each case using each of the Tresca and the von Mises yield conditions. Then upper and lower bounds on the maximum allowable limit spins where the body would globally fail assuming perfectly plastic behavior are derived. The elastic solutions at first yield always give a lower bound to that limit spin, but global failure generally does not occur until the spin is increased. A way to calculate an improved lower bound is illustrated. Upper bounds are found in a simple and new way. The method uses the fact that the volume-averaged stresses can be calculated directly from the loadings without the need for any actual stress solutions, and then it is proved that the use of those average stresses in the yield functions always gives an upper bound to the limit loads. That use of the statically determinate average stresses to obtain meaningful plastic upper bounds to limit loads is though to be a new method, and can be applied to any shape. Finally, several finite element calculations are used to determine the quantitative relations between the lower and upper bounds and the actual limit spins for ellipsoidal bodies.The results are of interest in the spin of planetary bodies, where they explain the nature of an average-stress approximate method, and in the analysis of spinning bodies in general. In addition, the approach gives a very interesting example of the utility of the limit analysis approaches of plasticity theories.     

Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

From the water wheel to turbines and hydroelectricity. Technological evolution and revolutions

Since its appearance in the first century BC, the water wheel has developed with increasing pre-industrial activities, and has been at the origin of the industrial revolution for metallurgy, textile mills, and paper mills. Since the nineteenth century, the water wheel has become highly efficient. The reaction turbine appeared by 1825, and continued to undergo technological development. The impulsion turbine appeared for high chutes, by 1880. Other turbines for low-head chutes were further designed. Turbine development was associated, after 1890, with the use of hydropower to generate electricity, both for industrial activities, and for the benefits of cities. A model “one city + one plant” was followed in the twentieth century by more complex and efficient schemes when electrical interconnection developed, together with pumped plants for energy storage.     

Plastic deformation modes in perforated sheets and their relation to yield and limit surfaces

Construction of mechanism-based plasticity theories for the homogenized response of heterogeneous materials requires identification of plastic deformation modes as a function of loading direction relative to the microstructural details. Herein, we employ an efficient homogenization theory to identify for the first time such deformation modes in plates under plane stress with hexagonal arrays of circular holes at small and intermmediate pore volume fractions, and establish their relation to the branches of initial and subsequent yield and limit surfaces. Newly introduced maps of the intrinsic geometric features of point-wise yield surfaces provide full-field picture of the investigated microstructures’ propensity for plastic strain initiation and localization. The identified characteristic plastic modes provide a rational explanation for the evolving geometric features of subsequent yield and limit surfaces whose branches represent different plastic flow mechanisms, as well as a basis for the construction of a mechanism-based homogenized plasticity theory for use in structural analysis algorithms. The results suggest the need for composite yield surfaces comprised of multiple branches in the construction of mechanism-based homogenized plasticity theory for the investigated class of porous materials.     

Predicting Noninertial Effects with Linear and Nonlinear Eddy-Viscosity, and Algebraic Stress Models

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coroiolis-modified eddy-viscosity model, a realizable nonlinear eddy-viscosity model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the validation of the turbulence models. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.     

Motion and shape control of soft robots and materials

Nonlinear Dynamics - A continuum-based approach for simultaneously controlling the motion and shape of soft robots and materials (SRM) is proposed. This approach allows for systematically computing...     

Digging and pushing lunar regolith: Classical soil mechanics and the forces needed for excavation and traction

There are many notional systems for excavating lunar regolith in NASA’s Exploration Vision. Quantitative system performance comparisons are scarce in the literature. This paper focuses on the required forces for excavation and traction as quantitative predictors of system feasibility. The rich history of terrestrial soil mechanics is adapted to extant lunar regolith parameters to calculate the forces. The soil mechanics literature often acknowledges the approximate results from the numerous excavation force models in use. An intent of this paper is to examine their variations in the lunar context. Six excavation models and one traction model are presented. The effects of soil properties are explored for each excavation model, for example, soil cohesion and friction, tool–soil adhesion, and soil density. Excavation operational parameters like digging depth, rake angle, gravity, and surcharge are examined. For the traction model, soil, operational, and machine design parameters are varied to probe choices. Mathematical anomalies are noted for several models. One conclusion is that the excavation models yield such disparate results that lunar-field testing is prudent. All the equations and graphs presented have been programmed for design use. Parameter ranges and units are included.     

Velocity and temperature dissimilarity in fully developed turbulent channel and plane Couette flows

The natural dissimilarity or decorrelation of stream-wise velocity and temperature fluctuations in fully developed turbulent channel and plane Couette flows was studied using direct numerical simulation (DNS). For both of the flow configurations, a Reynolds number of about 150 was used based on the friction velocity and half the distance between walls. Buoyancy effects were neglected, and only results with a molecular Prandtl number, Pr, equal to 1 are presented. The boundary conditions for the thermal field were a uniform source of energy in the domain and isothermal wall temperature for the channel and Couette flow, respectively. The importance of those events responsible for wall-normal turbulent fluxes in the generation of axial velocity and temperature dissimilarity was examined using conditional probability. It was found that the dissimilarity in the whole domain was higher in Couette than in channel flow. It was also found that for wall-normal turbulent fluxes (momentum and heat), the averaged dissimilarity in the whole domain was slightly more correlated with those events in the second or fourth quadrant, according to the quadrant analysis technique. For channel flow, the importance of both kinds of events was similar, while for Couette flow there was a predominance in the generation of dissimilarity by those events in the fourth quadrant. Also, for both flow configurations and throughout the wall-normal direction, it was found that in the buffer region there was a predominance of events in the fourth quadrant associated with dissimilarity for both wall-normal turbulent fluxes. In the frequency domain, the distribution of energy showed that there was a high-frequency shift experienced from the wall towards the centerline by the temperature spectrum with regards to the axial velocity spectrum, for which the action of the fluctuations of the wall-normal velocity was the main cause. In the central region of the flow, on the other hand, there was a global convergence of all spectra towards the pressure spectrum, with this convergence lower for Couette flow. Finally, it is shown that the dissimilarity in developed conditions is caused by the greater correlation existing for the temperature fluctuation with the instantaneous axial pressure gradient than for the velocity fluctuation with the instantaneous axial pressure gradient.     

Fatigue and static strength of notched and unnotched aluminum-alloy and steel specimens

This paper describes a method of presentation of fatigue data on three commonly used aircraft materials, 2024-T3 and 7075-T6 aluminum alloys and normalized SAE 4130 steel, such that variations in fatigue strength with stress-concentration factor can be shown. Comparisons of the fatigue strengths of 2024-T3 and 7075-T6 aluminum are made for the most useful range of stress-concentration factors. Static-strength results of notched and unnotched specimens of the three materials are presented to show how the strength varies with some parameters of the stress concentration. Comparison of the data with one theory for the strength of cracked specimens was made.     

材料的率相关与梯度二、四阶耦合模型的内尺度律与稳定性研究

在率相关与梯度塑性二阶耦合本构模型的基础上,提出了二、四阶率相关与梯度塑性耦合模型。采用简谐波的分析方法对材料的应变局部化及材料的稳定性进行了研究,得到了二、四阶耦合模型在一维情况下的内尺度律的变化及材料稳定性的关系,得到了波长变化的上下界及材料稳定性的条件;并对其进行了对比性研究,得到材料稳定点移动的规律。     

Accumulation and Precipitation of Salts during Groundwater Evaporation and Flow

A mathematical model for the precipitation of dissolved salts during groundwater evaporation is proposed. An asymptotic solution of the problem is obtained for a specified flow law. It is shown that there are two different regimes of solute precipitation which are determined by the evaporation front velocity and the groundwater flow rate. The dependence of the precipitated salt mass on the soil surface temperature, the atmospheric humidity, the initial solute concentration and the filtration rate is illustrated.