Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
See all formats and pricing
More options …
Volume 16, Issue 1


Volume 13 (2015)

A note on the uniqueness of 2D elastostatic problems formulated by different types of potential functions

José Luis Morales Guerrero
  • Department of Structures and Construction, Universidad Politécnica de Cartagena (UPCT), Campus Muralla del Mar, ETSII, 30202, Cartagena, Spain
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Manuel Cánovas Vidal
  • Corresponding author
  • Department of Metallurgical and Mining Engineering, Universidad Católica del Norte, Avda. Angamos 0610, Antofagasta, Chile
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ José Andrés Moreno Nicolás / Francisco Alhama López
Published Online: 2018-05-04 | DOI: https://doi.org/10.1515/phys-2018-0029


New additional conditions required for the uniqueness of the 2D elastostatic problems formulated in terms of potential functions for the derived Papkovich-Neuber representations, are studied. Two cases are considered, each of them formulated by the scalar potential function plus one of the rectangular non-zero components of the vector potential function. For these formulations, in addition to the original (physical) boundary conditions, two new additional conditions are required. In addition, for the complete Papkovich-Neuber formulation, expressed by the scalar potential plus two components of the vector potential, the additional conditions established previously for the three-dimensional case in z-convex domain can be applied. To show the usefulness of these new conditions in a numerical scheme two applications are numerically solved by the network method for the three cases of potential formulations.

Keywords: Papkovich-Neuber representation; uniqueness solution; network simulation method; linear elasticity; classical linear elasticity

PACS: 02.30.Em; 02.60.Lj; 02.70.-c; 04.20.Ex; 46.25.-y

1 Introduction

In linear elasticity, the equilibrium equation in terms of displacement, in absence of body forces, is that of Navier equation, whose analytical solution can be obtained by means of potential formulations which simplify the elastic problem by yielding uncoupled governing equations in terms of some potential unknowns. These alternative formulations have allowed to solve a great variety of interesting problems such as Boussinesq and Kelvin [1], among others. As regards numerical solutions, classical methods, such as finite-difference, derive solutions from practical problems based, for example, on the potential formulations of the Airy’s stress function – for 2D elasticity problems – or the Prandtl’s function – for torsion in prismatic bars. Wang et. al. [2] resume the set of potential formulations related to linear elasticity field and their applications. Recently, Weisz-Patrault et al. [3] have presented a new and compact representation of the elastic problem using only two quaternionic-valued monogenic potentials that can be applied to 3D problems. Among them, Papkovich-Neuber representation [4, 5] is extensively used. This work deals with singular aspects of this formulation referred to their uniqueness solutions, a complicated subject partially investigated in the last years [6,7,8]. The study of uniqueness of a potential solution can be useful from the point of view of a numerical analysis or if we want to connect analytical solutions applied to different parts of a domain.

For the case of Papkovich-Neuber (PN) representation, under certain conditions, the number of potential functions (one scalar φ0, and three components of the vector φx, φy φz.) of the general PN solution can be reduced to a derivate solution with no lack of completeness. Eubanks and Sternberg [6] studied the necessary conditions to delete either the scalar function or one of the rectangular components of the vector potential function. In addition, they considered the derived PN solution with the scalar and z-component of the vector potential for the axisymmetric case – named the Boussinesq solution [9].

The completeness of a derived PN solution ensures the suitability of the formulation for the problem under study but new additional conditions can be required to ensure the uniqueness in terms of potential functions as primary unknown. The terms completeness (i.e. the set of displacement functions in the derivate formulation is able to represent an arbitrary deformation of the elastic body) and uniqueness (i.e. the derived displacement potentials are unique for the given deformation) are used here in the sense given by Tran-Cong [8].

Stippes [7] was the first author to deal with the subject of uniqueness but always referring to the derived PN solution that results from deleting the scalar potential. This author noted that Eubanks and Sternberg [6] did not make any reference to this subject, even though the uniqueness condition was implicit in their demonstration of completeness. In a footnote, Stippes [7] wrote on the importance of the uniqueness conditions in reference to a numerical solution, but with arguments only valid for the case of deleting the scalar component of the potential representation. Twenty five years later, Tran-Cong [8], who did not apply his conclusions to particular problems as we do here, came back to the question and, based on the study of the null displacement field problem, proved the necessity of additional conditions for derived PN solutions other than those studied by Stippes [7]; more particularly, the case of deleting the rectangular z-component of the vector function when the problem is z-convex, as well as the case of Boussinesq solution. The uniqueness approach of Tran-Cong [8] completes the work of Eubanks and Sternberg [6] and Stippes [7]. Recently, Morales et al. [10] have successfully deduced results for axisymmetric problems in the Boussinesq solution, applying uniqueness conditions different from those deduced by Tran-Cong [8].

For the 2D elasticity case of plane strain in rectangular coordinates (the extension to plane stress is immediate), the complete PN solution is directly reduced to three potential functions. Choosing the xy-plane as the domain, and thus making the z-component of the vector potential zero, the complete PN solution is defined by three potential functions: φ0, φx, φy. The only possible derived PN functions (maintaining the scalar potential), are two: φ0, φx and φ0, φy. Golecki [11] suggested that in the case of multiply connected domains, the derived solutions can be multivalued, while Wang and Wang [12] demonstrated that the complete PN solution always remains single valued.

In this work, we investigate the uniqueness conditions for the 2D elasticity problems in rectangular coordinates, using as the main unknowns the potential functions of the complete and derived PN solutions. For the derived formulations, formed by two potentials, it is demonstrated that the solution is unique fixing two additional conditions while for the completed formulation, uniqueness additional conditions can be derived from the tridimensional z-convex case studied by Tran-Cong [8].

To demonstrate the suitability of the proposed conditions for the derived formulations, two applications are solved: a rectangular domain in which the displacement field is zero anywhere, and a domain whose boundary supports an arbitrary load distribution that gives rise to a nonzero displacement field. In addition, for the complete formulation, a particular solution for the above applications that provides the same elastic solutions is also shown. Applications are numerically solved by using the network method [13], a technique already successfully applied to elasticity by Morales et al. [14] and other engineering problems [15,16,17,18]. A detailed description of the application of the network model for the numerical solution of the elasticity problem formulated by potential functions can be found in Morales et al. [19]. Solutions are compared with those from the Navier formulation. The network models are generated by using the EPSNET_10 © software [20] and run in the code PSpice © [21].

2 Governing equations and boundary conditions

In term of displacements (u) Navier equation is written as:


whereλ and μ are the Lamé’s constants [1]. Using PN formulation, Eq. (1) is represented by four potential functions, named displacement potentials, in the form




where φ is a harmonic vector potential, φ0 a harmonic scalar potential, R - the position vector and ν - the Poisson ratio. Eqs. (2a, 2b, 2c) are a general solution of Eq. (1) and have proved to be complete for the general case [22].

In rectangular coordinates, the completed PN solutions for the 2D-elasticity plane strain case reduce the displacement field from Eq. (2a) to



with uz = ∂/∂z = 0, ux = ux(x,y) and uy = uy(x,y).

The corresponding governing Eq. (2b, 2c) can be written as




and finally, for the more general case, the mixed physical boundary conditions can be expressed in the form:


Su denotes the boundary surface points where the displacement, uib, is prescribed, while St refers to points where the traction values, tib, are imposed. σij is the stress field and nj the outer normal vector on St. Note that S = St+ Su represents the complete boundary surface, Figure 1.

Scheme of the boundary conditions: tractions and displacements [14]
Figure 1

Scheme of the boundary conditions: tractions and displacements [14]

The expressions that relate stress and potential functions, which are required for implementing the boundary condition (5), are [1]:




For each of the derived solutions, the resulting governing equations and boundary conditions are obtained by omitting one of the rectangular components, φx or φy, in Eqs. (3-6). Note that, since the physical boundary conditions for a given problem are unique, the problem can be formulated by different sets of PDEs, two governing equations for each derived solution or three governing equations for the complete PN solution.

3 Uniqueness of the potential PN solutions in 2D elasticity problems

For the following study, it is convenient to express the general PN solution (2a) as


with α = 4(1-ν). The new harmonic vector and scalar potentials, Ψ and φ, respectively, are related with φ and φ0 by the equations φ = 2μαΨ andφ0 = 2μφ.

According to Stippes [7] and Tran-Cong [8], the potential representation (7) is unique only with some additional conditions. One way to find them is to study the form of the potential functions corresponding to a null displacement field.

3.1 Uniqueness of derived PN solutions

Let us consider the derived solution of the completed PN representation reduced to the scalar potential and the y-component of the vector potential (φ, Ψy). Using the nomenclature of Eq. (7), the vector form of the displacement field is


For a zero displacement field (u = 0), representation (8) satisfies the equation


Applying the curl to this equation and taking into account the properties of this operator gives: ∂Ψy/∂z = ∂Ψy/∂x = 0, so that Ψy = Ψy(y). Also, applying the divergence operator and using the above result and the expressions ∇2 Ψ = 0 and ∇2 φ = 0, it is deduced that ∂Ψy/∂y = 0. In short, for solution (8) to represent a zero displacement field, the non-zero component of the vector potential must be an arbitrary constant,


Introducing this result in Eq. (9) yields


From the first and last addends of Eq. (11) we obtain ∂φ/∂x = ∂φ/∂z = 0; this implies that φ = g(y). Substituting this result in the second addend, the necessary condition (referring to the scalar potential φ) for the derived PN solution (8) to represent a zero displacement field is that φ is a plane parallel to the x-axe:


Fixing the arbitrary constants k1 and k2, the potential solution (8) is unique for a zero displacement field. By extension, it is also unique for any other displacement field represented by the derived PN solution (8). Figure 2 represents equations (10) and (12): two planes defined only by two constants, k1 and k2, since the slope of the z plane is determined by the expression k1 (α-1) = tanβ.

Derived PN solution (φ, Ψy) for a null displacement field
Figure 2

Derived PN solution (φ, Ψy) for a null displacement field

The uniqueness conditions can be expressed in the form of a Dirichlet condition


with a and b two arbitrary points of the domain, and two constants: c1 and c2, also arbitrarily chosen. Equation (13) is a suitable way of fixing constants kl and k2 and makes solutions (10) and (12) unique: k1 = c2 and k2 = c1-c2 (α-1) ay. Note that the notation of the potential functions introduced in (7), together with the arbitrability of choice constants in (13), let us use φ0(a) = c1 and φy(b) = c2 as additional conditions.

Other possibilities for establishing the uniqueness conditions are: i) to fix values of the scalar potential φ at two different points of the domain, φ(a) = c1 and φ(b) = c2, and ii) to fix a value of the scalar potential φ at one point and its derivative at the same point, φ (a) = cl (Dirichlet condition) and ∂φ/∂y|a = c2 (Neumann condition).

As regards the other derived solution (φ, Ψx), following similar reasoning, the general form of the potential solution for the null displacement is Ψx = k1 and φ = k1 (α-1)x + k2. The resulting uniqueness condition can be similarly expressed as: φ0(a) = c1 and φx(b) = c2, with a, b, c1 and c2 have the same meaning as in condition (13).

3.2 The complete PN solution

The complete PN solution for 2D elasticity, using the nomenclature of (7) is


Following the former subsection, the starting point for the uniqueness study is the zero displacement field:


This equation is equivalent to the null displacement field condition studied by Tran-Cong [8] for the three-dimensional case of z-convex domains. As a consequence, the uniqueness conditions are also the same. For the case of plane deformation, and by extension of plane tension, the uniqueness conditions according to Tran-Cong can be set in several ways, being the most suitable: i) fixing the value of Ψx at the 2D boundary, Figure 1, by means of a continuous arbitrary function, together with φ(a) = c1 and Ψy(b) = c2, or ii) fixing the value of Ψy at the 2D boundary by means of a continuous arbitrary function, together with φ(a) = c1 and Ψx(b) = c2. As in the former section, a and b are two arbitrary points of the domain and cl and c2 two constants, also arbitrarily chosen.

4 Applications

4.1 Rectangular plate with zero displacement field

Figure 3 shows the physical model (a) and the 10 × 10 grid (b) for a 2D domain under plane strain. Rigid body movement is restricted by setting zero the perpendicular displacements over the whole boundary unb = 0, Figure 3a; in addition, the tangential traction is zero in the whole boundary ttb = 0. The choice of this grid size is not relevant since we hope a null displacement field. Values of the geometry parameters are L = 1000 mm, H =1000 mm, while elastic constants are E (Young modulus) = 210 GPa and ν (Poisson ratio) = 0.3. Additional constants for the two derived solutions, section 3.1, are detailed in Table 1 using the grid reference of Figure 3b. In this figure additional conditions are plotted using the following code: circles, squares and diamonds in red for the potentials φ0, φx and φy, respectively. Note that the values of these additional conditions are not zero because we wish to show a general form of the potential surfaces – a zero value for these constants is a particular value that provides a null (trivial) solution for the potentials. In addition, Table 1 shows additional conditions related to the complete formulation needed to find a particular solution by means of the software EPSNET_10 © [20].

Table 1

Additional conditions for derived and completed PN solutions referring to cells in a 10 × 10 grid, Figure 3b

Physical model (a) and grid (b) for a rectangular plate under simple boundary conditions for a null displacement field
Figure 3

Physical model (a) and grid (b) for a rectangular plate under simple boundary conditions for a null displacement field

Simulation of the three models, one for each solution type, in the code PSpice © [21] using the software EPSNET_10© for preprocessing and postprocessing, gives the following results: Figure 4 shows the displacement field using a large scale to test the null displacement results while Figure 5 shows the potential solutions. Potentials of the derived solutions, Figures 5(a) and (b), agree with the theoretical shape of Figure 2 for the null displacement field. Potential results of the complete PN solution (see Figure 5c) are equivalent to those of derived solution (φ0, φy) due to the extra condition φx|∂S = 0. Changing this condition by other arbitrary continuous linear function, as indicated in the Table 2, new potential results are obtained (see Figure 6) being the displacement solution the same.

Solution of the displacement field by PSpice © and ESPNET_10 ©: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0. Vy) and c) completed PN solution (φ0, φx. φy)
Figure 4

Solution of the displacement field by PSpice © and ESPNET_10 ©: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0. Vy) and c) completed PN solution (φ0, φx. φy)

Potential solutions: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx, φy)
Figure 5

Potential solutions: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx, φy)

Potential solutions for the completed PN solution: a) φ0, b) φx, c) φy. The uniqueness conditions are shown in Table 2
Figure 6

Potential solutions for the completed PN solution: a) φ0, b) φx, c) φy. The uniqueness conditions are shown in Table 2

Table 2

Additional conditions for completed PN solutions referring to cells in a 10 × 10 grid (Figure 3b)

This application demonstrates that, in accordance with the uniqueness conditions mentioned in Section 3.1 and the results given in Figure 2, the solution for the potential functions using the formulations φ0, φx and φ0, φy has the form of constant planes for the vector component and inclined planes for the scalar components. As regards the formulation for φ0, φx and φy, a change at the boundary conditions for one of components of the vector potential function, retaining the fixed values of the other component as well as the scalar potential, provides solutions qualitatively different (with appreciable curvature in some regions) to the former, giving rise to the same solution for the displacements (Tables 1 and 2 and Figures 5c and 6).

4.2 Rectangular plate under arbitrary loads

The new boundary conditions referring to load (p = 100 MPa) and displacements are shown in Figure 7a. These are: zero normal displacement and zero tangential traction on the bottom and left boundaries, zero tangential and normal traction on the top boundary and zero tangential traction and specified normal traction on the right boundary. Geometry and elastic parameters of the material are the same as in the above problem. In this case, a more refined grid of 20 × 20 was used for a better resolution. In addition, all the additional conditions were set up to zero and placed as indicated in Figure 7b and Table 3. Figure 8 shows the displacement field of the domain for all the potential solutions as well as the Navier solution [14], while Figure 9 depicts the potential solutions whose forms depend on the physical boundary conditions (loads and displacements) and the arbitrary values of the additional conditions (for an only numerical solution of the potential functions). For each potential formulation, the harmonic functions are smooth surfaces which curve slightly at their boundaries; although these are very different potential surfaces, the displacement solutions for each one are the same.

Physical model (a) and grid (b) for a rectangular plate under arbitrary boundary conditions
Figure 7

Physical model (a) and grid (b) for a rectangular plate under arbitrary boundary conditions

Solutions of the non-zero displacement field obtained by ESPNET_10 ©: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx,φy) and d) Navier solution
Figure 8

Solutions of the non-zero displacement field obtained by ESPNET_10 ©: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx,φy) and d) Navier solution

Potential solutions for a non-zero displacement field: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx, φy)
Figure 9

Potential solutions for a non-zero displacement field: a) Derived PN solution (φ0, φx), b) derived PN solution (φ0, φy) and c) completed PN solution (φ0, φx, φy)

Table 3

Additional conditions for derived and completed PN solutions referring to cells in a 20 × 20 grid, Figure 7b

To see the influence of the grid size in the numerical results of the potential formulation solved by network method and to compare them with a common formulation of the elastic problem solved by a standard numerical code such as FEM, finite element method [23], Table 4 is presented. For the potential formulation (φ0, φx ), this table shows the maximum displacement located at the point of maximum load, right-bottom corner, as well as the comparisons with FEM solution for a same grid. As expected, deviations diminish as grid size increases and are always less than 1%.

Table 4

Displacements at the point with maximum load and comparison with FEM code solution for different grid size. Deviations are shown between brackets

5 Conclusions

Based on the study of the null displacement field, additional conditions must be imposed to obtain a unique solution for the 2D planar elasticity problem formulated by the two potential solutions derived from the Papkovich- Neuber representation. These are defined by a scalar potential and one of the components of the vector potential. For each solution, two additional conditions are required. One way for implementing these conditions is to fix an arbitrary value of scalar potential at one arbitrary point of the domain and the same for the component of the vector potential. Other possibilities for implementing these conditions (Dirichlet or Neumann types) are also proposed. For the complete Papkovich-Neuber formulation, formed by three unknown potential functions, it has been demonstrated that the additional conditions proposed by Tran-Cong [8] for the z-convex domains can also be applied. Again, there are many different forms for implementing these additional conditions. To demonstrate the reliability of these results two applications for each type of potential solutions were numerically solved by network method.


  • [1]

    Barber J.R., Elasticity, Solid Mechanics and its Applications 172, 3rd ed. Springer-Verlag, New York, 2010 Google Scholar

  • [2]

    Wang M.Z., Xu B.X., Gao C.F., Recent General Solutions in Linear Elasticity and Their Applications, Appl. Mech. Rev., 2008, 61, 030803-1/20 Web of ScienceGoogle Scholar

  • [3]

    Weisz-Patrault D., Bock S., Gürlebeck, K., Three-dimensional elasticity based on quaternion-valued potentials, Int. J. Solids Struct., 2014, 51 (19-20), 3422-3430. Web of ScienceCrossrefGoogle Scholar

  • [4]

    Papkovich P.F., Solution Générale des équations différentielles fondamentales ďélasticité exprimée par trois fonctions harmoniques, Compt. Rend. Acad. Sci., 1932, 195, 513-515. Google Scholar

  • [5]

    Neuber H., Ein neuer Ansatz zur Lösung räumlicher Probleme der Elastizitätstheorie, Z. Angew. Math. Mech., 1934, 14, 203-212. CrossrefGoogle Scholar

  • [6]

    Eubanks R.A., Sternberg E., On the completeness of the Boussinesq-Papkovich Stress Functions, J. Rat. Mech. Anal., 1956, 5, 735-746. Google Scholar

  • [7]

    Stippes M., Completeness of the Papkovich Potentials, Quart. Appl. Math., 1969, 26, 477-483. CrossrefGoogle Scholar

  • [8]

    Tran-Cong T., On the Completeness and Uniqueness of the Papkovich-Neuber and the Non-axisymmetric Boussinesq, Love, and Burgatti solutions in General Cylindrical Coordinates, J. Elast., 1995, 36, 227-255. CrossrefGoogle Scholar

  • [9]

    Boussinesq J., Application des potentiels a ľétude de ľéquilibre et du mouvement des solides élastiques, Gauthier-Villars, Paris, 1885 Google Scholar

  • [10]

    Morales J.L., Moreno J.A., Alhama F., New additional conditions for the numerical uniqueness of the Boussinesq and Timpe solutions of elasticity problems, Int. J. Comput. Math., 2012, 89, 1794-1807. Web of ScienceCrossrefGoogle Scholar

  • [11]

    Golecki J.J., On Many Valuedness of the Neuber-Papkovich in Two-Dimensional Elasticity, Mech. Res. Commun., 1974, 1, 85-74. CrossrefGoogle Scholar

  • [12]

    Wang M., Wang M.Z., On single-valuedness of the Papkovich-Neuber general solution in two-dimensional elasticity, Mech. Res. Commun., 1990, 17, 403-408. CrossrefGoogle Scholar

  • [13]

    González-Fernández C.F., Network Simulation Method, Horno Ed. Research Signpost, Trivandrum, India, 2002 Google Scholar

  • [14]

    Morales J.L., Moreno J.A., Alhama F., Numerical solutions of 2-D linear elastostatic problems by network method, CMES-Comput. Model. Eng. Sci., 2011, 76, 1-18. Google Scholar

  • [15]

    Alarcón M., Alhama F., González-Fernández C.F., Time-dependent heat transfer in a fin-wall assembly. New performance coefficient: thermal reverse admittance., Int. J. Therm. Sci., 2002, 41, 386-395. CrossrefGoogle Scholar

  • [16]

    Castro E., García-Hernández M.T., Gallego A., Transversal waves in beams via the network simulation method, J. Sound. Vibrat., 2005, 283, 997-1013. CrossrefGoogle Scholar

  • [17]

    Anwar Bég O., Zueco J., Takhar H.S., Unsteady magnetohydrodynamic Hartmann Couette flow and heat transfer in a Darcian channel with Hall current, ion slip, viscous and Joule heating effects: Network numerical solutions, Commun. Nonlinear Sci. Numer. Simul., 2009, 14, 1082-1097. CrossrefGoogle Scholar

  • [18]

    Cánovas M., Alhama I., Trigueros E., Alhama F., Numerical simulation of Nusselt-Rayleigh correlation in Bénard cells. A solution based on the network simulation method, Int. J. Numer. Methods Heat Fluid Flow, 2015, 25(5), 986-997. CrossrefGoogle Scholar

  • [19]

    Morales J.L., Moreno J.A., Alhama F., Numerical solution of 2D elastostatic problems formulated by potential functions, Appl. Math. Modell., 2013, 37 (9), 6339-6353. CrossrefGoogle Scholar

  • [20]

    EPSNET_10 ©, 2011. Network Simulation Group, Universidad Politécnica de Cartagena, Spain. Google Scholar

  • [21]

    PSpice 6.0 ©, 1994. Microsim Corporation, 20 Fairbanks, Irvine, California 92718, USA Google Scholar

  • [22]

    Mindlin R.D., Note on the Galerkin and Papkovich stress functions. Bull. Amer. Math. Soc., 1936, 42, 373-376. CrossrefGoogle Scholar

  • [23]

    ANSYS, Inc. (2004): ANSYS Release 9.0. ANSYS, Inc., South-pointe, 275 Technology Drive, Canonsburg, PA 15317, USA Google Scholar

About the article

Received: 2017-10-23

Accepted: 2018-02-06

Published Online: 2018-05-04

Citation Information: Open Physics, Volume 16, Issue 1, Pages 201–210, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0029.

Export Citation

© 2018 José Luis Morales Guerrero et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in