Search results

The partition of unity of the standard meshless Galerkin method is used as basis in expressing the discontinuity of the contact surface displacement, particularly by adding discontinuous terms into the displacement mode, and constructing the discontinuous meshless displacement field function. In this study the contact surface equation is aimed to derive from the improved Coulomb friction contact model.

In this paper based on the basic idea of meshless method, an improved moving least squares approximation function (expansion method based on out of unit division) is applied to the analysis of two-dimensional contact problems.

On the basis of this equation after discrete processing, it is combined with the discrete form of the virtual work equation with added contact conditions, and eventually transformed into a standard linear complementary problem. Moreover, it is solved by using the Lemke algorithm, and a corresponding example is provided in this research.

The proposed method can effectively control the mutual embedding of the contact surface, and the stress distribution that is the same as the actual situation can be obtained on the contact surface.

Nobody concerned with political economy can neglect the history of economic doctrines. Structural changes in the economy and society influence economic thinking and, conversely, innovative thought structures and attitudes have almost always forced economic institutions and modes of behaviour to adjust. We learn from the history of economic doctrines how a particular theory emerged and whether, and in which environment, it could take root. We can see how a school evolves out of a common methodological perception and similar techniques of analysis, and how it has to establish itself. The interaction between unresolved problems on the one hand, and the search for better solutions or explanations on the other, leads to a change in paradigma and to the formation of new lines of reasoning. As long as the real world is subject to progress and change scientific search for explanation must out of necessity continue.

This paper deals with the organizing of interactive product development. Developing products in interaction between firms may provide benefits in terms of specialization, increased innovation, and possibilities to perform development activities in parallel. However, the differentiation of product development among a number of firms also implies that various dependencies need to be dealt with across firm boundaries. How dependencies may be dealt with across firms is related to how product development is organized. The purpose of the paper is to explore dependencies and how interactive product development may be organized with regard to these dependencies.

The analytical framework is based on the industrial network approach, and deals with the development of products in terms of adaptation and combination of heterogeneous resources. There are dependencies between resources, that is, they are embedded, implying that no resource can be developed in isolation. The characteristics of and dependencies related to four main categories of resources (products, production facilities, business units and business relationships) provide a basis for analyzing the organizing of interactive product development.

Three in-depth case studies are used to explore the organizing of interactive product development with regard to dependencies. The first two cases are based on the development of the electrical system and the seats for Volvo’s large car platform (P2), performed in interaction with Delphi and Lear respectively. The third case is based on the interaction between Scania and Dayco/DFC Tech for the development of various pipes and hoses for a new truck model.

The analysis is focused on what different dependencies the firms considered and dealt with, and how product development was organized with regard to these dependencies. It is concluded that there is a complex and dynamic pattern of dependencies that reaches far beyond the developed product as well as beyond individual business units. To deal with these dependencies, development may be organized in teams where several business units are represented. This enables interaction between different business units’ resource collections, which is important for resource adaptation as well as for innovation. The delimiting and relating functions of the team boundary are elaborated upon and it is argued that also teams may be regarded as actors. It is also concluded that a modular product structure may entail a modular organization with regard to the teams, though, interaction between business units and teams is needed. A strong connection between the technical structure and the organizational structure is identified and it is concluded that policies regarding the technical structure (e.g. concerning “carry-over”) cannot be separated from the management of the organizational structure (e.g. the supplier structure). The organizing of product development is in itself a complex and dynamic task that needs to be subject to interaction between business units.

Natural convection in a square cavity with a centrally located partition is considered. While one of the side walls is fully active, the other is partly insulated. Numerical simulation, based on the finite element method, has been carried out for different lengths of the active surface. The results have been compared with the cases when the cavity is without partition as well as the case of a partitioned cavity with fully active side walls.

It has been well recognized that interface problems often contain strong singularities which make conventional numerical approaches such as uniform h‐ or p‐version of finite element methods (FEMs) inefficient. In this paper, the partition‐of‐unity finite element method (PUFEM) is applied to obtain solution for interface problems with severe singularities. In the present approach, asymptotical expansions of the analytical solutions near the interface singularities are employed to enhance the accuracy of the solution. Three different enrichment schemes for interface problems are presented, and their performances are studied. Compared to other numerical approaches such as h‐p version of FEM, the main advantages of the present method include: easy and simple formulation; highly flexible enrichment configurations; no special treatment needed for numerical integration and boundary conditions; and highly effective in terms of computational efficiency. Numerical examples are included to illustrate the robustness and performance of the three schemes in conjunction with uniform h‐ or p‐refinements. It shows that the present PUFEM formulations can significantly improve the accuracy of solution. Very often, improved convergence rate is obtained through enrichment in conjunction with p‐refinement.

It has been usual to prefer an enrichment pattern independent of the mesh when applying singular functions in the Generalized/eXtended finite element method (G/XFEM). This choice, when modeling crack tip singularities through extrinsic enrichment, has been understood as the only way to surpass the typical poor convergence rate obtained with the finite element method (FEM), on uniform or quasi-uniform meshes conforming to the crack. Then, the purpose of this study is to revisit the topological enrichment strategy in the light of a higher-order continuity obtained with a smooth partition of unity (PoU). Aiming to verify the smoothness' impacts on the blending phenomenon, a series of numerical experiments is conceived to compare the two GFEM versions: the conventional one, based on piecewise continuous PoU's, and another which considers PoU's with high-regularity.

The stress approximations right at the crack tip vicinity are qualified by focusing on crack severity parameters. For this purpose, the material forces method originated from the configurational mechanics is employed. Some attempts to improve solution using different polynomial enrichment schemes, besides the singular one, are discussed aiming to verify the transition/blending effects. A classical two-dimensional problem of the linear elastic fracture mechanics (LEFM) is solved, considering the pure mode I and the mixed-mode loadings.

The results reveal that, in the presence of smooth PoU's, the topological enrichment can still be considered as a suitable strategy for extrinsic enrichment. First, because such an enrichment pattern still can treat the crack independently of the mesh and deliver some advantage in terms of convergence rates, under certain conditions, when compared to the conventional FEM. Second, because the topological pattern demands fewer degrees of freedom and impacts conditioning less than the geometrical strategy.

Several outputs are presented, considering estimations for the J–integral and the angle of probable crack advance, this last computed from two different strategies to monitoring blending/transition effects, besides some comments about conditioning. Both h- and p-behaviors are displayed to allow a discussion from different points of view concerning the topological enrichment in smooth GFEM.

This paper aims to present a computational framework to generate numeric enrichment functions for two-dimensional problems dealing with single/multiple local phenomenon/phenomena. The two-scale generalized/extended finite element method (G/XFEM) approach used here is based on the solution decomposition, having global- and local-scale components. This strategy allows the use of a coarse mesh even when the problem produces complex local phenomena. For this purpose, local problems can be defined where these local phenomena are observed and are solved separately by using fine meshes. The results of the local problems are used to enrich the global one improving the approximate solution.

The implementation of the two-scale G/XFEM formulation follows the object-oriented approach presented by the authors in a previous work, where it is possible to combine different kinds of elements and analyses models with the partition of unity enrichment scheme. Beside the extension of the G/XFEM implementation to enclose the global–local strategy, the imposition of different boundary conditions is also generalized.

The generalization done for boundary conditions is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. The flexibility for the numerical analysis of the proposed framework is illustrated by several examples. Different analysis models, element formulations and enrichment functions are used, and the accuracy, robustness and computational efficiency are demonstrated.

This work shows a generalize imposition of different boundary conditions for global–local G/XFEM analysis through an object-oriented implementation. This generalization is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. Also, solving multiple local problems simultaneously and solving plate problems using global–local G/XFEM are other contributions of this work.

The effect of a horizontal partial porous partition on heat transfer and flow structure in a differentially heated square cavity is investigated. While the fluid flow is assumed to be governed by Navier—Stokes equations, fluid saturated porous media is assumed to be governed by Darcy’s equations. Standard Galerkin method of finite element formulation is applied for discretization of the system of equations. The non‐linearities in the discretized equations are treated with Newton‐Raphson scheme. The code developed is tested for validation for modified Rayleigh number Ra up to 400. The code is then applied to a differentially heated square cavity with a horizontal partial porous partition. While the thickness of the porous partition is found to have appreciable effect on heat transfer and flow field, width of the porous partition is found to have insignificant bearing on heat transfer except when the partition is very small and compatible to the thickness of the boundary layer developed. During the experimentation Darcy number and Rayleigh number are assumed to be constant at 10^–4 and 10⁶ respectively.

The purpose of this paper is to investigate propagation characteristics of seismic waves at the welded interface of an elastic solid and unsaturated poro-thermoelastic solid.

A theoretical formulation of partially saturated poro-thermoelastic solid is used in this study established by Zhou et al. (2019). The incidence of two primary waves (P and SV) is taken. The incident wave from the elastic solid induces two reflected waves and five refracted waves. Due to viscous pore fluids, partially saturated poro-thermoelastic solid behave dissipative, whereas elastic solid behaves non-dissipative. As a result, both reflected and incident waves are homogeneous. However, all the refracted waves are inhomogeneous. A non-singular system of linear equations is formed by the coefficients of reflection and refraction for a specified incident wave. The energy shares of various reflected and refracted waves are determined by using these reflection and refraction factors. Finally, a sensitivity analysis is performed, and the effect of critical variables on energy partitioning at the interface is observed. The numerical example shows that throughout the process of reflection/refraction, the energy of incidence is conserved at all angles of incidences.

This study demonstrated two refracted (homogeneous) and five refracted (inhomogeneous) waves due to the incident wave from elastic solid. The reflection and refraction coefficients and partitioning of incident energy are acquired as a part of diverse physical parameters of the partially saturated poro-thermoelastic media. The interference energies between unlike pairs of refracted waves have been discovered due to the dissipative behavior of unsaturated poro-thermoelastic solid.

The sensitivity of different energy shares to various aspects of the considered model is graphically analyzed for a specific numerical model. The energy balance is maintained by combining interaction energy and bulk wave energy shares.

The purpose of this paper is to develop and test an implicit scheme, accurate to the second order, for solving full Navier‐Stokes equations for three dimensional problems, using parallel algorithm.

Parallel solution to the 3‐D incompressible full Navier‐Stokes equations is presented, based on two fractional steps in time and finite element in space. The accuracy of the scheme is second order in both time and space domains. Large time‐step sizes, with Courant‐Friedrichs‐Lewy (CFL) numbers much larger than unity, are taken since the momentum equation is solved implicitly. A fourth order artificial viscosity term is added. In order to stabilize the numerical solution, fourth order artificial viscosity term is used for high Reynolds number flows. The domain decomposition technique is implemented for parallel solution to the problem with matching and non‐overlapping sub‐domains. It is aimed to study both a 3D free and mixed convection problems using the developed scheme. The segregate solution for temperature field is calibrated by a 3‐D free convection problem. Then the flow case where the forced convection is one order of magnitude higher than the free convection is studied.

It is observed that the long time solution to the flow field shows oscillatory behaviour as the Reynolds number of the flow doubled while keeping the ratio of the forced to free convection fixed. The solution using a parallel algorithm gives satisfactory results, in terms of computation time and accuracy, for the natural convection problem in cubic cavity, and, the forced cooling of a room with chilled ceiling having a parabolic geometry as presented at the end. It is observed that doubling the Reynolds number, while keeping all the parameters unchanged, varies the flow behaviour completely.

A code previously developed and published by the author only solved momentum equation and studied the velocity field. In this study, full Navier Stokes equation is solved and the code is calibrated with a well‐known 3D free‐convection for two different Rayleigh number cases and then 3D mixed convection problem is studied for two cases. Re=2000 case results, solved both by the scheme in this study and by commercial code, presented an interesting physics of the problem. For Re=2000 case, continuous cooling of the room is not possible. Doubling the Reynolds number, raising it from 1000 to 2000, while keeping all the parameters unchanged, varies the flow behaviour completely.