首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 234 毫秒
1.
Summary This note presents an exact solution for the stress and displacement field in an unbounded and transversely constrained elastic medium resulting from the motion of a plane heat source travelling through the medium at constant speed in the direction normal to the source plane.Nomenclature mass density - diffusivity - thermal conductivity - Q heat emitted by plane heat source per unit time per unit area - speed of propagation of plane heat source - shear modulus - Poisson's ratio - T temperature - x, y, z normal stress components - u x, uy, uz displacement components - c speed of irrotational waves - t time - x, y, z Cartesian coordinates - =x–vt moving coordinate  相似文献   

2.
The relative contributions of the viscous, elastic, and texture stresses in flows of modestly concentrated liquid crystalline solutions of poly(-benzyl glutamate) or PBG, at low shear rates are assessed by comparing published experimental rheological and light-scattering data to theoretical predictions. The data considered include ratios of Leslie-Ericksen viscosities, the shear viscosity jump at the transition from the isotropic to the liquid-crystalline state, the dependence of the shear viscosity on concentration, and the fraction of slowly relaxing stress after cessation of shearing flow. We quantify the relative contributions of viscous and elastic stresses by estimating from these data a value for the parameter ß* V str D * r /kT, where str , is the solvent-polymer friction coefficient and D r * is the rotary diffusivity in the isotropic phase at the concentration where the isotropic phase becomes unstable to liquid crystalline order. The experimental evidence suggests a value \* v 0.03, five times lower than the value, 0.15, for dilute solutions. The relative contributions of elastic and texture stresses can be quantified by a parameter in the mesoscopic theory for polydomain nematics; a value of 0.03 gives reasonable agreement with transient shearing experiments, and predicts that the shear viscosity will show an upturn at high polymer concentrations, in agreement with experiments.Dedicated to Prof. Dr. J. Meissner on the occasion of his retirement from the chair of Polymer Physics at the Eidgenössische Technische Hochschule (ETH) Zürich, Switzerland  相似文献   

3.
If a fluid enters an axially rotating pipe, it receives a tangential component of velocity from the moving wall, and the flow pattern change according to the rotational speed. A flow relaminarization is set up by an increase in the rotational speed of the pipe. It will be shown that the tangential- and the axial velocity distribution adopt a quite universal shape in the case of fully developed flow for a fixed value of a new defined rotation parameter. By taking into account the universal character of the velocity profiles, a formula is derived for describing the velocity distribution in an axially rotating pipe. The resulting velocity profiles are compared with measurements of Reich [10] and generally good agreement is found.Nomenclature b constant, equation (34) - D pipe diameter - l mixing length - l 0 mixing length in a non-rotating pipe - N rotation rate,N=Re /Re D - p pressure - R pipe radius - Re D flow-rate Reynolds number, - Re rotational Reynolds number, Re =v w D/ - Re* Reynolds number based on the friction velocity, Re*=v*R/ - (Re*)0 Reynolds number based on the friction velocity in a non-rotating pipe - Ri Richardson number, equation (10) - r coordinate in radial direction - dimensionless coordinate in radial direction, - v r ,v ,v z time mean velocity components - v r ,v ,v z velocity fluctations - v w tangential velocity of the pipe wall - v* friction velocity, - axial mean velocity - v ZM maximum axial velocity - dimensionless radial distance from pipe wall, - y + dimensionless radial distance from pipe wall - y 1 + constant - Z rotation parameter,Z =v w/v * =N Re D /2Re* - m eddy viscosity - ( m )0 eddy viscosity in a non-rotating pipe - coefficient of friction loss - von Karman constant - 1 constant, equation (31) - density - dynamic viscosity - kinematic viscosity  相似文献   

4.
A solution has been found to the problem of calculating the stress and displacement fields caused by a rectilinear dislocation in an anisotropic elastic plate. Special cases of anisotropy have been found with solutions represented by elementary functions.Certain problems in describing crystal plastic deformation phenomena make it vital to know the fields of the elastic stresses and displacements caused by an individual dislocation in a bounded crystal. It is interesting to study the effect of crystal boundaries on these fields with a simple model which approximates fairly closely to experimental conditions.The model selected is shown in Fig, 1. A dislocation with a Burgers vector (b1, b2, b3) is situated in an infinite elastic anisotropic plate of thickness 2h. The dislocation line is parallel to the plate boundaries. The following restriction is introduced in relation to the plate's elastic properties: the medium has a plane of elastic symmetry perpendicular to the dislocation line. The selection of the coordinate system and position of the dislocation are shown in Fig. 1. The requirement is to find the stresses and displacements at an arbitrary point in the plate.One limited special form of this problem has been solved by Kroupa [1]. The limitations which he introduced are as follows: the medium is isotropic, the dislocation is at the precise center of the band and the Burgers vector has only one component b2 differing from zero (the same coordinates were chosen in [1] as in Fig. 1).Thus Kroupa's results can be obtained from the results of the present work as a special case. Other special cases arising from this problem are those concerning the elastic stress and displacement fields caused by a dislocation in anisotropic semi-bounded [2] and bounded [3] media.It is immediately apparent that the problem is a plane one, in the sense that the fields to be found do not depend on coordinate z. Since the medium has a plane of elastic symmetry perpendicular to the dislocation line, it is clear from [4] that the system of stresses and strains in such a medium can be divided into two independent subsystems. The first of these is plane deformation with stress components xx, yy and xy differing from zero and displacement vector components ux and uy, the second is antiplane deformation with stress components xz and yz differing from zero and the displacement vector component uz.In the case under examination, the plane deformation is caused by the Burgers vector edge components bx and by and the antiplane deformation by the screw component bz. The solution is therefore divided into two stages, corresponding to edge and screw dislocations.In conclusion, I wish to thank A. M. Kosevich for his valuable advice and L. A. Pastur for his constant vigilance and assistance in the work.  相似文献   

5.
An infinite plate containing a finite through crack under tensile loading is analysed by Fourier transform based on the Kane-Mindlin kinematic assumptions for the quasi-three-dimensional deformation of plates in extension. The asymptotic expressions of stress and displacement fields near the crack tip, the variation of the stress intensity factor with the plate-thickness and the three-dimensional deformation zone near the crack tip are investigated. The results of the analysis show that, (a) the crack-tip stress and displacement fields accounting for the plate-thickness effects are different from the plane stress solutions and this is true even for extremely small parameter (=1–vh/6 a). In a very small region near the crack tip, plane strain solutions prevail; (b) the ratio of the stress intensity factor KI to the corresponding plane stress one KI, KI/K I o , approaches 1/(1–v2) as tends to zero; (c) plane stress solutions can give satisfactory results for points a distance from the crack tip greater than about three-fourths of the plate-thickness; (d) the linear elastic result for the zone of three-dimensional effects is approximately valid for an elasto-plastic material with linear strain-hardening when the plastic tangential mudulus Et is not very small.The Project Supported by National Natural Science Foundation of China.  相似文献   

6.
In this paper the slip phenomenon is considered as a stochastic process where the polymer segments (taken as Hookean springs) break off the wall due to the excessive tension imposed by the bulk fluid motion. The convection equation arising in network theories is solved for the special case of a polymer/wall interface to determine the time evolution of the configuration distribution function (Q, t). The stress tensor and the slip velocity are calculated by averaging the proper relations over a large number of polymer segments. Due to the fact that the model is probabilistic and time dependent, dynamic slip velocity calculations become possible for the first time and therefore some new insight is gained on the slip phenomenon. Finally, the model predictions are found to match macroscopic experimental data satisfactorily.Nomenclature rate of creation of polymer segments - g(Q) constant of rate of creation of polymer segments - rate of loss of polymer segments - h(Q) constant of rate of loss of polymer segments - h(Q) constant of rate of loss of polymer segments due to destruction of its B-link - H Hookean spring constant - k Boltzmann's constant - n unit vector normal to the polymer/wall interface - n 0 number density of polymer segments - n 0 surface number density of polymer segments - Q vector defining the size and orientation of a polymer segment - Q * critical length of a segment beyond which the tension may overcome the W adh - t time - t h howering time of broken polymer segments - T absolute temperature - W adh work of adhesion Greek Letters n nominal strain - strain - n nominal shear rate - shear rate - dimensionless constant in the expressions of h(Q), g(Q) - viscosity - T velocity gradient tensor - 0 time constant - standard deviation of vectors Q at equilibrium - w wall shear stress - stress tensor - 0 equilibrium configuration distribution function of Q - configuration distribution function of Q  相似文献   

7.
This paper deals with the stress distribution in a homogeneous isotropic elastic hemisphere embedded in a semi-infinite homogeneous isotropic elastic medium when a rigid circular disc is clamped to the plane face of the hemisphere and the stresses are caused by the rotation of the disc through an angle . The problem is reduced to the solution of a Fredholm integral equation of the second kind in the auxiliary function (t). An analytical expression for the torque T required to rotate the die through an angle is obtained in terms of (t). The Fredholm integral equation is solved numerically, and the numerical values of T are graphed.This work has been supported by the National Research Council of Canada through NRC-Research Grant No. A4177.  相似文献   

8.
The rheological properties of glass fibre-filled polypropylene melts have been investigated. A high pressure capillary rheometer has been used for the experimental study. The effect of shear rate, temperature, and fibre concentration on the melt viscosity and viscoelastic properties have been studied. An equation has been proposed to correlate the melt viscosity with shear rate, temperature and fibre content. A master curve relation on this basis has been brought out using the shift factora T . a T shift factor (=/ r ) - A i coefficients of the polynomical of eq. (1) (i = 0, 1, 2, ,n) - B constant in the AFE equation (eq. (2)) (Pa s) - B constant in eq. (3) - D extrudate diameter - d capillary diameter - activation energy at constant shear rate (kcal/mole) - E activation energy at constant shear stress (kcal/mole) - T melt temperature (K) - X fraction glass fibre by weight - shear rate (s–1) - shear viscosity (Pa s) - normal stress coefficient (Pa s2) - 1 2 first normal-stress difference (Pa) - shear stress (Pa) - r at reference temperature  相似文献   

9.
Summary This study deals with the geometrically nonlinear axisymmetric deflection characteristics of orthotropic thick shallow truncated conical springs. Marguerre type governing equations including the effect of shear deformation have been formulated in terms of normal displacement , shear rotation and stress function . Polynmial approximations are used for these field variables and the discretised equations are obtained using the orthogonal point collocation method. The influence of thickness-to-radius ratio, orthotropic parameter and edge conditions has been investigated. Significant influence of shear deformation on the response is found for moderately thick high-modulus composite clamped caps. Its effect is less pronounced for simply-supported caps.
Nichtlineares Verhalten orthotroper, dicker, flacher Tellerfedern
Übersicht In dieser Arbeit werden die nichtlinearen Kennlinien von orthotropen, dicken, flachen Tellerfedern untersucht. Die dem Problem zugrunde liegenden Gleichungen vom Marguerre-Typ schließen den Effekt der Schubdeformation ein und werden mit Hilfe der Auslenkung , des Schubwinkels und einer Spannungsfunktion formuliert. Diese Größen werden in Potenzreihen entwickelt, und mit dem Kollokationsverfahren ergeben sich diskretisierte Gleichungen. Der Einfluß des Verhältnisses von Dicke zu Radius sowie der Einfluß von Orthotropieparametern und Randbedingungen wird untersucht.
  相似文献   

10.
The problem of the interaction of a strongly underexpanded axisymmetric jet with an obstacle for which the normal to the surface makes an angle near /2 with the jet axis is rather laborious for numerical solution due to the high disequilibrium of the gas-dynamic parameters in the peripheral part of the jet and the three-dimensional nature of the flow in the interaction region. Therefore, the results at present available have mainly been obtained experimentally [1, 2]. Among the theoretical studies made in this direction, it is necessary to mention Ivanov and Nazarov's [3], which gives the results of numerical investigation of lateral interaction of a jet with obstacles of various shapes in the case of weakly underexpanded jets when the flow in the interaction region is everywhere supersonic. In the present paper, a study is made of the case when a jet exhausts into vacuum and in front of the obstacle there is a detached shock wave, behind which there is mixed subsonic and supersonic three-dimensional flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 49–54, November–December, 1982.We thank V. I. Uskov for assistance in the present work.  相似文献   

11.
An analytical solution of the E- model for the downstream evolution of a stationary and nearly homogeneous turbulent shear flow is presented. In case that the turbulent time scale has adjusted itself to the time scale imposed by the shear, an asymptotic solution can be derived from the full solution, which shows that both E and increase downstream exponentially. By comparing this asymptotic solution with experimental data a value for the unknown constant c l , in the -equation, is derived. Moreover, we find an expression for the downstream development of the variance of a scalar, which is also compared with experimental data. The analytical solution shows that a homogeneous shear flow with a uniform velocity gradient can only be obtained if the shear is sufficiently small. In the experiments this condition is not always satisfied. A discussion is given on how a nearly homogeneous shear flow can be obtained over a limited downstream interval by changing the initial conditions in E and , and a comparison is made with experimental data. Finally it is shown that better transverse homogeneity can be obtained by taking an exponential velocity profile instead of a linear profile.  相似文献   

12.
The inviscid transonic flow past a symmetric airfoil having a curvature minimum in the middle is numerically investigated. It is shown that at zero angle of attack both symmetric and asymmetric steady-state flow patterns can exist on a certain freestream Mach number range Mmin < Mmax. On this range, the asymmetric flows are stable against small perturbations, whereas the symmetric flows are stable only if M does not coincide with a singular Mach number at which small variations in M or can result in flow restructuring.  相似文献   

13.
The results of investigations of inviscid flow over inverted cones with nose consisting of a spherical segment were published for the first time in Soviet literature in [1–4]. In the present paper, a numerical solution to this problem is obtained using the improved algorithms of [5, 6], which have proved themselves well in problems of exterior flow over surfaces with positive angles of inclination to the oncoming flow. It is shown that the Mach number 2 M , equilibrium and nonequilibrium physicochemical transformations in air (H = 60 km, V = 7.4 km/sec, R0 = 1 m), and the angle of attack 0 40° influence the investigated pressure distributions. A comparison of the results of the calculations with drainage experiments for M = 6, = 0-25° confirms the extended region of applicability of the developed numerical methods. Also proposed is a simple correlation of the dependence on the Mach number in the range 1.5 M of the shape of the shock wave near a sphere in a stream of ideal gas with adiabatic exponent = 1.4.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 178–183, January–February, 1981.  相似文献   

14.
The dependence of the stressintensity factor near the tip of a growing crack in an SO120 acrylic plastic on the crackpropagation velocity KI(:v) within the range of 10–5 –300 m/sec is determined. Specific features of crack propagation associated with the shape of the curve KI(v), which has discontinuities and nonuniqueness intervals, are discussed.  相似文献   

15.
The stress distribution in a homogeneous beam subjected to three-point bending is investigated using the method of optical isodynes. The three stress components xx,yy and xy acting in the planes formed by the longitudinal and vertical axes of the beam are determined in three planes situated at different through the thickness locations with respect to the beam's midplane. The experimental results are subsequently correlated with the two-dimensional elasticity solution. It is illustrated that at locations sufficiently removed from the centrally applied concentrated load, good correlation between theory and experiment is obtained. In the regions where high stress gradients exist however, differences are observed in the in-plane stress distributions in the different planes. These differences are explained by the presence of the out of plane normal stress zz using the relations of optical isodynes. Greatest differences between theory and experiment are obtained for the in-plane shear stress component xy.  相似文献   

16.
Summary Previous work on the creeping flow of viscoelastic fluids past a sphere is reviewed. Theoretical analyses available in the literature were obtained for weakly elastic fluids and therefore they predict only a small influence of fluid elasticity on the drag. In this paper, an approximate theoretical analysis is given for the creeping flow past a rigid sphere in an unbounded medium. The analysis uses a variational principle to solve the equations of motion and continuity in conjunction with the Carreau constitutive equation. The theoretical results are presented in terms of a correction factor to the Newtonian drag coefficient. The correction factor is a function of the power law flow behaviour indexn, the ratio of limiting viscosities ( 0)/0 and a dimensionless time which reflects the elastic nature of the fluids. The results are presented in graphical form covering a realistic range of these dimensionless groups.In order to verify the theoretical predictions, the drag coefficient of a number of spheres was measured in a series of shear thinning elastic test fluids. The flow properties of the test fluids were independently measured with a Weissenberg Rheogoniometer. The power law index of the test fluids varied between 1.0 and 0.4. Particle Reynolds number based on 0 was in the range of 410–6 to 410–2. The difference between theoretically predicted values of drag coefficient and the experimentally measured values is less than ±7.5%. In addition, it is found that the Carreau viscosity equation can be used to predict the elastic parameter of primary normal stress difference with moderate to good accuracy for all the polymer solutions used in this work.
Zusammenfassung Einleitend wird ein Überblick über die früheren Untersuchungen betreffend die schleichende Strömung um eine Kugel gegeben. Die in der Literatur vorliegenden theoretischen Analysen sind auf schwach viskoelastische Flüssigkeiten beschränkt und sagen deshalb nur einen geringen Einfluß der Elastizität auf den Widerstand voraus. In dieser Veröffentlichung wird dagegen eine genäherte theoretische Analyse für die schleichende Strömung um eine starre Kugel in einem unendlich ausgedehnten Medium gegeben, bei welcher zur Lösung der Bewegungsgleichungen und der Kontinuitätsgleichung in Verbindung mit den rheologischen Stoffgleichungen vonCarreau ein Variationsprinzip verwendet wird. Die theoretischen Ergebnisse werden mittels eines Korrekturfaktors zum newtonschen Widerstandskoeffizienten beschrieben. Dieser Korrekturfaktor ist eine Funktion des Potenz-Gesetz-Exponentenn, des Verhältnisses der Grenzviskositäten ( 0)/0 und einer dimensionslosen Zeit, welche das elastische Verhalten kennzeichnet. Die Ergebnisse werden in graphischer Form unter Zugrundelegung eines realistischen Wertebereichs dieser dimensionslosen Gruppen dargestellt.Um diese theoretischen Voraussagen zu verifizieren, wurde der Widerstandskoeffizient für eine Anzahl von Kugeln in einer Reihe von Scherentzähung aufweisenden elastischen Probeflüssigkeiten gemessen. Die Fließeigenschaften dieser Flüssigkeiten wurden zusätzlich mit dem Weissenberg-Rheogoniometer bestimmt. Der Potenz-Gesetz-Exponent variierte dabei zwischen 1,0 und 0,4. Die auf den Kugeldurchmesser und die Nullviskosität bezogenen Reynolds-Zahlen lagen zwischen 410–6 und 410–2. Der Unterschied zwischen theoretisch vorausgesagten und experimentell bestimmten Widerstandskoeffizienten war kleiner als ±7,5%. Außerdem wurde noch gefunden, daß die Viskositätsgleichung vonCarreau dazu verwendet werden kann, den elastischen Parameter erste Normalspannungs-Differenz für alle in dieser Untersuchung verwendeten Polymerlösungen mit mäßiger bis guter Genauigkeit vorauszusagen.

Notation C d drag coefficient - d diameter of sphere - f external body forces in equation of motion [2] - F d drag force - g acceleration due to gravity - J integral defined in eq. [3] - n a parameter in the Carreau viscosity eq. [6] - p isotropic pressure term in equation of motion [2] - r,, spherical coordinates - R radius of sphere - Re 0, Re1 Reynolds numbers defined in eq. [16] - t time - u i ,u j velocities in equation of motion [2] - u r ,u r and components of velocity - V terminal velocity of sphere in unbounded medium - V volume, in eq. [3] - X correction factor to the drag force, eq. [14] - y,z dimensionless spherical coordinates, eq. [9] - ratio of two Reynolds numbers given by eq. [16] - shear rate - apparent viscosity - 0, zero shear rate and infinite shear rate viscosities respectively - a parameter in the Carreau viscosity eq. [6] - the dimensionless time, defined in eq. [11] - second invariant of the rate of deformation tensor - a parameter in the stream function, eq. [8] - stream function - p,f densities of sphere and fluid respectively With 7 figures and 1 table  相似文献   

17.
The results of experiments carried out in order to determine the principal characteristics of the process of displacement of one (nonwetting) fluid from inclined macroinhomogeneous porous media by another (wetting) fluid are presented. Irrespective of whether flow in inclined stratified nonhomogeneous formations or in zonally nonhomogeneous media (with a corresponding well distribution) is investigated, the term oblique stratification is used for describing these processes.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 125–131, November–December, 1992.  相似文献   

18.
Diffusion in anisotropic porous media   总被引:2,自引:0,他引:2  
An experimental system was constructed in order to measure the two distinct components of the effective diffusivity tensor in transversely isotropic, unconsolidated porous media. Measurements were made for porous media consisting of glass spheres, mica particles, and disks made from mylar sheets. Both the particle geometry and the void fraction of the porous media were determined experimentally, and theoretical calculations for the two components of the effective diffusivity tensor were carried out. The comparison between theory and experiment clearly indicates that the void fraction and particle geometry are insufficient to characterize the process of diffusion in anisotropic porous media. Roman Letters A interfacial area between - and -phases for the macroscopic system, m2 - A e area of entrances and exits of the -phase for the macroscopic system, m2 - A interfacial area contained within the averaging volume, m2 - a characteristic length of a particle, m - b average thickness of a particle, m - c A concentration of species A, moles/m3 - c o reference concentration of species A, moles/m3 - c A intrinsic phase average concentration of species A, moles/m3 - c a c Ac A, spatial deviation concentration of species A, moles/m3 - C c A/c 0, dimensionless concentration of species A - binary molecular diffusion coefficient, m2/s - D eff effective diffusivity tensor, m2/s - D xx component of the effective diffusivity tensor associated with diffusion parallel to the bedding plane, m2/s - D yy component of the effective diffusivity tensor associated with diffusion perpendicular to the bedding plane, m2/s - D eff effective diffusivity for isotropic systems, m2/s - f vector field that maps c A on to c a , m - h depth of the mixing chamber, m  相似文献   

19.
By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · xct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · xct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the n -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.  相似文献   

20.
By means of a cone and plate rheometer the relaxation of the shear stress and the first normal stress difference in polymer liquids upon cessation of a constant shear rate were examined. The experiments were conducted mostly in a high shear rate region of relevance for the processing of these materials. The relaxation behavior at these shear rates can only be measured accurately under extremely precise specifications of the rheometer. To determine under which conditions the integral normal thrust is a convenient measure for the relaxing local first normal stress difference the radial distribution of the pressure in the shear gap was measured. The shape of relaxation of both the shear stress and the first normal stress difference could be closely approximated for the entire measured shear rate and time range by a two parameter statistical function. In the range of measured shear rates, one of the parameters, the standard deviationS, is equal for the shear and the normal stress, and is independent of the shear rate within the limit of experimental error. The second parameter, the mean relaxation timet 50, of the shear stress andt 50, of the first normal stress difference, can be calculated approximately from the viscosity function and only a single relaxation experiment.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号