Based on B-spline wavelet on the interval (BSWI), two classes of truncated conical shell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conical shell element and BSWI moderately t...Based on B-spline wavelet on the interval (BSWI), two classes of truncated conical shell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conical shell element and BSWI moderately thick truncated conical shell element with independent slopedeformation interpolation. In the construction of wavelet-based element, instead of traditional polynomial interpolation, the scaling functions of BSWI were employed to form the shape functions through the constructed elemental transformation matrix, and then construct BSWI element via the variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkin method, the elemental displacement field represented by the coefficients of wavelets expansion was transformed into edges and internal modes via the constructed transformation matrix. BSWI element combines the accuracy of B-spline function approximation and various wavelet-based elements for structural analysis. Some static and dynamic numerical examples of conical shells were studied to demonstrate the present element with higher efficiency and precision than the traditional element.展开更多
An 8-noded locking-free degenerated isoparametric shell element is presented. A revised interpolation for shear strain terms was constructed in natural co-ordinate system such that all necessary modes (translation, ro...An 8-noded locking-free degenerated isoparametric shell element is presented. A revised interpolation for shear strain terms was constructed in natural co-ordinate system such that all necessary modes (translation, rotation and constant curvature) are preserved, which can be used to eliminate shear locking. A revised interpolation for membrane strains was produced in the local Cartesian co-ordinate system to overcome membrane locking behavior. The new 8-noded element has the proper rank, with the requisite number of zero eigenvalues each associated with a rigid mode. The element does not exhibit membrane or shear locking for large span-thickness ratio. The element does not form element mechanisms or extra spurious zero energy modes. Therefore, it can be used for both thin and thick shells.展开更多
The linear buckling problems of plates and shells were analysed using a recently developped quadrilateral,16-degrees of freedom flat shell element (called DKQ16).The geometrical stiffness matrix was established.Compar...The linear buckling problems of plates and shells were analysed using a recently developped quadrilateral,16-degrees of freedom flat shell element (called DKQ16).The geometrical stiffness matrix was established.Comparison of the numerical results for several typical problems shows that the DKQ16 element has a very good precision for the linear buckling problems of plates and shells.展开更多
In this work, a novel numerical method is developed for simulating arbitrary crack growth in pipes with the idea of enriched shape functions which can represent the discontinuity independent of the mesh. The concept o...In this work, a novel numerical method is developed for simulating arbitrary crack growth in pipes with the idea of enriched shape functions which can represent the discontinuity independent of the mesh. The concept of the enriched shape functions is introduced into the continuum-based (CB) shell element. Due to the advantage of CB shell element, the shell thickness varia- tion and surface connection can be concerned during the deformation. The stress intensity factors of the crack in the CB shell element are calculated by using the 'equivalent domain integral' method for 3D arbitrary non-planar crack. The maximum en- ergy release rate is used as a propagation criterion. This method is proved able to capture arbitrary crack growth path in pipes which is independent of the element mesh. Numerical examples of different fracture patterns in pipes are presented here.展开更多
The continuum-based(CB)shell theory is combined with the extended finite element method(X-FEM)in this paper to model crack propagation in shells under static and dynamic situations.Both jump function and asymptotic cr...The continuum-based(CB)shell theory is combined with the extended finite element method(X-FEM)in this paper to model crack propagation in shells under static and dynamic situations.Both jump function and asymptotic crack tip solution are adopted for describing the discontinuity and singularity of the crack in shells.Level set method(LSM)is used to represent the crack surface and define the enriched shape functions.Stress intensity factors(SIFs)are calculated by the displacement interpolation technique to prove the capability of the method and the maximum strain is applied for the fracture criterion.Also,an efficient integration scheme for the CB shell element with cracks is proposed.展开更多
A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis (MRA) framework is formulated out of a mutually nesting displacement s...A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis (MRA) framework is formulated out of a mutually nesting displacement subspace sequence, whose basis functions are constructed of scaling and shifting on the element domain of basic node shape functions. The basic node shape functions are constructed from shifting to other three quadrants around a specific node of a basic element in one quadrant and joining the corresponding node shape functions of four elements at the specific node. The MRA endows the proposed element with the resolution level (RL) to adjust the element node number, thus modulating structural analysis accuracy accordingly. The node shape functions of Kronecker delta property make the treatment of element boundary condition quite convenient and enable the stiffness matrix and the loading column vectors of the proposed element to be automatically acquired through quadraturing around nodes in RL adjusting. As a result, the traditional 4-node rectangular shell element is a mono-resolution one and also a special case of the proposed element. The accuracy of a structural analysis is actually determined by the RL, not by the mesh. The simplicity and clarity of node shape function formulation with the Kronecker delta property, and the rational MRA enable the proposed element method to be implemented more rationally, easily and efficiently than the conventional mono-resolution rectangular shell element method or other corresponding MRA methods.展开更多
A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpola...A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpolation where the bending and membrane stiffness matrices are calculated on the boundaries of the smoothing cells while the shear terms are approximated by independent interpolation functions in natural coordinates. The proposed element is robust, computationally inexpensive and free of locking. Since the integration is done on the element boundaries for the bending and membrane terms, the element is more accurate than the MITC4 element for distorted meshes. This will be demonstrated for several numerical examples.展开更多
Trace elemental associations and Sr-Nd isotopic compositions are of important to recognition of biogenic material from mixed marine sediments. The foraminifera shell from the Okinawa Trough strongly enriches Sr,P,Mn a...Trace elemental associations and Sr-Nd isotopic compositions are of important to recognition of biogenic material from mixed marine sediments. The foraminifera shell from the Okinawa Trough strongly enriches Sr,P,Mn and Ba, enriches Li,U,Th,Sc,Co,Cu,Pb,Zn,Cr,Rb,Y,Sb and light rare earth elements,slightly enriches V,Ga,Zr,Nb,Cd and middle rare earth elements,is short of Mo,In,Sn,Cs,Hf,Ta,W,Ti,Bi and heavy rare earth elements. The mechanism of elemental enrichment in forminifera is the concentrations of trace elements in sea water and selective absorption of trace elements during foraminifera living, as well as the geochemical affinity between major elements and trace elements.The REE (rare earth elements)partition pattern of foraminifera shell of the Okinawa Trough shows enrichment of middle rare earth elements with slightly negative Ce anomaly,which are different from those of foraminifera of the Pacific Ocean.The Sr,Nd isotopic ratios of the Okinawa Trough foraminifera are 0 709 769 and 0 512 162,respectively, which are different not only from those of oceanic water, but also from those of river water of China's Mainland, the former is slightly higher than those of oceanic water,but much lower than those of river water;the latter is slightly lower than those of oceanic water,but higher than those of river water,demonstrating that the Okinawa Trough sea water has been influenced by river water of China's Mainland.展开更多
This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc mod...This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger-Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid ele- ment, a hybrid-stress solid-shell element is formulated. The presented finite shell element is able to model arbitrary curved shell structures. Non-linear numerical examples demonstrate the ability of the proposed model to analyze nonlinear piezoelectric devices.展开更多
Based on the generalized vaxiational principle of magneto-thermo-elasticity of a ferromagnetic thin shell established (see, Analyses on nonlinear coupling of magneto-thermo- elasticity of ferromagnetic thin shell--I...Based on the generalized vaxiational principle of magneto-thermo-elasticity of a ferromagnetic thin shell established (see, Analyses on nonlinear coupling of magneto-thermo- elasticity of ferromagnetic thin shell--I), the present paper developed a finite element modeling for the mechanical-magneto-thermal multi-field coupling of a ferromagnetic thin shell. The numerical modeling composes of finite element equations for three sub-systems of magnetic, thermal and deformation fields, as well as iterative methods for nonlinearities of the geometrical large-deflection and the multi-field coupling of the ferromagnetic shell. As examples, the numerical simulations on magneto-elastic behaviors of a ferromagnetic cylindrical shell in an applied magnetic field, and magneto-thermo-elastic behaviors of the shell in applied magnetic and thermal fields are carried out. The results are in good agreement with the experimental ones.展开更多
In this paper, a kind of rationalism theory of shell is established which is of different mechanic characters in tension and in compression, and the finite element numerical analysis method is also described.
An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error b...An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error bounds in three-dimensional fracture mechanics simulation which involves enrichments to solve the singularity in crack tip. This enriched goal-oriented error estimation gives a chance to evaluate this continuum- based shell extended finite element method simulation. With comparisons of reliability to the stress intensity factor calculation in stretching and bending, the accuracy of the continuum-based shell extended finite element method simulation is evaluated, and the reason of error is discussed.展开更多
Projective-iterative version of finite element method has developed for numerical simulation of the stress-strain state of nonhomogeneous shell-type structures (shells with openings). Plastic deformation of the materi...Projective-iterative version of finite element method has developed for numerical simulation of the stress-strain state of nonhomogeneous shell-type structures (shells with openings). Plastic deformation of the material is taken into account when using the method of elastic solutions that reduce the solution of elastoplastic problems to solution of elastic problems. Developed PIV’s significant savings of computer calculation has been compared with the calculation on a fine mesh of traditional FEM. Designed scheme allows analysis of the mutual influence of openings. Analysis of the transformation zone of plastic deformation is developed. For definiteness, the cylindrical shell structures with several rectangular openings are considered.展开更多
The recently introduced even-odd rule has been shown to successfully represent chemical structures of ions and molecules. While comparing available drawings in the scientific literature with the list of compounds pred...The recently introduced even-odd rule has been shown to successfully represent chemical structures of ions and molecules. While comparing available drawings in the scientific literature with the list of compounds predicted by the even-odd rule, it became however obvious that existing compounds are fewer than expected. Several predicted compounds involving many covalent bonds have apparently never been experimentally observed. Neutral oxygen for instance is expected to have 6 valence electrons, whereas oxygen can only build a maximum of two bonds, as in water. This specificity is observed for elements in the top-right corner of the periodic table. For compounds to contain only single covalent bonds, and thus follow the even-odd rule, further explanations are necessary. The present paper proposes that those specific elements experience a transfer of electrons from the valence shell into the inner shell, making them unavailable for further bonding. These elements will be described as organic, hereby providing a clear and hopefully unifying definition of the term. In opposition, inorganic elements have a constant inner shell no matter their electrical state or the number of bonds they maintain. More than 70 compounds involving 11 elements of the main group are studied, revealing a progression from fully inorganic elements at the left of the periodic table to fully organic elements. The transition between inorganic or organic elements is made of few elements that take an organic form when negatively charged;they are labelled semi-organic. The article concludes that the fully organic elements of the main group are Oxygen and Fluorine, whereas semi-organic elements are more numerous: C, N, S, Cl, Se, Br and I. Thus, the even-odd rule becomes fully compatible with scientific knowledge of compounds in liquid or gaseous phase.展开更多
The response of random plate and shell construction is analyzed with the stochastic finite element method (SFEM). Random material properties and geometric dimensions of construction are involved in this paper. A simpl...The response of random plate and shell construction is analyzed with the stochastic finite element method (SFEM). Random material properties and geometric dimensions of construction are involved in this paper. A simplified isoparametric local average model is used to describe the random field. Numerical results of the examples indicate that the approach presented herein is an economical and efficient solution for such an analysis compared with Monte Carlo simulation (MCS).展开更多
In this paper,a boundary element scheme for arbitrary elastic thin shells is elaborated,Based on BEM of 3D linear elasticity and Kirchhoff's hypothesis,boundary integral equations for shells are deduced. As a resu...In this paper,a boundary element scheme for arbitrary elastic thin shells is elaborated,Based on BEM of 3D linear elasticity and Kirchhoff's hypothesis,boundary integral equations for shells are deduced. As a result,only Kelvin's solution is used,the difficulty,in finding a fundamental solution of arbitrary shells is successfully avoided.展开更多
基金Project supported by the National Natural Science Foundation of China (Nos. 50335030, 50505033 and 50575171)National Basic Research Program of China (No. 2005CB724106)Doctoral Program Foundation of University of China(No. 20040698026)
文摘Based on B-spline wavelet on the interval (BSWI), two classes of truncated conical shell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conical shell element and BSWI moderately thick truncated conical shell element with independent slopedeformation interpolation. In the construction of wavelet-based element, instead of traditional polynomial interpolation, the scaling functions of BSWI were employed to form the shape functions through the constructed elemental transformation matrix, and then construct BSWI element via the variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkin method, the elemental displacement field represented by the coefficients of wavelets expansion was transformed into edges and internal modes via the constructed transformation matrix. BSWI element combines the accuracy of B-spline function approximation and various wavelet-based elements for structural analysis. Some static and dynamic numerical examples of conical shells were studied to demonstrate the present element with higher efficiency and precision than the traditional element.
文摘An 8-noded locking-free degenerated isoparametric shell element is presented. A revised interpolation for shear strain terms was constructed in natural co-ordinate system such that all necessary modes (translation, rotation and constant curvature) are preserved, which can be used to eliminate shear locking. A revised interpolation for membrane strains was produced in the local Cartesian co-ordinate system to overcome membrane locking behavior. The new 8-noded element has the proper rank, with the requisite number of zero eigenvalues each associated with a rigid mode. The element does not exhibit membrane or shear locking for large span-thickness ratio. The element does not form element mechanisms or extra spurious zero energy modes. Therefore, it can be used for both thin and thick shells.
文摘The linear buckling problems of plates and shells were analysed using a recently developped quadrilateral,16-degrees of freedom flat shell element (called DKQ16).The geometrical stiffness matrix was established.Comparison of the numerical results for several typical problems shows that the DKQ16 element has a very good precision for the linear buckling problems of plates and shells.
基金supported by the National Natural Science Foundation of China (Grant No. 11011140335)
文摘In this work, a novel numerical method is developed for simulating arbitrary crack growth in pipes with the idea of enriched shape functions which can represent the discontinuity independent of the mesh. The concept of the enriched shape functions is introduced into the continuum-based (CB) shell element. Due to the advantage of CB shell element, the shell thickness varia- tion and surface connection can be concerned during the deformation. The stress intensity factors of the crack in the CB shell element are calculated by using the 'equivalent domain integral' method for 3D arbitrary non-planar crack. The maximum en- ergy release rate is used as a propagation criterion. This method is proved able to capture arbitrary crack growth path in pipes which is independent of the element mesh. Numerical examples of different fracture patterns in pipes are presented here.
基金supported by the National Natural Science Foundation of China(Grant No.11372157)
文摘The continuum-based(CB)shell theory is combined with the extended finite element method(X-FEM)in this paper to model crack propagation in shells under static and dynamic situations.Both jump function and asymptotic crack tip solution are adopted for describing the discontinuity and singularity of the crack in shells.Level set method(LSM)is used to represent the crack surface and define the enriched shape functions.Stress intensity factors(SIFs)are calculated by the displacement interpolation technique to prove the capability of the method and the maximum strain is applied for the fracture criterion.Also,an efficient integration scheme for the CB shell element with cracks is proposed.
基金financial support by the Open Foundation of Chongqing Key Laboratory of Geomechanics and Geoenvironment Protection(Logistical Engineering University)(No.GKLGGP 2013-02)
文摘A multi-resolution rectangular shell element with membrane-bending based on the Kirchhoff-Love theory is proposed. The multi-resolution analysis (MRA) framework is formulated out of a mutually nesting displacement subspace sequence, whose basis functions are constructed of scaling and shifting on the element domain of basic node shape functions. The basic node shape functions are constructed from shifting to other three quadrants around a specific node of a basic element in one quadrant and joining the corresponding node shape functions of four elements at the specific node. The MRA endows the proposed element with the resolution level (RL) to adjust the element node number, thus modulating structural analysis accuracy accordingly. The node shape functions of Kronecker delta property make the treatment of element boundary condition quite convenient and enable the stiffness matrix and the loading column vectors of the proposed element to be automatically acquired through quadraturing around nodes in RL adjusting. As a result, the traditional 4-node rectangular shell element is a mono-resolution one and also a special case of the proposed element. The accuracy of a structural analysis is actually determined by the RL, not by the mesh. The simplicity and clarity of node shape function formulation with the Kronecker delta property, and the rational MRA enable the proposed element method to be implemented more rationally, easily and efficiently than the conventional mono-resolution rectangular shell element method or other corresponding MRA methods.
文摘A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpolation where the bending and membrane stiffness matrices are calculated on the boundaries of the smoothing cells while the shear terms are approximated by independent interpolation functions in natural coordinates. The proposed element is robust, computationally inexpensive and free of locking. Since the integration is done on the element boundaries for the bending and membrane terms, the element is more accurate than the MITC4 element for distorted meshes. This will be demonstrated for several numerical examples.
文摘Trace elemental associations and Sr-Nd isotopic compositions are of important to recognition of biogenic material from mixed marine sediments. The foraminifera shell from the Okinawa Trough strongly enriches Sr,P,Mn and Ba, enriches Li,U,Th,Sc,Co,Cu,Pb,Zn,Cr,Rb,Y,Sb and light rare earth elements,slightly enriches V,Ga,Zr,Nb,Cd and middle rare earth elements,is short of Mo,In,Sn,Cs,Hf,Ta,W,Ti,Bi and heavy rare earth elements. The mechanism of elemental enrichment in forminifera is the concentrations of trace elements in sea water and selective absorption of trace elements during foraminifera living, as well as the geochemical affinity between major elements and trace elements.The REE (rare earth elements)partition pattern of foraminifera shell of the Okinawa Trough shows enrichment of middle rare earth elements with slightly negative Ce anomaly,which are different from those of foraminifera of the Pacific Ocean.The Sr,Nd isotopic ratios of the Okinawa Trough foraminifera are 0 709 769 and 0 512 162,respectively, which are different not only from those of oceanic water, but also from those of river water of China's Mainland, the former is slightly higher than those of oceanic water,but much lower than those of river water;the latter is slightly lower than those of oceanic water,but higher than those of river water,demonstrating that the Okinawa Trough sea water has been influenced by river water of China's Mainland.
基金Supported by the National Natural Science Foundation of China (Grant No. 10672111)the Major Project of the Natural Science Foundation of Jiangsu Province (Grant No. BK2006725)
文摘This paper presents eight-node solid-shell elements for geometric non-linear analyze of piezoelectric structures. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger-Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid ele- ment, a hybrid-stress solid-shell element is formulated. The presented finite shell element is able to model arbitrary curved shell structures. Non-linear numerical examples demonstrate the ability of the proposed model to analyze nonlinear piezoelectric devices.
基金supported by he National Natural Science Foundation of China (No.10872081)the Fok Ying-Tong Education Foundation for Young Teachers in the Higher Education Institutions of China (No.111005)
文摘Based on the generalized vaxiational principle of magneto-thermo-elasticity of a ferromagnetic thin shell established (see, Analyses on nonlinear coupling of magneto-thermo- elasticity of ferromagnetic thin shell--I), the present paper developed a finite element modeling for the mechanical-magneto-thermal multi-field coupling of a ferromagnetic thin shell. The numerical modeling composes of finite element equations for three sub-systems of magnetic, thermal and deformation fields, as well as iterative methods for nonlinearities of the geometrical large-deflection and the multi-field coupling of the ferromagnetic shell. As examples, the numerical simulations on magneto-elastic behaviors of a ferromagnetic cylindrical shell in an applied magnetic field, and magneto-thermo-elastic behaviors of the shell in applied magnetic and thermal fields are carried out. The results are in good agreement with the experimental ones.
文摘In this paper, a kind of rationalism theory of shell is established which is of different mechanic characters in tension and in compression, and the finite element numerical analysis method is also described.
基金Project supported by the National Natural Science Foundation of China(No.10876100)
文摘An enriched goal-oriented error estimation method with extended degrees of freedom is developed to estimate the error in the continuum-based shell extended finite element method. It leads to high quality local error bounds in three-dimensional fracture mechanics simulation which involves enrichments to solve the singularity in crack tip. This enriched goal-oriented error estimation gives a chance to evaluate this continuum- based shell extended finite element method simulation. With comparisons of reliability to the stress intensity factor calculation in stretching and bending, the accuracy of the continuum-based shell extended finite element method simulation is evaluated, and the reason of error is discussed.
文摘Projective-iterative version of finite element method has developed for numerical simulation of the stress-strain state of nonhomogeneous shell-type structures (shells with openings). Plastic deformation of the material is taken into account when using the method of elastic solutions that reduce the solution of elastoplastic problems to solution of elastic problems. Developed PIV’s significant savings of computer calculation has been compared with the calculation on a fine mesh of traditional FEM. Designed scheme allows analysis of the mutual influence of openings. Analysis of the transformation zone of plastic deformation is developed. For definiteness, the cylindrical shell structures with several rectangular openings are considered.
文摘The recently introduced even-odd rule has been shown to successfully represent chemical structures of ions and molecules. While comparing available drawings in the scientific literature with the list of compounds predicted by the even-odd rule, it became however obvious that existing compounds are fewer than expected. Several predicted compounds involving many covalent bonds have apparently never been experimentally observed. Neutral oxygen for instance is expected to have 6 valence electrons, whereas oxygen can only build a maximum of two bonds, as in water. This specificity is observed for elements in the top-right corner of the periodic table. For compounds to contain only single covalent bonds, and thus follow the even-odd rule, further explanations are necessary. The present paper proposes that those specific elements experience a transfer of electrons from the valence shell into the inner shell, making them unavailable for further bonding. These elements will be described as organic, hereby providing a clear and hopefully unifying definition of the term. In opposition, inorganic elements have a constant inner shell no matter their electrical state or the number of bonds they maintain. More than 70 compounds involving 11 elements of the main group are studied, revealing a progression from fully inorganic elements at the left of the periodic table to fully organic elements. The transition between inorganic or organic elements is made of few elements that take an organic form when negatively charged;they are labelled semi-organic. The article concludes that the fully organic elements of the main group are Oxygen and Fluorine, whereas semi-organic elements are more numerous: C, N, S, Cl, Se, Br and I. Thus, the even-odd rule becomes fully compatible with scientific knowledge of compounds in liquid or gaseous phase.
基金SupportedbytheNationalNaturalScienceFoundationofChina (No .5 96 780 39) .
文摘The response of random plate and shell construction is analyzed with the stochastic finite element method (SFEM). Random material properties and geometric dimensions of construction are involved in this paper. A simplified isoparametric local average model is used to describe the random field. Numerical results of the examples indicate that the approach presented herein is an economical and efficient solution for such an analysis compared with Monte Carlo simulation (MCS).
基金The project supported by National Natural Science Foundation of China
文摘In this paper,a boundary element scheme for arbitrary elastic thin shells is elaborated,Based on BEM of 3D linear elasticity and Kirchhoff's hypothesis,boundary integral equations for shells are deduced. As a result,only Kelvin's solution is used,the difficulty,in finding a fundamental solution of arbitrary shells is successfully avoided.