热粘塑性介质中的温度率相关理论

本文提出了一种热粘塑性介质中热波的温度率相关理论，利用Ｃｌａｕｓｉｕｓ－Ｄｕｈｅｍ不等式和广义的Ｈｅｌｍｈｏｌｔｚ自由能，在一些合理的假设下，严格推导了温度率相关介质中的本构模型，即用Ｍｕｌｌｅｒ的冷度函数Ｔ代替经典热力学的绝对温度后，经典热力学中的本构理论可以方便地推广到温度率相关介质中。     

热粘塑性体的积分—微分型本构关系

应用ИЛЪЮШИН关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史泛；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分型本构方程这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价。对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实     

一种基于位错机制的动态应变时效模型

动态应变时效是由位错与溶质原子的相互作用引起的，只考虑位错与位错芯内的溶质原子（位错芯气团）的相互作用，在位错热激活运动机制的Zerilli-Armstrong热粘塑性本构模型的基础上，加以改进，并加入位错和位错芯片团的相互作用的影响，建立了一种可定量描写动态变变时效现象的本构模型，所得到的本构模型以Zerilli-Armstrong模型为基础，不仅可以描写动态应变时效现象，还可以描写金属在很大温度（77K－1000K）和应变率（10^-4-10^4s^-1)范围内的力学行为，本构模型对钽的拟合和预测与实验结果有较好的吻合。     

高温合金材料循环相关热机械疲劳寿命预测

在变温非线性运动强化规律所描述的高温合金材料热机械寿命应力－应变循环特性的基础上，讨论了应变控制的循环相关热机械疲劳寿命预测技术，所建模型采用了由应变以密度表示的损伤参数，并且引入了温度损伤系数，考虑了温度变化范围以及温度循环和应变循环相位关系对疲劳寿命的影响，在确定模型的一些参数，采用等温力学试验和疲劳试验的数据，为了把等温疲劳研究成果推广到变温疲劳分析领域，开辟了新的途径。     

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

含损伤材料的热粘塑性本构关系和柱壳破裂研究

以含内变量的本构关系理论为基础 ,结合材料损伤演化方程 ,并考虑了温度和损伤对材料参数的影响 ,得到了增量形式的热粘塑性本构关系的普适显式表达式。然后使用Bodner幂函数型粘塑性模型 ,具体推导了其增量形式的热粘塑性本构方程。接着结合在实践中有重要意义的内部爆炸载荷作用下的柱壳破裂问题 ,建立了含损伤热粘塑性柱壳破裂问题的完备方程组 ,使用有限差分方法 ,完成了对问题的数值模拟 ,并对结果进行了分析。计算结果与实验结果符合良好。     

气固两相流风机的三元计算与二元计算的比较

在单相耦合模型的基础上，采用合适的粒子反弹模型和磨损模型，分析、推导、编制了三元粒子运动、碰撞的程序，并与二元粒子程序进行了比较，验证了该三元程序的优越性，为气固两相流离心风机叶轮的研制提供了一个行之有效的理论工具。     

热-力联合作用下柱壳结构变形的计算分析

本文研究了激光连续加热和恒内压作用下柱壳结构的变形规律。以4340钢材料作为研究对象，其本构方程选用热粘塑性本构模型（Johnson-Cook模型）；选取某一应变率为临界值（本文中临界应变率取为1s^-1)。考察了激光功率密度，预载荷大小，激光作用时间等对结构变形规律的影响。主要结论有：结构失稳是导致激光辐照下充内压柱壳破坏的重要原因；激光功率密度和内压越高，结构破坏时间越短；根据辐照时间的长短，预测到了结构急速破坏，延迟破坏和不破坏三种模式。本工作对于深入认识激光作用下预载结构的热－力联合破坏有一定意义。     

一个考虑循环应变幅值历史效应的粘塑性本构模型

本文提出了一个考虑材料应变幅值历史效应的粘塑性本构模型。在该模型中，引入了三个具有不同演化速率的背应力演化方程；建立了非弹性应变幅值历史记忆模型，对各向同性变形阻力，引入了具有先前加载历史记忆的演化方程。将本文模型用于１Ｃｒ１８Ｎｉ９Ｔｉ不锈钢循环变形行为描述中，其预言结果与实验结果吻合得很好，表明该模型能很好地描述材料的循环应变幅值历史下的循环变形行为。     

多晶金属弹粘塑性的取向元模型

基于“三维组集式本构模型”［１～３］，运用功共轭原理和场平均方法，发展了一种“取向元”概念——多晶聚集体内具有相同取向的滑移系集的一个体积平均意义上的代表性元素；由一维的率（粘性）敏感“取向元”，通过在取向空间里的积分，获得了三维的弹粘塑性本构方程式，并模拟了蠕变和松弛现象．与以往的粘塑性模型相比，本文模型不仅能同时考虑多晶金属材料的率相关性和路径相关性，而且由于引入了物理机制，因而具有较强的预测能力．数值算例也展示了该模型的合理性     

Modeling of thermally activated dislocation glide and plastic flow through local obstacles

A unified phenomenological model is developed to study the dislocation glide through weak obstacles during the first stage of plastic deformation in metals. This model takes into account both the dynamical responses of dislocations during the flight process and thermal activations while dislocations are bound by obstacle arrays. The average thermal activation rate is estimated using an analytical model based on the generalized Friedel relations. Then, the average flight velocity after an activation event is obtained numerically by discrete dislocation dynamics (DD). To simulate the dynamical dislocation behavior, the inertia term is implemented into the equation of dislocation motion within the DD code. The results from the DD simulations, coupled with the analytical model, determine the total dislocation velocity as a function of the stress and temperatures. By choosing parameters typical of the face centered cubic metals, the model reproduces both obstacle control and drag control motion in low and high velocity regimes, respectively. As expected by other string models, dislocation overshoots of obstacles caused by the dislocation inertia at the collisions are enhanced as temperature goes down.     

A micromechanical model of hardening, rate sensitivity and thermal softening in BCC single crystals

The present paper is concerned with the development of a micromechanical model of the hardening, rate-sensitivity and thermal softening of bcc crystals. In formulating the model, we specifically consider the following unit processes: double-kink formation and thermally activated motion of kinks; the close-range interactions between primary and forest dislocations, leading to the formation of jogs; the percolation motion of dislocations through a random array of forest dislocations introducing short-range obstacles of different strengths; dislocation multiplication due to breeding by double cross-slip; and dislocation pair annihilation. The model is found to capture salient features of the behavior of Ta crystals such as: the dependence of the initial yield point on temperature and strain rate; the presence of a marked stage I of easy glide, specially at low temperatures and high strain rates; the sharp onset of stage II hardening and its tendency to shift towards lower strains, and eventually disappear, as the temperature increases or the strain rate decreases; the parabolic stage II hardening at low strain rates or high temperatures; the stage II softening at high strain rates or low temperatures; the trend towards saturation at high strains; the temperature and strain-rate dependence of the saturation stress; and the orientation dependence of the hardening rate.     

A single theory for some quasi-static,supersonic, atomic,and tectonic scale applications of dislocations

We describe a model based on continuum mechanics that reduces the study of a significant class of problems of discrete dislocation dynamics to questions of the modern theory of continuum plasticity. As applications, we explore the questions of the existence of a Peierls stress in a continuum theory, dislocation annihilation, dislocation dissociation, finite-speed-of-propagation effects of elastic waves vis-a-vis dynamic dislocation fields, supersonic dislocation motion, and short-slip duration in rupture dynamics.     

Effect of solutes on dislocation motion —a phase-field simulation

Based on recent advances in phase-field models for integrating phase and defect microstructures as well as dislocation dynamics, the interactions between diffusional solutes and moving dislocations under applied stresses are studied in three dimensions. A new functional form for describing the eigenstrains of dislocations is proposed, eliminating the dependence of the magnitude of the dislocation Burgers vector on the applied stress and providing correct stress fields of dislocations. A relationship between the velocity of the dislocation and the applied stress is obtained by theoretical analysis and numerical simulations. The operation of Frank–Read sources in the presence of diffusional solutes, the effect of chemical interactions in solid solution on the equilibrium distribution of Cottrell atmosphere, and the drag effect of Cottrell atmosphere on dislocation motion are examined. The results demonstrate that the phase-field model correctly describes the long-range elastic interactions and short-range chemical interactions that control dislocation motion.     

Tensile response of passivated films with climb-assisted dislocation glide

The tensile response of single crystal films passivated on two sides is analysed using climb enabled discrete dislocation plasticity. Plastic deformation is modelled through the motion of edge dislocations in an elastic solid with a lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation incorporated through a set of constitutive rules. The dislocation motion in the films is by glide-only or by climb-assisted glide whereas in the surface passivation layers dislocation motion occurs by glide-only and penalized by a friction stress. For realistic values of the friction stress, the size dependence of the flow strength of the oxidised films was mainly a geometrical effect resulting from the fact that the ratio of the oxide layer thickness to film thickness increases with decreasing film thickness. However, if the passivation layer was modelled as impenetrable, i.e. an infinite friction stress, the plastic hardening rate of the films increases with decreasing film thickness even for geometrically self-similar specimens. This size dependence is an intrinsic material size effect that occurs because the dislocation pile-up lengths become on the order of the film thickness. Counter-intuitively, the films have a higher flow strength when dislocation motion is driven by climb-assisted glide compared to the case when dislocation motion is glide-only. This occurs because dislocation climb breaks up the dislocation pile-ups that aid dislocations to penetrate the passivation layers. The results also show that the Bauschinger effect in passivated thin films is stronger when dislocation motion is climb-assisted compared to films wherein dislocation motion is by glide-only.     

关于工程地震实践若干问题

本文介绍和讨论了有关工程地震实践中的某些成果和问题:(1)考虑近源地震动饱和的衰减场模式及其转换；(2)核电厂厂址设计基准地面运动和不同核法规对人工设计时程合成的技术要求对比；(3)水利水电工程场地地震安全性评价；(4)以金沙江中下游地区为例介绍地震设防区划和设防区划图系；(5)从科学性和工程可接受性, 以工程法规为基础提出工程活动断裂的定义、断裂活动性工程分类和地震断错形变工程评价。     

Microstructural effects on shear instability and shear band spacing

A constitutive relation that accounts for the thermally activated dislocation motion and microstructure interaction is used to study the stability of a homogeneous solution of equations governing the simple shearing deformations of a thermoviscoplastic body. An instability criterion and an upper bound for the growth rate of the infinitesimal deformations superimposed on the homogeneous solution are derived. By adopting Wright and Ockendon's postulate, i.e., the wavelength of the dominant instability mode with the maximum growth rate determines the minimum spacing between shear bands, the shear band spacing is computed. The effect of the initial dislocation density, the nominal strain-rate, and parameters describing the initial thermal activation and the initial microstructure interaction on the shear band spacing are delineated.     

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.