首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
2.
    
Linear elasticity problems posed on cracked domains, or domains with re‐entrant corners, yield singular solutions that deteriorate the optimality of convergence of finite element methods. In this work, we propose an optimally convergent finite element method for this class of problems. The method is based on approximating a much smoother function obtained by locally reparameterizing the solution around the singularities. This reparameterized solution can be approximated using standard finite element procedures yielding optimal convergence rates for any order of interpolating polynomials, without additional degrees of freedom or special shape functions. Hence, the method provides optimally convergent solutions for the same computational complexity of standard finite element methods. Furthermore, the sparsity and the conditioning of the resulting system is preserved. The method handles body forces and crack‐face tractions, as well as multiple crack tips and re‐entrant corners. The advantages of the method are showcased for four different problems: a straight crack with loaded faces, a circular arc crack, an L‐shaped domain undergoing anti‐plane deformation, and lastly a crack along a bimaterial interface. Optimality in convergence is observed for all the examples. A proof of optimal convergence is accomplished mainly by proving the regularity of the reparameterized solution. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

3.
    
This work presents a hybrid element formulation for the three‐dimensional penalty finite element analysis of incompressible Newtonian fluids. The formulation is based on a mixed variational statement in which velocity and stresses are treated as independent field variables. The main advantage of this formulation is that it bypasses the use of ad hoc techniques such as selective reduced integration that are commonly used in penalty‐based finite element formulations, and directly yields high accuracy for the velocity and stress fields without the need to carry out smoothing. In addition, since the stress degrees of freedom are condensed out at an element level, the cost of solving for the global degrees of freedom is the same as in a standard penalty finite element method, although the gain in accuracy for both the velocity and stress (including the pressure) fields is quite significant. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

4.
    
A numerical technique that is based on the integration of the asymptotic solution in the numerical framework for computing the local singular behavior of Stokes flow near a sharp corner is presented. Moffat's asymptotic solution is used, and special enriched shape functions are developed and integrated in the extended finite element method (X‐FEM) framework to solve the Navier–Stokes equations. The no‐slip boundary condition on the walls of the corner is enforced via the use of Lagrange multipliers. Flows around corners with different angles are simulated, and the results are compared with both those of the known analytic solution and the X‐FEM with no special enrichment near the corner. The results of the present technique are shown to greatly reduce the error made in computing the pressure and velocity fields near a corner tip without the need for mesh refinement near the corner. The method is then applied to the estimation of the permeability of a network of fibers, where it is shown that the local small‐scale pressure singularities have a large impact on the large‐scale network permeability. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
    
The damage caused by impulse waves generated in water bodies by fast landslides can be very high in terms of human lives and economic losses. The complex phenomena taking place in this highly unsteady process are difficult to model because three interacting phases: air, water and soil are involved. Solutions currently available are based on either closed form equations supported experimentally or the depth integrated Navier–Stokes equations. The latter, although of more general applicability, requires knowledge of the evolution of the bathimetry and slide drag forces and their applicability may be restricted by the steep slopes existing in most real cases. To avoid these limitations, the authors propose the solution of the full Navier–Stokes equations, using indicator functions to assign the material properties to each spatial point in the domain. The method performance is illustrated by comparison against the experimental results obtained in a physical model of an actual case. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

6.
The viscous incompressible flow in a wedge between a rigid plane and a surface of constant shear stress is calculated by use of the Mellin transform. For wedge angles below a critical value the asymptotic solution near the vertex is given by a local similarity solution. The respective stream function grows quadratically with the distance from the origin. For supercritical wedge angles the similarity solution breaks down and the leading order solution for the stream function grows with a power law having an exponent less than two. At the critical angle logarithmic terms appear in the stream function. The asymptotic dependence of the stream function found here is the same as for the 'hinged plate' problem. It is shown that the validity of the Stokes flow assumption is restricted to a vanishingly small distance from the vertex when the wedge angle is above critical and when the region of nonzero constant shear stress is extended to infinity. The relevance of the present result for technical flow systems is pointed out by comparison with the numerically calculated flow in a thermocapillary liquid bridge.  相似文献   

7.
The dynamics and heat transfer of a droplet in a mixture of gas with fine particles (dusty gas) are investigated in the presence of deposition of fine particles on the droplet (dust catching) and of phase transformations (evaporation, condensation). Separately treated is a case in which only phase transformations occur. Analytical solutions of the problem are found in the limiting cases of Stokes and Newtonian flow past a droplet (corresponding to low and high values of the Reynolds number). The combined effect of dust catching and phase transformations on the droplet motion is investigated numerically.  相似文献   

8.
In this paper, a cooled composite interface corner consisting of two bonded dissimilar materials is considered as a plane problem. With the complex variable method, the thermal residual stress field is studied analytically. It is found that the regular stress term possesses the singularity either of lnr or ln2r. The exact expressions for the corresponding singular stress field are presented.  相似文献   

9.
    
This paper presents the application of weight function method for the calculation of stress intensity factors (K) and T‐stress for surface semi‐elliptical crack in finite thickness plates subjected to arbitrary two‐dimensional stress fields. New general mathematical forms of point load weight functions for K and T have been formulated by taking advantage of the knowledge of a few specific weight functions for two‐dimensional planar cracks available in the literature and certain properties of weight function in general. The existence of the generalised forms of the weight functions simplifies the determination of specific weight functions for specific crack configurations. The determination of a specific weight function is reduced to the determination of the parameters of the generalised weight function expression. These unknown parameters can be determined from reference stress intensity factor and T‐stress solutions. This method is used to derive the weight functions for both K and T for semi‐elliptical surface cracks in finite thickness plates, covering a wide range of crack aspect ratio (a/c) and relative depth (a/t) at any point along the crack front. The derived weight functions are then validated against stress intensity factor and T‐stress solutions for several linear and nonlinear two‐dimensional stress distributions. These derived weight functions are particularly useful for the development of two‐parameter fracture and fatigue models for surface cracks subjected to fluctuating nonlinear stress fields, such as these resulting from surface treatment (shot peening), stress concentration or welding (residual stress).  相似文献   

10.
    
A unified framework of dual‐primal finite element tearing and interconnecting (FETI‐DP) algorithms is proposed for solving the system of linear equations arising from the mixed finite element approximation of incompressible Stokes equations. A distinctive feature of this framework is that it allows using both continuous and discontinuous pressures in the algorithm, whereas previous FETI‐DP methods only apply to discontinuous pressures. A preconditioned conjugate gradient method is used in the algorithm with either a lumped or a Dirichlet preconditioner, and scalable convergence rates are proved. This framework is also used to describe several previously developed FETI‐DP algorithms and greatly simplifies their analysis. Numerical experiments of solving a two‐dimensional incompressible Stokes problem demonstrate the performances of the discussed FETI‐DP algorithms represented under the same framework.Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
    
The motion of a two‐dimensional glacier is considered. At each time step, given the shape of the glacier, ice is modelled as an incompressible non‐Newtonian fluid and a non‐linear elliptic problem has to be solved to obtain the horizontal velocity field. Then, the upper surface of the glacier is updated by solving a transport equation. Finite element techniques are used to compute the velocity field whereas the transport equation is solved using a Lax–Wendroff scheme. Numerical results are compared to experiments on Gries glacier (Wallis, Switzerland) between 1961 and 1991. Then, a predition for 2021 is proposed. Copyright 2004 John Wiley & Sons, Ltd.  相似文献   

12.
13.
14.
A boundary integral equation (BIE) formulation for elasticity problems with mixed boundary conditions, proposed by Parton and Perlin (Mathematical Methods of the Theory of Elasticity, Mir, Moscow, 1984), is implemented in this paper using quadratic boundary elements. The formulation is specialised to Stokes flow problems by setting the Poisson ratio to 0·5 in the relevant kernels. The implementation is used to analyse non-trivial three dimensional problems in elasticity and Stokes flows. The results compare well with those obtained by a direct boundary element method. An outline of the extension of the formulation to non-linear problems is also given.  相似文献   

15.
    
This investigation focuses on an alternative approach to topology optimization problems involving incompressible materials using the P1‐nonconforming finite element. Instead of using the mixed displacement‐pressure formulation, a pure displacement‐based approach can be employed for finite element formulation owing to the Poisson locking‐free property of the P1‐nonconforming element. Moreover, because the P1‐nonconforming element has linear shape functions that are defined at element vertices, it has considerably fewer degrees of freedom than other quadrilateral nonconforming elements and its implementation is as simple as that of the conforming bilinear element. Various problems dealing with incompressible materials and pressure‐loaded structures found in published works are solved to verify the applicability of the proposed method. The application of the method is extended to the optimal design of fluid channels in the Stokes flow. This is done by expressing pressure in terms of volumetric strain rates and developing a velocity‐field‐only finite element formulation. The optimization results obtained from all the problems considered in this study are in close agreement with those found in the literature. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

16.
    
In this work, a mixed variational formulation to simulate quasi‐incompressible electro‐active or magneto‐active polymers immersed in the surrounding free space is presented. A novel domain decomposition is used to disconnect the primary coupled problem and the arbitrary free‐space mesh update problem. Exploiting this decomposition, we describe a block‐iterative approach to solving the linearised multiphysics problem, and a physically and geometrically based, three‐parameter method to update the free space mesh. Several application‐driven example problems are implemented to demonstrate the robustness of the mixed formulation for both electro‐elastic and magneto‐elastic problems involving both finite deformations and quasi‐incompressible media. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

17.
The main goal of this contribution is to provide an improved mixed finite element for quasi‐incompressible linear elasticity. Based on a classical least‐squares formulation, a modified weak form with displacements and stresses as process variables is derived. This weak form is the basis for a finite element with an advanced fulfillment of the momentum balance and therefore with a better performance. For the continuous approximation of stresses and displacements on the triangular and tetrahedral elements, lowest‐order Raviart–Thomas and linear standard Lagrange interpolations can be used. It is shown that coercivity and continuity of the resulting asymmetric bilinear form could be established with respect to appropriate norms. Further on, details about the implementation of the least‐squares mixed finite elements are given and some numerical examples are presented in order to demonstrate the performance of the proposed formulation. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

18.
This study aimed to investigate crack behaviour at the internal and external surfaces of the cement layer in total hip replacement. A three‐dimensional model of the femur with the cemented prosthesis was developed and analysed. Cracks were placed on the internal, external and both internal and external surfaces of the cement layer. Stress intensity factors were measured during gait. Results revealed that the stress intensity factors modes I and III were the most dominant in the crack propagation in the cement layer. The domain of mode I was the medial and lateral sides of the cement layer. Meanwhile, the domain of mode III was the anterior and posterior sides of the cement layer. The stress intensity factor and distance from the distal end indicated an inverse relationship. The internal and external cracks had no significant interaction. Moreover, stress intensity factors at the external surface of the cement layer were higher than those on the internal surface.  相似文献   

19.
    
The possibility of using free‐slip conditions within the context of the particle finite element method (PFEM) is investigated. For high Reynolds number engineering applications in which tangential effects at the fluid–solid boundaries are not of primary interest, the use of free‐slip conditions can alleviate the need for very fine boundary layer meshes. Two novel ways for the imposition of free‐slip conditions in the framework of the PFEM are presented. The proposed approach emphasizes robustness and simplicity, while retaining a sufficient level of generality. These two methods are then tested in the case of dam break and sloshing problems, and their respective advantages and drawbacks are discussed. It is also shown how the use of free‐slip conditions can indirectly improve mass conservation properties of the PFEM, even when coarse meshes are employed. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

20.
This paper is concerned with applying the mixture theory of two chemically inert incompressible Newtonian fluids to some simple unsteady flows in the annular region between two infinitely long coaxial cylinders. The equations governing the motion of the binary mixture under discussion are reduced to a system of coupled partial differential equations. With the help of finite Hankel transforms, the exact solutions of these equations are obtained in series form for the following three problems: (i) unsteady axial Couette flow in an annulus, (ii) unsteady Poiseuille flow in an annulus, (iii) unsteady circular Couette flow in an annulus.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号