期刊文献+
共找到251篇文章
< 1 2 13 >
每页显示 20 50 100
PERTURBATION METHODS IN SEMILINEAR ELLIPTIC PROBLEMS INVOLVING CRITICAL HARDY-SOBOLEV EXPONENT 被引量:4
1
作者 蓝永艺 唐春雷 《Acta Mathematica Scientia》 SCIE CSCD 2014年第3期703-712,共10页
In this article, we deal with a class of semilinear elliptic equations which are perturbations of the problems with the critical Hardy-Sobolev exponent. Some existence results are given via an abstract perturbation me... In this article, we deal with a class of semilinear elliptic equations which are perturbations of the problems with the critical Hardy-Sobolev exponent. Some existence results are given via an abstract perturbation method in critical point theory. 展开更多
关键词 Critical Hardy-Sobolev exponent semilinear elliptic equation perturbation methods positive radial solution
在线阅读 下载PDF
Numerical Simulation via Homotopy Perturbation Approach of a Dissipative Squeezed Carreau Fluid Flow Due to a Sensor Surface
2
作者 Sara I.Abdelsalam W.Abbas +2 位作者 Ahmed M.Megahed Hassan M.H.Sadek M.S.Emam 《Frontiers in Heat and Mass Transfer》 2025年第5期1511-1527,共17页
This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro... This study rigorously examines the interplay between viscous dissipation,magnetic effects,and thermal radiation on the flow behavior of a non-Newtonian Carreau squeezed fluid passing by a sensor surface within a micro cantilever channel,aiming to deepen our understanding of heat transport processes in complex fluid dynamics scenarios.The primary objective is to elucidate how physical operational parameters influence both the velocity of fluid flow and its temperature distribution,utilizing a comprehensive numerical approach.Employing a combination of mathematical modeling techniques,including similarity transformation,this investigation transforms complex partial differential equations into more manageable ordinary ones,subsequently solving them using the homotopy perturbation method.By analyzing the obtained solutions and presenting them graphically,alongside detailed analysis,the study sheds light on the pivotal role of significant parameters in shaping fluid movement and energy distribution.Noteworthy observations reveal a substantial increase in fluid velocity with escalating magnetic parameters,while conversely,a contrasting trend emerges in the temperature distribution,highlighting the intricate relationship between magnetic effects,flow dynamics,and thermal behavior in non-Newtonian fluids.Further,the suction velocity enhance both the local skin friction and Nusselt numbers,whereas theWeissenberg number reduces them,opposite to the effect of the power-law index. 展开更多
关键词 Homotopy perturbation method squeezed flow Carreau fluid sensor surface magnetic field viscous dissipation
在线阅读 下载PDF
FREE VIBRATION ANALYSIS OF SIMPLY SUPPORTED BEAM WITH BREATHING CRACK USING PERTURBATION METHOD 被引量:20
3
作者 Mousa Rezaee Reza Hassannejad 《Acta Mechanica Solida Sinica》 SCIE EI 2010年第5期459-470,共12页
In this paper, a new analytical method for vibration analysis of a cracked simply supported beam is investigated. By considering a nonlinear model for the fatigue crack, the governing equation of motion of the cracked... In this paper, a new analytical method for vibration analysis of a cracked simply supported beam is investigated. By considering a nonlinear model for the fatigue crack, the governing equation of motion of the cracked beam is solved using perturbation method. The solution of the governing equation reveals the superhaxmonics of the fundamental frequency due to the nonlinear effects in the dynamic response of the cracked beam. Furthermore, considering such a solution, an explicit expression is also derived for the system damping changes due to the changes in the crack parameters, geometric dimensions and mechanical properties of the cracked beam. The results show that an increase in the crack severity and approaching the crack location to the middle of the beam increase the system damping. In order to validate the results, changes in the fundamental frequency ratios against the fatigue crack severities are compared with those of experimental results available in the literature. Also, a comparison is made between the free response of the cracked beam with a given crack depth and location obtained by the proposed analytical solution and that of the numerical method. The results of the proposed method agree with the experimental and numerical results. 展开更多
关键词 free vibration cracked beam breathing crack damping factor perturbation method
原文传递
Method of Semi Analytic Perturbation Weighted Residuals for Nonlinear Bending of Shallow Shells
4
作者 李云飞 黄怡筠 《Journal of Beijing Institute of Technology》 EI CAS 2000年第1期15-19,共5页
The semi? analytic perturbation weighted residuals method was used to solve the nonlinear bending problem of shallow shells, and the fifth order B spline was taken as trial function to seek an efficient method for n... The semi? analytic perturbation weighted residuals method was used to solve the nonlinear bending problem of shallow shells, and the fifth order B spline was taken as trial function to seek an efficient method for nonlinear bending problem of shallow shells. The results from the present method are in good agreement with those derived from other methods. The present method is of higher accuracy, lower computing time and wider adaptability. In addition, the design of computer program is simple and it is easy to be programmed. 展开更多
关键词 shallow shell nonlinear bending perturbation method weighted residuals method spline function
在线阅读 下载PDF
Perturbation method of studying the EI Nifio oscillation with two parameters by using the delay sea-air oscillator model 被引量:4
5
作者 杜增吉 林万涛 莫嘉琪 《Chinese Physics B》 SCIE EI CAS CSCD 2012年第9期32-36,共5页
The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation usin... The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation using the ENSO model. Based on a class of the oscillator of the ENSO model, a approximate solution of the corresponding problem is studied employing the perturbation method. 展开更多
关键词 NONLINEAR perturbation method E1 Nino-southern oscillation model
原文传递
Asymptotical solutions of coupled nonlinear Schrodinger equations with perturbations 被引量:2
6
作者 程雪苹 林机 叶丽军 《Chinese Physics B》 SCIE EI CAS CSCD 2007年第9期2503-2509,共7页
In this paper Lou's direct perturbation method is applied to the perturbed coupled nonlinear Schrodinger equations to obtain their asymptotical solutions, which include not only the zero-order solutions but also the ... In this paper Lou's direct perturbation method is applied to the perturbed coupled nonlinear Schrodinger equations to obtain their asymptotical solutions, which include not only the zero-order solutions but also the first-order modifications. Based on the asymptotical solutions, the effects of perturbations on soliton parameters and the collision between two solitons are then discussed in brief. Furthermore, we directly simulate the perturbed coupled nonlinear SchrSdinger equations by split-step Fourier method to check the validity of the direct perturbation method. It turns out that our analytical results are well supported by the numerical calculations. 展开更多
关键词 direct perturbation method perturbed coupled nonlinear Schrodinger equations soli- tons asymptotical solutions
原文传递
HOMOTOPY PERTURBATION SOLUTION AND PERIODICITY ANALYSIS OF NONLINEAR VIBRATION OF THIN RECTANGULAR FUNCTIONALLY GRADED PLATES 被引量:2
7
作者 A.Allahverdizadeh R.Oftadeh +1 位作者 M.J.Mahjoob M.H.Naei 《Acta Mechanica Solida Sinica》 SCIE EI CSCD 2014年第2期210-220,共11页
In this paper nonlinear analysis of a thin rectangular functionally graded piate is formulated in terms of von-Karman's dynamic equations. Functionaily Graded Material (FGM) properties vary through the constant thi... In this paper nonlinear analysis of a thin rectangular functionally graded piate is formulated in terms of von-Karman's dynamic equations. Functionaily Graded Material (FGM) properties vary through the constant thickness of the plate at ambient temperature. By expansion of the solution as a series of mode functions, we reduce the governing equations of motion to a Duffing's equation. The homotopy perturbation solution of generated Duffing's equation is also obtained and compared with numerical solutions. The sufficient conditions for the existence of periodic oscillatory behavior of the plate are established by using Green's function and Schauder's fixed point theorem. 展开更多
关键词 nonlinear vibration FGM rectangular plate Schauder's fixed point theorem homotopy perturbation method
原文传递
Ground-state energy of beryllium atom with parameter perturbation method 被引量:2
8
作者 Feng Wu LijuanMeng 《Chinese Physics B》 SCIE EI CAS CSCD 2018年第9期237-240,共4页
We present a perturbation study of the ground-state energy of the beryllium atom by incorporating double parameters in the atom's Hamiltonian. The eigenvalue of the Hamiltonian is then solved with a double-fold pertu... We present a perturbation study of the ground-state energy of the beryllium atom by incorporating double parameters in the atom's Hamiltonian. The eigenvalue of the Hamiltonian is then solved with a double-fold perturbation scheme,where the spin-spin interaction of electrons from different shells of the atom is also considered. Calculations show that the obtained ground-state energy is in satisfactory agreement with experiment. It is found that the Coulomb repulsion of the inner-shell electrons enhances the effective nuclear charge seen by the outer-shell electrons, and the shielding effect of the outer-shell electrons to the nucleus is also notable compared with that of the inner-shell electrons. 展开更多
关键词 parameter perturbation method double-fold perturbation scheme ground-state energy
原文传递
Thermoelastic Analysis of Non-uniform Pressurized Functionally Graded Cylinder with Variable Thickness Using First Order Shear Deformation Theory(FSDT) and Perturbation Method 被引量:1
9
作者 KHOSHGOFTAR M J MIRZAALI M J RAHIMI G H 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 2015年第6期1149-1156,共8页
Recently application of functionally graded materials(FGMs) have attracted a great deal of interest. These materials are composed of various materials with different micro-structures which can vary spatially in FGMs... Recently application of functionally graded materials(FGMs) have attracted a great deal of interest. These materials are composed of various materials with different micro-structures which can vary spatially in FGMs. Such composites with varying thickness and non-uniform pressure can be used in the aerospace engineering. Therefore, analysis of such composite is of high importance in engineering problems. Thermoelastic analysis of functionally graded cylinder with variable thickness under non-uniform pressure is considered. First order shear deformation theory and total potential energy approach is applied to obtain the governing equations of non-homogeneous cylinder. Considering the inner and outer solutions, perturbation series are applied to solve the governing equations. Outer solution for out of boundaries and more sensitive variable in inner solution at the boundaries are considered. Combining of inner and outer solution for near and far points from boundaries leads to high accurate displacement field distribution. The main aim of this paper is to show the capability of matched asymptotic solution for different non-homogeneous cylinders with different shapes and different non-uniform pressures. The results can be used to design the optimum thickness of the cylinder and also some properties such as high temperature residence by applying non-homogeneous material. 展开更多
关键词 non-homogenous cylinder First order Shear Deformation Theory matched asymptotic method perturbation method functionally graded material
在线阅读 下载PDF
Regular perturbation solution of Couette flow(non-Newtonian)between two parallel porous plates: a numerical analysis with irreversibility 被引量:1
10
作者 M.NAZEER M.I.KHAN +5 位作者 S.KADRY Yuming CHU F.AHMAD W.ALI M.IRFAN M.SHAHEEN 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI CSCD 2021年第1期127-142,共16页
The unavailability of wasted energy due to the irreversibility in the process is called the entropy generation.An irreversible process is a process in which the entropy of the system is increased.The second law of the... The unavailability of wasted energy due to the irreversibility in the process is called the entropy generation.An irreversible process is a process in which the entropy of the system is increased.The second law of thermodynamics is used to define whether the given system is reversible or irreversible.Here,our focus is how to reduce the entropy of the system and maximize the capability of the system.There are many methods for maximizing the capacity of heat transport.The constant pressure gradient or motion of the wall can be used to increase the heat transfer rate and minimize the entropy.The objective of this study is to analyze the heat and mass transfer of an Eyring-Powell fluid in a porous channel.For this,we choose two different fluid models,namely,the plane and generalized Couette flows.The flow is generated in the channel due to a pressure gradient or with the moving of the upper lid.The present analysis shows the effects of the fluid parameters on the velocity,the temperature,the entropy generation,and the Bejan number.The nonlinear boundary value problem of the flow problem is solved with the help of the regular perturbation method.To validate the perturbation solution,a numerical solution is also obtained with the help of the built-in command NDSolve of MATHEMATICA 11.0.The velocity profile shows the shear thickening behavior via first-order Eyring-Powell parameters.It is also observed that the profile of the Bejan number has a decreasing trend against the Brinkman number.Whenηi→0(i=1,2,3),the Eyring-Powell fluid is transformed into a Newtonian fluid. 展开更多
关键词 Couette flow Eyring-Powell fluid entropy generation perturbation method Bejan number
在线阅读 下载PDF
INTERPOLATION PERTURBATION METHOD FORSOLVING NONLINEAR PROBLEMS 被引量:1
11
作者 袁镒吾 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 1997年第11期0-0,0-0+0-0+0-0+0,共9页
In this paper, using the interpolation perturbation method. the author seeks tosolve several nonlinear problems. Numerical examples show that the method Df thispaper has good accuracy.
关键词 INTERPOLATION singular perturbation method NONLINEAR
在线阅读 下载PDF
PERTURBATIONAL FINITE VOLUME METHOD FOR THE SOLUTION OF 2-D NAVIER-STOKES EQUATIONS ON UNSTRUCTURED AND STRUCTURED COLOCATED MESHES 被引量:1
12
作者 高智 代民果 +1 位作者 李桂波 柏威 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2005年第2期242-251,共10页
Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were de... Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4. 展开更多
关键词 colocated grid structured grid unstructured grid perturbation finite volume method incompressible fluid NS equations SIMPLEC algorithm MSIMPLEC algorithm SIMPLER algorithm
在线阅读 下载PDF
He’s Homotopy Perturbation Method and Fractional Complex Transform for Analysis Time Fractional Fornberg-Whitham Equation 被引量:1
13
作者 Yanni Zhang Jing Pang 《Sound & Vibration》 EI 2021年第4期295-303,共9页
In this article,time fractional Fornberg-Whitham equation of He’s fractional derivative is studied.To transform the fractional model into its equivalent differential equation,the fractional complex transform is used ... In this article,time fractional Fornberg-Whitham equation of He’s fractional derivative is studied.To transform the fractional model into its equivalent differential equation,the fractional complex transform is used and He’s homotopy perturbation method is implemented to get the approximate analytical solutions of the fractional-order problems.The graphs are plotted to analysis the fractional-order mathematical modeling. 展开更多
关键词 Time fractional Fornberg-Whitham equation fractional complex transform He’s homotopy perturbation method
在线阅读 下载PDF
A New Approximation to the Linear Matrix Equation AX = B by Modification of He’s Homotopy Perturbation Method 被引量:1
14
作者 Amir Sadeghi 《Advances in Linear Algebra & Matrix Theory》 2016年第2期23-30,共8页
It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obt... It is well known that the matrix equations play a significant role in engineering and applicable sciences. In this research article, a new modification of the homotopy perturbation method (HPM) will be proposed to obtain the approximated solution of the matrix equation in the form AX = B. Moreover, the conditions are deduced to check the convergence of the homotopy series. Numerical implementations are adapted to illustrate the properties of the modified method. 展开更多
关键词 Matrix Equation Homotopy perturbation Method CONVERGENCE Diagonally Dominant Matrix Regular Splitting
在线阅读 下载PDF
PERTURBATION SOLUTION TO 3-D NONLINEAR SUPERCAVITATING FLOW PROBLEMS
15
作者 叶取源 何友声 朱世权 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 1993年第1期13-21,共9页
A 3-D nonlinear problem of supercavitating flow past an axisymmetric body at a small angle of attack is investigated by means of the perturbation method and Fourier-cosine-expansion method. The first three order pertu... A 3-D nonlinear problem of supercavitating flow past an axisymmetric body at a small angle of attack is investigated by means of the perturbation method and Fourier-cosine-expansion method. The first three order perturbation equations are derived in detail and solved numerically using the boundary integral equation method and iterative techniques. Computational results of the hydrodynamic characteristics and cavity shapes of each order are presented for nonaxisymmetric supercavitating flow past cones with various apex-angles at differ- ent cavitation numbers. The numerical results are found in good agreement with experimental data. 展开更多
关键词 supercavitating flow perturbation method nonlinear solution boundary integral equation method Fourier expansion elliptic integrals
在线阅读 下载PDF
High order symplectic conservative perturbation method for time-varying Hamiltonian system 被引量:1
16
作者 Ming-Hui Fu Ke-Lang Lu Lin-Hua Lan 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2012年第3期885-890,共6页
This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order... This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly. 展开更多
关键词 Time-varying Hamiltonian system High ordermultiplicative perturbation Symplectic conservation expo-nential matrix Precise time integration method
在线阅读 下载PDF
Perturbation solutions for asymmetric laminar flow in porous channel with expanding and contracting walls
17
作者 张燕 林平 司新辉 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2014年第2期203-220,共18页
The cases of large Reynolds number and small expansion ratio for the asym- metric laminar flow through a two-dimensional porous expanding channel are considered. The Navier-Stokes equations are reduced to a nonlinear ... The cases of large Reynolds number and small expansion ratio for the asym- metric laminar flow through a two-dimensional porous expanding channel are considered. The Navier-Stokes equations are reduced to a nonlinear fourth-order ordinary differential equation by introducing a time and space similar transformation. A singular perturbation method is used for the large suction Reynolds case to obtain an asymptotic solution by matching outer and inner solutions. For the case of small expansion ratios, we are able to obtain asymptotic solutions by double parameter expansion in either a small Reynolds number or a small asymmetric parameter. The asymptotic solutions indicate that the Reynolds number and expansion ratio play an important role in the flow behavior. Nu- merical methods are also designed to confirm the correctness of the present asymptotic solutions. 展开更多
关键词 singular perturbation method regular perturbation method porousexpanding channel expansion ratio
在线阅读 下载PDF
A Significant Look at the Effects of Persian Gulf Environmental Conditions on Sound Scattering Based on Small Perturbation Method
18
作者 Parviz Ghadimi Alireza Bolghasi +1 位作者 Mohammad A. Feizi Chekab Rahim Zamanian 《Journal of Marine Science and Application》 CSCD 2015年第4期413-424,共12页
The main goal of this paper is to investigate sound scattering from the sea surface, by Kuo's small perturbation method (SPM), in the Persian Gulf's environmental conditions. Accordingly the SPM method is reviewed... The main goal of this paper is to investigate sound scattering from the sea surface, by Kuo's small perturbation method (SPM), in the Persian Gulf's environmental conditions. Accordingly the SPM method is reviewed, then it is demonstrated how it can accurately model sound scattering from the sea surface. Since in Kuo's approach, the effects of surface roughness and sub-surface bubble plumes on incident sounds can be studied separately, it is possible to investigate the importance of each mechanism in various scattering regimes. To conduct this study, wind and wave information reported by Arzanah station as well as some numerical atmospheric models for the Persian Gulf are presented and applied to examine sound scattering from the sea surface in the Persian Gulf region. Plots of scattering strength by Kuo's SPM method versus grazing angle for various frequencies, wave heights, and wind speeds are presented. The calculated scattering strength by the SPM method for various frequencies and wind speeds are compared against the results of critical sea tests 7 (CST-7). The favorable agreement achieved for sound scattering in the Persian Gulf region is indicative of the fact that the SPM method can quite accurately model and predict sound scattering from the sea surface. 展开更多
关键词 sea surface Persian Gulf small perturbation method (SPM) wind rose plots wave rose plots sound scattering surface roughness sub-surface bubble plumes
在线阅读 下载PDF
A NEW MATRIX PERTURBATION METHOD FOR ANALYTICAL SOLUTION OF THE COMPLEX MODAL EIGENVALUE PROBLEM OF VISCOUSLY DAMPED LINEAR VIBRATION SYSTEMS
19
作者 吕振华 冯振东 方传流 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 1991年第8期767-776,共10页
A new matrix perturbation analysis method is presented for efficient approximate solution of the complex modal quadratic generalized eigenvalue problem of viscously damped linear vibration systems. First, the damping ... A new matrix perturbation analysis method is presented for efficient approximate solution of the complex modal quadratic generalized eigenvalue problem of viscously damped linear vibration systems. First, the damping matrix is decomposed into the sum of a proportional-and a nonproportional-damping parts, and the solutions of the real modal eigenproblem with the proportional dampings are determined, which are a set of initial approximate solutions of the complex modal eigenproblem. Second, by taking the nonproportional-damping part as a small modification to the proportional one and using the matrix perturbation analysis method, a set of approximate solutions of the complex modal eigenvalue problem can be obtained analytically. The result is quite simple. The new method is applicable to the systems with viscous dampings-which do not deviate far away from the proportional-damping case. It is particularly important that the solution technique be also effective to the systems with heavy, but not over, dampings. The solution formulas of complex modal eigenvlaues and eigenvectors are derived up to second-order perturbation terms. The effectiveness of the perturbation algorithm is illustrated by an exemplar numerical problem with heavy dampings. In addition, the practicability of approximately estimating the complex modal eigenvalues, under the proportional-damping hypothesis, of damped vibration systems is discussed by several numerical examples. 展开更多
关键词 vibration analysis damped system complex mode quadratic eigenvalue problem matrix perturbation method
在线阅读 下载PDF
A NOTE ON DELTA-PERTURBATION EXPANSION METHOD
20
作者 HE Ji-huan(何吉欢) 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2002年第6期634-638,共5页
The Delta-perturbation expansion method, a kind of new perturbation technique depending upon an artificial parameter Delta was studied. The study reveals that the method exits some advantages, but also exits some limi... The Delta-perturbation expansion method, a kind of new perturbation technique depending upon an artificial parameter Delta was studied. The study reveals that the method exits some advantages, but also exits some limitations. To overcome the limitations, the so-called linearized perturbation method proposed by HE Ji-huan can be powerfully applied. 展开更多
关键词 perturbation method artificial parameter nonlinear equation homotropy
在线阅读 下载PDF
上一页 1 2 13 下一页 到第
使用帮助 返回顶部