Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

A solving method for a system of partial differential equations with an application to the bending problem of a thick plate

Ｓｏｍｅｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓｃａｎｂｅｅｘｐｒｅｓｓｅｄｂｙｓｕｃｈａｓｙｓｔｅｍｏｆｌｉｎｅａｒｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓａｓＬ１［ｕ１（ｑ）］＝Ｆ１（ｑ）,Ｌ２［ｕ２（ｑ）］＝Ｆ２（ｑ）＋ｆ２［ｕ１（ｑ）］,…Ｌｎ［ｕｎ（ｑ）］＝Ｆｎ（ｑ）＋ｆｎ［ｕ１（ｑ）,ｕ２（ｑ）,…,ｕｎ－１（ｑ）］,（１）ａｎｄｔｈｅｉｒｉｎｉｔｉａｌａｎｄｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎｓ,ｗｈｅｒｅｑｉｓｔｈｅｓｅｔｏｆａｌｌａｒｇｕｍｅｎｔｓｗ…     

Numerical simulation of ground water flow via a new approach to the local radial point interpolation meshless method

A novel approach to local radial point interpolation meshless (LRPIM) method is introduced to investigate the influence of leakage on tidal response in a coastal leaky confined aquifer system, based on a local weighted residual method with the Heaviside step function as the weighting function over a local sub-domain. The present approach is a truly meshless method based only on a number of randomly located nodes. In this approach, neither global background integration mesh nor domain integration is needed. Radial basis functions (RBFs) interpolation is employed in shape function and its derivatives construction for evaluating the local weak form integrals. Due to satisfaction of kronecker delta property in RBF interpolation, no special treatment is needed to impose the essential boundary conditions. In order to obtain the optimum parameters, shape parameters of multiquadrics (MQ)-RBF are tuned and studied. The leakage has a significant impact on the tidal behaviour of the confined aquifer. The numerical results of this research indicate that both tidal amplitude of groundwater head in the aquifer and the distance over which the aquifer can be disturbed by the tide are considerably reduced by leakage. The novelty of the approach is the use of a local Heaviside weight function in the LRPIM which does not need local domain integration and only integrations on the boundary of the local domains are needed. Therefore, in this research a new local Heaviside weight function has been proposed. Numerical results are presented and compared with the results of analytical solution. It is observed that the obtained results agreed very well with the results of analytical solution. The numerical results show that the use of a local Heaviside weight function in the LRPIM is highly accurate, fast and robust. It is also noticed that this novel meshless approach using MQ radial basis is very stable.     

Applications of the dual integral formulation in conjunction with fast multipole method to the oblique incident wave problem

In this paper, the dual integral formulation is derived for the modified Helmholtz equation in the propagation of oblique incident wave passing a thin barrier (zero thickness) by employing the concept of fast multipole method (FMM) to accelerate the construction of an influence matrix. By adopting the addition theorem, the four kernels in the dual formulation are expanded into degenerate kernels that separate the field point and the source point. The source point matrices decomposed in the four influence matrices are similar to each other or only to some combinations. There are many zeros or the same influence coefficients in the field point matrices decomposed in the four influence matrices, which can avoid calculating the same terms repeatedly. The separable technique reduces the number of floating‐point operations from O((N)²) to O(N log^a(N)), where N is the number of elements and a is a small constant independent of N. Finally, the FMM is shown to reduce the CPU time and memory requirement, thus enabling us to apply boundary element method (BEM) to solve water scattering problems efficiently. Two‐moment FMM formulation was found to be sufficient for convergence in the singular equation. The results are compared well with those of conventional BEM and analytical solutions and show the accuracy and efficiency of the FMM. Copyright © 2008 John Wiley & Sons, Ltd.     

A generalized constitutive elasticity law for GLPD micromorphic materials,with application to the problem of a spherical shell subjected to axisymmetric loading conditions

In this work we propose to replace the GLPD hypo-elasticity law by a more rigorous generalized Hooke's law based on classical material symmetry characterization assumptions. This law introduces in addition to the two well-known Lame's moduli, five constitutive constants. An analytical solution is derived for the problem of a spherical shell subjected to axisymmetric loading conditions to illustrate the potential of the proposed generalized Hooke's law.     

Uniaxial extension and shear in combination with cyclic straining as a promising method for studying the viscoelasticity of polymers over a wide range of stresses and strains atT > T g

Simple relations are found between the fracture strains, the long-term durability and the initial values of parameters characterizing polymer rheology at deformation rates and stresses approaching zero. Fracture regimes are determined by two groups of parameters. One includes critical values of stresses. They are invariant with repect to temperature and molecular weight of the polymers; the values of critical stresses for different polymer compositions differ by a factor of 10 to 20. The second group of critical parameters includes the rates of deformation determined by the initial viscosity. The latter may vary by many orders of magnitude. There exists a universal critical value determining polymer fracture independent of linear macromolecule composition, its molecular weight, the temperature and the way of attaining a given state. This value is the recoverable strain and is equal to 0.5 according to Hencky. There exists a relation between the maximum value of recoverable strain in the transition region from the rubbery to the leathery state and the extensibility of macromolecules for polymers with various molecular weights. Quenching of the polymer near the maximum recoverable strain makes it possible to obtain high strength samples. Overspurt regimes for polymer flow have also been studied. It has been shown that this causes polymer static electrification. Simple and unique dependences of the charge density on temperature and polymer molecular weight have been established.Presented at the 28th IUPAC International Symposium on Macromolecules, June 1982, Amherst, Mass. (USA)Dedicated to Prof. Dr. Josef Schurz on the occassion of his 60th birthday     

Perturbation solution to the nonlinear problem of oblique water exit of an axisymmetric body with a large exit-angle

In this paper,a nonlinear,unsteady3-D free surface problem of the oblique water exitof an axisymmetric body with a large water exit-angle was investigated by means of theperturbation method in which the complementary angle a of the water exit angle waschosen as a small parameter.The original3-D problem was solved by expanding it into apower series of a and reduced to a number of2-D problems.The integral expressions forthe first three order solutions were given in terms of the complete elliptic functions of thefirst and second kinds.The zeroth-order solution didn‘t turn out to be a linear problem asusual but a nonlinear one corresponding to the vertical water exit for the same body.Computational results were presented for the free surface shapes and the forces exerted upto the second order during the oblique water exit of a series of ellipsoids with various ratiosof length to diameter at different Froude numbers.