High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit (IC) packaging equipment. In this paper, we derive the transient nonlinear dynamicres...High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit (IC) packaging equipment. In this paper, we derive the transient nonlinear dynamicresponse equations of high-acceleration mechanisms, which reveal that stiffness, frequency, damping, and driving frequency are the primary factors. Therefore, we propose a new structural optimization and velocity-planning method for the precision positioning of a high-acceleration mechanism based on optimal spatial and temporal distribution of inertial energy. For structural optimization, we first reviewed the commonly flexible multibody dynamic optimization using equivalent static loads method (ESLM), and then we selected the modified ESLM for optimal spatial distribution of inertial energy; hence, not only the stiffness but also the inertia and frequency of the real modal shapes are considered. For velocity planning, we developed a new velocity-planning method based on nonlinear dynamic-response optimization with varying motion conditions. Our method was verified on a high-acceleration die bonder. The amplitude of residual vibration could be decreased by more than 20% via structural optimization and the positioning time could be reduced by more than 40% via asymmetric variable velocity planning. This method provides an effective theoretical support for the precision positioning of high-acceleration low-load mechanisms.展开更多
Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control.The fragment spatial distribution of prismatic casings is more complex than that of traditional c...Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control.The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings.In this study,numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted.A new numerical method,which adds the Lagrangian marker points to the Eulerian grid,was proposed to track the multi-material interfaces and material dynamic fractures.Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship.Thereafter,the fragment spatial distributions of the prismatic casing with different fragment sizes,fragment shapes,and casing geometries were obtained using the numerical method.Moreover,fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes,and the experimental data were compared with the numerical results.The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data.Finally,a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.展开更多
Analytical method for the distributions of axial-load and stress is based on elastic assumption, but the threaded connections are often in plastic deformation stage in practice. Meanwhile the strain in the threaded co...Analytical method for the distributions of axial-load and stress is based on elastic assumption, but the threaded connections are often in plastic deformation stage in practice. Meanwhile the strain in the threaded connection is difficult to measure. So it is necessary to study the reliable numerical method. At present neither the convergence analysis of the computational results nor the elastic-plastic analysis in the loading-unloading process are studied. In this paper, von Mises plasticity and kinematic hardening model is used to describe the material response. A new convergence criterion for nonlinear finite element analysis of the loading-unloading process is proposed. An axisymmetric finite element model according to the proposed convergence criterion is developed and used to analyze the distributions of axial-load and stress. It can be conclude that the stress distribution analysis is more dependent on the mesh density than the axial-load distribution analysis. The stress distribution result indicates that with increasing of applied load, the engaged threads close to the nut-bearing surface become plastic firstly. The axial-load distribution result reveals that the load percentage carried by single thread depends on the position of thread and load intensity. When the load is relatively small, the applied load is mainly carried by the engaged threads near the nut-bearing surface, when the load is larger, the differences of percentages for all threads become small. The proposed convergence analyzing procedure is applicable for other nonlinear analyses. The obtained distributions of axial-load and stress can be a reference of engineering application.展开更多
In this paper, primary manufacturing and assembling errors of three-ring gear reducer (TRGR) are analyzed. TRGR is a new transmission type whose eccentric phase difference between middle ring plate and side ring pla...In this paper, primary manufacturing and assembling errors of three-ring gear reducer (TRGR) are analyzed. TRGR is a new transmission type whose eccentric phase difference between middle ring plate and side ring plates is 120°, Its mass of middle ring plate is equal to that of side ring plate or 180°, and its inass of middle ring plate is twice of that of side ring plate, which affects load distribution between ring plates. The primary manufacturing and assembling errors include eccentric error of eccentric sheath E111, internal gear plate E1 and output external gear E11. A new theoretical method is presented in this paper, which converts load on ring plates into the dedendum bending stress of ring plate to calculate load distribution coefficient ( LDC ), by means of gap element method (GEM), one of finite element method (FEM). The theoretical calculation and experimental study, which measures ring plate dedendum bending stress by means of sticking strain gauges on the dedendum of middle ring plate internal gear and side ring plate internal gears, are presented. The theoretical calculation and comparison with experiment result of LDC are implemented an two kinds of three-ring gear reducers whose eccentric phase difference between eccentric sheaths is 120° and 180°respectively. The research indicates that the result of theoretical calculation is consistent with that of experimental study. That is to say, the theoretical calculation method is feasible.展开更多
The mathematical model of current methods for determining load distribution along contact lines of the instantaneously engaged teeth of gears has a defect that its numerical solution is not unique. In order to overcom...The mathematical model of current methods for determining load distribution along contact lines of the instantaneously engaged teeth of gears has a defect that its numerical solution is not unique. In order to overcome this defect, a new method based on three-dimensional finite element compliance matrix method improved by smoothing method is presented. Calculated examples of helical gears show that the new method can calculate tooth load distribution accurately and agree well with the results of experiments.展开更多
The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition exp...The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing anal an integral inequality. Examples of part of these singular sources such as Kelvin's point force, Point-Ring-Couple (PRC) etc. are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.展开更多
In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.Thi...In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.展开更多
The analytical method based on "Hertz theory on normal contact of elastic solids" and the numerical method based on finite element method (FEM) calculating the contact stress of face-gear drive with spur inv...The analytical method based on "Hertz theory on normal contact of elastic solids" and the numerical method based on finite element method (FEM) calculating the contact stress of face-gear drive with spur involute pinion were introduced, and their relative errors are below 10%, except edge contact, which turns out that these two methods can compute contact stress of face-gear drive correctly and effectively. An agreement of the localized bearing contact stress is gotten for these two methods, making sure that the calculation results of FEM are reliable. The loaded meshing simulations of multi-tooth FEM model were developed, and the determination of the transmission error and the maximal load distribution factor of face-gear drive under torques were given. A formula for the maximal load distribution factor was proposed. By introducing the maximal load distribution factor in multi-tooth contact zone, a method for calculating the maximal contact stress in multi-tooth contact can be given. Compared to FEM, the results of these formulae are proved to be reliable, and the relative errors are below 10%.展开更多
Jacket-type offshore platforms are widely used for oil, gas field, and energy development in shallow water. The design of a jacket structure is highly dependent on target environmental variables. This study focuses on...Jacket-type offshore platforms are widely used for oil, gas field, and energy development in shallow water. The design of a jacket structure is highly dependent on target environmental variables. This study focuses on a strategy to estimate design loads for offshore jacket structures based on an environmental contour approach. In addition to the popular conditional distribution model, various classes of bivariate copulas are adopted to construct joint distributions of environmental variables. Analytical formulations of environmental contours based on various models are presented and discussed in this study. The design loads are examined by dynamic response analysis of jacket platform. Results suggest that the conditional model is not recommended for use in estimating design loads in sampling locations due to poor fitting results. Independent copula produces conservative design loads and the extreme response obtained using the conditional model are smaller than those determined by copulas. The suitability of a model for contour construction varies with the origin of wave data. This study provides a reference for the design load estimation of jacket structures and offers an alternative procedure to determine the design criteria for offshore structures.展开更多
Ring footings are suitable for the structures like tall transmission towers, chimneys, silos and oil storages.These types of structures are susceptible to horizontal loads(wind load) in addition to their dead weight.I...Ring footings are suitable for the structures like tall transmission towers, chimneys, silos and oil storages.These types of structures are susceptible to horizontal loads(wind load) in addition to their dead weight.In the literature, very little or no effort has been made to study the effect of ring footing resting on reinforced sand when subjected to eccentric, inclined and/or eccentric-inclined loadings. This paper aims to study the behavior of ring footing resting on loose sand and/or compacted randomly distributed fiberreinforced sand(RDFS) when subjected to eccentric(0 B, 0.05 B and 0.1 B, where B is the outer diameter of ring footing), inclined(0°,5°,10°, 15°,-5°,-10° and-15°)and eccentric-inclined loadings by using a finite element(FE) software PLAXIS 3 D. The behavior of ring footing is studied by using a dimensionless factor called reduction factor(RF). The numerical model used in the PLAXIS 3 D has been validated by conducting model plate load tests. Moreover, an empirical expression using regression analysis has been presented which will be helpful in plotting a load-settlement curve for the ring footing.展开更多
Load and resistance factors are generally obtained using the first order reliability method(FORM)in which the design point should be determined and derivative-based iterations used.In this article,the thirdmoment reli...Load and resistance factors are generally obtained using the first order reliability method(FORM)in which the design point should be determined and derivative-based iterations used.In this article,the thirdmoment reliability index,based on the three-parameter lognormal(3P-lognormal)distribution,is investigated.A simple method based on the third-moment method for estimating load and resistance factors is then proposed,and a simple formula for the target mean resistance is also presented to avoid iterative computations.Unlike the currently used method,the proposed method can be used to determine load and resistance factors,even when the probability density functions(PDFs)of the basic random variables are not available.Moreover,the proposed method does not require the iterative computation of derivatives or any design points.Thus,the method provides a more convenient and effective way to estimate load and resistance factors in practical engineering applications.Numerical examples are presented to demonstrate the advantages of the proposed third moment method for determining load and resistance factors.展开更多
Clustering approaches are one of the probabilistic load flow(PLF)methods in distribution networks that can be used to obtain output random variables,with much less computation burden and time than the Monte Carlo simu...Clustering approaches are one of the probabilistic load flow(PLF)methods in distribution networks that can be used to obtain output random variables,with much less computation burden and time than the Monte Carlo simulation(MCS)method.However,a challenge of the clustering methods is that the statistical characteristics of the output random variables are obtained with low accuracy.This paper presents a hybrid approach based on clustering and Point estimate methods.In the proposed approach,first,the sample points are clustered based on the𝑙-means method and the optimal agent of each cluster is determined.Then,for each member of the population of agents,the deterministic load flow calculations are performed,and the output variables are calculated.Afterward,a Point estimate-based PLF is performed and the mean and the standard deviation of the output variables are obtained.Finally,the statistical data of each output random variable are modified using the Point estimate method.The use of the proposed method makes it possible to obtain the statistical properties of output random variables such as mean,standard deviation and probabilistic functions,with high accuracy and without significantly increasing the burden of calculations.In order to confirm the consistency and efficiency of the proposed method,the 10-,33-,69-,85-,and 118-bus standard distribution networks have been simulated using coding in Python®programming language.In simulation studies,the results of the proposed method have been compared with the results obtained from the clustering method as well as the MCS method,as a criterion.展开更多
基金supported by the National Key Basic Research Program of China (2011CB013104)National Natural Science Foundation of China (U1134004)+2 种基金Guangdong Provincial Natural Science Foundation (2015A030312008)Science and Technology Program of Guangzhou (201510010281)Guangdong Provincial Science and Technology Plan (2013B010402014)
文摘High-speed and precision positioning are fundamental requirements for high-acceleration low-load mechanisms in integrated circuit (IC) packaging equipment. In this paper, we derive the transient nonlinear dynamicresponse equations of high-acceleration mechanisms, which reveal that stiffness, frequency, damping, and driving frequency are the primary factors. Therefore, we propose a new structural optimization and velocity-planning method for the precision positioning of a high-acceleration mechanism based on optimal spatial and temporal distribution of inertial energy. For structural optimization, we first reviewed the commonly flexible multibody dynamic optimization using equivalent static loads method (ESLM), and then we selected the modified ESLM for optimal spatial distribution of inertial energy; hence, not only the stiffness but also the inertia and frequency of the real modal shapes are considered. For velocity planning, we developed a new velocity-planning method based on nonlinear dynamic-response optimization with varying motion conditions. Our method was verified on a high-acceleration die bonder. The amplitude of residual vibration could be decreased by more than 20% via structural optimization and the positioning time could be reduced by more than 40% via asymmetric variable velocity planning. This method provides an effective theoretical support for the precision positioning of high-acceleration low-load mechanisms.
基金supported by the National Natural Science Foundation of China(Grant No.11822203and 11702026)。
文摘Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control.The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings.In this study,numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted.A new numerical method,which adds the Lagrangian marker points to the Eulerian grid,was proposed to track the multi-material interfaces and material dynamic fractures.Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship.Thereafter,the fragment spatial distributions of the prismatic casing with different fragment sizes,fragment shapes,and casing geometries were obtained using the numerical method.Moreover,fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes,and the experimental data were compared with the numerical results.The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data.Finally,a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.
基金supported by Vehicular Diesel Engine Development Program of China (Grant No. DEDP0202)
文摘Analytical method for the distributions of axial-load and stress is based on elastic assumption, but the threaded connections are often in plastic deformation stage in practice. Meanwhile the strain in the threaded connection is difficult to measure. So it is necessary to study the reliable numerical method. At present neither the convergence analysis of the computational results nor the elastic-plastic analysis in the loading-unloading process are studied. In this paper, von Mises plasticity and kinematic hardening model is used to describe the material response. A new convergence criterion for nonlinear finite element analysis of the loading-unloading process is proposed. An axisymmetric finite element model according to the proposed convergence criterion is developed and used to analyze the distributions of axial-load and stress. It can be conclude that the stress distribution analysis is more dependent on the mesh density than the axial-load distribution analysis. The stress distribution result indicates that with increasing of applied load, the engaged threads close to the nut-bearing surface become plastic firstly. The axial-load distribution result reveals that the load percentage carried by single thread depends on the position of thread and load intensity. When the load is relatively small, the applied load is mainly carried by the engaged threads near the nut-bearing surface, when the load is larger, the differences of percentages for all threads become small. The proposed convergence analyzing procedure is applicable for other nonlinear analyses. The obtained distributions of axial-load and stress can be a reference of engineering application.
基金Sponsored by the National Natural Science Foundation of China(Grant No.59575007).
文摘In this paper, primary manufacturing and assembling errors of three-ring gear reducer (TRGR) are analyzed. TRGR is a new transmission type whose eccentric phase difference between middle ring plate and side ring plates is 120°, Its mass of middle ring plate is equal to that of side ring plate or 180°, and its inass of middle ring plate is twice of that of side ring plate, which affects load distribution between ring plates. The primary manufacturing and assembling errors include eccentric error of eccentric sheath E111, internal gear plate E1 and output external gear E11. A new theoretical method is presented in this paper, which converts load on ring plates into the dedendum bending stress of ring plate to calculate load distribution coefficient ( LDC ), by means of gap element method (GEM), one of finite element method (FEM). The theoretical calculation and experimental study, which measures ring plate dedendum bending stress by means of sticking strain gauges on the dedendum of middle ring plate internal gear and side ring plate internal gears, are presented. The theoretical calculation and comparison with experiment result of LDC are implemented an two kinds of three-ring gear reducers whose eccentric phase difference between eccentric sheaths is 120° and 180°respectively. The research indicates that the result of theoretical calculation is consistent with that of experimental study. That is to say, the theoretical calculation method is feasible.
文摘The mathematical model of current methods for determining load distribution along contact lines of the instantaneously engaged teeth of gears has a defect that its numerical solution is not unique. In order to overcome this defect, a new method based on three-dimensional finite element compliance matrix method improved by smoothing method is presented. Calculated examples of helical gears show that the new method can calculate tooth load distribution accurately and agree well with the results of experiments.
文摘The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing anal an integral inequality. Examples of part of these singular sources such as Kelvin's point force, Point-Ring-Couple (PRC) etc. are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.
基金funded by the Funds from Joint National-Local Engineering Research Center for Safe and Precise Coal Mining(Grant No.EC2021004).
文摘In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.
基金Project(50875263) supported by the National Natural Science Foundation of ChinaProject(2011CB706800) supported by the National Basic Research Program of ChinaProject(2010ssxt172) supported by the Natural Science Foundation of Hunan Province,China
文摘The analytical method based on "Hertz theory on normal contact of elastic solids" and the numerical method based on finite element method (FEM) calculating the contact stress of face-gear drive with spur involute pinion were introduced, and their relative errors are below 10%, except edge contact, which turns out that these two methods can compute contact stress of face-gear drive correctly and effectively. An agreement of the localized bearing contact stress is gotten for these two methods, making sure that the calculation results of FEM are reliable. The loaded meshing simulations of multi-tooth FEM model were developed, and the determination of the transmission error and the maximal load distribution factor of face-gear drive under torques were given. A formula for the maximal load distribution factor was proposed. By introducing the maximal load distribution factor in multi-tooth contact zone, a method for calculating the maximal contact stress in multi-tooth contact can be given. Compared to FEM, the results of these formulae are proved to be reliable, and the relative errors are below 10%.
基金supported by the National Key Research and Development Program (No. 2016YFC0303401)the National Natural Science Foundation of China (No. 51779236)the National Natural Science Foundation of China–Shandong Joint Fund Project (No. U1706226)。
文摘Jacket-type offshore platforms are widely used for oil, gas field, and energy development in shallow water. The design of a jacket structure is highly dependent on target environmental variables. This study focuses on a strategy to estimate design loads for offshore jacket structures based on an environmental contour approach. In addition to the popular conditional distribution model, various classes of bivariate copulas are adopted to construct joint distributions of environmental variables. Analytical formulations of environmental contours based on various models are presented and discussed in this study. The design loads are examined by dynamic response analysis of jacket platform. Results suggest that the conditional model is not recommended for use in estimating design loads in sampling locations due to poor fitting results. Independent copula produces conservative design loads and the extreme response obtained using the conditional model are smaller than those determined by copulas. The suitability of a model for contour construction varies with the origin of wave data. This study provides a reference for the design load estimation of jacket structures and offers an alternative procedure to determine the design criteria for offshore structures.
文摘Ring footings are suitable for the structures like tall transmission towers, chimneys, silos and oil storages.These types of structures are susceptible to horizontal loads(wind load) in addition to their dead weight.In the literature, very little or no effort has been made to study the effect of ring footing resting on reinforced sand when subjected to eccentric, inclined and/or eccentric-inclined loadings. This paper aims to study the behavior of ring footing resting on loose sand and/or compacted randomly distributed fiberreinforced sand(RDFS) when subjected to eccentric(0 B, 0.05 B and 0.1 B, where B is the outer diameter of ring footing), inclined(0°,5°,10°, 15°,-5°,-10° and-15°)and eccentric-inclined loadings by using a finite element(FE) software PLAXIS 3 D. The behavior of ring footing is studied by using a dimensionless factor called reduction factor(RF). The numerical model used in the PLAXIS 3 D has been validated by conducting model plate load tests. Moreover, an empirical expression using regression analysis has been presented which will be helpful in plotting a load-settlement curve for the ring footing.
基金This study was supported by the National Natural Science Foundation of China(Grant No.51008313)the Sheng-hua Lie-ying Program of Central South University,a grant from the National High Technology Research and Development Program of China(863 Program,No.2009AA11Z101)the Joint Research Fund for Overseas Chinese,Hong Kong and Macao Young Scholars(No.50828801)from the National Natural Science Foundation of China。
文摘Load and resistance factors are generally obtained using the first order reliability method(FORM)in which the design point should be determined and derivative-based iterations used.In this article,the thirdmoment reliability index,based on the three-parameter lognormal(3P-lognormal)distribution,is investigated.A simple method based on the third-moment method for estimating load and resistance factors is then proposed,and a simple formula for the target mean resistance is also presented to avoid iterative computations.Unlike the currently used method,the proposed method can be used to determine load and resistance factors,even when the probability density functions(PDFs)of the basic random variables are not available.Moreover,the proposed method does not require the iterative computation of derivatives or any design points.Thus,the method provides a more convenient and effective way to estimate load and resistance factors in practical engineering applications.Numerical examples are presented to demonstrate the advantages of the proposed third moment method for determining load and resistance factors.
文摘Clustering approaches are one of the probabilistic load flow(PLF)methods in distribution networks that can be used to obtain output random variables,with much less computation burden and time than the Monte Carlo simulation(MCS)method.However,a challenge of the clustering methods is that the statistical characteristics of the output random variables are obtained with low accuracy.This paper presents a hybrid approach based on clustering and Point estimate methods.In the proposed approach,first,the sample points are clustered based on the𝑙-means method and the optimal agent of each cluster is determined.Then,for each member of the population of agents,the deterministic load flow calculations are performed,and the output variables are calculated.Afterward,a Point estimate-based PLF is performed and the mean and the standard deviation of the output variables are obtained.Finally,the statistical data of each output random variable are modified using the Point estimate method.The use of the proposed method makes it possible to obtain the statistical properties of output random variables such as mean,standard deviation and probabilistic functions,with high accuracy and without significantly increasing the burden of calculations.In order to confirm the consistency and efficiency of the proposed method,the 10-,33-,69-,85-,and 118-bus standard distribution networks have been simulated using coding in Python®programming language.In simulation studies,the results of the proposed method have been compared with the results obtained from the clustering method as well as the MCS method,as a criterion.