The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the effi...The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm,which hinders their application to high-dimensional engineering problems.To address these issues,this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps,the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm (SAGA) and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A'nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
Reasonable degree of the landscape pattern spatial distribution directly influences the sustainable use of regional land resources. Aiming at the unreasonable distribution of agricultural ecological landscape pattern,...Reasonable degree of the landscape pattern spatial distribution directly influences the sustainable use of regional land resources. Aiming at the unreasonable distribution of agricultural ecological landscape pattern, the deterioration of ecological environment, the Cellular Automaton (CA) principles were used to establish the optimized rules, such as landscape suitability rules, landscape prior rules and restraint conditions; and the standardization of the spatial data was realized by the competence of GIS spatial data handling and data spatial analysis, and finally, landscape pattern spatial optimization model was established with the support of MATLAB platform, i.e. LPSO Model. The spatial pattern optimization of agricultural landscape in west Jilin Province has been realized, which also laid a theoretical foundation for the proper spatial distribution of landscape pattern in west Jilin Province and realizing the sustainable agricultural development.展开更多
Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary a...Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms (EAs) and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.展开更多
In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model f...In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.展开更多
基于MATLAB矢量化的物质点法(material point method,MPM)框架,分析车身前防撞梁的碰撞冲击问题。MPM在每一迭代步将物理参数在物质点和背景网格间相互映射,使用MATLAB矢量化框架可以使用户在快速入门的同时保证求解效率,其应力更新采...基于MATLAB矢量化的物质点法(material point method,MPM)框架,分析车身前防撞梁的碰撞冲击问题。MPM在每一迭代步将物理参数在物质点和背景网格间相互映射,使用MATLAB矢量化框架可以使用户在快速入门的同时保证求解效率,其应力更新采用车身结构材料的弹塑性本构模型。前防撞梁碰撞冲击数值算例结果表明,MPM可以保证求解精度,同时矢量化技术可以大幅提高求解效率。展开更多
Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and ther...Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.展开更多
This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method fo...This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.展开更多
A topology optimization method based on the solid isotropic material with penalization interpolation scheme is utilized for designing gradient coils for use in magnetic resonance microscopy.Unlike the popular stream f...A topology optimization method based on the solid isotropic material with penalization interpolation scheme is utilized for designing gradient coils for use in magnetic resonance microscopy.Unlike the popular stream function method,the proposed method has design variables that are the distribution of conductive material.A voltage-driven transverse gradient coil is proposed to be used as micro-scale magnetic resonance imaging(MRI)gradient coils,thus avoiding introducing a coil-winding pattern and simplifying the coil configuration.The proposed method avoids post-processing errors that occur when the continuous current density is approximated by discrete wires in the stream function approach.The feasibility and accuracy of the method are verified through designing the z-gradient and y-gradient coils on a cylindrical surface.Numerical design results show that the proposed method can provide a new coil layout in a compact design space.展开更多
Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindric...Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindrical Bessel equation eigenvalue problem on general interval is given, second, by mean of compared method, we obtaining roots of characteristic equation with Matlab program is discussed.展开更多
This research extensively evaluates three leading mathematical software packages: Python, MATLAB, and Scilab, in the context of solving nonlinear systems of equations with five unknown variables. The study’s core obj...This research extensively evaluates three leading mathematical software packages: Python, MATLAB, and Scilab, in the context of solving nonlinear systems of equations with five unknown variables. The study’s core objectives include comparing software performance using standardized benchmarks, employing key performance metrics for quantitative assessment, and examining the influence of varying hardware specifications on software efficiency across HP ProBook, HP EliteBook, Dell Inspiron, and Dell Latitude laptops. Results from this investigation reveal insights into the capabilities of these software tools in diverse computing environments. On the HP ProBook, Python consistently outperforms MATLAB in terms of computational time. Python also exhibits a lower robustness index for problems 3 and 5 but matches or surpasses MATLAB for problem 1, for some initial guess values. In contrast, on the HP EliteBook, MATLAB consistently exhibits shorter computational times than Python across all benchmark problems. However, Python maintains a lower robustness index for most problems, except for problem 3, where MATLAB performs better. A notable challenge is Python’s failure to converge for problem 4 with certain initial guess values, while MATLAB succeeds in producing results. Analysis on the Dell Inspiron reveals a split in strengths. Python demonstrates superior computational efficiency for some problems, while MATLAB excels in handling others. This pattern extends to the robustness index, with Python showing lower values for some problems, and MATLAB achieving the lowest indices for other problems. In conclusion, this research offers valuable insights into the comparative performance of Python, MATLAB, and Scilab in solving nonlinear systems of equations. It underscores the importance of considering both software and hardware specifications in real-world applications. The choice between Python and MATLAB can yield distinct advantages depending on the specific problem and computational environment, providing guidance for researchers and practitioners in selecting tools for their unique challenges.展开更多
This paper presents a MATLAB implementation of the material-field series-expansion(MFSE)topology optimization method.The MFSE method uses a bounded material field with specified spatial correlation to represent the st...This paper presents a MATLAB implementation of the material-field series-expansion(MFSE)topology optimization method.The MFSE method uses a bounded material field with specified spatial correlation to represent the structural topology.With the series-expansion method for bounded fields,this material field is described with the characteristic base functions and the corresponding coefficients.Compared with the conventional density-based method,the MFSE method decouples the topological description and the finite element discretization,and greatly reduces the number of design variables after dimensionality reduction.Other features of this method include inherent control on structural topological complexity,crisp structural boundary description,mesh independence,and being free from the checkerboard pattern.With the focus on the implementation of the MFSE method,the present MATLAB code uses the maximum stiffness optimization problems solved with a gradientbased optimizer as examples.The MATLAB code consists of three parts,namely,the main program and two subroutines(one for aggregating the optimization constraints and the other about the method of moving asymptotes optimizer).The implementation of the code and its extensions to topology optimization problems with multiple load cases and passive elements are discussed in detail.The code is intended for researchers who are interested in this method and want to get started with it quickly.It can also be used as a basis for handling complex engineering optimization problems by combining the MFSE topology optimization method with non-gradient optimization algorithms without sensitivity information because only a few design variables are required to describe relatively complex structural topology and smooth structural boundaries using the MFSE method.展开更多
High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis mode...High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function(RBF) based on fuzzy clustering(ARFC) is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration(SLE) is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method(ARSM) proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.展开更多
Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is stron...Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is strongly perceived and may be a serious source of driver annoyance. Thus, many studies regarding the cavity noise mechanism and its reduction have already been conducted. In this work, a vibro-acoustic coupled analysis was conducted between a tire structure and air cavity. Using this analysis, we can more accurately simulate the tire noise performance in the region of the cavity resonance frequency. An analysis of the effects of variation of tire contour design factors was conducted, using design-of-experiments methods. Finally, a multi-objective optimization was performed using in-house codes to reduce the cavity noise level while minimizing the loss of other performances, such as diminished ride comfort and handling caused by the variations of contour. As a result of this optimization, an optimized contour shape was derived, which satisfied the multi-objective performances.展开更多
The metallic antenna design problem can be treated as a problem to find the optimal distribution of conductive material in a certain domain. Although this problem is well suited for topology optimization method, the v...The metallic antenna design problem can be treated as a problem to find the optimal distribution of conductive material in a certain domain. Although this problem is well suited for topology optimization method, the volumetric distribution of conductive material based on 3D finite element method (FEM) has been known to cause numerical bottlenecks such as the skin depth issue, meshed 'air regions' and other numerical problems. In this paper a topology optimization method based on the method of moments (MoM) for configuration design of planar metallic antenna was proposed. The candidate structure of the planar metallic antenna was approximately considered as a resistance sheet with position-dependent impedance. In this way, the electromagnetic property of the antenna can be analyzed easily by using the MoM to solve the radiation problem of the resistance sheet in a finite domain. The topology of the antenna was depicted with the distribution of the impedance related to the design parameters or relative densities. The conductive material (metal) was assumed to have zero impedance, whereas the non-conductive material was simulated as a material with a finite but large enough impedance. The interpolation function of the impedance between conductive material and non-conductive material was taken as a tangential function. The design of planar metallic antenna was optimized for maximizing the efficiency at the target frequency. The results illustrated the effectiveness of the method.展开更多
基金supported by the Innovation Fund Project of the Gansu Education Department(Grant No.2021B-099).
文摘The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm,which hinders their application to high-dimensional engineering problems.To address these issues,this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps,the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm (SAGA) and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A'nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
文摘Reasonable degree of the landscape pattern spatial distribution directly influences the sustainable use of regional land resources. Aiming at the unreasonable distribution of agricultural ecological landscape pattern, the deterioration of ecological environment, the Cellular Automaton (CA) principles were used to establish the optimized rules, such as landscape suitability rules, landscape prior rules and restraint conditions; and the standardization of the spatial data was realized by the competence of GIS spatial data handling and data spatial analysis, and finally, landscape pattern spatial optimization model was established with the support of MATLAB platform, i.e. LPSO Model. The spatial pattern optimization of agricultural landscape in west Jilin Province has been realized, which also laid a theoretical foundation for the proper spatial distribution of landscape pattern in west Jilin Province and realizing the sustainable agricultural development.
文摘Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms (EAs) and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.
基金This work was supported by the National Natural Science Foundation of China(10071037)
文摘In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.
文摘基于MATLAB矢量化的物质点法(material point method,MPM)框架,分析车身前防撞梁的碰撞冲击问题。MPM在每一迭代步将物理参数在物质点和背景网格间相互映射,使用MATLAB矢量化框架可以使用户在快速入门的同时保证求解效率,其应力更新采用车身结构材料的弹塑性本构模型。前防撞梁碰撞冲击数值算例结果表明,MPM可以保证求解精度,同时矢量化技术可以大幅提高求解效率。
基金University Putra Malaysia under Putra Grant No.9531200。
文摘Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.
基金supported by the National Natural Science Foundation of China(12171106)the Natural Science Foundation of Guangxi Province(2020GXNSFDA238017 and 2018GXNSFFA281007)the Shanghai Sailing Program(21YF1430300)。
文摘This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51675506 and 51275504)the German Research Foundation(DFG)(Grant Nos.#ZA 422/5-1 and#ZA 422/6-1)
文摘A topology optimization method based on the solid isotropic material with penalization interpolation scheme is utilized for designing gradient coils for use in magnetic resonance microscopy.Unlike the popular stream function method,the proposed method has design variables that are the distribution of conductive material.A voltage-driven transverse gradient coil is proposed to be used as micro-scale magnetic resonance imaging(MRI)gradient coils,thus avoiding introducing a coil-winding pattern and simplifying the coil configuration.The proposed method avoids post-processing errors that occur when the continuous current density is approximated by discrete wires in the stream function approach.The feasibility and accuracy of the method are verified through designing the z-gradient and y-gradient coils on a cylindrical surface.Numerical design results show that the proposed method can provide a new coil layout in a compact design space.
文摘Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindrical Bessel equation eigenvalue problem on general interval is given, second, by mean of compared method, we obtaining roots of characteristic equation with Matlab program is discussed.
文摘This research extensively evaluates three leading mathematical software packages: Python, MATLAB, and Scilab, in the context of solving nonlinear systems of equations with five unknown variables. The study’s core objectives include comparing software performance using standardized benchmarks, employing key performance metrics for quantitative assessment, and examining the influence of varying hardware specifications on software efficiency across HP ProBook, HP EliteBook, Dell Inspiron, and Dell Latitude laptops. Results from this investigation reveal insights into the capabilities of these software tools in diverse computing environments. On the HP ProBook, Python consistently outperforms MATLAB in terms of computational time. Python also exhibits a lower robustness index for problems 3 and 5 but matches or surpasses MATLAB for problem 1, for some initial guess values. In contrast, on the HP EliteBook, MATLAB consistently exhibits shorter computational times than Python across all benchmark problems. However, Python maintains a lower robustness index for most problems, except for problem 3, where MATLAB performs better. A notable challenge is Python’s failure to converge for problem 4 with certain initial guess values, while MATLAB succeeds in producing results. Analysis on the Dell Inspiron reveals a split in strengths. Python demonstrates superior computational efficiency for some problems, while MATLAB excels in handling others. This pattern extends to the robustness index, with Python showing lower values for some problems, and MATLAB achieving the lowest indices for other problems. In conclusion, this research offers valuable insights into the comparative performance of Python, MATLAB, and Scilab in solving nonlinear systems of equations. It underscores the importance of considering both software and hardware specifications in real-world applications. The choice between Python and MATLAB can yield distinct advantages depending on the specific problem and computational environment, providing guidance for researchers and practitioners in selecting tools for their unique challenges.
基金The authors acknowledge the support of the National Key R&D Program of China(Grant No.2017YFB0203604)the National Natural Science Foundation of China(Grant Nos.11902064 and 11772077)the Liaoning Revitalization Talents Program,China(Grant No.XLYC1807187).
文摘This paper presents a MATLAB implementation of the material-field series-expansion(MFSE)topology optimization method.The MFSE method uses a bounded material field with specified spatial correlation to represent the structural topology.With the series-expansion method for bounded fields,this material field is described with the characteristic base functions and the corresponding coefficients.Compared with the conventional density-based method,the MFSE method decouples the topological description and the finite element discretization,and greatly reduces the number of design variables after dimensionality reduction.Other features of this method include inherent control on structural topological complexity,crisp structural boundary description,mesh independence,and being free from the checkerboard pattern.With the focus on the implementation of the MFSE method,the present MATLAB code uses the maximum stiffness optimization problems solved with a gradientbased optimizer as examples.The MATLAB code consists of three parts,namely,the main program and two subroutines(one for aggregating the optimization constraints and the other about the method of moving asymptotes optimizer).The implementation of the code and its extensions to topology optimization problems with multiple load cases and passive elements are discussed in detail.The code is intended for researchers who are interested in this method and want to get started with it quickly.It can also be used as a basis for handling complex engineering optimization problems by combining the MFSE topology optimization method with non-gradient optimization algorithms without sensitivity information because only a few design variables are required to describe relatively complex structural topology and smooth structural boundaries using the MFSE method.
基金supported by National Natural Science Foundation of China (Grant Nos. 50875024,51105040)Excellent Young Scholars Research Fund of Beijing Institute of Technology,China (Grant No.2010Y0102)Defense Creative Research Group Foundation of China(Grant No. GFTD0803)
文摘High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function(RBF) based on fuzzy clustering(ARFC) is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration(SLE) is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method(ARSM) proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.
文摘Cavity resonance noise of passenger car tires is generated by interacting excitation between a tire structure and the fill gas (air), and generally lies in a frequency range of 200?250 Hz. As such, this noise is strongly perceived and may be a serious source of driver annoyance. Thus, many studies regarding the cavity noise mechanism and its reduction have already been conducted. In this work, a vibro-acoustic coupled analysis was conducted between a tire structure and air cavity. Using this analysis, we can more accurately simulate the tire noise performance in the region of the cavity resonance frequency. An analysis of the effects of variation of tire contour design factors was conducted, using design-of-experiments methods. Finally, a multi-objective optimization was performed using in-house codes to reduce the cavity noise level while minimizing the loss of other performances, such as diminished ride comfort and handling caused by the variations of contour. As a result of this optimization, an optimized contour shape was derived, which satisfied the multi-objective performances.
基金supported by the National Natural Science Foundation of China (Grants 11332004, 11372063, and 11572073)111 Project (Grant B14013)the Fundamental Research Funds for the Central Universities (Grant DUT15ZD101)
文摘The metallic antenna design problem can be treated as a problem to find the optimal distribution of conductive material in a certain domain. Although this problem is well suited for topology optimization method, the volumetric distribution of conductive material based on 3D finite element method (FEM) has been known to cause numerical bottlenecks such as the skin depth issue, meshed 'air regions' and other numerical problems. In this paper a topology optimization method based on the method of moments (MoM) for configuration design of planar metallic antenna was proposed. The candidate structure of the planar metallic antenna was approximately considered as a resistance sheet with position-dependent impedance. In this way, the electromagnetic property of the antenna can be analyzed easily by using the MoM to solve the radiation problem of the resistance sheet in a finite domain. The topology of the antenna was depicted with the distribution of the impedance related to the design parameters or relative densities. The conductive material (metal) was assumed to have zero impedance, whereas the non-conductive material was simulated as a material with a finite but large enough impedance. The interpolation function of the impedance between conductive material and non-conductive material was taken as a tangential function. The design of planar metallic antenna was optimized for maximizing the efficiency at the target frequency. The results illustrated the effectiveness of the method.