期刊文献+
共找到2,896篇文章
< 1 2 145 >
每页显示 20 50 100
A Deep Learning Approach to Shape Optimization Problems for Flexoelectric Materials Using the Isogeometric Finite Element Method
1
作者 Yu Cheng Yajun Huang +3 位作者 Shuai Li Zhongbin Zhou Xiaohui Yuan Yanming Xu 《Computer Modeling in Engineering & Sciences》 SCIE EI 2024年第5期1935-1960,共26页
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization... A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization. 展开更多
关键词 Shape optimization deep learning flexoelectric structure finite element method isogeometric
下载PDF
A Cell-Based Linear Smoothed Finite Element Method for Polygonal Topology Optimization
2
作者 Changkye Lee Sundararajan Natarajan +1 位作者 Seong-Hoon Kee Jurng-Jae Yee 《Computer Modeling in Engineering & Sciences》 SCIE EI 2022年第6期1615-1634,共20页
The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basi... The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy. 展开更多
关键词 Smoothed finite element method linear smoothing function topology optimization solid isotropic material with penalization(SIMP) polygonal finite element cell
下载PDF
Gradient Recovery Based Two-Grid Finite Element Method for Parabolic Integro-Differential Optimal Control Problems
3
作者 Miao Yang 《Journal of Applied Mathematics and Physics》 2024年第8期2849-2865,共17页
In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and ... In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results. 展开更多
关键词 Optimal Control Problem Gradient Recovery Two-Grid finite element method
下载PDF
SHAPE OPTIMIZATION DESIGN OF CRUSH TOOTH USING FINITE ELEMENT METHOD
4
作者 Gong Kaijun Yang Qingfu Zhang Junxia (Taiyuan University of Technology) Zhu Wennong (Jilin University of Technology) 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 1994年第2期130-134,共17页
The load condition of the mineral sizer is discussed. Acoording to three typical load cases presented, FEM and optimization models are established. Computed the model problems by use of ANSYS software and sequential u... The load condition of the mineral sizer is discussed. Acoording to three typical load cases presented, FEM and optimization models are established. Computed the model problems by use of ANSYS software and sequential unconstrained minimization technique, the stress distribution on c rush tooth and it's optimal shape are aquired. After analyzing the results, the relation between stress distribution and structural parameters is presented. 展开更多
关键词 CRUSH Optimal design finite element method
全文增补中
OPTIMIZATION OF AUTOBODY PANEL STAMPING PROCESS BASED ON DYNAMIC EXPLICIT FINITE ELEMENT METHOD 被引量:6
5
作者 X. G. Bao D. N. He +2 位作者 D. Lu C. X. Li J. L. Cheng and J. Y. Jiang( 1) National Mold and Dies CAD Engineering Research Center, Shanghai 200030, China 2) Shanghai Volkswagen Automotive Company Ltd.(SVW), Shanghai 201805, China) 《Acta Metallurgica Sinica(English Letters)》 SCIE EI CAS CSCD 2000年第1期387-393,共7页
Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this... Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained. 展开更多
关键词 dynamic explicit finite element method autobody panel optimization stamping process
下载PDF
FORMING MICHELL TRUSS IN THREE-DIMENSIONS BY FINITE ELEMENT METHOD
6
作者 周克民 李俊峰 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2005年第3期381-388,共8页
The finite element method to form Michell truss in three_dimensions is presented.The orthotropic composite with fiber_reinforcement is employed as the material model to simulate Michell truss.The orientation and densi... The finite element method to form Michell truss in three_dimensions is presented.The orthotropic composite with fiber_reinforcement is employed as the material model to simulate Michell truss.The orientation and densities of fibers at nodes are taken as basic design variables.The stresses and strains at nodes are calculated by finite element method.An iteration scheme is suggested to adjust the orientations of fibers to be along the orientations of principal stresses, and the densities of fibers according to the strains in the orientations of fibers.The strain field satisfying Michell criteria and truss_like continuum are achieved after several iterations. Lastly, the Michell truss is showed by continuous lines, which are formed according to the orientations of fibers at nodes. Several examples are used to demonstrate the efficiency of the presented approach. 展开更多
关键词 structural optimization finite element method topology optimization Michell truss stress constraint
下载PDF
Analytical Model of Nonlinear Semi-rigid Frames Using Finite Element Method
7
作者 Shahrin Mohammad Ahmad Baharuddin Abd Rahman +1 位作者 Cher Siang Tan Yeong Huei Lee 《Journal of Architectural Environment & Structural Engineering Research》 2020年第4期22-27,共6页
Performance-based design for a constructional steel frame in nonlinear-plastic region requires an improvement in order to achieve a reliable structural analysis.The need to explicitly consider the nonlinear behaviour ... Performance-based design for a constructional steel frame in nonlinear-plastic region requires an improvement in order to achieve a reliable structural analysis.The need to explicitly consider the nonlinear behaviour of structures makes the numerical modelling approach much more favourable than expensive and potentially dangerous experimental work.The parameters considered in the analysis are not limited to the linear change of geometry and material yielding,but also include the effect of large deformations,geometrical imperfections,load eccentricities,residual stresses,strain-unloading,and the nonlinear boundary conditions.Such analysis requires the use of accurate mathematical modelling and effective numerical procedures for solving equations of equilibrium.With that in mind,this paper presents the mathematical formulations and finite element procedures of nonlinear inelastic steel frame analysis with quasi-static semi-rigid connections.Verification and validation of the developed analytical procedures are conducted and good agreements are obtained.It is an approach that enables the structural behaviour of constructional steel frames to be traced throughout the entire range of loading until failure.It also provides information on the derivation of the structural analysis by using finite element method. 展开更多
关键词 finite element method NONLINEARITY Steel frame Semi-rigid connection
下载PDF
Error estimates of H^1-Galerkin mixed finite element method for Schrdinger equation 被引量:28
8
作者 LIU Yang LI Hong WANG Jin-feng 《Applied Mathematics(A Journal of Chinese Universities)》 SCIE CSCD 2009年第1期83-89,共7页
An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t... An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition. 展开更多
关键词 H1-Galerkin mixed finite element method Schrdinger equation LBB condition optimal error estimates
下载PDF
A Priori and A Posteriori Error Estimates of Streamline Diffusion Finite Element Method for Optimal Control Problem Governed by Convection Dominated Diffusion Equation 被引量:5
9
作者 Ningning Yan Zhaojie Zhou 《Numerical Mathematics(Theory,Methods and Applications)》 SCIE 2008年第3期297-320,共24页
In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existenc... In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results. 展开更多
关键词 Constrained optimal control problem convection dominated diffusion equation stream-line diffusion finite element method a priori error estimate a posteriori error estimate.
下载PDF
ADAPTIVE FINITE ELEMENT METHOD BASED ON OPTIMAL ERROR ESTIMATES FOR LINEAR ELLIPTIC PROBLEMS ON CONCAVE CORNER DOMAINS 被引量:1
10
作者 汤雁 《Transactions of Tianjin University》 EI CAS 2001年第1期64-67,共4页
An adaptive finite element procedure designed for specific computational goals is presented,using mesh refinement strategies based on optimal or nearly optimal a priori error estimates for the finite element method an... An adaptive finite element procedure designed for specific computational goals is presented,using mesh refinement strategies based on optimal or nearly optimal a priori error estimates for the finite element method and using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions.The proposed procedure is analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that the method can generate the correct type of refinements and lead to the desired control under consideration. 展开更多
关键词 adaptive finite element method priori error estimate optimal mesh design
全文增补中
Optimal design of butterfly-shaped linear ultrasonic motor using finite element method and response surface methodology 被引量:9
11
作者 时运来 陈超 赵淳生 《Journal of Central South University》 SCIE EI CAS 2013年第2期393-404,共12页
A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to ... A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8. 展开更多
关键词 linear ultrasonic motor PIEZOELECTRIC optimal design response surface methodology finite element method
下载PDF
Self-adaptive strategy for one-dimensional finite element method based on EEP method with optimal super-convergence order 被引量:4
12
作者 袁驷 邢沁妍 +1 位作者 王旭 叶康生 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2008年第5期591-602,共12页
Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite ele... Based on the newly-developed element energy projection (EEP) method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method (FEM) is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach. 展开更多
关键词 finite element method (FEM) self-adaptive solution super-convergence optimal convergence order element energy projection condensed shape functions
下载PDF
Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients 被引量:2
13
作者 Andrew Lundberg Pengtao Sun +1 位作者 Cheng Wang Chen-song Zhang 《Computer Modeling in Engineering & Sciences》 SCIE EI 2019年第4期35-62,共28页
The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully disc... The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated. 展开更多
关键词 TRANSIENT STOKES interface problem JUMP COEFFICIENTS DISTRIBUTED LAGRANGE multiplier fictitious domain method mixed finite element an optimal error estimate stability
下载PDF
Dynamic analysis and optimization of car's body-in-white based on superelement approach 被引量:3
14
作者 张建润 孙庆鸿 陈南 《Journal of Southeast University(English Edition)》 EI CAS 2004年第3期324-327,共4页
An approach based on the superelement theory is proposed, and it is applied to model the car's body-in-white as well as to dynamic simulation and optimization. This approach can improve the calculation speed and d... An approach based on the superelement theory is proposed, and it is applied to model the car's body-in-white as well as to dynamic simulation and optimization. This approach can improve the calculation speed and do the dynamic optimization among substructures respectively in the car's body design. To meet the car's design of harshness, a dynamic optimal design model, based on the mean square of vertical displacement response at two points of the car floor, is proposed. Satisfactory results are achieved. 展开更多
关键词 finite element method Mathematical models optimization Structural design
下载PDF
Independent continuous mapping for topological optimization of frame structures 被引量:10
15
作者 Yunkang Sui Jiazheng Du Yingqiao Guo 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2006年第6期611-619,共9页
Based on the Independent Continuous Mapping method (ICM), a topological optimization model with continuous topological variables is built by introducing three filter functions for element weight, element allowable s... Based on the Independent Continuous Mapping method (ICM), a topological optimization model with continuous topological variables is built by introducing three filter functions for element weight, element allowable stress and element stiffness, which transform the 0-1 type discrete topological variables into continuous topological variables between 0 and 1. Two methods for the filter functions are adopted to avoid the structural singularity and recover falsely deleted elements: the weak material element method and the tiny section element method. Three criteria (no structural singularity, no violated constraints and no change of structural weight) are introduced to judge iteration convergence. These criteria allow finding an appropriate threshold by adjusting a discount factor in the iteration procedure. To improve the efficiency, the original optimization model is transformed into a dual problem according to the dual theory and solved in its dual space. By using MSC/Nastran as the structural solver and MSC/Patran as the developing platform, a topological optimization software of frame structures is accomplished. Numerical examples show that the ICM method is very efficient for the topological optimization of frame structures. 展开更多
关键词 frame structures Topological optimization ICM method Filter functions element elimination
下载PDF
Permanent Magnet Shape Optimization Method for PMSM Air Gap Flux Density Harmonics Reduction 被引量:6
16
作者 Chengsi Liu Yongxiang Xu +2 位作者 Jibin Zou Guodong Yu Liang Zhuo 《CES Transactions on Electrical Machines and Systems》 CSCD 2021年第4期284-290,共7页
This paper proposed a permanent magnet optimization method to suppress the air gap flux density harmonic of permanent magnet synchronous motor(PMSM).The method corrected the effective air gap length of the motor,calcu... This paper proposed a permanent magnet optimization method to suppress the air gap flux density harmonic of permanent magnet synchronous motor(PMSM).The method corrected the effective air gap length of the motor,calculated the magnetization length of the permanent in the case of parallel magnetization,and took the influence of the permanent magnet relative permeability into consideration.Based on these works,for a given sinusoidal air gap flux density waveform,the corresponding structural parameters can be calculated,so as to achieve the optimization of the permanent magnet.By using this method to optimize the shape of the magnet,the fundamental wave of the air gap flux density can be retained to the greatest extent,so as to eliminate harmonics and maintain the output capacity at the same time.The feasibility and accuracy of the method have been verified by finite element analysis(FEA)and prototype machine experiment.This method is simple and time-saving,and has a satisfactory accuracy,which provides a reference method for permanent magnet optimization of PMSM. 展开更多
关键词 Permanent magnet synchronous motor(PMSM) analytical method finite element analysis(FEA) magnet shape optimization.
下载PDF
A Study on the Multi-Objective Optimization Method of Brackets in Ship Structures
17
作者 LIU Fan HU Yu-meng +2 位作者 FENG Guo-qing ZHAO Wei-dong ZHANG Ming 《China Ocean Engineering》 SCIE EI CSCD 2022年第2期208-222,共15页
The shape and size optimization of brackets in hull structures was conducted to achieve the simultaneous reduction of mass and high stress,where the parametric finite element model was built based on Patran Command La... The shape and size optimization of brackets in hull structures was conducted to achieve the simultaneous reduction of mass and high stress,where the parametric finite element model was built based on Patran Command Language codes.The optimization procedure was executed on Isight platform,on which the linear dimensionless method was introduced to establish the weighted multi-objective function.The extreme processing method was applied and proved effective to normalize the objectives.The bracket was optimized under the typical single loads and design waves,accompanied by the different proportions of weights in the objective function,in which the safety factor function was further established,including yielding,buckling,and fatigue strength,and the weight minimization and safety maximization of the bracket were obtained.The findings of this study illustrate that the dimensionless objectives share equal contributions to the multi-objective function,which enhances the role of weights in the optimization. 展开更多
关键词 BRACKETS parametric finite element model multi-objective optimization extreme processing method safety factor function weighted multi-objective function
下载PDF
A Simple and Efficient Structural Topology Optimization Implementation Using Open-Source Software for All Steps of the Algorithm:Modeling, Sensitivity Analysis and Optimization
18
作者 Rafael Marin Ferro Renato Pavanello 《Computer Modeling in Engineering & Sciences》 SCIE EI 2023年第8期1371-1397,共27页
This work analyzes the implementation of a continuous method of structural topology optimization(STO)using open-source software for all stages of the topology optimization problem:modeling,sensitivity analysis and opt... This work analyzes the implementation of a continuous method of structural topology optimization(STO)using open-source software for all stages of the topology optimization problem:modeling,sensitivity analysis and optimization.Its implementation involves three main components:numerical analysis using the Finite Element Method(FEM),sensitivity analysis using an Adjoint method and an optimization solver.In order to allow the automated numerical solution of Partial Differential Equations(PDEs)and perform a sensitivity analysis,FEniCS and Dolfin Adjoint software are used as tools,which are open-source code.For the optimization process,Ipopt(Interior Point OPTimizer)is used,which is a software package for nonlinear optimization scale designed to find(local)solutions of mathematical optimization problems.The topological optimization method used is based on the SIMP-Solid IsotropicMaterial with Penalization interpolation.The considered problem is the minimization of compliance/maximization of stiffness,considering the examples of recurrent structures in the literature in 2D and 3D.A density filtering algorithm based on Helmholtz formulation is used.The complete code involves 51 lines of programming and is presented and commented in detail in this article. 展开更多
关键词 Structural topology optimization finite element method FEniCS Dolfin Adjoint
下载PDF
Topology optimization design of the main structure of 6-DOF manipulator based on the variable density method
19
作者 Yujian RUI Yongming LIU +4 位作者 Zhen ZHANG Zhijian TU Jian WU Neng WEI Zhuanzhe ZHAO 《Mechanical Engineering Science》 2022年第1期14-21,共8页
In order to solve the problems of too large mass,too complex structure and poor flexibility of the 6 DOF manipulator,the topological optimization theory based on variable density method is applied to the 6-DOF manipul... In order to solve the problems of too large mass,too complex structure and poor flexibility of the 6 DOF manipulator,the topological optimization theory based on variable density method is applied to the 6-DOF manipulator,the topology optimization of the main structural components of the manipulator is carried out with the help of the finite element software ANSYS,and the optimized structure is simplified according to the density distribution of the units and the requirements of manufacturability.the results are compared and analysed by static mechanics.It shows that the whole mass of the 6-DOF manipulator is reduced by 47.23%without changing the original mechanical properties after topological optimization,and the optimized model can meet the requirements of manufacturability,the optimization effect is signifcant,which can be used as a reference for the structure optimization of the 6-DOF manipulator. 展开更多
关键词 Tariable density method topology optimization 6-DOF manipulator finite element analysis
下载PDF
Cavity noise sensitivity analysis of tire contour design factors and application of contour optimization methodology 被引量:7
20
作者 KIM Seong-rae SUNG Ki-deug +1 位作者 LEE Dong-woo HUH Sun-chul 《Journal of Central South University》 SCIE EI CAS 2012年第8期2386-2393,共8页
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. 展开更多
关键词 tire cavity noise vibro-acoustic coupled analysis finite element method boundary element method (BEM) sensitivity optimization
下载PDF
上一页 1 2 145 下一页 到第
使用帮助 返回顶部