Determining the form of a body resulting in minimum heat flux,with laminar flow in the boundary layer

The exact solution of the problem of determining the optimal body shape for which the total thermal flux will be minimal for high supersonic flow about the body involves both computational and theoretical difficulties. Therefore, at the present time wide use is made of the inverse method, based on comparing the thermal fluxes for bodies of various specified form [1, 2]. The results of such calculations cannot always replace the solution of the direct variational problem. Therefore it is advisable to consider the direct variational problem of determining the form of a body with minimal thermal flux by using the approximate Newton formula for finding the gasdynamic parameters at the edge of the boundary layer. This approach has been used in finding the form of the body of minimal drag in an ideal fluid [3–5] arid with account for friction [6], and also for determining the form of a thin two-dimensional profile with minimal thermal flux for given aerodynamic characteristics [7].     

摩擦约束塑性力学变分不等原理的半反推法

带摩擦约束的弹塑性接触问题,由于摩擦约束条件是一种判别性的条件,它的变分问题的逆问题的研究比较困难。本文对弹塑性接触力学中的变分不等问题的逆问题进行了研究,改进了半反推法并将其应用到弹塑性变分不等原理的研究中,导出了摩擦约束弹塑性增量广义变分不等原理中的能量泛函,消除了用拉氏乘子法可能产生的临界变分现象,在证明中,巧妙地处理了增量表示的接触摩擦边界条件,避免了使用非线性泛函分析和凸分析,简化了证明。     

一类Signorini问题的边界变分不等式

本文讨论了一类简化的Signorini问题。首先将原问题和一个边值问题建立联系，其次将原问题的解分解为不带不等边界条件的变分方程的解和一个变分不等式的解。然后利用边值问题的边界积分方程将变分不等式等价地化解为边界变分不等式。这样原求区域上的第一类椭圆变分不等式问题化解为求一个区域上的变分方程和一个边界变分不等式。最后说明了边界变分不等式解的存在唯一性。文末计算了柱面和半无限刚性基础的摩擦接触问题。结论表明文中方法具有较好的精度。     

应用变域变分有限元法解平面叶栅反命题

反命题作为一种可变(未知)边界问题近年来得到了广泛的研究。本文给出了亚声速平面叶栅反命题计算的势函数变域变分有限元解法。变域变分通过把可变边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而使可变边界求解和流场分析可以完全耦合进行。本文针对亚声速平面叶栅的反命题，根据泛函的驻值必要条件，介绍了变域变分有限元方法的求解过程，最后给出了两个数值算例。这两个算例均采用四节点矩形单元的插值基函数，第一个算例用于检验程序的可靠性，第二个算例设计了一个给定叶面马赫数分布的叶型，并与试验结果进行对照。     

Boundary mixed variational inequality in friction problem

IntroductionThestaticffichonProblemdiscussedinthespaperisakindofnonlinearunilateralboundalvaluePIDbleth.Itiswell-knownthatordinarylinearelliphcboundaryvalueProblemoftencormspondstOCertainlinearvariahonalequahon,theexistenceanduniquenessofitssolutioncanbeobtainedbyusingthein-MilgramtheoremorBabuskatheorem.HoweverthisunilateralboundalvallueProblemcormSPOndstoakindof"dxedvariahonalinequality['J.Sinceitcontainsthenonlinearindifferentialfunchonal,theusuallinearizationmethodsareuseless.Thereby…     

Variational principles for hydrodynamic impact problems

We first establish the rigorous field equations of the two continuous stages before and after entering water. Then correspondently, we obtain the specific variational principles, bounded theorems, and boundary integral equations of the second stage problems. The existence of solutions are proved and the scheme of solving the solutions are provided. Finally, as a numerical example, the ship's wave resistence problem is used to demonstrate the specific application of the second stage problems and its accuracy. Then we provide a rigorous and sound theoretical basis of variational finite element method and boundary element method for calculating the accurately fundamental equations.     

CONTACT PROBLEMS AND DUAL VARIATIONAL INEQUALITY OF 2－D ELASTOPLASTIC BEAM THEORY

ＣＯＮＴＡＣＴＰＲＯＢＬＥＭＳＡＮＤＤＵＡＬＶＡＲＩＡＴＩＯＮＡＬＩＮＥＱＵＡＬＩＴＹＯＦ２－ＤＥＬＡＳＴＯＰＬＡＳＴＩＣＢＥＡＭＴＨＥＯＲＹＤａｖｉｄＹａｎｇＧａｏ（高扬）（ＤｅｐａｒｔｍｅｎｔｏｆＭａｔｈｅｍａｔｉｃｓ．ＶｉｒｇｉｎｉａＰｏｌｙｔ...     

Specific impulse losses due to friction and dispersion in a gas-film cooled liquid rocket engine nozzle

Conclusions The proposed method and program for calculating the compressible turbulent boundary layer in rocket engine nozzles with gas film cooling make it possible to determine the specific impulse losses due to friction, the heat fluxes and other characteristics of the flow. The calculations are based on the numerical solution of the equations of gas dynamics in the boundary layer approximation using a three-parameter differential turbulence model.The calculations for nozzles without film cooling showed that the contours occupying a narrow interval between the families of contours with uniform and variational characteristics have the minimum impulse losses due to friction and dispersion. In contrast to the known results, the loss minimum is displaced relative to nozzles with a variational characteristic (Rao nozzles) towards truncated nozzles with a uniform characteristic.The dependence of the maximum heat transfer to the wall in the critical throat section of the nozzle on the rate of flow of fuel into the film has been determined for nozzles with film cooling. It is shown that as the film flow rate increases, the friction losses decrease, and the minimum of the impulse losses due to friction and dispersion is shifted towards the contours with a variational characteristic, which have the minimum dispersion losses. The total impulse losses, which take into account the change in the fuel component ratio in the flow core due to the diversion of part of the fuel into the film, increase with increase in the film flow rate.The results of our numerical investigation of the effect of the contour shape and film flow rate indicate that the contour with a variational characteristic, which has near-minimum specific impulse losses due to friction and dispersion, should be used as the optimum contour for LRE nozzles.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 82–93, May–June, 1993.deceased.The authors wish to thank their colleagues at the Énergomash NPO L. P. Vereshchak and L. K. Danilyuk for assisting with the calculations, the participants in G. A. Lyubimov's seminar for discussing the results obtained, and D. A. Mel'nikov, U. G. Pirumov, and A. A. Sergienko for valuable advice.     

调和函数边界积分方程的充要条件

本文以调和函数的边值问题为例,探讨了边界积分方程的充要条件.文中首次提出了超定问题的概念,并建立了超定问题有解的一个充要条件,它也就是直接变量边界积分方程的一个充要条件.文中首次阐明了边界积分方程与变分原理的内在的联系,还指出了间接变量与直接变量两类边界积分方程之间存在着一一对应的关系.文中的慨念、思路和论点不难用于其它有变分原理的问题的边界积分方程.     

Analysis of a quasistatic contact problem with adhesion and nonlocal friction for viscoelastic materials

A mathematical model is established to describe a contact problem between a deformable body and a foundation. The contact is bilateral and modelled with a nonlocal friction law, in which adhesion is taken into account. Evolution of the bonding field is described by a first-order differential equation. The materials behavior is modelled with a nonlinear viscoelastic constitutive law. A variational formulation of the mechanical problem is derived, and the existence and uniqueness of the weak solution can be proven if the coefficient of friction is sufficiently small. The proof is based on arguments of time-dependent variational inequalities, differential equations, and the Banach fixed-point theorem.     

Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.     

Development of a solution method for frictional contact problems with complex loading

In this work, solution methods for frictional contact problems are extended to the case of moving punches and to the external loading history-dependent system states. To solve the frictional contact problems in the contact area, an iterative method is developed and implemented. Solutions of two-dimensional problems are constructed using the boundary element method. Numerical analysis is aimed at the quantitative study of effects such as the interaction of contact pressure and friction forces, estimates of the friction force differences due to the differences in the choice of local basis for the calculation of normal pressure and friction forces, and evaluation of the effects of complex loading (rotation of the rigid punch after its preliminary penetration into the solid). We find that, for the same definition of the friction force, different initial approximations lead to the same solution. At the same time, the friction forces defined either as projections onto the common tangent plane or as projections onto the plane tangent to the punch can differ quite substantially. Similar conclusions are derived for the solutions corresponding to single or multiple loading steps. The work relies on the variational principle for the solution of contact problems and numerical algorithms developed for the problems with one-sided constraints. The variational principle was first applied by Signorini [1] to the determination of the stress-strain state in a linearly deformed body in a rigid smooth shell. The modern view of the problem and its generalizations to the frictional problems and some other problems involving unilateral constraints in given in the monograph [2]. Finite difference and finite element methods in application to the problems with unilateral constraints are described in [3]. Analytical solution methods are developed in the monographs [4–6].     

Axisymmetrical bodies with minimal resistance and with minimal flow toward the surface of the body,with different characters of the flow in the boundary layer

Approximate formulas are obtained which permit calculating the distribution of the friction stress and of the local heat flux over the surface of a body of arbitrary form, with a given pressure distribution, both with laminar and with turbulent flow conditions in the boundary layer. These approximate equations were used to solve the variational problems of determination of the form of axisymmetrical bodies with a minimal resistance or with a minimal total flow of heat to the surface (in the class of bodies consisting of a flat leading edge and a lateral surface) in a hypersonic flow of viscous gas. In the solution of variational problems of determination of the optimal form of a body from the conditions of minimal resistance or of a minimal total heat flux toward the surface, we must be able to determine the distribution of the pressure, the friction stress, and the local heat flux along the surface of a body of arbitrary form. At large Reynolds numbers, the problem of determining the pressure distribution comes down to solving the Euler equation with corresponding boundary conditions. However, at the present time there are no effective methods for solving this problem (at least from the point of view of using the methods for solution of variational problems); in the solution of variational problems, to determine the pressure distribution this forces us to use various approximate methods (for example, the method of tangential wedges or cones, the Newton method, etc.). The use of such approximate formulas renders unfeasible on exact solution of the equations of the boundary layer, for which the distribution of the gas-dynamic parameters at the outer limit of the boundary layer (including also the pressure distribution) must be known previously. This makes it necessary to construct approximate formulas to determine the friction stress, _w, and the local heat flux, q_w, whose accuracy in an arbitrary body will be determined by the accuracy of the assignment of the gas-dynamic parameters at the outer limit of the boundary layer. We give below the derivation of such dependences for laminar and turbulent flow in a boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 94–102, March–April, 1971.     

A class of variational difference schemes for a singular perturbation problem

In this paper, a singularly perturbed boundary value problem for second order self-adjoint ordinary differential equation is discussed. A class of variational difference schemes is constructed by the finite element method. Uniform convergence about small parameter is proved under a weaker smooth condition with respect to the coefficients of the equation. The schemes studied in refs. [1], [3], [4] and [5] belong to the class.     

平移断层的倾角对地震产生的影响

利用边界单元法和滑移弱化摩擦本构关系分析了平移断层上地震的产生。从依赖于速率和状态的摩擦本构关系出发,通过忽略速度的影响得到了滑移弱化摩擦本构关系。建立了两种本构关系之间的联系,使得在两种模型中可以使用共同的参数。通过将地球表面模拟成一个包含在无穷大弹性介质中的无穷大裂纹,已有的边界积分方法可以直接用来分析断层的滑移。由于断层上的摩擦本构关系的非线性,得到的方程也是非线性的,采用牛顿迭代法进行求解。通过数值计算得到了平移断层上滑移位移、速度及摩擦力的分布规律。考察了平移断层的倾角对地震产生的影响,计算结果表明断层的倾角越小,地震产生的位置离地球表面越近且地震产生所需滑移的时间越短。     

Invariant Boundary‐Value Problems of an Optimally Controlled Boundary Layer

The group properties of equations of the variational problem on finding a continuous law for velocity distribution of liquid injection into an incompressible laminar boundary layer in a planar case, which ensures the minimum friction force acting on the airfoil, are considered. It is shown that the optimal injection velocity on a wedge with x=0 is finite.     

Variational analysis of thermomechanically coupled steady-state rolling problem

A steady-state, rigid-plastic rolling problem for temperature and strain-rate dependent materials with nonlocal friction is considered. A variational formulation is derived, coupling a nonlinear variational inequality for the velocity and a nonlinear vari- ational equation for the temperature. The existence and uniqueness results are obtained by a proposed fixed point method.     

Analytical solution for the problem of maximum exit velocity under Coulomb friction in gravity flow discharge chutes

In this paper, an analytical solution for the problem of finding profiles of gravity flow discharge chutes required to achieve maximum exit velocity under Coulomb friction is obtained by application of variational calculus. The model of a particle which moves down a rough curve in a uniform gravitational field is used to obtain a solution of the problem for various boundary conditions. The projection sign of the normal reaction force of the rough curve onto the normal to the curve and the restriction requiring that the tangential acceleration be non-negative are introduced as the additional constraints in the form of inequalities. These inequalities are transformed into equalities by introducing new state variables. Although this is fundamentally a constrained variational problem, by further introducing a new functional with an expanded set of unknown functions, it is transformed into an unconstrained problem where broken extremals appear. The obtained equations of the chute profiles contain a certain number of unknown constants which are determined from a corresponding system of nonlinear algebraic equations. The obtained results are compared with the known results from the literature.     

Une formulation variationnelle du problème de contact avec frottement de Coulomb

A variational relationship is proposed as the weak form of the large deformation contact problem with Coulomb friction. It is a mixed relationship involving both the displacements and the multipliers; the weighting functions are the virtual displacements and the virtual multipliers. It is shown that the proposed weak form is equivalent to the strong form of the initial/boundary value contact problem and the multipliers are equal to the contact tractions. To cite this article: A. Le van, T.H.T. Nguyen, C. R. Mecanique 336 (2008).     

Problem of equilibrium of the timoshenko plate containing a crack on the boundary of an elastic inclusion with an infinite shear rigidity

A problem of equilibrium of a composite plate consisting of a matrix and an elastic inclusion with a through crack along the boundary of this inclusion is studied. The matrix deformation is described by the Timoshenko model, and the elastic inclusion deformation is described by the Kirchhoff-Love model. Conditions of mutual non-penetration of the crack edges are imposed on the curve that describes the crack. Unique solvability of the variational problem is proved. A system of boundary conditions on the curve bounding (in the mid-plane) the elastic inclusion is obtained. A differential formulation of the problem equivalent to the initial variational formulation is given.