首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G N . Two time constants, s and 1, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios 1/ s and k / s were determined solely by N for each polymer species. The ratio 1/ s was proportional to N 4.5 and N 3.5 for PMMA and PS solutions, respectively. The ratio k / s was approximately proportional to N 2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants s and G N are known.  相似文献   

2.
Die swell of filled polymer melts   总被引:1,自引:0,他引:1  
The Barus effect in polypropylene and polystyrene blended with a variety of fillers at various concentrations was investigated using a capillary extrusion rheometer. If the die swell is defined as the square of the ratio of the extrudate diameterd to the die diameterD, it is found to depend on the apparent shear stress W . Below a certain value of w the relation =B B A applies. The die swell, M , of a filled polymer depends on the type, size and volume fraction of the filler. In particular,A increases as the volume fraction increases and is largest for powders, smaller for flakes and smallest for fibres, whereasB shows the opposite trend but to a lesser extent.  相似文献   

3.
Transients in melt spinning of isothermal power law and Newtonian fluids were found to be governed by an extremely simple partial differential equation 2 ( 1/n )/() = 0 in Lagrangian coordinates where is the cross-sectional area,n the power law exponent, the time and the the time at which a fluid molecule constituting the spinline left the spinneret. The general integral 1/n =f() +g () of the above governing equation containing two arbitrary functions represents physically attainable spinline transients. Hitherto unknown analytical transient solutions of the above governing equation were obtained for the response of isothermal constant tension spinlines to a stepwise change in tension, spinneret hole area, extrusion speed or extrusion viscosity and for the starting transient in gravitational spinning. Linearized perturbation solutions and the stability limit of the spinline derived from the above new found nonlinear solutions were in agreement with previous findings and the above nonlinear response of the spinline to a step increase in the spinneret hole area was found to be equivalent to Orowan's tandem cylinder model of dent growth in filament stretching.  相似文献   

4.
The pseudoplastic flow of suspensions, alumina or styrene-acrylamide copolymer particles in water or an aqueous solution of glycerin has been studied by the step-shear-rate method. The relation between the shear rate,D, and the shear stress,, in the step-shear-rate measurements, where the state of dispersion was considered to be constant, was expressed as = AD 1/2 +CD. The effective solid volume fraction,ø F, andA were dependent on the shear rate and expressed byø F =aD b andA = D . Combining the above relations, the steady flow curve was expressed by = D 1/2 + + 0 (1 – a D b/0.74)–1.85 D, where 0 is the viscosity of the medium.With an increase in solid volume fraction and a decreases in the absolute value of the-potential, the flow behavior of the suspensions changed from Newtonian ( = = b = 0), slightly pseudoplastic ( = b = 0), pseudoplastic ( = 0) to a Bingham-like behavior.The change in viscosity of the medium had an effect on the change in the effective volume fraction.  相似文献   

5.
F. Durst  R. Haas 《Rheologica Acta》1981,20(2):179-192
Zusammenfassung Es werden theoretische Überlegungen zusammenfassend dargestellt, welche die Streckung und Ausrichtung von flexiblen Makromolekülen in stationären einfachen Dehnströmungen beschreiben. Die Makromoleküle werden hierbei als EDNE-(endlich dehnbare, nichtlinear elastische) Hanteln modelliert. Für den Fall niedriger bzw. hoher Dehnungsraten werden Dehnviskositätsgleichungen für Strömungen mit verdünnten Polymerlösungen angegeben.Die Arbeit vergleicht die abgeleiteten theoretischen Gleichungen mit experimentellen Ergebnissen, welche für Porenströmungen erhalten wurden; Porenströmungen weisen Dehnströmungen auf. Anhand der durchgeführten experimentellen Untersuchungen, in denen alle die den Druckverlust maßgebend beeinflussenden strömungsmechanischen und physikalisch-chemischen Parameter variiert wurden, kann gezeigt werden, daß sich die aufgezeigten theoretischen Zusammenhänge quantitativ bestätigen lassen.Schlüsselwörter Dehnströmung, Makromolekülmodell, Porenströmung, EDNE-Hantelmodell, Polymerlösung
Summary The present paper summarizes theoretical considerations regarding the elongation of flexible macromolecules in simple steady elongational flows. The macromolecules are treated as FENE(finite extensible, nonlinear elastic)-dumbbells. Equations for extensional viscosity are given for flows of dilute polymer solutions applicable at low and high elongation rates.The present paper compares the derived theoretical relationships with experimental results. These results were obtained in porous media flows, which exhibit strong elongational rates. It can be shown on the basis of the experimental investigations, that all fluid mechanic and physico-chemical parameters that influence the measured pressure losses responded as predicted by the theory.

a Mark-Houwink-Exponent - A Avogadro-Konstante - b Verhältnis von Molekülzeitkonstanten - c Polymergewichtskonzentration - d Kugeldurchmesser der Schüttung - D Diffusionskonstante - De Deborahzahl - f Reibungsbeiwert der Porenströmung - F Kraftvektor des Hantelmodells - g Erdbeschleunigung - H Hookesche Federkonstante des Makromoleküls - k Boltzmann-Konstante - k 1,2,3 empirische Konstanten - K Mark-Houwink-Konstante - l 0 Länge des Monomeren - L Länge des statistischen Fadenelementes - L 0 Maximallänge des gestreckten Polymermoleküls - L Bezugslänge für den Druckverlust der Porenströmung - m Masse des statistischen Fadenelementes - m 0 Masse des Monomeren - Molarität - M Molekulargewicht des Polymeren - n Porosität der Kugelschüttung - n 0 Hantelkonzentration - N Anzahl der statistischen Fadenelemente - p Druckverlust der Porenströmung - P Polymerisationsgrad - R Endpunktabstand des Makromoleküls - R 0 maximaler Endpunktabstand des gestreckten Moleküls - mittlerer Endpunktabstand des Moleküls - Orientierungsvektor des Hantelmodells - Re Reynoldszahl der Porenströmung - t Zeit - T Temperatur - mittlere Filtergeschwindigkeit der Porenströmung - v Strömungsfeld - Aufweitungsparameter - Bindungswinkel zweier Kohlenstoffatome - Dehnungsrate - Stokesscher Reibungsfaktor - dynamische Viskosität - * reduzierte Viskosität - [] Grenzviskositätszahl - Dehnviskosität - * reduzierte Dehnviskosität - Widerstandskennzahl der Porenströmung - v kinematische Viskosität - Dichte des Fluids - H Hookesche Relaxationszeit des EDNE-Hantelmodells - H,e Hookesche Relaxationszeit des linear elastischen Hantelmodells - R Relaxationszeit des starren Hantelmodells - zz , yy Normalspannungen - Volumenkonzentration - fl. dimensionsloser Faktor des Strömungsfeldes - 0 Konstante der Flory-Fox-Gleichung - Verteilungsfunktion des Hantelmodells - eq. Gleichgewichtsverteilungsfunktion - a aufgeweitet - e effektiv - max maximal - p polymer - s solvent, Lösungsmittel - Theta-Zustand Mit 12 Abbildungen und 2 Tabellen  相似文献   

6.
Bifurcation phenomena from standing pulse solutions of the problem is considered. (>0) is a sufficiently small parameter and is a positive one. It is shown that there exist two types of destabilization of standing pulse solutions when decreases. One is the appearance of travelling pulse solutions via the static bifurcation and the other is that of in-phase breathers via the Hopf bifurcation. Furthermore which type of destabilization occurs first with decreasing is discussed for the piecewise linear nonlinearities f and g.  相似文献   

7.
An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of *.Here the mixture component temperature difference relaxation time T /, where is the relaxation time for the light component. If 1, 1, so that T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.  相似文献   

8.
Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U 2, dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u 1 velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - w wall shear stress  相似文献   

9.
Zusammenfassung Der Wärmeübergang bei turbulenter Film kondensation strömenden Dampfes an einer waagerechten ebenen Platte wurde mit Hilfe der Analogie zwischen Impuls-und Wärmeaustausch untersucht. Zur Beschreibung des Impulsaustausches im Film wurde ein Vierbereichmodell vorgestellt. Nach diesem Modell wird die wellige Phasengrenze als starre rauhe Wand angesehen. Die Abhängigkeit einer Schubspannungs-Nusseltzahl von der Film-Reynoldszahl und Prandtlzahl wurde berechnet und dargestellt.
A model for turbulent film condensation of flowing vapour
The heat transfer in turbulent film condensation of flowing vapour on a horizontal flat plate was investigated by means of the analogy between momentum and heat transfer. To describe the momentum transfer in the film a four-region model was presented. With this model the wavy interfacial surface is treated as a stiff rough wall. A shear Nusselt number has been calculated and represented as a function of film Reynolds number and Prandtl number.

Formelzeichen a Temperaturleitkoeffizient - k Mischungswegkonstante - k s äquivalente Sandkornrauhigkeit - Nu x lokale Schubspannungs-Nusseltzahl,Nu x=xxv/uw - Pr Prandtlzahl,Pr=v/a - Pr t turbulente Prandtlzahl,Pr t =m/q - q Wärmestromdichte q - R Wärmeübergangswiderstand - Rf Wärmeübergangswiderstand des Films - Re F Reynoldszahl der Filmströmung - T Temperatur - U, V Geschwindigkeitskomponenten des Dampfes in waagerechter und senkrechter Richtung - u, Geschwindigkeitskomponenten des Kondensats in waagerechter und senkrechter Richtung - V Querschwankungsgeschwindigkeit des Kondensats und des Dampfes - u /gtD Schubspannungsgeschwindigkeit an der Phasengrenze für die Dampfgrenzschicht, uD =(/)1/2 - u F Schubspannungsgeschwindigkeit an der Phasengrenze für den Kondensatfilm,u F =(/)1/2 - u w Schubspannungsgeschwindigkeit an der Wand der Kühlplatte,u w =(w/)1/2 - y Wandabstand - x Wärmeübergangskoeffizient - gemittelte Kondensatfilmdicke - s Dicke der zähen Schicht der Filmströmung an der welligen Phasengrenze - 4 Dicke der zähen Schicht der Filmströmung an der gemittelten glatten Phasengrenze - Wärmeleitzahl - dynamische Viskosität - v kinematische Viskosität - Dichte - Oberflächenspannung - w Wandschubspannung - Schubspannung an der Phasengrenzfläche - m turbulente Impulsaustauschgröße - q turbulente Wärmeaustauschgröße Indizes d Wert des Dampfes - w Wert an der Wand - x lokaler Wert inx - Wert an der Phasengrenze Stoffgrößen ohne Index gelten für das Kondensat  相似文献   

10.
In this paper, the compactness of quasi-conforming element spaces and the—convergence of quasi-conforming element method are discussed. The well-known Rellich compactness theorem is generalized to the sequences of quasi-conforming element spaces with certain properties, and the generalized Poincare inequality. The generalized Friedrichs inequality and the generalized inequality of Poincare-Friedrichs are proved true for them. The error estimates are also given. It is shown that the quasi-conforming element method is convergent if the quasi-conforming element spaces have the approximability and the strong continuity, and satisfy the rank condition of element and pass the test IPT. As practical examples, 6-parameter, 9-paramenter, 12-paramenter, 15-parameter, 18-parameter and 21-paramenter quasi-conforming elements are shown to be convergent, and their L22()-errors are O(h), O(h), O(h 2 ), O(h 2 ), O(h ), and O(h 4 ) respectively.  相似文献   

11.
A technique is described which employs automated image processing of hydrogen-bubble flow visualization pictures to establish local, instantaneous velocity profile information. Hydrogen bubble flow visualization sequences are recorded using a high-speed video system and then digitized, stored, and evaluated by a VAX 11/780 computer. Employing special smoothing and gradient detection algorithms, individual bubble-lines are computer identified, which allows local velocity profiles to be constructed using time-of-flight techniques. It is demonstrated how this techniques may be used to 1) determine local velocity behavior as a function of position and time, 2) evaluate time-averaged turbulence properties, and 3) correlate probe-type turbulent burst detection techniques with the corresponding visualization data.List of symbols Re Reynolds number based on momentum thickness, u / - t + nondimensional time tu 2 / - T VITA variance averaging time period - u shear velocity = - u local instantaneous streamwise velocity,x-direction - u local fluctuating streamwise velocity,x-direction - u + nondimensional streamwise velocity, /u - local normal velocity,y-direction - w local spanwise velocity,z-direction - x + nondimensional coordinate in streamwise direction xu /v - y + nondimensional coordinate normal to wall, yu /v Greek momentum thickness, - kinematic viscosity - w wall shear stress This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984  相似文献   

12.
A new procedure for the reduction of Preston tube data is introduced, based on the van Driest transformation. It appears to give results agreeing with the better calibration experiments, although a significant assumption in its derivation is violated.List of Symbols M s Mach number sensed by Preston tube - M Friction Mach number (=u/wall sound speed) - R Gas constant - T w Wall temperature - d Diameter of Preston tube - h Height of effective centre of Preston tube - p Preston tube pressure difference reading - p i Equivalent incompressible Preston tube reading - p w Wall pressure - r Recovery factor (=0.896) - u Friction velocity (=[w/wall density]1/2) - Empirical constant allowing for departure from Crocco temperature-velocity correlation (=0.975) - Specific heat ratio - Fluid kinematic viscosity - w Wall shear-stress  相似文献   

13.
Summary Similarity laws for the mean flow and scaling laws for the turbulent motion are used in an attempt to obtain a general expression for the eddy viscosity of equilibrium layers. It is found that =0.09 w 2 /w*, in which w 2 is a Reynolds stress representative for the region of overlap between the law of the wall and the velocity-defect law, while w* is the logarithmic slope of the mean velocity profile in that region. The distinction between w and w* is related to the strong inhomogeneity of the mean rate of strain in the inner layer. The results of the theory agree with experimental evidence obtained from transpired equilibrium layers.  相似文献   

14.
Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< –1; = linear viscoelastic relaxation time) and shear-thickening for larger (>–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and (t, ). In the shear-thinning regime, the and dependence of +(t,) and (t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson  相似文献   

15.
The behavior of the neutral stability curves is investigated for various values of the particle relaxation time and mass concentration 0 100 and 0 f 0.1. It is shown that as increases from zero the flow is at first destabilized and then at >6 becomes stable, while at >40 the stabilizing effect of the dispersed phase grows weaker. It is found that there is a certain interval 10< <40 on which the flow is most stable.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 46–53, January–February, 1986.  相似文献   

16.
We establish the saddle-point property of the system of functional differential equations (t) = Ax(t) + Bx((t)) + C ((t)) + f (x(t), x((t))), (0) = 0.Translated from Neliniini Kolyvannya, Vol. 7, No. 3, pp. 302–310, July–September, 2004.  相似文献   

17.
An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/y - c = 2µ - d = 0 2–2cy - f() = (–0)2+2c(–y) - r Radius - r b Bob radius - r c Cup radius - r y Yield radius - r 0 Stationary surface radius - r Rotating Stationary radius - Y 0 Shear rate parameter = /µ Greek letters Shear rate - = (r y /r b )2– 1 - µ Bingham viscosity - µ0 Initial differential viscosity - µ µ0 - Rotation speed - Angular velocity - Shear stress - b Bob shear stress - B Bingham stress - y (True) yield stress - 0 Stress parameter = B Y 0 - B - y   相似文献   

18.
McAdams  J. E.  Williams  M. C. 《Rheologica Acta》1986,25(2):102-109
Theta solvents for polystyrene are prepared from high-viscosity blends of styrene and low-molecular-weight polystyrene, and then used to make dilute solutions with monodisperse polystyrene solutes of high-M = 2.3, 6.0, 9.0, 18.0 · 105. A Weissenberg rheogoniometer is used to measure the non-Newtonian viscosity as a function of shear stress, for low values, and also the complex viscosity components and as functions of frequency. A capillary viscometer is used for high- measurements of(). Viscometric properties, at room temperature, are analyzed as functions of high-molecular-weight solute concentrationc with parameters of constant or to obtain [()], [ ()], and [ ()]. Such a collection of data has apparently not previously been available for polymers in theta solvents (in which Gaussian chain statistics prevail). Also unique is the achievement of high stress ( = 2 104 Pa) at low shear rate, by virtue of high solvent viscosity which is not characteristic of other known theta solvents.  相似文献   

19.
Correlations for corrections to hot-wire data for the effects of wall proximity within the viscous sublayer are usually presented in the form u/u = F (y u /). The application of such correlations requires a prior knowledge of the wall shear stress; alternatively, the correlation must be used in an iterative fashion. It is shown in the present note that any such correlation may be recast with no loss of generality in the explicit form u/u m = f (y u m/), which is more convenient for use.List of symbols u difference between measured and actual velocities, u mu - u m measured velocity - u shear velocity, - u + on-dimensional velocity, u/u - y distance from wall - y + non-dimensional distance from wall, y u / - fluid density - fluid kinematic viscosity - s wall shear stress  相似文献   

20.
The effects of finite measuring volume length on laser velocimetry measurements of turbulent boundary layers were studied. Four different effective measuring volume lengths, ranging in spanwise extent from 7 to 44 viscous units, were used in a low Reynolds number (Re=1440) turbulent boundary layer with high data density. Reynolds shear stress profiles in the near-wall region show that u v strongly depends on the measuring volume length; at a given y-position, u v decreases with increasing measuring volume length. This dependence was attributed to simultaneous validations on the U and V channels of Doppler bursts coming from different particles within the measuring volume. Moments of the streamwise velocity showed a slight dependence on measuring volume length, indicating that spatial averaging effects well known for hot-films and hot-wires can occur in laser velocimetry measurements when the data density is high.List of symbols time-averaged quantity - u wall friction velocity, ( w /)1/2 - v kinematic viscosity - d p pinhole diameter - l eff spanwise extent of LDV measuring volume viewed by photomultiplier - l + non-dimensional length of measuring volume, l eff u /v - y + non-dimensional coordinate in spanwise direction, y u /v - z + non-dimensional coordinate in spanwise direction, z u /v - U + non-dimensional mean velocity, /u - u instantaneous streamwise velocity fluctuation, U &#x2329;U - v instantaneous normal velocity fluctuation, V–V - u RMS streamwise velocity fluctuation, u 21/2 - v RMS normal velocity fluctuation, v 21/2 - Re Reynolds number based on momentum thickness, U 0/v - R uv cross-correlation coefficient, u v/u v - R12(0, 0, z) two point correlation between u and v with z-separation, <u(0, 0, 0) v (0, 0, z)>/<u(0, 0, 0) v (0, 0, 0)> - N rate at which bursts are validated by counter processor - T Taylor time microscale, u (dv/dt2)–1/2  相似文献   

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

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