The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology ...The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
The existing analytical solutions are extended to obtain the stress fields and the stress intensity factors(SIFs) of two unequal aligned cracks emanating from an elliptical hole in an infinite isotropic plane. A confo...The existing analytical solutions are extended to obtain the stress fields and the stress intensity factors(SIFs) of two unequal aligned cracks emanating from an elliptical hole in an infinite isotropic plane. A conformal mapping is proposed and combined with the complex variable method. Due to some difficulties in the calculation of the stress function, the mapping function is approximated and simplified via the applications of the series expansion. To validate the obtained solution, several examples are analyzed with the proposed method, the finite element method, etc. In addition, the effects of the lengths of the cracks and the ratio of the semi-axes of the elliptical hole(a/b) on the SIFs are studied. The results show that the present analytical solution is applicable to the SIFs for small cracks.展开更多
The strength of composite plate with different hole-shapes is always one of the most important but complicated issues in the application of the composite material.The holes will lead to mutations and discontinuity to ...The strength of composite plate with different hole-shapes is always one of the most important but complicated issues in the application of the composite material.The holes will lead to mutations and discontinuity to the structure.So the hole-edge stress concentration is always a serious phenomenon.And the phenomenon makes the structure strength decrease very quickly to form dangerous weak points.Most partial damage begins from these weak points.According to the complex variable functions theory,the accurate boundary condition of composite plate with different hole-shapes is founded by conformal mapping method to settle the boundary condition problem of complex hole-shapes.Composite plate with commonly hole-shapes in engineering is studied by several complex variable stress fimction.The boundary integral equations are founded based on exact boundary conditions.Then the exact hole-edge stress analytic solution of composite plate with rectangle holes and wing manholes is resolved.Both of offset axis loadings and its influences on the stress concentration coefficient of the hole-edge are discussed.And comparisons of different loads along various offset axis on the hole-edge stress distribution of orthotropic plate with rectangle hole or wing manhole are made.It can be concluded that hole-edge with continuous variable curvatures might help to decrease the stress concentration coefficient;and smaller angle of outer load and fiber can decrease the stress peak value.展开更多
By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem...By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.展开更多
In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied....In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.展开更多
文摘The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
文摘The existing analytical solutions are extended to obtain the stress fields and the stress intensity factors(SIFs) of two unequal aligned cracks emanating from an elliptical hole in an infinite isotropic plane. A conformal mapping is proposed and combined with the complex variable method. Due to some difficulties in the calculation of the stress function, the mapping function is approximated and simplified via the applications of the series expansion. To validate the obtained solution, several examples are analyzed with the proposed method, the finite element method, etc. In addition, the effects of the lengths of the cracks and the ratio of the semi-axes of the elliptical hole(a/b) on the SIFs are studied. The results show that the present analytical solution is applicable to the SIFs for small cracks.
基金supported by National Natural Science Foundation of China(No.50675209)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education of China(No.200724).
文摘The strength of composite plate with different hole-shapes is always one of the most important but complicated issues in the application of the composite material.The holes will lead to mutations and discontinuity to the structure.So the hole-edge stress concentration is always a serious phenomenon.And the phenomenon makes the structure strength decrease very quickly to form dangerous weak points.Most partial damage begins from these weak points.According to the complex variable functions theory,the accurate boundary condition of composite plate with different hole-shapes is founded by conformal mapping method to settle the boundary condition problem of complex hole-shapes.Composite plate with commonly hole-shapes in engineering is studied by several complex variable stress fimction.The boundary integral equations are founded based on exact boundary conditions.Then the exact hole-edge stress analytic solution of composite plate with rectangle holes and wing manholes is resolved.Both of offset axis loadings and its influences on the stress concentration coefficient of the hole-edge are discussed.And comparisons of different loads along various offset axis on the hole-edge stress distribution of orthotropic plate with rectangle hole or wing manhole are made.It can be concluded that hole-edge with continuous variable curvatures might help to decrease the stress concentration coefficient;and smaller angle of outer load and fiber can decrease the stress peak value.
文摘By using the complex variable method and conformal mapping, the diffraction of flexural waves and dynamic stress concentrations in thick plates with a cavity have been studied. A general solution of the stress problem of the thick plate satisfying the boundary conditions on the contour of an arbitrary cavity is obtained. By employing the orthogonal function expansion technique, the dynamic stress problem can be reduced to the solution of an infinite algebraic equation series. As an example, the numerical results for the dynamic stress concentration factor in thick plates with a circular, elliptic cavity are graphically presented. The numerical results are discussed.
文摘In this paper, based on complex variables and conformal mapping methods, using the refined dynamic equation of plates, elastic wave scattering and dynamic stress concentrations in plates with two cutouts were studied. Applying the orthogonal function expansion method, the problem to be solved can be reduced into the solution of a set of infinite algebraic equations. According to free boundary conditions, numerical results of dynamic moment concentration factors in thick plates with two circular cutouts analyze that: there will be more complex interaction changes between two-cutout situation than single cutout situation. In the case of low frequency or high frequency and thin plate, the hole-spacing in the absence of coupling interactions was larger or smaller. The numerical results and method can be used to analyze the dynamics and strength of plate-like structures.
基金The National Natural Science Foundation of China(11362018,11261045)the Specialized Research Fund for the Doctoral Program of Higher Education of China(20116401110002)the Natural Science Foundation of Ningxia University(ZR15026).