期刊文献+
共找到136篇文章
< 1 2 7 >
每页显示 20 50 100
Topology Optimization of Orthotropic Materials Using the Improved Element-Free Galerkin (IEFG) Method
1
作者 Wenna He Yichen Yang +1 位作者 Dongqiong Liang Heng Cheng 《Computers, Materials & Continua》 2025年第4期1415-1437,共23页
In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s... In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials. 展开更多
关键词 Solid isotropic microstructures with penalization method variable density method sensitivity analysis improved element-free Galerkin method meshless method
在线阅读 下载PDF
Advances in the Improved Element-Free Galerkin Methods:A Comprehensive Review
2
作者 Heng Cheng YichenYang Yumin Cheng 《Computer Modeling in Engineering & Sciences》 2025年第12期2853-2894,共42页
The element-free Galerkin(EFG)method,which constructs shape functions via moving least squares(MLS)approximation,represents a fundamental and widely studied meshless method in numerical computation.Although it achieve... The element-free Galerkin(EFG)method,which constructs shape functions via moving least squares(MLS)approximation,represents a fundamental and widely studied meshless method in numerical computation.Although it achieves high computational accuracy,the shape functions are more complex than those in the conventional finite element method(FEM),resulting in great computational requirements.Therefore,improving the computational efficiency of the EFG method represents an important research direction.This paper systematically reviews significant contributions fromdomestic and international scholars in advancing the EFGmethod.Including the improved element-free Galerkin(IEFG)method,various interpolating EFG methods,four distinct complex variable EFG methods,and a series of dimension splitting meshless methods.In the numerical examples,the effectiveness and efficiency of the three methods are validated by analyzing the solutions of the IEFG method for 3D steadystate anisotropic heat conduction,3D elastoplasticity,and large deformation problems,as well as the performance of two-dimensional splitting meshless methods in solving the 3D Helmholtz equation. 展开更多
关键词 Meshless method improved element-free Galerkin method singular weight function nonsingular weight function interpolating element-free Galerkin method complex variable element-free Galerkin method dimension splitting method dimension splitting meshless method
在线阅读 下载PDF
GLOBAL CONVERGENCE OF A CLASS OF OPTIMALLY CONDITIONED SSVM METHODS
3
作者 杨正方 夏爱生 +1 位作者 韩立兴 刘光辉 《Transactions of Tianjin University》 EI CAS 1997年第1期73-76,共4页
This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are glob... This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions. 展开更多
关键词 optimally conditioned self scaling variable metric methods global convergence unconstrained optimization
在线阅读 下载PDF
Improved finite difference method for pressure distribution of aerostatic bearing 被引量:4
4
作者 郑书飞 蒋书运 《Journal of Southeast University(English Edition)》 EI CAS 2009年第4期501-505,共5页
An improved finite difference method (FDM)is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aero... An improved finite difference method (FDM)is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses (less than 2 Win)are considered, the improved calculation method still yields a result and converges fast. 展开更多
关键词 aerostatic bearing: pressure distribution: Reynolds equation: finite difference method variable step size
在线阅读 下载PDF
On the feasibility of variable separation method based on Hamiltonian system for plane magnetoelectroelastic solids 被引量:5
5
作者 侯国林 阿拉坦仓 《Chinese Physics B》 SCIE EI CAS CSCD 2008年第8期2753-2758,共6页
The eigenvalue problem of an infinite-dimensional Hamiltonian operator appearing in the isotropic plane magnetoelectroelastic solids is studied. First, all the eigenvalues and their eigenfunctions in a rectangular dom... The eigenvalue problem of an infinite-dimensional Hamiltonian operator appearing in the isotropic plane magnetoelectroelastic solids is studied. First, all the eigenvalues and their eigenfunctions in a rectangular domain are solved directly. Then the completeness of the eigenfunction system is proved, which offers a theoretic guarantee of the feasibility of variable separation method based on a Hamiltonian system for isotropic plane magnetoelectroelastic solids. Finally, the general solution for the equation in the rectangular domain is obtained by using the symplectic Fourier expansion method. 展开更多
关键词 magnetoelectroelastic solid variable separation method COMPLETENESS general solution
原文传递
Physical Parameter Identification Method Based on Modal Analysis for Two-axis On-road Vehicles:Theory and Simulation 被引量:3
6
作者 ZHENG Minyi ZHANG Bangji +1 位作者 ZHANG Jie ZHANG Nong 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 2016年第4期756-764,共9页
Physical parameters are very important for vehicle dynamic modeling and analysis.However,most of physical parameter identification methods are assuming some physical parameters of vehicle are known,and the other unkno... Physical parameters are very important for vehicle dynamic modeling and analysis.However,most of physical parameter identification methods are assuming some physical parameters of vehicle are known,and the other unknown parameters can be identified.In order to identify physical parameters of vehicle in the case that all physical parameters are unknown,a methodology based on the State Variable Method(SVM)for physical parameter identification of two-axis on-road vehicle is presented.The modal parameters of the vehicle are identified by the SVM,furthermore,the physical parameters of the vehicle are estimated by least squares method.In numerical simulations,physical parameters of Ford Granada are chosen as parameters of vehicle model,and half-sine bump function is chosen to simulate tire stimulated by impulse excitation.The first numerical simulation shows that the present method can identify all of the physical parameters and the largest absolute value of percentage error of the identified physical parameter is 0.205%;and the effect of the errors of additional mass,structural parameter and measurement noise are discussed in the following simulations,the results shows that when signal contains 30 d B noise,the largest absolute value of percentage error of the identification is 3.78%.These simulations verify that the presented method is effective and accurate for physical parameter identification of two-axis on-road vehicles.The proposed methodology can identify all physical parameters of 7-DOF vehicle model by using free-decay responses of vehicle without need to assume some physical parameters are known. 展开更多
关键词 Parameter identification free-decay response state variable method modal parameter physical parameter
在线阅读 下载PDF
A new complex variable meshless method for transient heat conduction problems 被引量:5
7
作者 王健菲 程玉民 《Chinese Physics B》 SCIE EI CAS CSCD 2012年第12期42-50,共9页
In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is pres... In this paper, based on the improved complex variable moving least-square (ICVMLS) approximation, a new complex variable meshless method (CVMM) for two-dimensional (2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper. 展开更多
关键词 meshless method improved complex variable moving least-square approximation com-plex variable meshless method transient heat conduction problem
原文传递
A fast integration method for translating-pulsating source Green's function in Bessho form 被引量:5
8
作者 Chao-bang YAO Wen-cai DONG 《Journal of Zhejiang University-Science A(Applied Physics & Engineering)》 SCIE EI CAS CSCD 2014年第2期108-119,共12页
The singularities and oscillatory performance of translating-pulsating source Green's function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitut... The singularities and oscillatory performance of translating-pulsating source Green's function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method (VSM), and steepest descent integration method (SDIM) were used to evaluate this type of Green's function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves. 展开更多
关键词 Translating-pulsating source Green's function Oscillatory performance False singularities point Steepest descentintegration method (SDIM) Variable substitution method (VSM)
原文传递
Solving mKdV-sinh-Gordon equation by a modified variable separated ordinary differential equation method 被引量:4
9
作者 谢元喜 《Chinese Physics B》 SCIE EI CAS CSCD 2009年第12期5123-5132,共10页
By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact sol... By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact solutions including some new formal solutions are successfully picked up for the mKdV-sinh-Gordon equation by this approach. 展开更多
关键词 modified variable separated ODE method mKdV-sinh-Gordon equation explicit andexact solution
原文传递
An improved complex variable element-free Galerkin method for two-dimensional elasticity problems 被引量:4
10
作者 Bai Fu-Nong Li Dong-Ming +1 位作者 Wang Jian-Fei Cheng Yu-Min 《Chinese Physics B》 SCIE EI CAS CSCD 2012年第2期56-65,共10页
In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar... In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares (CVMLS) approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin (EFG) method and the ICVMLS approximation, the improved complex variable element-free Galerkin (ICVEFG) method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method. 展开更多
关键词 meshless method improved complex variable moving least-squares approximation improved complex variable element-free Galerkin method ELASTICITY
原文传递
A Combined Shape and Topology Optimization Based on Isogeometric Boundary Element Method for 3D Acoustics 被引量:2
11
作者 Jie Wang Fuhang Jiang +1 位作者 Wenchang Zhao Haibo Chen 《Computer Modeling in Engineering & Sciences》 SCIE EI 2021年第5期645-681,共37页
A combined shape and topology optimization algorithm based on isogeometric boundary element method for 3D acoustics is developed in this study.The key treatment involves using adjoint variable method in shape sensitiv... A combined shape and topology optimization algorithm based on isogeometric boundary element method for 3D acoustics is developed in this study.The key treatment involves using adjoint variable method in shape sensitivity analysis with respect to non-uniform rational basis splines control points,and in topology sensitivity analysis with respect to the artificial densities of sound absorption material.OpenMP tool in Fortran code is adopted to improve the efficiency of analysis.To consider the features and efficiencies of the two types of optimization methods,this study adopts a combined iteration scheme for the optimization process to investigate the simultaneous change of geometry shape and distribution of material to achieve better noise control.Numerical examples,such as sound barrier,simple tank,and BeTSSi submarine,are performed to validate the advantage of combined optimization in noise reduction,and to demonstrate the potential of the proposed method for engineering problems. 展开更多
关键词 Combined shape and topology optimization isogeometric boundary element method shape sensitivity analysis topology sensitivity analysis adjoint variable method sound absorption material
在线阅读 下载PDF
A new complex variable element-free Galerkin method for two-dimensional potential problems 被引量:4
12
作者 程玉民 王健菲 白福浓 《Chinese Physics B》 SCIE EI CAS CSCD 2012年第9期43-52,共10页
In this paper, based on the element-free Galerkin (EFG) method and the improved complex variable moving least- square (ICVMLS) approximation, a new meshless method, which is the improved complex variable element-f... In this paper, based on the element-free Galerkin (EFG) method and the improved complex variable moving least- square (ICVMLS) approximation, a new meshless method, which is the improved complex variable element-free Galerkin (ICVEFG) method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square (CVMLS) approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method. 展开更多
关键词 meshless method improved complex variable moving least-square approximation im- proved complex variable element-free Galerkin method potential problem
原文传递
The shape optimization of the arterial graft design by level set methods 被引量:1
13
作者 JIANG Dong HAN Dan-fu HU Xian-liang 《Applied Mathematics(A Journal of Chinese Universities)》 SCIE CSCD 2016年第2期205-218,共14页
The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is on... The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is one of the important aspects in the arterial graft design problem since it affects the flow of the blood significantly. As an attractive design tool for this problem, level set methods are quite efficient for obtaining better shape of the graft. In this paper, a cubic spline level set method and a radial basis function level set method are designed to solve the arterial graft design problem. In both approaches, the shape of the arterial graft is implicitly tracked by the zero-level contour of a level set function and a high level of smoothness of the graft is achieved. Numerical results show the efficiency of the algorithms in the arterial graft design. 展开更多
关键词 Arterial graft Shape optimization Level set method Adjoint variable method.
在线阅读 下载PDF
The complex variable meshless local Petrov-Galerkin method of solving two-dimensional potential problems 被引量:1
14
作者 杨秀丽 戴保东 张伟伟 《Chinese Physics B》 SCIE EI CAS CSCD 2012年第10期49-55,共7页
Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential proble... Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square(MLS) approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin(MLPG) method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method. 展开更多
关键词 meshless method complex variable moving least-square method complex variable meshless local Petrov-Galerkin method potential problems
原文传递
The complex variable reproducing kernel particle method for two-dimensional elastodynamics 被引量:2
15
作者 陈丽 程玉民 《Chinese Physics B》 SCIE EI CAS CSCD 2010年第9期59-70,共12页
On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is present... On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM. 展开更多
关键词 meshless method reproducing kernel particle method complex variable reproducing kernel particle method elastodvnamics
原文传递
New complex variable meshless method for advection-diffusion problems 被引量:1
16
作者 王健菲 程玉民 《Chinese Physics B》 SCIE EI CAS CSCD 2013年第3期92-98,共7页
In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equi... In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency. 展开更多
关键词 meshless method improved complex variable moving least-square approximation improved complex variable meshless method advection-diffusion problem
原文传递
Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems 被引量:2
17
作者 陈丽 马和平 程玉民 《Chinese Physics B》 SCIE EI CAS CSCD 2013年第5期67-74,共8页
In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE metho... In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method. 展开更多
关键词 complex variable reproducing kernel particle method finite element method combined method transient heat conduction
原文传递
Fuzzy norm method for evaluating random vibration of airborne platform from limited PSD data 被引量:7
18
作者 Wang Zhongyu Wang Yanqing +1 位作者 Wang Qian Zhang Jianjun 《Chinese Journal of Aeronautics》 SCIE EI CAS CSCD 2014年第6期1442-1450,共9页
For random vibration of airborne platform, the accurate evaluation is a key indicator to ensure normal operation of airborne equipment in flight. However, only limited power spectral density(PSD) data can be obtaine... For random vibration of airborne platform, the accurate evaluation is a key indicator to ensure normal operation of airborne equipment in flight. However, only limited power spectral density(PSD) data can be obtained at the stage of flight test. Thus, those conventional evaluation methods cannot be employed when the distribution characteristics and priori information are unknown. In this paper, the fuzzy norm method(FNM) is proposed which combines the advantages of fuzzy theory and norm theory. The proposed method can deeply dig system information from limited data, which probability distribution is not taken into account. Firstly, the FNM is employed to evaluate variable interval and expanded uncertainty from limited PSD data, and the performance of FNM is demonstrated by confidence level, reliability and computing accuracy of expanded uncertainty. In addition, the optimal fuzzy parameters are discussed to meet the requirements of aviation standards and metrological practice. Finally, computer simulation is used to prove the adaptability of FNM. Compared with statistical methods, FNM has superiority for evaluating expanded uncertainty from limited data. The results show that the reliability of calculation and evaluation is superior to 95%. 展开更多
关键词 Expanded uncertainty Fuzzy norm method Limited PSD data Random vibration Reliability Variable interval
原文传递
Solving (2+1)-dimensional sine-Poisson equation by a modified variable separated ordinary differential equation method 被引量:1
19
作者 苏卡林 谢元喜 《Chinese Physics B》 SCIE EI CAS CSCD 2010年第10期40-48,共9页
By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ordinary differential equation method is presented for solving the (2 + 1)-dimensional sine-Poisson equa... By introducing a more general auxiliary ordinary differential equation (ODE), a modified variable separated ordinary differential equation method is presented for solving the (2 + 1)-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the (2 + 1)-dimensional sine-Poisson equation are derived in a simple manner by this technique. 展开更多
关键词 modified variable separated ODE method (2 1)-dimensional sine-Poisson equation explicit and exact solution
原文传递
DYNAMIC RESPONSE OF ELASTIC RECTANGULAR PLATES BY SPLINE STATE VARIABLE METHOD
20
作者 陈荣毅 沈小璞 沈鹏程 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2000年第6期691-698,共8页
Application of spline element and state space method for analysis of dynamic response of elastic rectangular plates is presented. The spline element method is used for space domain and the state space method in contro... Application of spline element and state space method for analysis of dynamic response of elastic rectangular plates is presented. The spline element method is used for space domain and the state space method in control theory of system is used for time domain. A state variable recursive scheme is developed, then the dynamic response of structure can he calculated directly. Several numerical examples are given. The results which are presented to demonstrate the accuracy and efficiency of the present method are quite satisfactory. 展开更多
关键词 spline element state space method dynamic response recursive scheme spline state variable method
在线阅读 下载PDF
上一页 1 2 7 下一页 到第
使用帮助 返回顶部