Some Approaches to the Solution of Problems on Thin Shells with Variable Geometrical and Mechanical Parameters

A number of approaches to the solution of stress problems for anisotropic inhomogeneous shells in the classical formulation are discussed. A review is made of approaches to the solution of one- and two-dimensional static problems for thin shells with variable parameters and to the solution of stress–strain problems for anisotropic shells of revolution under axisymmetric and non-axisymmetric loading, shallow convexo-convex shells, noncircular cylindrical shells, plates of various shapes, and shells of complex geometry     

Asymmetrical dynamic propagation problems on mode Ⅲ interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on mode Ⅲ interface crack are studied. The universal representations of analytical solutions are obtained by the approaches of serf-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.     

Asymmetrical dynamic propagation problems on mode III interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on modeⅢinterface crack are studied. The universal representations of analytical solutions are obtained by the approaches of self-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.     

一维六方压电准晶中唇形匀速移动裂纹问题的解析解

运用复变函数方法,通过构造保角映射,研究了在电非渗透型边界条件下一维六方压电准晶中唇形快速传播裂纹问题的反平面剪切问题,获得了III型裂纹动态的应力强度因子和电位移强度因子的解析解．研究结果表示,当裂纹传播速度趋于零时,动力学问题就还原为了静力学问题的解．当唇形裂纹的高度趋于零时,所得结果可以退化为Griffth裂纹问题．     

Application of boundary integral equations for analyzing the dynamics of elastic,viscoelastic, and poroelastic bodies

Two approaches (classical and nonclassical) of the boundary integral equation method for solving three-dimensional dynamical boundary value problems of elasticity, viscoelasticity, and poroelasticity are considered. The boundary integral equation model is used for porous materials. The Kelvin–Voigt model and the weakly singular hereditary Abel kernel are used to describe the viscoelastic properties. Both approaches permit solving the dynamic problems exactly not only in the isotropic but also in the anisotropic case. The boundary integral equation solution scheme is constructed on the basis of the boundary element technique. The numerical results obtained by the classical and nonclassical approaches are compared.     

High-order discontinuous Galerkin solver on hybrid anisotropic meshes for laminar and turbulent simulations

Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax （BL） model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory （HWENO） limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.     

Bluff-body Thermal Property and Initial State Effects on a Laminar Premixed Flame Anchoring Pattern

Bluff-body stabilized laminar flames remain at the root of many industrial applications. Such a simple flame arrangement although steady results from complex chemical, flow mixing as well as solid body thermal interactions that are still today misunderstood. Numerically, accurate predictions of such non linear problems require Conjugate Heat Transfer (CHT) approaches that are seldom because of the need for complex fluid flow solvers as well as multi-physics coupling strategies that are computationally expensive and difficult to master. Such numerical tools however provide access to fundamental elements otherwise inaccessible. Relying on Direct Numerical Simulation (DNS) CHT based predictions, the following work underlines several key features of importance to predict and understand square bluff-body stabilized flames. In the case of fluid only predictions, where the bluff-body wall temperature is fixed and assumed constant, three possible flame topologies are obtained and respectively qualified as anchored, lifted and bowed flames. Out of these three stable flow solutions, only two topologies are found physically possible whenever computed in a CHT context. Furthermore, depending on the solid material and the initial solution, the non linear CHT problem exhibits multiple solutions highlighting the complex coupling that can arise. As evidenced by these simple flame problems, such a behavior higlights the potential difficulties of predicting flame wall interaction problems where coupling schemes and turbulent closures / modeling will be required.     

裂纹体分析的权函数理论与应用: 回顾和展望

断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题.      

A new method for the non-linear deflection analysis of an infinite beam resting on a non-linear elastic foundation

The aim of this paper is to develop a new method of analyzing the non-linear deflection behavior of an infinite beam on a non-linear elastic foundation. Non-linear beam problems have traditionally been dealt with by semi-analytical approaches that involve small perturbations or by numerical methods, such as the non-linear finite element method. In this paper, in contrast, a transformed non-linear integral equation that governs non-linear beam deflection behavior is formulated to develop a new method for non-linear solutions. The proposed method requires an iteration to solve non-linear problems, but is fairly simple and straightforward to apply. It also converges quickly, whereas traditional non-linear solution procedures are generally quite complex in application. Mathematical analysis of the proposed method is performed. In addition, illustrative examples are presented to demonstrate the validity of the method developed in the present study.     

On the two basic problems of complex damping theory

In this paper, two important basic problems relating to the complex damping theory are discussed. The two problems are: (1) the essence of complex damping, and (2) the dynamic solution of structures with complex damping in time domain. To investigate the above problems, the frequency-dependent viscous damping theory is established and the correct and applicable solution of forced vibration for a complex damping oscillator is given beforehand.Supported by the Joint Earthquake Foundation     

An implicit pseudo‐spectral multi‐domain method for the simulation of incompressible flows

A pseudo‐spectral method for the solution of incompressible flow problems based on an iterative solver involving an implicit treatment of linearized convective terms is presented. The method allows the treatment of moderately complex geometries by means of a multi‐domain approach and it is able to cope with non‐constant fluid properties and non‐orthogonal problem domains. In addition, the fully implicit scheme yields improved stability properties as opposed to semi‐implicit schemes commonly employed. Key components of the method are a Chebyshev collocation discretization, a special pressure–correction scheme, and a restarted GMRES method with a preconditioner derived from a fast direct solver. The performance of the proposed method is investigated by considering several numerical examples of different complexity, and also includes comparisons to alternative solution approaches based on finite‐volume discretizations. Copyright © 2003 John Wiley & Sons, Ltd.     

关于复阻尼理论的两个基本问题

本文讨论有关复阻尼理论的两个重要的基本问题:复阻尼的实质问题和复阻尼结构动力响应的时域解法问题.为探究这两个问题,先建立了频率相关粘性阻尼理论,然后再给出复阻尼振子强迫振动响应的正确和适用的解答.     

Using discrete fourier series to solve boundary-value stress problems for elastic bodies with complex geometry and structure

The paper deals with some approaches to solving linear and nonlinear boundary-value stress problems for elastic bodies with complex geometry and structure. The problems are described by partial differential equations solved using discrete Fourier series. The results obtained are presented in the form of plots and tables     

On the theory of plane inhomogeneous waves in anisotropic elastic media

In the theory of plane inhomogeneous elastic waves, the complex wave vector constituted by two real vectors in a given plane may be described with the aid of two complex scalar parameters. Either of those parameters may be taken as a free one in the characteristic condition assigned to the wave equation. This alternative underlies the two fundamental approaches in the theory, namely, one associated with the Stroh eigenvalue problem and the other with the generalized Christoffel eigenvalue problem. The two approaches are identical insofar as a partial nondegenerate wave solution (partial mode) is concerned, but they differ in the fundamental solution (wave packet) assembling, and their dissimilarity is also revealed in the presence of degeneracies, which may involve either of the two governing parameters or both of them. Therefore, use of both approaches is essential for studying the degeneracy phenomenon in the theory of inhomogeneous waves. The criteria for different types of degeneracy, related to a double eigenvalue of the Stroh matrix or the Christoffel matrix and at the same time to a repeated root of the characteristic condition, are formulated by appeal to the matrix algebra and to the theory of polynomial equations. On this basis, dimensions of the manifolds, associated with degeneracy of different types in the space of variables, are established for elastic media of unrestricted anisotropy. The relation to the boundary-value problems is discussed.     

Minimizing sequences of variational problems with small parameters

A class of variational problems with small parameters is studied. Their zeroth-order asymptotic solutions are constructed. It is shown that the zeroth-order asymptotic solution is just the minimizing sequence of variational problems as the small parameter approaches to zero.     

Reliability Index and Failure Probability

Abstract

Numerical algorithms for the solution of nonlinear algebraic equation systems are discussed. Special application to the mechanism and multibody system kinematic analysis, as well as to the problems of constraint stabilization during dynamics simulation is regarded. Special attention is paid to the approaches of a separate solution of the differential equations and constraint stabilization. Numerical procedures that are effective additions to the well-known algorithms based on the Newton-Raphson method are presented. The problems of loss of precision and achievement of large unreal increments of the varying parameters are discussed. The traditional Newton-Raphson method is modified by applying a step reduction procedure that is developed numerically for the symbolic form of kinematic and dynamic equations. An optimization method for stabilization of constraints using the mass matrix of dynamic equations is suggested. According to the objective function defined the stabilization procedure provides minimal deviations of the parameters and their velocities with respect to the solution of the differential equations. No generalized coordinate partitioning is required either for solution of the dynamic equations or for stabilization of the constraints. Several examples of kinematic analysis of single and four contour plane mechanisms and constraint stabilization are solved, and the results are compared. The advantages of the algorithms developed are tested with a high-degree of initial deviation from the real solution. It is also shown that the step correction algorithm could provide admissible solution even when, in many cases, the classical approaches are not reliable. An example of the direct and inverse kinematic problem solutions of the four-degrees-of-freedom spatial platform is presented.     

Shakedown Analysis of Elastoplastic Arches

Abstract

The use of mathematical programming methods proves to be of interest in the formulation and the solution of some rather tedious problems of structural plasticity. This paper extends the application of linear programming techniques to elastoplastic arches under variable repeated loadings by adopting linearized models for the arch shape and the yield criteria of its sections. The concept of M and M + N mechanisms proves convenient in formulating the equilibrium equations by the static approach. With nonlinear yield conditions the formulations presented will still hold, but the solution of corresponding nonlinear programming problems will likely become more complex. In addition to their practical applications, the linear programming approaches suggested and illustrated have the advantage of allowing the incremental collapse load to be derived by a direct, algorithmic procedure, rather than by the trial and error procedures suggested in earlier investigations.     

多层正交异性矩形薄板弯曲振动稳定解析解

构造了带有补充项的双重正弦傅里叶级数通解来求解各种边界条件的多层正交各向异性矩形薄板的弯曲、振动和稳定问题.将坐标轴取在中性面上,求出用挠度表示的应力表达式,然后由横截面上每单位宽度的应力合成板的内力;再将层合板的内力代入板的平衡方程中得到板的控制方程,将多层板的物理参数折算为等价的单层板物理参数;最后联立控制方程与边界条件,求得未知量的系数并代入本文的通解中.本文的通解不需要叠加即可求解各种边界条件的板的弯曲、振动和稳定问题;现有的对于单层板的研究都可以用本文的方法拓展到多层板领域;对于复杂边界条件的板,也可以使用该通解分析.     

A typhoon-like vortex solution of incompressible 3D inviscid flow

This letter obtains a typhoon-like vortex solution by solving the Bragg-Hawthorne equation. The solution describes spiral paths of fluid material element on the Bernoulli surface, whereas some new exact solutions are obtained which are bounded in the whole region. The first one is a continued umbrella vortex solution, which is a typhoon-like vortex. The second one is a multi-planar solution, which is periodic in z-coordinate. Within each layer, there is an umbrella vortex solution similar to the first one. The above explicit solutions can be applied to the study of radial structure of typhoon. Both the solutions and the approaches used in the present work can also be applied to other complex flows.     

Analytical solutions of dynamic mode I crack under the conditions of displacement boundary

By the approaches of the theory of complex variable functions, the problems of dynamic mode I crack under the condition of displacement boundary are investigated. For this kind of dynamic crack extension problems with arbitrary index of self-similarity, the universal representations of analytical solutions are facilely deduced by the methods of self-similar functions. Analytical solutions of the stresses, displacements and stress intensity factors are readily acquired using the methods of self-similar functions. The problems studied can be very easily translated into Riemann–Hilbert problems and their closed solutions are gained rather straightforward in terms of this technique. According to corresponding material properties, the mutative rule of stress intensity factor was illustrated very well. Using those solutions and superposition theorem, the solutions of arbitrarily complex problems can be attained.