Project supported by the Natural Science Foundation of Shanghai, China (Grant No.13ZR1415900).
Abstract
A hybrid natural element method (HNEM) for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation system of the HNEM for viscoelasticity problems is obtained using the Hellinger–Reissner variational principle. In contrast to the natural element method (NEM), the HNEM can directly obtain the nodal stresses, which have higher precisions than those obtained using the moving least-square (MLS) approximation. Some numerical examples are given to demonstrate the validity and superiority of this HNEM.
The meshless method is a numerical method that is an alternative to the finite element method (FEM). The trial function of the meshless method is constructed from a set of nodes, which allows the meshless method not to generate a mesh, so it will save more preliminary workload than the FEM. It has been given great attention by many researchers since the element-free Galerkin (EFG) method was presented by Belytschko.[1] In addition to the EFG method, [2, 3] there are many other meshless methods, such as the improved element-free Galerkin method, [4– 7] the complex variable meshless method, [8– 10] the complex variable element-free Galerkin (CVEFG) method, [11– 17] the interpolating element-free Galerkin (IEFG) method, [18, 19] the radial basis function (RBF) method, [20] the finite point method (FPM), [21, 22] the meshless local Petrov– Galerkin (MLPG) method, [23] the reproducing kernel particle method (RKPM), [24– 28] the meshless manifold method, [29, 30] the boundary element-free method (BEFM), [31– 40] and the improved boundary element-free method (IBEFM).[41– 44]
The natural element method (NEM) is a meshless method, first presented by Braun and Sambridge; [45] this method has been successfully applied in solid mechanics and fluid mechanics. Sukumar et al.[46] successfully solved partial differential equations of the elliptic type in solid mechanics by using this method. Bueche et al.[47] researched the NEM for elastokinetics problems. Cueto et al.[48] extended the NEM based on the concept of an a -shape and realized a three-dimensional analysis. Matinez et al.[49] applied the NEM in the analysis of fluids. Calvo et al.[50] studied the NEM for large deformation problems. Daneshmand and Niroomandi discussed the vibration modes of fluid-structure systems using the NEM.[51] Daneshmand et al.[52] used the NEM to calculate the two-dimensional flow under sluice gates. Shahrokhabadi et al.[53] used the NEM to analyze seepage problems of underground water. The natural element approximation of the Reissner– Mindlin plate for locking-free numerical analysis of plate-like thin elastic structures was established by Cho et al.[54]
Because the test functions of the NEM are smooth (C∞ )except at the nodes where they are C0, [55– 57] we cannot obtain the nodal stresses. To overcome this disadvantage, we can use the MLS approximation to calculate the nodal stresses by the nodal displacements obtained using the NEM. This method was first used in the FEM by Tabbara et al.[58] Except for this method, the hybrid natural element method (HNEM) can also obtain the nodal stresses, which has been successfully applied in elasticity problems and elastoplasticity problems.[59, 60] The HNEM adopts displacement interpolation and stress interpolation respectively using natural neighbor interpolation, [61, 62] and the corresponding formulae are obtained according to the Hellinger– Reissner variational principle.[63, 64]
Viscoelasticity problems are important nonlinearity problems in material engineering. Viscoelastic materials have the characteristics of both an elastic solid and a viscous fluid, so the stress– strain relationship is difficult to model. Peng and Xu studied the constitutive relationship of permanent deformation in asphalt pavements.[65] The constitutive relationship of viscoelastic materials correlates with time. Thus the effect of time on the material should be considered in the calculation. That the time parameter is important in the process of calculation brings great difficulty in solving the viscoelasticity mechanics problems, so analytical solutions can be obtained for only a few simple viscoelastic problems.[66] The solutions of complex viscoelastic problems in engineering mainly rely on numerical calculations, for which the meshless methods have been widely used. Yang and Liu combined the element-free Galerkin method and a precise algorithm in the time domain to study two-dimensional viscoelasticity problems.[67] Using the meshless local Petrov– Galerkin method, Sladek et al.[68, 69] analyzed anisotropic and linear viscoelastic solids, and discussed quasistatic and transient dynamic problems in two-dimensional linear viscoelastic media. Based on the rigid/viscoplastic element-free Galerkin method, the massive metal forming process was discussed by Guan et al.[70] Canelas and Sensale studied a boundary knot method for harmonic elastic and viscoelastic problems.[71] A complex variable element-free Galerkin method for linear viscoelasticity problems under the creep condition was analyzed by Cheng et al.[72] Three-dimensional viscoelasticity problems were analyzed via the improved element-free Galerkin method by Peng et al.[73]
Based on the natural neighbor interpolation, the HNEM for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. According to the linear viscoelastic Hellinger– Reissner variational principle, the equation system of the viscoelastic hybrid natural element method under the creep condition is established. Finally, the effectiveness and superiority of this method are proved through numerical examples in this paper.
2. The shape functions for the HNEM
The natural neighbor interpolation (NNI)[61, 62] in the NEM is adopted to construct the shape functions for the HNEM. Consider a bounded domain Ω in two dimensions that is described by a set N of m scattered nodes: N = {N1, N2, . . . , Nm}, m ≥ 3; their coordinates are xi, i = 1, 2, . . . , m. The Voronoi diagram V ( N) of the set N is a subdivision of the domain into regions V ( Ni) , where each region V ( Ni) is associated with a node Ni, such that any point in V ( Ni) is closer to Ni than to any other node Nj ( j ≠ i) . The Voronoi diagram in essence partitions space into closest-point regions. The region V ( Ni) (Voronoi cell) for a node Ni is a polygon
where d( xi, xj) is the Euclidean distance between xi and xj.
A part of the Voronoi diagram V ( N) is shown in Fig. 1. Suppose x to be any point within the region, the natural neighbor nodes around x are found according to the empty circumcircle criterion of Delaunay.[74] The natural neighbor nodes of x are 1, 5, 6, and 7, as shown in Fig. 2.
Fig. 2. The 1st-order and the 2nd-order Voronoi cells about x.
After the natural neighbor nodes of x are determined, we can build a local interpolation format as follows:
where f( x) is the interpolation physical quantity for node x; i is the serial number of natural neighbor nodes around x; n is the total number of natural neighbor nodes around x; fi is the physical quantity value of node i; and ϕ i( x) is the shape function of node xi. The shape function can be chosen between Sibsonian and non-Sibsonian interpolations.[62]
The Sibsonian interpolation can be written as
where Aj( x) is the area of overlap of Voronoi cells between node j and x, and A( x) is the area of the Voronoi cell of x. The four regions shown in Fig.2 are the second-order cells, while their union (closed polygon abcd) is the first-order Voronoi cell, and the shape function ϕ 1( x) is given by
The derivatives of ϕ i( x) can be obtained by differentiating Eq. (3)
(6)
For the traditional mechanical problems, we can take Eq. (3) as the displacement interpolation function
where ϕ i( x) is the shape function associated with node i, and ui is the displacement of node i.
By definition of the shape function given in Eq. (3), we have the following relations:
Referring to Fig. 2, if x is to coincide with any node, for instance node 1, then it is readily seen that ϕ 1( x) = 1 and ϕ i( x) = 0, i ≠ 1, i.e.,
which implies that the natural neighbor interpolation passes through the node values, and the node unknowns ui are the nodal displacements, which is in contrast to most meshless approximations, where the nodal parameters ui are not the nodal displacements. Thus the essential boundary conditions can be imposed easily by using the natural neighbor interpolation, as in the FEM.
3. The HNEM for two-dimensional linear viscoelasticity problems
3.1. Constitutive relationships of linear viscoelastic materials
Under the condition of small deformation, the stress tensor of an isotropic material can be described as
where
The strain tensor can be similarly written as
where
The constitutive relationship of linear viscoelastic materials can be written as
where J1 ( t) and J2 ( t) are two creep compliances. The J1 ( t) varies for different models. For the case of the Maxwell model,
for the case of the Kelvin model,
and for the case of the three-parameter model,
where G, G1, and G2 are the shear elastic moduli, and η is the viscosity constant of each model. The J2 ( t) is
where K is the volumetric modulus of elasticity, and
with E being the elasticity modulus, and ν being the Poisson ratio.
3.2. The HNEM for two-dimensional viscoelasticity problems
Under the creep condition, the stress σ is constant, while the strain ε ( t) and the displacement u ( t) depend on time, as shown in Eq. (17). For the two-dimensional problems, by assuming that the body force in the problem domain Ω , the traction on the natural boundary Γ σ , and the displacement ū of the essential boundary Γ u are all known, the basic equations of viscoelasticity problems can be expressed as follows.
The equilibrium equation can be written as
where
The strain– displacement relationship can be written as
where
The stress– strain relationship can be written as
and the constitutive matrix C is
where J1 and K can be obtained from Eqs. (18)– (20) and (22).
The boundary conditions can be written as
where
and n1 and n2 are the units outward normal to boundary Γ σ , respectively.
Suppose that there are N nodes in problem domain Ω , the displacement and the stress at an arbitrary point x can be expressed as
where n is the total number of natural neighbor nodes around x, and ϕ i is the interpolation function of displacement u and stress σ , which can be constructed by Eq. (3).
Equations (38) and (39) can be written in matrix form as
where Φ ( x) and P( x) are matrices of shape functions, q and β are the nodal displacement vector and the nodal stress vector, respectively,
The properties of linear viscoelastic materials are related to time only and have nothing to do with the loading history. So the stress and the displacement at each point in time can be calculated directly by the Hellinger– Reissner variational principle. At each point in time, the Hellinger– Reissner variational principle can be expressed as
Based on the stationary condition of the variation principle, i.e.,
the equations of the HNEM for solving the generalized displacement can be found as follows:
We can obtain the nodal displacements at any time by solving Eq. (64), and the nodal stresses can be obtained by substituting the nodal displacements in Eq. (60). Then the displacement and the stress at any point in the domain can be gained using Eqs. (40) and (41). Thus the HNEM for two-dimensional viscoelasticity problems under the creep condition is established.
4. Numerical examples
Four examples are discussed to demonstrate the applicability of the HNEM for two-dimensional linear viscoelasticity problems. The results obtained using the HNEM are compared with the analytical solutions and the ones of the NEM and the FEM. The nodal stresses obtained from the NEM are calculated with the MLS approximation.
The Sibson interpolation is used as the shape function in the numerical examples presented in this section. Delaunay triangles are used as the integration cells, and in each integration cell, three Gauss points are used for the Gaussian quadrature.
4.1. Cylinder under a distributed inner pressure
A cylinder with uniform pressure p on the inner surface is shown in Fig.3, the outer radius of the cylinder is b = 5 m, and the inner radius is a = 1 m. The inner pressure is p = 30 kPa and the volume deformation is elastic, elasticity modulus E = 1.0 × 106 Pa, Poisson ratio ν = 0.25. For the rheological properties of shear deformation, we respectively adopt the Maxwell model and the Kelvin model, and the parameters are
Fig. 3. Cylinder under a distributed inner pressure.
Regardless of the body force, the cylinder can be studied as a plane strain problem.
In the polar coordinates, the analytical solution of the radial displacement of a Maxwell material is given by[66]
where
Equation (65) shows that there is a flow phenomenon in the radial displacement of the Maxwell material, so it can be applied only in a short time; otherwise, the deformation is too large, and the solution of elastic displacement based on the small strain is not correct.
In the polar coordinates, the analytical solution of the radial displacement of a Kelvin material is given by[66]
where
Equation (67) shows that the displacement of the Kelvin material is zero at the loading instant, and it tends to be steady when the time is long enough.
Because of the symmetry, only one quarter of the cylinder needs to be considered in the analysis, as shown in Fig. 4. A total of 11 × 11 nodes are arranged in the domain as shown in Fig. 5.
fig. 5. Nodal distribution and Delaunay triangles of a quarter of the cylinder.
The radial displacements of nodes at x2 = 0 when t = 30 s using the two models are given in Figs. 6 and 8. The displacements of node (1, 0) depending on time using the two models are shown in Figs. 7 and 9. From Figs. 6– 9, we can see that the results obtained by the HNEM have precision similar with that of the NEM and are in accord with the analytical solution.
Fig. 9. Time variation of radial displacement at point (1, 0) obtained using the Kelvin model.
4.2. Cantilever beam with a distributed load
A cantilever beam with a distributed load is shown in Figs. 10. The length of the beam is L = 8 m and the width is H = 2 m. The distributed load is q = 30 kPa, and the volume deformation is elastic with E = 1.0× 108 Pa, ν = 0.25. The rheological property of shear deformation satisfies the Kelvin model, and the corresponding parameters are G = 2× 108 Pa and η = 6 × 108 Pa · s. Regardless of the body force, the cantilever beam can be discussed in terms of a plane stress problem.
An array of 13 × 5 nodes is distributed on the beam, as shown in Figs. 11. The vertical and the horizontal displacements of the nodes at the left side of the cantilever beam are set as zero when calculated. Under the same nodal distribution, the numerical results obtained using the HNEM and the NEM are shown in Figs. 12– 15.From Figs. 12– 14, we can see that the nodal displacements obtained using the HNEM are in good agreement with those obtained using the NEM and the FEM; while the nodal stresses obtained using the HNEM are more precise than those from the NEM, as shown in Figs. 15.
Fig. 15. Horizontal normal stress at x1 = 2 m when t = 20 s.
4.3. Beam subjected to simple bending
Consider a beam subjected to simple bending, as shown in Figs. 16. The length is L = 10 m and the width is H = 4 m. The volume deformation is elastic and E = 1.0 × 106 Pa, ν = 0.3. The three-parameter model is adopted to describe the rheological property of the shear deformation, and the parameters are G1 = 5.0 × 105 Pa, G2 = 1.0 × 106 Pa, and η = 2.0 × 106 Pa · s. This example is discussed as a plane stress problem without considering the body force.
By taking advantage of the symmetry of the model, only a quarter of the model is considered in this analysis. At axes of x1 and x2, the horizontal displacements of the nodes are zero, and the vertical displacements of the nodes at the origin are set to zero for the restriction of rigid body displacements. As shown in Figs. 17, 13 × 7 nodes are arranged on a quarter of the model.
Fig. 17. Nodal distribution and Delaunay triangles in a quarter of the beam.
Under the same nodal distribution, the nodal displacements obtained using the HNEM and the NEM are described in Figs. 18– 20, and the nodal stresses are shown in Figs. 21. From Figs. 18– 20, we can see that the nodal displacements obtained using the HNEM are in good agreement with those obtained using the NEM and the FEM. The nodal stresses from the HNEM agree well with those from the FEM, while the precision of the nodal stresses is higher than that from the NEM, as shown in Figs. 21.
Fig. 21. Horizontal normal stress at x1 = 2.5 m when t = 30 s.
4.4. Plate with a central hole under an axial uniform tension
A plate with a central hole subjected to an axial uniform tension is shown in Figs. 22. The length of the plate is L = 10 m, the width is H = 4 m, and the radius of the hole is R = 1 m. The tensile load is q = 1 kPa, and the volume deformation is elastic with elasticity modulus E = 1.0× 108 Pa, and Poisson ratio ν = 0.25. Using the three-parameter viscoelastic model, the corresponding parameters are G1 = 8.0 × 107 Pa, G2 = 1.0 × 108 Pa, and η = 2.0 × 108 Pa · s. Regardless of the body force, the plate can be discussed as a plane stress problem.
Fig. 22. Plate with a central hole under an axial uniform tension.
There are 200 nodes distributed on the plate, as shown in Figs. 23. The vertical and the horizontal movements of the nodes at the left side of the plate are restricted in the procedure. The nodes near the hole are numerous so as to improve the calculation precision. The numerical results obtained by the HNEM and the NEM are shown in Figs. 24– 26. From Figs. 24 and 25, we can see that the nodal displacements solved by the HNEM fit well with those from the NEM and the FEM. In Figs. 26, it is shown that the nodal stresses gained by the HNEM are close to those by the FEM and have higher precision than those by the NEM.
Fig. 26. Horizontal normal stresses of nodes at x1 = − 2.5 m when t = 30 s.
5. Conclusion
In this paper, according to the properties of linear viscoelastic materials, the HNEM for linear viscoelastic materials under the creep condition is presented. The HNEM solves the problem that the nodal stresses cannot be obtained directly by the NEM.
Numerical examples show that the HNEM can obtain precise nodal displacements of linear viscoelastic materials, which illustrates the effectiveness of the method proposed in this paper. Under the same nodal distribution, the HNEM has a precision of nodal displacements similar to that of the NEM, while the precision of nodal stresses using the HNEM is higher than that using the NEM by the MLS approximation, which shows the superiority of this method.
ShahrokhabadiS, ToufighM M and GholizadehR2010GeoShanghai International Conference-Geoenvironmental Engineering and GeotechnicsJune 3–5, 2010Shanghai, China245DOI:10.1061/41105(378)34[Cited within:1]
SibsonR1980Mathematical Proceedings of the Cambridge Philosophical SocietyJanuary, 1980Cambridge, England 151DOI:10.1017/S0305004100056589[Cited within:1]
BrownJ L1994An international conference on curves and surfaces on Wavelets, images, and surface fitting1994Chamonix-Mont-Blanc, France67[Cited within:1]
Abstract In this paper, the element-free Galerkin (EFG) method and improved moving least-squares (IMLS) approximation are combined. An improved FEG (IEFG) method for two-dimensional elasticity is discussed, and the coupling of the IEFG method and the boundary element method (BEM) is presented. In the IMLS approximation, an orthogonal function system with a weight function is used as the basis function. The IMLS approximation has greater computational efficiency and precision than the existing moving least-squares (MLS) approximation, and does not lead to an ill-conditioned system of equations. There are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method that is formed with the IMLS approximation fewer nodes are selected in the entire domain than are selected using the conventional EFG method. Hence, the IEFG method should result in a higher computing speed. Based on the IMLS approximation and the IEFG method, a direct coupling of the IEFG method and the BEM is discussed for two-dimensional elasticity problems, and the corresponding formulae of the coupled method are obtained. The coupled method does not need a new sub-domain between the IEFG method and the BEM sub-domains. Selected numerical examples are solved using the coupled method.
... [1] In addition to the EFG method,[2,3] there are many other meshless methods, such as the improved element-free Galerkin method,[4#cod#x2013 ...
Abstract This paper presents an improved moving least-squares (IMLS) approximation in which the orthogonal function system with a weight function is used as the basis function. The IMLS approximation has greater computational efficiency and precision than the existing moving least-squares (MLS) approximation, and does not lead to an ill-conditioned system of equations. By combining the element-free Galerkin (EFG) method and the IMLS approximation, an improved element-free Galerkin (IEFG) method for two-dimensional elasticity is derived. There are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method that is formed with the IMLS approximation fewer nodes are selected in the entire domain than are selected in the conventional EFG method. Hence, the IEFG method should result in a higher computing speed. For two-dimensional fracture problems, the enriched basis function is used at the tip of the crack to give an enriched IEFG method. When the enriched IEFG method is used, the singularity of the stresses at the tip of the crack can be shown better than that in the IEFG method. To provide a demonstration, numerical examples are solved using the IEFG method and the enriched IEFG method.
In this paper, we derive an improved element-free Galerkin (IEFG) method for two-dimensional linear elastodynamics by employing the improved moving least-squares (IMLS) approximation. In comparison with the conventional moving least-squares (MLS) approximation function, the algebraic equation system in IMLS approximation is well-conditioned. It can be solved without having to derive the inverse matrix. Thus the IEFG method may result in a higher computing speed. In the IEFG method for two-dimensional linear elastodynamics, we employed the Galerkin weak form to derive the discretized system equations, and the Newmark time integration method for the time history analyses. In the modeling process, the penalty method is used to impose the essential boundary conditions to obtain the corresponding formulae of the IEFG method for two-dimensional elastodynamics. The numerical studies illustrated that the IEFG method is efficient by comparing it with the analytical method and the finite element method. (C) 2013 Elsevier Ltd. All rights reserved.
Zhang, Zan 1 ;Hao, S. Y. 2 ;Liew, K. M. 3 ;Cheng, Y. M. 4 ;
1
2013
0.0
0.0
... 7] the complex variable meshless method,[8#cod#x2013 ...
The moving least-square approximation with complex variables (MLSCV) is developed on the basis of moving least-square approximation. The advantages of MLSCV are that the approximation function of a 2-D problem is formed with 1-D basis function, and the meshless method obtained has greater computational efficiency. A meshless method with complex variables for 2-D elasticity is then presented using MLSCV, and the formulae of the meshless method with complex variables are obtained. Compared with the conventional meshless method, the meshless method with complex variables introduced in this paper has greater precision and computational efficiency. Some examples are given.
Effects of an ultra-strong magnetic field on electron capture rates for Co-55 are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of a magnetar are given. The results show that electron capture rates on Co-55 are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future study of the evolution of magnetars.
Du Jun 1 ;Li Ping-Ping 1 ;Luo Xia 1,2 ;
Effects of ultra-strong magnetic field on electron capture rates for 55 Co are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of magnetar are given. The results show that electron capture rates on 55 Co are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future studying the evolution of magnetar.
... 10] the complex variable element-free Galerkin (CVEFG) method,[11#cod#x2013 ...
1
2009
1.483
0.0
... 10] the complex variable element-free Galerkin (CVEFG) method,[11#cod#x2013 ...
Abstract Based on the complex variable moving least-squares (CVMLS) approximation and element-free Galerkin (EFG) method, the complex variable element-free Galerkin (CVEFG) method for two-dimensional elasto-plasticity problems is presented in this paper. The CVMLS approximation is an approximation method for a vector function. Under the same node distribution the meshless method based on the CVMLS approximation has higher precision than the one based on the moving least-squares (MLS) approximation. For two-dimensional elasto-plasticity problems, the Galerkin weak form is employed to obtain the equations system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional elasto-plasticity problems are obtained. Compared with the EFG method, the CVEFG method can obtain greater precision. For the purposes of demonstration, some selected numerical examples are solved using the CVEFG method.
1
2012
1.148
1.2429
1
2012
1.148
1.2429
1
2012
0.0
0.0
1
2012
1.483
0.0
1
2014
0.0
0.0
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
1
2011
1.483
0.0
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
1
2012
1.596
0.0
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
Combining the interpolation function, which has the delta function property and is constructed on the basis of radial basis functions and polynomial functions, using the local boundary integral equation method (LBIE), the local boundary integral equation method based on radial basis functions is presented for potential problem in this paper. The corresponding discrete equations are obtained. Comparing with the other meshless boundary integral equation methods, the present method has simpler numerical procedures, lower computation cost and higher accuracy. In addition, the essential boundary conditions can be implemented directly. Some numerical results to show the efficiency of the present method are given.
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
In this paper, the finite point method is used to obtain the solution of a one-d imensional inverse heat conduction problem with a source parameter, and the corr esponding discrete equations are obtained. Compared with the numerical methods b ased on mesh, the finite point method only needs the scattered nodes instead of meshing the domain of the problem. The finite point method is a meshless method in which the moving least-square approximation is used to form the meshless appr oximation functions. And the collocation method is used to discretize the govern ing partial differential equations. The finite point method has the advantages o f simpler numerical procedures, lower computation cost and arbitrary nodes. The result of a numerical example is presented to show the method is effective.
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
1
2011
1.148
1.2429
... 17] the interpolating element-free Galerkin (IEFG) method,[18, 19] the radial basis function (RBF) method,[20] the finite point method (FPM),[21, 22] the meshless local Petrov#cod#x2013 ...
On the basis of reproducing kernel particle method(RKPM), the complex variable reproducing kernel particle method (CVRKPM) is discussed. The advantage of the CVRKPM is that the correction function of a 2-D problem is formed with 1-D basis function when the shape function is obtained. Then, we apply the complex variable method to two-dimensional transient heat conduction problems. In combination with the Galerkin weak form of transient heat conduction problems, the penalty method is employed to enforce the essential boundary conditions, the CVRKPM for transient heat conduction problems is investigated and the corresponding formulae are obtained. Compared with the conventional RKPM, the CVRKPM introduced in this paper has a higher precision and a lower computation cost. Some examples given in this paper verify the effectivity of the proposed method.
Effects of an ultra-strong magnetic field on electron capture rates for Co-55 are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of a magnetar are given. The results show that electron capture rates on Co-55 are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future study of the evolution of magnetars.
Du Jun 1 ;Li Ping-Ping 1 ;Luo Xia 1,2 ;
Effects of ultra-strong magnetic field on electron capture rates for 55 Co are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of magnetar are given. The results show that electron capture rates on 55 Co are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future studying the evolution of magnetar.
Effects of an ultra-strong magnetic field on electron capture rates for Co-55 are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of a magnetar are given. The results show that electron capture rates on Co-55 are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future study of the evolution of magnetars.
Du Jun 1 ;Li Ping-Ping 1 ;Luo Xia 1,2 ;
Effects of ultra-strong magnetic field on electron capture rates for 55 Co are analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, and the electron capture rates on 10 abundant iron group nuclei at the surface of magnetar are given. The results show that electron capture rates on 55 Co are increased greatly in the ultra-strong magnetic field, by about 3 orders of magnitude generally. These conclusions play an important role in future studying the evolution of magnetar.
... 28] the meshless manifold method,[29, 30] the boundary element-free method (BEFM),[31#cod#x2013 ...
In the paper, the meshless manifold method (MMM) is utilized to analyze transient deformations in dynamic fracture. The MMM is based on the partition of unity method and the finite coverage approximation which provides a unified framework for solving problems involving both continuums and dis-continuums. The method can treat crack problem easily because the shape function is not affected by the discontinuity in the domain. For localization problems at the tip of the discontinuity, these shape functions are more effective than those used in other numerical methods. The method avoids the disadvantages of other meshless methods in which the tip of a discontinuous crack is not considered. In meshless manifold method, the finite coverage approximation is used to construct the shape functions that overcome influences of the interior discontinuities in the displacement. Consequently, the meshless manifold method has some advantages in solving the discontinuity problems when the discontinuities are complex. When the dynamic fracture mechanics is analyzed by the MMM, the weak formulation of the partial differential equation for elastic dynamics is derived from the method of weighted residuals (MWR). The discrete space of the domain is used for the MMM. The Newmark family of methods is used for the time integration scheme. At last, the validity and accuracy of the MMM are illustrated by two numerical examples of which the numerical results agree with the analytical solution.
By combining the reproducing kernel particle method (RKPM) and boundary integral equation method for elasticity, the reproducing kernel particle boundary element-free (RKP-BEF) method is presented in this paper. Formulae for the RKP-BEF method are derived. The discrete boundary integral equations of the RKP-BEF method are established by considering the numerical integral schemes and the treatment of singular integrals, and the formulae of the displacement and stress of inner points for the RKP-BEF method are given. The RKP-BEF method has a higher precision as the smoothness of the shape function of RKP method is the same as that of the kernel function, and the values of polynomials at interpolating points can be exactly reconstructed. Numerical examples are given for verifying the effectiveness and correctness of the RKP-BEF method presented in this paper.
Abstract Combining the boundary integral equation (BIE) method and improved moving least-squares (IMLS) approximation, a direct meshless BIE method, which is called the boundary element-free method (BEFM), for two-dimensional potential problems is discussed in this paper. In the IMLS approximation, the weighted orthogonal functions are used as the basis functions; then the algebra equation system is not ill-conditioned and can be solved without obtaining the inverse matrix. Based on the IMLS approximation and the BIE for two-dimensional potential problems, the formulae of the BEFM are given. The BEFM is a direct numerical method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily; thus, it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
... 40] and the improved boundary element-free method (IBEFM) ...
The interpolating moving least-squares (IMLS) method is discussedfirst in this paper. And the formulae of the IMLS method obtained byLancaster are revised. Then on the basis of the boundaryelement-free method (BEFM), combining the boundary integral equation(BIE) method with the IMLS method, the improved boundaryelement-free method (IBEFM) for two-dimensional potential problemsis presented, and the corresponding formulae of the IBEFM areobtained. In the BEFM, boundary conditions are applied directly, butthe shape function in the MLS does not satisfy the property ofthe Kronecker δ function. This is a problem of the BEFM, andmust be solved theoretically. In the IMLS method, when the shape functionsatisfies the property of the Kronecker δ function, then theboundary conditions, in the meshless method based on the IMLSmethod, can be applied directly. Then the IBEFM, based on the IMLSmethod, is a direct meshless boundary integral equation method inwhich the basic unknown quantity is the real solution of the nodalvariables, and the boundary conditions can be applied directly andeasily, thus it gives a greater computational precision. Somenumerical examples are presented to demonstrate the method.
School of Computer Engineering and Science, Shanghai University, Shanghai 200072, China
The interpolating moving least-squares (IMLS) method is discussedfirst in this paper. And the formulae of the IMLS method obtained byLancaster are revised. Then on the basis of the boundaryelement-free method (BEFM), combining the boundary integral equation(BIE) method with the IMLS method, the improved boundaryelement-free method (IBEFM) for two-dimensional potential problemsis presented, and the corresponding formulae of the IBEFM areobtained. In the BEFM, boundary conditions are applied directly, butthe shape function in the MLS does not satisfy the property ofthe Kronecker δ function. This is a problem of the BEFM, andmust be solved theoretically. In the IMLS method, when the shape functionsatisfies the property of the Kronecker δ function, then theboundary conditions, in the meshless method based on the IMLSmethod, can be applied directly. Then the IBEFM, based on the IMLSmethod, is a direct meshless boundary integral equation method inwhich the basic unknown quantity is the real solution of the nodalvariables, and the boundary conditions can be applied directly andeasily, thus it gives a greater computational precision. Somenumerical examples are presented to demonstrate the method.
An effective algorithm for the finite-horizon linear quadratic continuous terminal control is proposed. It is the combination of existing continuous soft and hard terminal control. We apply the algorithm to the automatic landing control of OH-6A helicopters. Numerical demonstration shows that, whether noise exists or not, the algorithm has less computation time and less feedback gains than existing hard terminal control while generally achieving the same terminal accuracy. The optimization problem which represents the hard terminal control can be addressed by sweep method and transit matrix method. It is also discovered that transit matrix method is a crucial point for improving terminal accuracy.
Xia, X. 1 ;Xu, Z. 1 ;
... 44] ...
1
1995
38.597
0.0
... [45] this method has been successfully applied in solid mechanics and fluid mechanics ...
1
1998
2.056
0.0
... [46] successfully solved partial differential equations of the elliptic type in solid mechanics by using this method ...
1
2000
2.432
0.0
... [47] researched the NEM for elastokinetics problems ...
1
2002
2.056
0.0
... [48] extended the NEM based on the concept of an #cod#x0061 ...
1
2003
0.0
0.0
... [49] applied the NEM in the analysis of fluids ...
1
2005
2.056
0.0
... [50] studied the NEM for large deformation problems ...
1
2007
0.0
0.0
... [51] Daneshmand et al ...
1
2010
0.0
0.0
... [52] used the NEM to calculate the two-dimensional flow under sluice gates ...
1
2010
0.0
0.0
... [53] used the NEM to analyze seepage problems of underground water ...
1
2013
0.0
0.0
... [54] ...
1
1980
0.0
0.0
... )except at the nodes where they are C0,[55#cod#x2013 ...
1
1990
0.0
0.0
1
1994
0.0
0.0
... 57] we cannot obtain the nodal stresses ...
1
1994
0.0
0.0
... [58] Except for this method, the hybrid natural element method (HNEM) can also obtain the nodal stresses, which has been successfully applied in elasticity problems and elastoplasticity problems ...
1
2012
0.0
0.558
... [59,60] The HNEM adopts displacement interpolation and stress interpolation respectively using natural neighbor interpolation,[61,62] and the corresponding formulae are obtained according to the Hellinger#cod#x2013 ...
1
2014
1.383
0.0
... [59,60] The HNEM adopts displacement interpolation and stress interpolation respectively using natural neighbor interpolation,[61,62] and the corresponding formulae are obtained according to the Hellinger#cod#x2013 ...
2
1997
0.408
0.0
... [59,60] The HNEM adopts displacement interpolation and stress interpolation respectively using natural neighbor interpolation,[61,62] and the corresponding formulae are obtained according to the Hellinger#cod#x2013 ...
... The shape functions for the HNEMThe natural neighbor interpolation (NNI)[61,62] in the NEM is adopted to construct the shape functions for the HNEM ...
3
2001
0.0
0.0
... [59,60] The HNEM adopts displacement interpolation and stress interpolation respectively using natural neighbor interpolation,[61,62] and the corresponding formulae are obtained according to the Hellinger#cod#x2013 ...
... The shape functions for the HNEMThe natural neighbor interpolation (NNI)[61,62] in the NEM is adopted to construct the shape functions for the HNEM ...
... [62] ...
1
2011
0.0
0.0
... [63,64] ...
1
1982
2.056
0.0
... [63,64] ...
1
2006
0.0
0.0
... [65] The constitutive relationship of viscoelastic materials correlates with time ...
3
1975
0.0
0.0
... [66] The solutions of complex viscoelastic problems in engineering mainly rely on numerical calculations, for which the meshless methods have been widely used ...
... In the polar coordinates, the analytical solution of the radial displacement of a Maxwell material is given by[66] ...
... In the polar coordinates, the analytical solution of the radial displacement of a Kelvin material is given by[66](67)
Abstract In this paper, element free Galerkin (EFG) method is combined with a precise algorithm in the time domain for solving viscoelasticity problems. By expanding variables at discretized time intervals, an initial boundary value problem is converted into a series of recurrent boundary value problems which can be conveniently solved by EFG/EFG-FE with a self-adaptive computing process. There is no requirement of iteration for the solution of non-linear cases. Satisfactory numerical results are obtained for both static and dynamical viscoelasticity problems.
... [67] Using the meshless local Petrov#cod#x2013 ...
Abstract A meshless method based on the local Petrov–Galerkin approach is proposed for stress analysis in two-dimensional (2D), anisotropic and linear elastic/viscoelastic solids with continuously varying material properties. The correspondence principle is applied for non-homogeneous, anisotropic and linear viscoelastic solids where the relaxation moduli are separable in space and time. The inertial dynamic term in the governing equations is considered too. A unit step function is used as the test functions in the local weak-form. It leads to local boundary integral equations (LBIEs). The analyzed domain is divided into small subdomains with a circular shape. The moving least squares (MLS) method is adopted for approximating the physical quantities in the LBIEs. For time-dependent problems, the Laplace-transform technique is utilized. Several numerical examples are given to verify the accuracy and the efficiency of the proposed method.
... [68, 69] analyzed anisotropic and linear viscoelastic solids, and discussed quasistatic and transient dynamic problems in two-dimensional linear viscoelastic media ...
A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of quasi-static and transient dynamic problems in two-dimensional (2-D) nonhomogeneous linear viscoelastic media. A unit step function is used as the test functions in the local weak form. It is leading to local boundary integral equations (LBIEs) involving only a domain-integral in the case of transient dynamic problems. The correspondence principle is applied to such nonhomogeneous linear viscoelastic solids where relaxation moduli are separable in space and time variables. Then, the LBIEs are formulated for the Laplace-transformed viscoelastic problem. The analyzed domain is covered by small subdomains with a simple geometry such as circles in 2-D problems. The moving least squares (MLS) method is used for approximation of physical quantities in LBIEs.
1.Institute of Construction and Architecture, Slovak Academy of Sciences, 84503 Bratislava, Slovakia 2.Department of Civil Engineering, University of Siegen, D-57068 Siegen, Germany 3.Institute of Applied Mechanics, Graz University of Technology, A-8010 Graz, Austria
... [68, 69] analyzed anisotropic and linear viscoelastic solids, and discussed quasistatic and transient dynamic problems in two-dimensional linear viscoelastic media ...
1
2007
0.0
0.0
... [70] Canelas and Sensale studied a boundary knot method for harmonic elastic and viscoelastic problems ...
Abstract The boundary knot method is a promising meshfree, integration-free, boundary-type technique for the solution of partial differential equations. It looks for an approximation of the solution in the linear span of a set of specialized radial basis functions that satisfy the governing equation of the problem. The boundary conditions are taken into account through the collocation technique. The specialized radial basis function for harmonic elastic and viscoelastic problems is derived, and a boundary knot method for the solution of these problems is proposed. The completeness issue regarding the proposed set of radial basis functions is discussed, and a formal proof of incompleteness for the circular ring problem is presented. In order to address the numerical performance of the proposed method, some numerical examples considering simple and complex domains are solved.
... [71] A complex variable element-free Galerkin method for linear viscoelasticity problems under the creep condition was analyzed by Cheng et al ...
1
2012
1.148
1.2429
... [72] Three-dimensional viscoelasticity problems were analyzed via the improved element-free Galerkin method by Peng et al ...
1
2014
1.596
0.0
... [73] ...
1
1977
0.0
0.0
... [74] The natural neighbor nodes of x are 1, 5, 6, and 7, as shown in Fig ...