The variation of new Gurtin-type region-wise variational principles results in continuous conditions, boundary conditions, all equations and relations in linear thermopiezoelectric elastodynamics. Gurtin-type region-w...The variation of new Gurtin-type region-wise variational principles results in continuous conditions, boundary conditions, all equations and relations in linear thermopiezoelectric elastodynamics. Gurtin-type region-wise variational principles comprise very important parts of linear thermopiezoelectric elastodynamics, and can fully characterize the initial-boundary-value problem in linear thermopiezoelectric elastodynamics.展开更多
The finite element method (FEM) plays a valuable role in computer modeling and is beneficial to the mechanicaldesign of various structural parts. However, the elements produced by conventional FEM are easily inaccurat...The finite element method (FEM) plays a valuable role in computer modeling and is beneficial to the mechanicaldesign of various structural parts. However, the elements produced by conventional FEM are easily inaccurate andunstable when applied. Therefore, developing new elements within the framework of the generalized variationalprinciple is of great significance. In this paper, an 8-node plane hybrid finite element with 15 parameters (PHQ8-15β) is developed for structural mechanics problems based on the Hellinger-Reissner variational principle.According to the design principle of Pian, 15 unknown parameters are adopted in the selection of stress modes toavoid the zero energy modes.Meanwhile, the stress functions within each element satisfy both the equilibrium andthe compatibility relations of plane stress problems. Subsequently, numerical examples are presented to illustrate theeffectiveness and robustness of the proposed finite element. Numerical results show that various common lockingbehaviors of plane elements can be overcome. The PH-Q8-15β element has excellent performance in all benchmarkproblems, especially for structures with varying cross sections. Furthermore, in bending problems, the reasonablemesh shape of the new element for curved edge structures is analyzed in detail, which can be a useful means toimprove numerical accuracy.展开更多
The difference of constitutive character and large deformation as to soil mass are basic questions to analyze deformational feature. According to the description method of limited deformation, the large deformation co...The difference of constitutive character and large deformation as to soil mass are basic questions to analyze deformational feature. According to the description method of limited deformation, the large deformation consolidation equations of soil mass were created and its variational principles were rigorously testified. The regionwise variational principles of consolidation theory were deduced using sub-structure continuous condition of region-wise. Quoting the method of Lagrangian multiplier operator, generalized variational principles of region-wise of large deformation consolidation in the nonconstrained condition were created and approved.展开更多
Let(X,φ) be a nonautonomous dynamical system.In this paper,we introduce the notions of packing topological entropy and measure-theoretical upper entropy for nonautonomous dynamical systems.Moreover,we establish the v...Let(X,φ) be a nonautonomous dynamical system.In this paper,we introduce the notions of packing topological entropy and measure-theoretical upper entropy for nonautonomous dynamical systems.Moreover,we establish the variational principle between the packing topological entropy and the measure-theoretical upper entropy.展开更多
This paper summarized the recent development on Herglotz’s generalized variational principle and its symmetries and conserved quantities for nonconservative dynamical systems.Taking Lagrangian mechanics,Hamiltonian m...This paper summarized the recent development on Herglotz’s generalized variational principle and its symmetries and conserved quantities for nonconservative dynamical systems.Taking Lagrangian mechanics,Hamiltonian mechanics and Birkhoffian mechanics as three research frames,we introduce Herglotz’s generalized variational principle,dynamical equations of Herglotz type,Noether symmetry and conserved quantities,and their generalization to time-delay dynamics,fractional dynamics and time-scale dynamics,and put forward some problems as suggestions for future research.展开更多
From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given....From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and the initial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.展开更多
An expression of the generalized principle of virtual work for the boundary value problem of the linear and anisotropic electromagnetic field is given. Using Chien's method, a pair of generalized variational principl...An expression of the generalized principle of virtual work for the boundary value problem of the linear and anisotropic electromagnetic field is given. Using Chien's method, a pair of generalized variational principles (GVPs) are established, which directly leads to all four Maxwell's equations, two intensity-potential equations, two constitutive equations, and eight boundary conditions. A family of constrained variational principles is derived sequentially. As additional verifications, two degenerated forms are obtained, equivalent to two known variational principles. Two modified GVPs are given to provide the hybrid finite element models for the present problem.展开更多
The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisi...The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.展开更多
In this paper, the variational principles of hydrodynamic problems for the incompressible and compressible viscous fluids are established. These principles are principles of maximum power losses. Their generalized var...In this paper, the variational principles of hydrodynamic problems for the incompressible and compressible viscous fluids are established. These principles are principles of maximum power losses. Their generalized variational principles are also discussed on the basis of Lagrangian multiplier methods.展开更多
The relations of all generalized variational principles in elasticity are studied by employing the invariance theorem of field theory. The infinitesimal scale transformation in field theory was employed to investigate...The relations of all generalized variational principles in elasticity are studied by employing the invariance theorem of field theory. The infinitesimal scale transformation in field theory was employed to investigate the equivalent theorem. Among the results found particularly interesting are those related to that all generalized variational principles in elasticity are equal to each other. Also studied result is that only two variables are independent in the functional and the stress-strain relation is the variational constraint condition for all generalized variational principles in elasticity. This work has proven again the conclusion of Prof. Chien Wei-zang.展开更多
Naturally curved and twisted closed thin-walled slender beams of composite material undergoing small strains, large displacements and rotations have been investigated, and an incomplete generalized variational functio...Naturally curved and twisted closed thin-walled slender beams of composite material undergoing small strains, large displacements and rotations have been investigated, and an incomplete generalized variational function on theory of elasticity with finite displacement is established far these beams with complete constrained boundaries at two ends. The balance equations as well as all boundary conditions concerned have been deduced from functional stationary value condition. The above-mentioned method can also be extended to other various incomplete constrained boundaries conveniently. In addition, the fundamental equations and concerned formulas in the small displacement theory of the beams can be derived by using above results.展开更多
According to generalized variational principles suitable for linear elastic incompatible displacement elements given by Professor Chien Wei-zang, using crack tip singular element and isoparametric surrounding element ...According to generalized variational principles suitable for linear elastic incompatible displacement elements given by Professor Chien Wei-zang, using crack tip singular element and isoparametric surrounding element given by the author of this paper, we will study the St. Venant's torsional bar with a radial vertical crack and compare the present computed results with the results of reference [2], The present computed results show that, using the method provided in this paper, satisfactory convergent solution can be obtained under lower degree of freedom.展开更多
The paper has proved that Hellinger-Reissuer and Hu-Washizu variational principles are but equivalent principles in elasticity by following three ways: 1) Lagrange multiplier method. The paper points out that only a n...The paper has proved that Hellinger-Reissuer and Hu-Washizu variational principles are but equivalent principles in elasticity by following three ways: 1) Lagrange multiplier method. The paper points out that only a new independent variable can be introduced when one constraint equation has been eliminated by one Lagrange multiplier, which must be expressed as a function of the original variable(s) and/or the new introduced variable after identification. In using Lagrange multiplier method to deduce Hu-Washizu principle from the minimum potential energy principle, which has only one kind of independent variable namely displacement, by eliminating the constraint equations of stress-displacement relations, one can only obtain a principle with two kinds of variables namely displacement and stress; 2) involutory transformation, with such method Hu-Washizu variational principle can be deduce directly from the Hellinger-Reissner variational principle under the same variational constraints of Stress-strain relation, and vice verse; 3)semi-inverse method, by which both of the above variational principles can be deduced from the minimum potential energy principle with tile same variational constraints. So the three kinds of variational functions in Hu-Washizu variational principle are not independent to each other,the stress-strain relationships are still its constraint conditions.展开更多
Based on generalized variational principles, an element called MR-12 was constructed for the static and dynamic analysis of thin plates with orthogonal anisotropy. Numerical results showed that this incompatible eleme...Based on generalized variational principles, an element called MR-12 was constructed for the static and dynamic analysis of thin plates with orthogonal anisotropy. Numerical results showed that this incompatible element converges very rapidly and has good accuracy. It was demonstrated that generalized varialional principles arc useful and effective in founding incompatible clement.Moreover, element MR-12 is easy for implementation since it does not differ very much from the common rectangular element R-12 of thin plate.展开更多
This paper discusses the generalized variational principles founded by the technique of Lagrangian multipliers in structural mechanics and analyzes the nonlinear statically indeterminate structures. It is assumed that...This paper discusses the generalized variational principles founded by the technique of Lagrangian multipliers in structural mechanics and analyzes the nonlinear statically indeterminate structures. It is assumed that the stress-strain relationship of the materials of structures has the form of namely, the physical equations of structures have the shape of exponential functions. Several examples are given to illustrate the statically indeterminate structures such as the trusses, beams, frames and torsional bars.展开更多
4 semi-analytical approach for the dynamic response of general thin plates whichemployes finite element discretization in space domain and a series of representation intime demain is developed on the basis of Curtin v...4 semi-analytical approach for the dynamic response of general thin plates whichemployes finite element discretization in space domain and a series of representation intime demain is developed on the basis of Curtin variational principles.The formulationof time series is also investigated so that the dynamic response of plates with arbitraryshape and boundary constraints can be achieved with adequate accuracy.展开更多
The generalized variational principles of isothermal quasi-static fluid full-filled elastic solids are established by using Variational Integral Method. Then by introducing constraints, several kinds of variational pr...The generalized variational principles of isothermal quasi-static fluid full-filled elastic solids are established by using Variational Integral Method. Then by introducing constraints, several kinds of variational principles are worked out, including five-field variable, four-field variable, three-field variable and two-field variable formulations. Some new variational principles are presented besides the principles noted in the previous works. Based on variational principles, finite element models can be set up.展开更多
Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical...Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fuid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small de- formation of the solid phase, small velocity of the fuid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fuid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles, especially Hu-Washizu type variational principles, for the initial boundary value problems of dy- namic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fuid-saturated elastic porous media, which have been discussed previously.展开更多
For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process...For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.展开更多
Based upon the Hellinger-Reissner (H-R) mixed variational principle for three-dimensional elastic bodies, the modified H-R mixed variational theorem for magnetoelectroelastic bodies was established. The state-vector e...Based upon the Hellinger-Reissner (H-R) mixed variational principle for three-dimensional elastic bodies, the modified H-R mixed variational theorem for magnetoelectroelastic bodies was established. The state-vector equation of magnetoelectroelastic plates was derived from the proposed theorem by performing the variational operations. To lay a theoretical basis of the semi-analytical solution applied with the magnetoelectroelastic plates, the state-vector equation for the discrete element in plane was proposed through the use of the proposed principle. Finally, it is pointed out that the modified H-R mixed variational principle for pure elastic, single piezoelectric or single piezomagnetic bodies are the special cases of the present variational theorem.展开更多
文摘The variation of new Gurtin-type region-wise variational principles results in continuous conditions, boundary conditions, all equations and relations in linear thermopiezoelectric elastodynamics. Gurtin-type region-wise variational principles comprise very important parts of linear thermopiezoelectric elastodynamics, and can fully characterize the initial-boundary-value problem in linear thermopiezoelectric elastodynamics.
基金the National Natural Science Foundation of China(No.11572210).
文摘The finite element method (FEM) plays a valuable role in computer modeling and is beneficial to the mechanicaldesign of various structural parts. However, the elements produced by conventional FEM are easily inaccurate andunstable when applied. Therefore, developing new elements within the framework of the generalized variationalprinciple is of great significance. In this paper, an 8-node plane hybrid finite element with 15 parameters (PHQ8-15β) is developed for structural mechanics problems based on the Hellinger-Reissner variational principle.According to the design principle of Pian, 15 unknown parameters are adopted in the selection of stress modes toavoid the zero energy modes.Meanwhile, the stress functions within each element satisfy both the equilibrium andthe compatibility relations of plane stress problems. Subsequently, numerical examples are presented to illustrate theeffectiveness and robustness of the proposed finite element. Numerical results show that various common lockingbehaviors of plane elements can be overcome. The PH-Q8-15β element has excellent performance in all benchmarkproblems, especially for structures with varying cross sections. Furthermore, in bending problems, the reasonablemesh shape of the new element for curved edge structures is analyzed in detail, which can be a useful means toimprove numerical accuracy.
文摘The difference of constitutive character and large deformation as to soil mass are basic questions to analyze deformational feature. According to the description method of limited deformation, the large deformation consolidation equations of soil mass were created and its variational principles were rigorously testified. The regionwise variational principles of consolidation theory were deduced using sub-structure continuous condition of region-wise. Quoting the method of Lagrangian multiplier operator, generalized variational principles of region-wise of large deformation consolidation in the nonconstrained condition were created and approved.
基金the National Natural Science Foundation of China (11871188, 12031019)。
文摘Let(X,φ) be a nonautonomous dynamical system.In this paper,we introduce the notions of packing topological entropy and measure-theoretical upper entropy for nonautonomous dynamical systems.Moreover,we establish the variational principle between the packing topological entropy and the measure-theoretical upper entropy.
基金supported by the National Natural Science Foundations of China (Nos. 11972241,11572212,11272227)the Natural Science Foundation of Jiangsu Province(No. BK20191454).
文摘This paper summarized the recent development on Herglotz’s generalized variational principle and its symmetries and conserved quantities for nonconservative dynamical systems.Taking Lagrangian mechanics,Hamiltonian mechanics and Birkhoffian mechanics as three research frames,we introduce Herglotz’s generalized variational principle,dynamical equations of Herglotz type,Noether symmetry and conserved quantities,and their generalization to time-delay dynamics,fractional dynamics and time-scale dynamics,and put forward some problems as suggestions for future research.
文摘From the Boltzmann's constitutive law of viscoelastic materials and the linear theory of elastic materials with voids, a constitutive model of generalized force fields for viscoelastic solids with voids was given. By using the variational integral method, the convolution-type functional was given and the corresponding generalized variational principles and potential energy principle of viscoelastic solids with voids were presented. It can be shown that the variational principles correspond to the differential equations and the initial and boundary conditions of viscoelastic body with voids. As an application, a generalized variational principle of viscoelastic Timoshenko beams with damage was obtained which corresponds to the differential equations of generalized motion and the initial and boundary conditions of beams. The variational principles provide a way for solving problems of viscoelastic solids with voids.
基金Project supported by the National Natural Science Foundation of China (No. 60304009) and the Natural Science Foundation of Hebei Province of China (No. F2005000385)
文摘An expression of the generalized principle of virtual work for the boundary value problem of the linear and anisotropic electromagnetic field is given. Using Chien's method, a pair of generalized variational principles (GVPs) are established, which directly leads to all four Maxwell's equations, two intensity-potential equations, two constitutive equations, and eight boundary conditions. A family of constrained variational principles is derived sequentially. As additional verifications, two degenerated forms are obtained, equivalent to two known variational principles. Two modified GVPs are given to provide the hybrid finite element models for the present problem.
文摘The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.
文摘In this paper, the variational principles of hydrodynamic problems for the incompressible and compressible viscous fluids are established. These principles are principles of maximum power losses. Their generalized variational principles are also discussed on the basis of Lagrangian multiplier methods.
文摘The relations of all generalized variational principles in elasticity are studied by employing the invariance theorem of field theory. The infinitesimal scale transformation in field theory was employed to investigate the equivalent theorem. Among the results found particularly interesting are those related to that all generalized variational principles in elasticity are equal to each other. Also studied result is that only two variables are independent in the functional and the stress-strain relation is the variational constraint condition for all generalized variational principles in elasticity. This work has proven again the conclusion of Prof. Chien Wei-zang.
文摘Naturally curved and twisted closed thin-walled slender beams of composite material undergoing small strains, large displacements and rotations have been investigated, and an incomplete generalized variational function on theory of elasticity with finite displacement is established far these beams with complete constrained boundaries at two ends. The balance equations as well as all boundary conditions concerned have been deduced from functional stationary value condition. The above-mentioned method can also be extended to other various incomplete constrained boundaries conveniently. In addition, the fundamental equations and concerned formulas in the small displacement theory of the beams can be derived by using above results.
文摘According to generalized variational principles suitable for linear elastic incompatible displacement elements given by Professor Chien Wei-zang, using crack tip singular element and isoparametric surrounding element given by the author of this paper, we will study the St. Venant's torsional bar with a radial vertical crack and compare the present computed results with the results of reference [2], The present computed results show that, using the method provided in this paper, satisfactory convergent solution can be obtained under lower degree of freedom.
文摘The paper has proved that Hellinger-Reissuer and Hu-Washizu variational principles are but equivalent principles in elasticity by following three ways: 1) Lagrange multiplier method. The paper points out that only a new independent variable can be introduced when one constraint equation has been eliminated by one Lagrange multiplier, which must be expressed as a function of the original variable(s) and/or the new introduced variable after identification. In using Lagrange multiplier method to deduce Hu-Washizu principle from the minimum potential energy principle, which has only one kind of independent variable namely displacement, by eliminating the constraint equations of stress-displacement relations, one can only obtain a principle with two kinds of variables namely displacement and stress; 2) involutory transformation, with such method Hu-Washizu variational principle can be deduce directly from the Hellinger-Reissner variational principle under the same variational constraints of Stress-strain relation, and vice verse; 3)semi-inverse method, by which both of the above variational principles can be deduced from the minimum potential energy principle with tile same variational constraints. So the three kinds of variational functions in Hu-Washizu variational principle are not independent to each other,the stress-strain relationships are still its constraint conditions.
文摘Based on generalized variational principles, an element called MR-12 was constructed for the static and dynamic analysis of thin plates with orthogonal anisotropy. Numerical results showed that this incompatible element converges very rapidly and has good accuracy. It was demonstrated that generalized varialional principles arc useful and effective in founding incompatible clement.Moreover, element MR-12 is easy for implementation since it does not differ very much from the common rectangular element R-12 of thin plate.
文摘This paper discusses the generalized variational principles founded by the technique of Lagrangian multipliers in structural mechanics and analyzes the nonlinear statically indeterminate structures. It is assumed that the stress-strain relationship of the materials of structures has the form of namely, the physical equations of structures have the shape of exponential functions. Several examples are given to illustrate the statically indeterminate structures such as the trusses, beams, frames and torsional bars.
文摘4 semi-analytical approach for the dynamic response of general thin plates whichemployes finite element discretization in space domain and a series of representation intime demain is developed on the basis of Curtin variational principles.The formulationof time series is also investigated so that the dynamic response of plates with arbitraryshape and boundary constraints can be achieved with adequate accuracy.
文摘The generalized variational principles of isothermal quasi-static fluid full-filled elastic solids are established by using Variational Integral Method. Then by introducing constraints, several kinds of variational principles are worked out, including five-field variable, four-field variable, three-field variable and two-field variable formulations. Some new variational principles are presented besides the principles noted in the previous works. Based on variational principles, finite element models can be set up.
基金Project supported by the National Natural Science Foundation of China(No.10272070)and the Development Foun-dation of the Education Commission of Shanghai,China.
文摘Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fuid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small de- formation of the solid phase, small velocity of the fuid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fuid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles, especially Hu-Washizu type variational principles, for the initial boundary value problems of dy- namic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fuid-saturated elastic porous media, which have been discussed previously.
基金supported by the National Natural Science Foundations of China (Grant 11502286)
文摘For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.
基金Project supported by the National Natural Science Foundation of China (No. 10072038)the Special Fund for PhD Program of Education Ministry of China (No. 2000005616)
文摘Based upon the Hellinger-Reissner (H-R) mixed variational principle for three-dimensional elastic bodies, the modified H-R mixed variational theorem for magnetoelectroelastic bodies was established. The state-vector equation of magnetoelectroelastic plates was derived from the proposed theorem by performing the variational operations. To lay a theoretical basis of the semi-analytical solution applied with the magnetoelectroelastic plates, the state-vector equation for the discrete element in plane was proposed through the use of the proposed principle. Finally, it is pointed out that the modified H-R mixed variational principle for pure elastic, single piezoelectric or single piezomagnetic bodies are the special cases of the present variational theorem.