# The calculation of stress-strain state of anisotropic composite finite-element area with different boundary conditions on the surface

## The numerical-analytic approach to solve the stress-strain state problems of an arbitrary finite element simple-connected composite object with different surface mounting types. The main problem is solvedby using the elasticity theory spatial model.

 Рубрика Физика и энергетика Вид контрольная работа Язык английский Дата добавления 27.01.2014 Размер файла 320,5 K

### Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

The calculation of stress-strain state of anisotropic composite finite-element area with different boundary conditions on the surface

Bergulyov A.S.

The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.

Composite finite-element areas have widespread application in many fields of modern technology, such as shipbuilding, space industry, construction building etc. An important issue in terms of providing obdurability when exploiting relevant construction elements, is the matter of obtaining information about their stress-strain state.The pursuance of the research on the basis of the elasticity theory spatial model of the finite element composite areas of arbitrary shape, is connected with calculation problems. That's why one can findlimited number of scientific works on the subject.

Here we suggest an effective numerical analytic research approach of stress-strain state of the layered anisotropic finite-element areas on the basis of the elasticity theory. The approach is based upon using the spline approximation and collocation methods, that help us to reduce our space boundary problem for differential equation in partial derivatives systems to the corresponding problem for ordinary differential equations. The latter is solved by using the stable numerical technique of discrete orthogonalization.

Basic equations

Let K denote the anisotropic plate composed of conjoined cubic plates. (Fig.1)

Fig.1

Let us consider the stress-strain state problem of such fixed thickness plate in a rectangular coordinate system (X, Y, Z). According to the linear elastic theory, original equations will appear as follows:

Equilibrium equations:

(1)

Where follows the pairing tangent stress law:

Cauchy equations:

(2)

Physical equations expressing Hooke's law for anisotropic material:

(3)

As we see, twenty one independent components occur in formula (3). Generally in practice, the objects under consideration have only one plane of elastic symmetry tangent to the coordinate surface z=const. In this situation, the number of explanatory variables is equal to thirteen, and Hooke's law for such objects changes to:

(4)

Coefficients of this system are determined from:

, , , , ,

, , , , (5)

, , , .

In these evaluations:

, , - elasticity modulus in x, y, z-directions respectively;

- shear modulus for the planes parallel to coordinate planes x=const, y=const, z=const respectively;

, ,,,, - Poisson's ratios characterizing lateral contraction when stretched along the coordinate axes; , - ratios characterizing shears in the planes parallel to the coordinate planes, resulting from shearing stresses that the planes tangent to other coordinate planes experience. , - coefficients of mutual influence characterizing shears in coordinate planes, due to normal stresses;

, , - coefficients of mutual influence characterizing expansions in length due to shear stress.

Further, by applying elementary transformations we obtain resolving system of three second-order partial differential equationsdescribing the stress-strain state of the plate from the systems (1), (2), (4):

(6)

Coefficients in system (6), are determined by the mechanical characteristics of the material, taking into account ratios (5)

In practice, the formulas for the anisotropic object with a single plane elastic symmetry, are developed from the formulas for the orthotropic object by rotating the coordinate system by an angle about axis OZ

In this connection, the coefficients of the elastic compliance matrixes иare in both cases interrelated by the following formulas:

(7)

By this means, knowing coefficients for an orthotropic object and the angle of rotation we can formulate a well-posed problem for the anisotropic object with a single plane of elastic symmetry.

The solution of (6) will be sought in the form of (8)

(8)

Here: N,M - the dimension of the collocation grid for each cubical plate. K - the number of cubic subplates of the plate [1].

Statement of the boundary conditions

Fig. 2

Let's consider cubic component of the plate [1] (fig. 2). Given normal and sheer stresses on faces z=0, z=1, that will determine distribution of stresses, strains and shifts within a thick-walled plate.

One can express from the set of equations (2), (4) that:

In matrix shape, the system is of the form:

(9)

Shear stresses may be defined as follows:

Thus, boundary conditions on the faces z=0, z=c are described by the system:

(10)

Approximation of solutions of the system (10), is performed with the help of functions (8).

Boundary conditions of three types, can be set on the lateral faces: conditions, corresponding to rigidly fixed contour, conditions, corresponding to simply supported contour, conditions for free contour. Approximations of these boundary conditions is achieved by fitting spline coefficients, which determine functions.functions, will be sought in the form of:

(11)

(12)

Let's consider the cubic component of the plate as an example [1] (fig. 2) and boundary conditions on x=0 face in case of rigidly fixed and simply supported contours.

1. Rigidly fixed contour.

This case corresponds to the conditions:

(13)

Let's define from the formula (9)

As , it follows that:

. Then, for fulfillment condition, it's necessary and sufficient for condition to be met. (14)

It follows that:

The latter equation follows from . Therefore, it is sufficient for at least one of the values not to be equal to zero simultaneously, to satisfy the condition (14).

I.e., the conditions (13), can be replaced by the conditions:

(15)

For the conditions (15) fulfillment is sufficient to satisfy the system:

(16)

For this, we choose the coefficients:

(17)

Taking into consideration that:

(18)

We obtain:

As we can see, selected in the form of (17), the coefficients satisfy the system (16), which automatically causes the conditions (13) fulfillment. Thus, for rigidly fixed contours x=0, functions take on form:

(19)

(20)

2. Simply supported contour.

This case satisfies the conditions:

(21)

From the preceding considerations it is seen, that conditions can be replaced by condition. For this case, we set the coefficients:

(22)

(23)

For specified coefficients according to (22), (23), we have:

As one can see, the coefficients can to satisfy the boundary conditions (21) if selected in the form of (22), (23). Therefore, for the simply supported x=0 contour, functions take on the form:

(24)

(25)

Andfunctions:

(26)

(27)

The boundary conditions for the free contour are defined similarly:

The case of non-homogenous orthotropic material

In the case of non-homogenous Z -coordinate orthotropic material, stiffness properties characterizing stress-strain state of such plate, will take the form:

-elasticity modulus in x, y, z-directions respectively;

- elasticity modulus for the planes, parallel to planes respectively.

In this case, the coefficients of the elastic compliance matrix assume the form:

(28)

After that, we have the system of three partial differential equations of the second order, which describes the stress-strain state of thick-walled rectangular orthotropic non-homogenous in Z-coordinate plate:

(29)

coefficients are defined by the coefficients of the elastic compliance matrix, described by the formulas (28)

If the material is non-homogenous in three spatial coordinates, the coefficients of the system (29) are the coefficients of all three spatial coordinates:

Statement of the problem. Let us consider a simple-connected body of arbitrary shape with an arbitrary mounting on the subsurfaces. Let T denote it. We will consider it in the Cartesian XYZ reference system. The problem is to find the stress-strain state of the current object, consisting of an arbitrary material composition, with an arbitrary mounting on the surface,at an arbitrary load on it (steady or point load) and with arbitrarily given volumetric forces (weight, thermal expansion etc.).

Fig. 3

Solution procedure

Let us arrange the considered object T in the Cartesian reference system in such a way to make XOY, XOZ, YOZ planes tangent to the object, and each point of the object satisfy the condition [Fig. 3]. Consider the case of a composite layered object. We divide the object into horizontal layers, with the lower and upper bounds that are parallel to XOY plane. Each of these layers is approximated by cubic subplates [Fig. 1]. Stress-strain state of K-area, thus obtained, will be described by the formulas (6), (10), (13), (21). The boundary conditions at the joining point will be predetermined by the formulas:

(30)

After the approximation of the required displacement functions in the form of (8), we obtain the system of ordinary differential equations of the following order:

(31)

Here: R - the number of layers of the simple-connected T-object, - the number of cubic subplates in the layer, - the number of the spline collocation points per direction. The problem can be solvedby using the stable numerical technique of discrete orthogonalization (1, 2).

Conclusions

anisotropic composite elasticity

This paper presents the numerical-analytic approach to solve the stress-strain state problems of an arbitrary finite element simple-connected composite object with different surface mounting types. The problem is solvedby using the elasticity theory spatial model. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. End problem of ordinary higher-order differential equation system is solved by using the stable numerical technique of discrete orthogonalization (1, 2).

Reference list

1. The solution of the spatial boundary value problems about rectangular plates bending/ Grigorenko A.Y., Bergulev A.S., Jaremchenko / / Proceedings of the National Academy of Sciences of Ukraine. - 2010. - № 10 - P.44-51.

2. Third-dimensional boundary value problem solution for the bend of rectangular disks / Grigorenko, A. Ya., Bergulev А.S., Yaremchenko // Reports of the National academy of sciences of Ukraine. - 2010. - №10 - P.44-51.

Размещено на Allbest.ru

### Подобные документы

• The chiral model of graphene based on the order parameter is suggested in the long-wave approximation, the ideal graphene plane being determined by the kink-like solution. Corrugation of the graphene surface is described in the form of ripple and rings.

статья [211,7 K], добавлен 23.05.2012

• The danger of cavitation and surface elements spillway structures in vertical spillway. Method of calculation capacity for vortex weirs with different geometry swirling device, the hydraulic resistance and changes in specific energy swirling flow.

статья [170,4 K], добавлен 22.06.2015

• A cosmological model to explain the origins of matter, energy, space, time the Big Bang theory asserts that the universe began at a certain point in the distant past. Pre-twentieth century ideas of Universe’s origins. Confirmation of the Big Bang theory.

реферат [37,2 K], добавлен 25.06.2010

• Stress in beams. Thin walled beams. Mechanical beam quality depends on several of its characteristics. The size and shape of its cross-section. Determining the size and shape of the cross section peppered. Сlosed or open cross sections of a beam.

презентация [100,6 K], добавлен 30.11.2013

• The principles of nonlinear multi-mode coupling. Consider a natural quasi-linear mechanical system with distributed parameters. Parametric approach, the theory of normal forms, according to a method of normal forms. Resonance in multi-frequency systems.

реферат [234,3 K], добавлен 14.02.2010

• Consideration of the need to apply nanotechnology in agriculture to improve nutrition in the soil, management of toxic elements in the hydrosphere, monitoring the ecological state of land, spraying of mineral substances, purifying water surfaces.

реферат [12,3 M], добавлен 25.06.2010

• The Rational Dynamics. The Classification of Shannon Isomorphisms. Problems in Parabolic Dynamics. Fundamental Properties of Hulls. An Application to the Invertibility of Ultra-Continuously Meager Random Variables. Fundamental Properties of Invariant.

диссертация [1,6 M], добавлен 24.10.2012

• The problems as the types of sentences in English, their classification, the problem of composite sentences. Sentences with only one predication and with more than one predication: simple and composite sentence. Types of sentences according to structure.

курсовая работа [25,5 K], добавлен 07.07.2009

• English stress is as a phenomenon. The nature of word stress and prominence. The placement of word stress. The questions of typology of accentual structure. Degrees of stress and rhythmical tendency. Practical analysis showing the types of stress.

курсовая работа [48,8 K], добавлен 03.05.2015

• Design of composite insulators. Structure of composite insulators. Typically a composite insulator. Designing composite insulators. The design requirements for composite insulators for 154-kV service. Еlectrical and mechanical aspects.

реферат [136,2 K], добавлен 08.01.2008

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.