The distribution network exhibits complex structural characteristics,which makes fault localization a challenging task.Especially when a branch of the multi-branch distribution network fails,the traditional multi-bran...The distribution network exhibits complex structural characteristics,which makes fault localization a challenging task.Especially when a branch of the multi-branch distribution network fails,the traditional multi-branch fault location algorithm makes it difficult to meet the demands of high-precision fault localization in the multi-branch distribution network system.In this paper,the multi-branch mainline is decomposed into single branch lines,transforming the complex multi-branch fault location problem into a double-ended fault location problem.Based on the different transmission characteristics of the fault-traveling wave in fault lines and non-fault lines,the endpoint reference time difference matrix S and the fault time difference matrix G were established.The time variation rule of the fault-traveling wave arriving at each endpoint before and after a fault was comprehensively utilized.To realize the fault segment location,the least square method was introduced.It was used to find the first-order fitting relation that satisfies the matching relationship between the corresponding row vector and the first-order function in the two matrices,to realize the fault segment location.Then,the time difference matrix is used to determine the traveling wave velocity,which,combined with the double-ended traveling wave location,enables accurate fault location.展开更多
Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-freque...Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.展开更多
As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initiall...As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initially proposed by Jiang et al.(Computational and Applied Mathematics,2021,40:174),through the utilization of a convex combination technique.And this improvement allows for an adaptive search direction by integrating a newly constructed spectral gradient-type restart strategy.Then,we develop a new spectral CGM by employing an inexact line search to determine the step size.With the application of the weak Wolfe line search,we establish the sufficient descent property of the proposed search direction.Moreover,under general assumptions,including the employment of the strong Wolfe line search for step size calculation,we demonstrate the global convergence of our new algorithm.Finally,the given unconstrained optimization test results show that the new algorithm is effective.展开更多
Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms(PKMs), but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method fo...Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms(PKMs), but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom(DOFs) and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R(T: translational DOF; R: rotational DOF), are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs' type synthesis. Moreover, this method can provide theoretical guidance for other PKMs' type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods' problems in terms of complicated, not intuitive and unsuitable for practical application.展开更多
This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two node...This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.展开更多
In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Comb...In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.展开更多
In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs...In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.展开更多
The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake d...In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.展开更多
In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-differenc...In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
The Finite Element Method of Lines (FEMOL) is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to...The Finite Element Method of Lines (FEMOL) is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor (SIF) and crack opening displacement (COD) are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.展开更多
In this work, we applied two electromagnetic models for the characterization of a planar structure including a flat, thick copper conductor. Indeed the first model is consisted by modeling two metal ribbons without bu...In this work, we applied two electromagnetic models for the characterization of a planar structure including a flat, thick copper conductor. Indeed the first model is consisted by modeling two metal ribbons without bulkiness, placed one above the other at a distance of h<sub>2</sub> equal to the thickness of the thick conductor. This approach has been implemented and tested by the iterative method. The results of simulations have been compared with those calculated by the Ansoft HFSS software, and they are in good concordance, validating the method of analysis used. The second model is based on the calculation of the effective permittivity of the medium containing the thick conductor. This medium consists of a metallic region of complex relative permittivity , the rest of this medium is filled with air e<sub>r</sub><sub>2</sub> = 1. The effective permittivity e<sub>eff</sub> calculated from these two relative permittivity e<sub>r</sub><sub>2</sub> and . Comparing the simulation results of this new formulation of the iterative method with those calculated by the software Ansoft HFSS shows that they are in good matching which validates the second model.展开更多
At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line...At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle, so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory, the theory of slip line field is put forward in this paper, and then the ultimate beating capacity of strip footing is obtained based on the associated flow rule and the non-associated flow nile individually. These two results are identical since the ultimate bearing capacity is independent of flow role. On the contrary, the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associat- ed flow rule is incorrect.展开更多
We present a time domain hybrid method to realize the fast coupling analysis of transmission lines excited by space electromagnetic fields, in which parallel finite-difference time-domain (FDTD) method, interpolation ...We present a time domain hybrid method to realize the fast coupling analysis of transmission lines excited by space electromagnetic fields, in which parallel finite-difference time-domain (FDTD) method, interpolation scheme, and Agrawal model-based transmission line (TL) equations are organically integrated together. Specifically, the Agrawal model is employed to establish the TL equations to describe the coupling effects of space electromagnetic fields on transmission lines. Then, the excitation fields functioning as distribution sources in TL equations are calculated by the parallel FDTD method through using the message passing interface (MPI) library scheme and interpolation scheme. Finally, the TL equations are discretized by the central difference scheme of FDTD and assigned to multiple processors to obtain the transient responses on the terminal loads of these lines. The significant feature of the presented method is embodied in its parallel and synchronous calculations of the space electromagnetic fields and transient responses on the lines. Numerical simulations of ambient wave acting on multi-conductor transmission lines (MTLs), which are located on the PEC ground and in the shielded cavity respectively, are implemented to verify the accuracy and efficiency of the presented method.展开更多
In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wol...In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wolfe line search conditions.Numerical results show that the new method is efficient and stationary by comparing with PRP+ method,so it can be widely used in scientific computation.展开更多
The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the meth...The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the method of lines (MOL) is introduced to solve the temperature field of FGM. The basic idea of the method is to semi-discretize the governing equation into a system of ordinary differential equations (ODEs) defined on discrete lines by means of the finite difference method. The temperature field of FGM can be obtained by solving the ODEs. The functions of thermal properties are directly embodied in these equations and these properties are not discretized in the domain. Thus, difficulty of FEM and BEM is overcome by the method. As a numerical example, the temperature field of a plane problem is analyzed for FGMs through varying thermal conductivity coefficient by the MOL.展开更多
In this paper we consider the global convergence of any conjugate gradient method of the form d1=-g1,dk+1=-gk+1+βkdk(k≥1)with any βk satisfying sume conditions,and with the strong wolfe line search conditions.Under...In this paper we consider the global convergence of any conjugate gradient method of the form d1=-g1,dk+1=-gk+1+βkdk(k≥1)with any βk satisfying sume conditions,and with the strong wolfe line search conditions.Under the convex assumption on the objective function,we preve the descenf property and the global convergence of this method.展开更多
基金This work was funded by the project of State Grid Hunan Electric Power Research Institute(No.SGHNDK00PWJS2210033).
文摘The distribution network exhibits complex structural characteristics,which makes fault localization a challenging task.Especially when a branch of the multi-branch distribution network fails,the traditional multi-branch fault location algorithm makes it difficult to meet the demands of high-precision fault localization in the multi-branch distribution network system.In this paper,the multi-branch mainline is decomposed into single branch lines,transforming the complex multi-branch fault location problem into a double-ended fault location problem.Based on the different transmission characteristics of the fault-traveling wave in fault lines and non-fault lines,the endpoint reference time difference matrix S and the fault time difference matrix G were established.The time variation rule of the fault-traveling wave arriving at each endpoint before and after a fault was comprehensively utilized.To realize the fault segment location,the least square method was introduced.It was used to find the first-order fitting relation that satisfies the matching relationship between the corresponding row vector and the first-order function in the two matrices,to realize the fault segment location.Then,the time difference matrix is used to determine the traveling wave velocity,which,combined with the double-ended traveling wave location,enables accurate fault location.
文摘Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards.
基金supported by the National Natural Science Foundation of China (No.72071202)the Key Laboratory of Mathematics and Engineering ApplicationsMinistry of Education。
文摘As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initially proposed by Jiang et al.(Computational and Applied Mathematics,2021,40:174),through the utilization of a convex combination technique.And this improvement allows for an adaptive search direction by integrating a newly constructed spectral gradient-type restart strategy.Then,we develop a new spectral CGM by employing an inexact line search to determine the step size.With the application of the weak Wolfe line search,we establish the sufficient descent property of the proposed search direction.Moreover,under general assumptions,including the employment of the strong Wolfe line search for step size calculation,we demonstrate the global convergence of our new algorithm.Finally,the given unconstrained optimization test results show that the new algorithm is effective.
基金supported by National Natural Science Foundation of China(Grant No.51135008)National Basic Research Program of China(973 Program,Grant No.2013CB035400)China Postdoctoral Science Foundation(Grant Nos.2012M520256,2013T60107)
文摘Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms(PKMs), but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom(DOFs) and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R(T: translational DOF; R: rotational DOF), are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs' type synthesis. Moreover, this method can provide theoretical guidance for other PKMs' type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods' problems in terms of complicated, not intuitive and unsuitable for practical application.
基金supported by the National Natural Science Foundation of China (Grant No.11072052)the National High Technology Research and Development Program of China (863 Program,Grant No.2006AA09A109-3)
文摘This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.
文摘In this paper, a new class of three term memory gradient method with non-monotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.
文摘In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.
文摘The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
基金Project supported by the National Key Research and Development Program of China(Nos.2019YFE0192600,2017YFE0132000,and 2019YFB1503700)the National Natural Science Foundation of China(Nos.51761135012 and 11872248)。
文摘In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.
文摘In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
文摘The Finite Element Method of Lines (FEMOL) is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor (SIF) and crack opening displacement (COD) are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.
文摘In this work, we applied two electromagnetic models for the characterization of a planar structure including a flat, thick copper conductor. Indeed the first model is consisted by modeling two metal ribbons without bulkiness, placed one above the other at a distance of h<sub>2</sub> equal to the thickness of the thick conductor. This approach has been implemented and tested by the iterative method. The results of simulations have been compared with those calculated by the Ansoft HFSS software, and they are in good concordance, validating the method of analysis used. The second model is based on the calculation of the effective permittivity of the medium containing the thick conductor. This medium consists of a metallic region of complex relative permittivity , the rest of this medium is filled with air e<sub>r</sub><sub>2</sub> = 1. The effective permittivity e<sub>eff</sub> calculated from these two relative permittivity e<sub>r</sub><sub>2</sub> and . Comparing the simulation results of this new formulation of the iterative method with those calculated by the software Ansoft HFSS shows that they are in good matching which validates the second model.
文摘At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle, so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory, the theory of slip line field is put forward in this paper, and then the ultimate beating capacity of strip footing is obtained based on the associated flow rule and the non-associated flow nile individually. These two results are identical since the ultimate bearing capacity is independent of flow role. On the contrary, the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associat- ed flow rule is incorrect.
基金Project supported by the National Natural Science Foundation of China(Grant No.61701057)the Chongqing Research Program of Basic Research and Frontier Technology,China(Grant No.cstc2017jcyjAX0345).
文摘We present a time domain hybrid method to realize the fast coupling analysis of transmission lines excited by space electromagnetic fields, in which parallel finite-difference time-domain (FDTD) method, interpolation scheme, and Agrawal model-based transmission line (TL) equations are organically integrated together. Specifically, the Agrawal model is employed to establish the TL equations to describe the coupling effects of space electromagnetic fields on transmission lines. Then, the excitation fields functioning as distribution sources in TL equations are calculated by the parallel FDTD method through using the message passing interface (MPI) library scheme and interpolation scheme. Finally, the TL equations are discretized by the central difference scheme of FDTD and assigned to multiple processors to obtain the transient responses on the terminal loads of these lines. The significant feature of the presented method is embodied in its parallel and synchronous calculations of the space electromagnetic fields and transient responses on the lines. Numerical simulations of ambient wave acting on multi-conductor transmission lines (MTLs), which are located on the PEC ground and in the shielded cavity respectively, are implemented to verify the accuracy and efficiency of the presented method.
基金Supported by the Fund of Chongqing Education Committee(KJ091104)
文摘In this paper,an efficient conjugate gradient method is given to solve the general unconstrained optimization problems,which can guarantee the sufficient descent property and the global convergence with the strong Wolfe line search conditions.Numerical results show that the new method is efficient and stationary by comparing with PRP+ method,so it can be widely used in scientific computation.
文摘The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the method of lines (MOL) is introduced to solve the temperature field of FGM. The basic idea of the method is to semi-discretize the governing equation into a system of ordinary differential equations (ODEs) defined on discrete lines by means of the finite difference method. The temperature field of FGM can be obtained by solving the ODEs. The functions of thermal properties are directly embodied in these equations and these properties are not discretized in the domain. Thus, difficulty of FEM and BEM is overcome by the method. As a numerical example, the temperature field of a plane problem is analyzed for FGMs through varying thermal conductivity coefficient by the MOL.
基金This work is supported by the National Natural Science Foundation of China
文摘In this paper we consider the global convergence of any conjugate gradient method of the form d1=-g1,dk+1=-gk+1+βkdk(k≥1)with any βk satisfying sume conditions,and with the strong wolfe line search conditions.Under the convex assumption on the objective function,we preve the descenf property and the global convergence of this method.