首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The gravitational interaction between two objects on similar orbits can effect noticeable changes in the orbital evolution even if the ratio of their masses to that of the central body is vanishingly small. Christou (Icarus 174:215–229, 2005) observed an occasional resonant lock in the differential node \(\varDelta \varOmega \) between two members in the Himalia irregular satellite group of Jupiter in the N-body simulations (corresponding mass ratio \(\sim 10^{-9}\)). Using a semianalytical approach, we have reproduced this phenomenon. We also demonstrate the existence of two additional types of resonance, involving angle differences \(\varDelta \omega \) and \(\varDelta (\varOmega +\varpi )\) between two group members. These resonances cause secular oscillations in eccentricity and/or inclination on timescales \(\sim \)1 Myr. We locate these resonances in (aei) space and analyse their topological structure. In subsequent N-body simulations, we confirm these three resonances and find a fourth one involving \(\varDelta \varpi \). In addition, we study the occurrence rates and the stability of the four resonances from a statistical perspective by integrating 1000 test particles for 100 Myr. We find \(\sim \)10 to 30 librators for each of the resonances. Particularly, the nodal resonance found by Christou is the most stable: 2 particles are observed to stay in libration for the entire integration.  相似文献   

2.
Jupiter’s satellites are subject to strong tidal forces which result in variations of the gravitational potential and deformations of the satellites’ surfaces on the diurnal tidal cycle. Such variations are described by the Love numbers \(k_2\) and \(h_2\) for the tide-induced potential variation due to internal mass redistribution and the radial surface displacement, respectively. The phase-lags \( \phi _{k_2}\) and \( \phi _{h_2}\) of these complex numbers contain information about the rheological and dissipative states of the satellites. Starting from interior structure models and assuming a Maxwell rheology to compute the tidal deformation, we calculate the phase-lags in application to Ganymede and Europa. For both satellites we assume a decoupling of the outer ice-shell from the deep interior by a liquid subsurface water ocean. We show that, in this case, the phase-lag difference \(\varDelta \phi = \phi _{k_2}- \phi _{h_2}\) can provide information on the rheological and thermal state of the deep interiors if the viscosities of the deeper layers are small. In case of Ganymede, phase-lag differences can reach values of a few degrees for high-pressure ice viscosities \({<}10^{14}\) Pa s and would indicate a highly dissipative state of the deep interior. In this case \(\varDelta \phi \) is dominated by dissipation in the high-pressure ice layer rather than dissipation within the ice-I shell. These phase lags would be detectable from spacecraft in orbit around the satellite. For Europa \(\varDelta \phi \) could reach values exceeding \(20^\circ \) and phase-lag measurements could help distinguish between (1) a hot dissipative silicate mantle which would in thermal equilibrium correspond to a very thin outer ice-I shell and (2) a cold deep interior implying that dissipation would mainly occur in a thick (several tens of km) outer ice-I shell. These measurements are highly relevant for ESA’s Jupiter Icy Moons Explorer (JUICE) and NASA’s Europa Multiple Flyby Mission, both targeted for the Jupiter system.  相似文献   

3.
Pulsation period changes in Mira type variables are investigated using the stellar evolution and nonlinear stellar pulsation calculations. We considered the evolutionary sequence of stellar models with initial mass \({M_{ZAMS}} = \;3{M_ \odot }\) and population I composition. Pulsations of stars in the early stage of the asymptotic giant branch are shown to be due to instability of the fundamental mode. In the later stage of evolution when the helium shell source becomes thermally unstable the stellar oscillations occur in either the fundamental mode (for the stellar luminosuty \(L < 5.4 \times {10^3}{L_ \odot }\)) or the first overtone (\(L > 7 \times {10^3}{L_ \odot }\)). Excitation of pulsations is due to the κ-mechanism in the hydrogen ionization zone. Stars with intermediate luminosities \(5.4 \times {10^3}{L_ \odot } < L < 7 \times {10^3}{L_ \odot }\) were found to be stable against radial oscillations. The pulsation period was determined as a function of evolutionary time and period change rates \(\dot \Pi \) were evaluated for the first ten helium flashes. The period change rate becomes the largest in absolute value \((\dot \Pi /\Pi \approx - {10^{ - 2}}y{r^{ - 1}})\) between the helium flash and the maximum of the stellar luminosity. Period changes with rate \(\left| {\dot \Pi /\Pi } \right| \geqslant - {10^{ - 3}}y{r^{ - 1}}\) take place during ≈500 yr, that is nearly one hundredth of the interval between helium flashes.  相似文献   

4.
Many models of eruptive flares or coronal mass ejections (CMEs) involve formation of a current sheet connecting the ejecting CME flux rope with a magnetic loop arcade. However, there is very limited observational information on the properties and evolution of these structures, hindering progress in understanding eruptive activity from the Sun. In white-light images, narrow coaxial rays trailing the outward-moving CME have been interpreted as current sheets. Here, we undertake the most comprehensive statistical study of CME-rays to date. We use SOHO/LASCO data, which have a higher cadence, larger field of view, and better sensitivity than any previous coronagraph. We compare our results to a previous study of Solar Maximum Mission (SMM) CMEs, in 1984?–?1989, having candidate magnetic disconnection features at the CME base, about half of which were followed by coaxial bright rays. We examine all LASCO CMEs during two periods of minimum and maximum activity in Solar Cycle 23, resulting in many more events, \(\sim130\) CME-rays, than during SMM. Important results include: The occurrence rate of the rays is \(\sim11~\%\) of all CMEs during solar minimum, but decreases to \(\sim7~\%\) at solar maximum; this is most likely related to the more complex coronal background. The rays appear on average 3?–?4 hours after the CME core, and are typically visible for three-fourths of a day. The mean observed current sheet length over the ray lifetime is \(\sim12~R_{\odot}\), with the longest current sheet of \(18.5~R_{\odot}\). The mean CS growth rates are \(188~\mbox{km}\,\mathrm{s}^{-1}\) at minimum and \(324~\mbox{km}\,\mathrm{s}^{-1}\) at maximum. Outward-moving blobs within several rays, which are indicative of reconnection outflows, have average velocities of \(\sim350~\mbox{km}\,\mathrm{s}^{-1}\) with small positive accelerations. A pre-existing streamer is blown out in most of the CME-ray events, but half of these are observed to reform within \(\sim1\) day. The long lifetime and long lengths of the CME-rays challenge our current understanding of the evolution of the magnetic field in the aftermath of CMEs.  相似文献   

5.
Hongqi Zhang 《Solar physics》2016,291(12):3501-3517
We present the photospheric energy density of magnetic fields in two solar active regions (one of them recurrent) inferred from observational vector magnetograms, and compare it with other available differently defined energy parameters of magnetic fields in the photosphere. We analyze the magnetic fields in Active Regions NOAA 6580-6619-6659 and 11158. The quantity \(\frac{1}{4\pi}{\mathbf{B}}_{n}\cdot{\mathbf{B}}_{p}\) is an important energy parameter that reflects the contribution of magnetic shear to the difference between the potential (\(\mathbf{B}_{p}\)) and the non-potential magnetic field (\(\mathbf{B}_{n}\)), and also the contribution to the free magnetic energy near the magnetic neutral lines in the active regions. It is found that the photospheric mean magnetic energy density shows clear changes before the powerful solar flares in Active Region NOAA 11158, which is consistent with the change in magnetic fields in the flaring lower atmosphere.  相似文献   

6.
The determination of analytical expressions which, including the main perturbative effects, allow the retrieval of the orbit elements of a probe represents an important requirement in designing science trajectories. One of these perturbations is given by the third body attraction. The case in which the perturbing body moves on a plane coincident with the equatorial plane of the primary body has been investigated in previous studies and equations able to provide the temporal evolution of the orbit elements have been determined and applied to the main moons of the Solar System. In this paper an extension of this topic has been carried out and equations which allow the determination of the orbit evolution have been analytically retrieved in the general case in which one or more perturbing bodies describe elliptical and inclined orbits with respect to the equatorial plane of the primary. Then, introducing these equations into the periodicity condition for the probe ground track, and considering the \(J_{2}\) and \(J_{4}\) effects coming from the primary body, an equation able to provide repeating ground track orbits has been determined.  相似文献   

7.
Analysis of the radial velocities based on spectra of high (near the H α line) and moderate (4420–4960 Å) resolutions supplemented by the published radial velocities has revealed the binarity of a bright member of the young open star cluster χ Per, the star V622 Per. The derived orbital elements of the binary show that the lines of both components are seen in its spectrum, the orbital period is 5.2 days, and the binary is in the phase of active mass exchange. The photometric variability of the star is caused by the ellipsoidal shape of its components. Analysis of the spectroscopic and photometric variabilities has allowed the absolute parameters of the binary’s orbit and its components to be found. V622 Per is shown to be a classical Algol with moderate mass exchange in the binary. Mass transfer occurs from the less massive (\({M_1} = 9.1 \pm 2.7{M_ \odot }\)) but brighter (\(\log {L_1} = 4.52 \pm 0.10{L_ \odot }\)) component onto the more massive (\({M_2} = 13.0 \pm 3.5{M_ \odot }\)) and less bright (\(\log {L_2} = 3.96 \pm 0.10{L_ \odot }\)) component. Analysis of the spectra has confirmed an appreciable overabundance of CNO-cycle products in the atmosphere of the primary component. Comparison of the positions of the binary’s components on the T eff–log g diagram with the age of the cluster χ Per points to a possible delay in the evolution of the primary component due to mass loss by no more than 1–2Myr.  相似文献   

8.
Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.  相似文献   

9.
We consider the three-dimensional bounded motion of a test particle around razor-thin disk configurations, by focusing on the adiabatic invariance of the vertical action associated with disk-crossing orbits. We find that it leads to an approximate third integral of motion predicting envelopes of the form \(Z(R)\propto [\varSigma (R)]^{-1/3}\), where R is the radial galactocentric coordinate, Z is the z-amplitude (vertical amplitude) of the orbit and \(\varSigma \) represents the surface mass density of the thin disk. This third integral, which was previously formulated for the case of flattened 3D configurations, is tested for a variety of trajectories in different thin-disk models.  相似文献   

10.
More than 80 giant planets are known by mass and radius. Their interior structure in terms of core mass, number of layers, and composition however is still poorly known. An overview is presented about the core mass M core and envelope mass of metals M Z in Jupiter as predicted by various equations of state. It is argued that the uncertainty about the true H/He EOS in a pressure regime where the gravitational moments J 2 and J 4 are most sensitive, i.e. between 0.5 and 4 Mbar, is in part responsible for the broad range \(M_{\mathit{core}}=0{-}18\:M_{\oplus }\), \(M_{Z}=0{-}38\:M_{\oplus }\), and \(M_{\mathit{core}}+M_{Z}=14{-}38\:M_{\oplus }\) currently offered for Jupiter. We then compare the Jupiter models obtained when we only match J 2 with the range of solutions for the exoplanet \(\mathrm{GJ}\:436\mathrm{b}\), when we match an assumed tidal Love number k 2 value.  相似文献   

11.
This paper studies the existence and stability of equilibrium points under the influence of small perturbations in the Coriolis and the centrifugal forces, together with the non-sphericity of the primaries. The problem is generalized in the sense that the bigger and smaller primaries are respectively triaxial and oblate spheroidal bodies. It is found that the locations of equilibrium points are affected by the non-sphericity of the bodies and the change in the centrifugal force. It is also seen that the triangular points are stable for 0<μ<μ c and unstable for \(\mu_{c}\le\mu <\frac{1}{2}\), where μ c is the critical mass parameter depending on the above perturbations, triaxiality and oblateness. It is further observed that collinear points remain unstable.  相似文献   

12.
The kinetic Alfven waves in the presence of homogeneous magnetic field plasma with multi-ions effect are investigated. The dispersion relation and normalised damping rate are derived for low-\(\beta\) plasma using kinetic theory. The effect of density variation of \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions is observed on frequency and damping rate of the wave. The variation of frequency (\(\omega\)) and normalised damping rate (\(\gamma / \varOmega_{H^{ +}} \)) of the wave are studied with respect to \(k_{ \bot} \rho_{j}\), where \(k_{ \bot} \) is the perpendicular wave number, \(\rho_{j}\) is the ion gyroradius and \(j \) denotes \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions. The variation with \(k_{ \bot} \rho_{j}\) is considered over wide range. The parameters appropriate to cusp region are used for the explanation of results. It is found that with hydrogen and helium ions gyration, the frequency of wave is influenced by the density variation of \(\text{H}^{+}\) and \(\text{He}^{+}\) ions but remains insensitive to the change in density of \(\text{O}^{+}\) ions. For oxygen ion gyration, the frequency of wave varies over a short range only for \(\text{O}^{+}\) ion density variation. The wave shows damping at lower altitude due to variation in density of lighter \(\text{H}^{+}\) and \(\text{He}^{+}\) ions whereas at higher altitude only heavy \(\text{O}^{+}\) ions contribute in wave damping. The damping of wave may be due to landau damping or energy transfer from wave to particles. The present study signifies that the both lighter and heavier ions dominate differently to change the characteristics of kinetic Alfven wave and density variation is also an important parameter to understand wave phenomena in cusp region.  相似文献   

13.
We use the data for the \(\text{H}\beta\) emission-line, far-ultraviolet (FUV) and mid-infrared 22 μm continuum luminosities to estimate star formation rates \(\langle \mbox{SFR} \rangle \) averaged over the galaxy lifetime for a sample of about 14000 bursting compact star-forming galaxies (CSFGs) selected from the Data Release 12 (DR12) of the Sloan Digital Sky Survey (SDSS). The average coefficient linking \(\langle \mbox{SFR} \rangle \) and the star formation rate \(\mbox{SFR}_{0}\) derived from the \(\text{H}\beta\) luminosity at zero starburst age is found to be 0.04. We compare \(\langle \mbox{SFR} \rangle \mbox{s}\) with some commonly used SFRs which are derived adopting a continuous star formation during a period of \({\sim}\,100~\mbox{Myr}\), and find that the latter ones are 2–3 times higher. It is shown that the relations between SFRs derived using a geometric mean of two star-formation indicators in the UV and IR ranges and reduced to zero starburst age have considerably lower dispersion compared to those with single star-formation indicators. We suggest that our relations for \(\langle \mbox{SFR} \rangle \) determination are more appropriate for CSFGs because they take into account a proper temporal evolution of their luminosities. On the other hand, we show that commonly used SFR relations can be applied for approximate estimation within a factor of \({\sim}\,2\) of the \(\langle \mbox{SFR} \rangle \) averaged over the lifetime of the bursting compact galaxy.  相似文献   

14.
The most used method to calculate the coronal electron temperature [\(T_{\mathrm{e}} (r)\)] from a coronal density distribution [\(n_{\mathrm{e}} (r)\)] is the scale-height method (SHM). We introduce a novel method that is a generalization of a method introduced by Alfvén (Ark. Mat. Astron. Fys. 27, 1, 1941) to calculate \(T_{\mathrm{e}}(r)\) for a corona in hydrostatic equilibrium: the “HST” method. All of the methods discussed here require given electron-density distributions [\(n_{\mathrm{e}} (r)\)] which can be derived from white-light (WL) eclipse observations. The new “DYN” method determines the unique solution of \(T_{\mathrm{e}}(r)\) for which \(T_{\mathrm{e}}(r \rightarrow \infty) \rightarrow 0\) when the solar corona expands radially as realized in hydrodynamical solar-wind models. The applications of the SHM method and DYN method give comparable distributions for \(T_{\mathrm{e}}(r)\). Both have a maximum [\(T_{\max}\)] whose value ranges between 1?–?3 MK. However, the peak of temperature is located at a different altitude in both cases. Close to the Sun where the expansion velocity is subsonic (\(r < 1.3\,\mathrm{R}_{\odot}\)) the DYN method gives the same results as the HST method. The effects of the other free parameters on the DYN temperature distribution are presented in the last part of this study. Our DYN method is a new tool to evaluate the range of altitudes where the heating rate is maximum in the solar corona when the electron-density distribution is obtained from WL coronal observations.  相似文献   

15.
We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.  相似文献   

16.
This paper deals with the photo-gravitational restricted four-body problem (PR4BP) with variable mass. Following the procedure given by Gascheau (C. R. 16:393–394, 1843) and Routh (Proc. Lond. Math. Soc. 6:86–97, 1875), the conditions of linear stability of Lagrange triangle solution in the PR4BP are determined. The three radiating primaries having masses \(m_{1}\), \(m_{2}\) and \(m_{3}\) in an equilateral triangle with \(m_{2}=m_{3}\) will be stable as long as they satisfy the linear stability condition of the Lagrangian triangle solution. We have derived the equations of motion of the mentioned problem and observed that there exist eight libration points for a fixed value of parameters \(\gamma (\frac{m \ \text{at time} \ t}{m \ \text{at initial time}}, 0<\gamma\leq1 )\), \(\alpha\) (the proportionality constant in Jeans’ law (Astronomy and Cosmogony, Cambridge University Press, Cambridge, 1928), \(0\leq\alpha\leq2.2\)), the mass parameter \(\mu=0.005\) and radiation parameters \(q_{i}, (0< q_{i}\leq1, i=1, 2, 3)\). All the libration points are non-collinear if \(q_{2}\neq q_{3}\). It has been observed that the collinear and out-of-plane libration points also exist for \(q_{2}=q_{3}\). In all the cases, each libration point is found to be unstable. Further, zero velocity curves (ZVCs) and Newton–Raphson basins of attraction are also discussed.  相似文献   

17.
The Diffuse Galactic Syncrotron Emission (DGSE) is the most important diffuse foreground component for future cosmological 21-cm observations. The DGSE is also an important probe of the cosmic ray electron and magnetic field distributions in the turbulent interstellar medium (ISM) of our galaxy. In this paper we briefly review the Tapered Gridded Estimator (TGE) which can be used to quantify the angular power spectrum C ? of the sky signal directly from the visibilities measured in radio-interferometric observations. The salient features of the TGE are: (1) it deals with the gridded data which makes it computationally very fast, (2) it avoids a positive noise bias which normally arises from the system noise inherent to the visibility data, and (3) it allows us to taper the sky response and thereby suppresses the contribution from unsubtracted point sources in the outer parts and the side lobes of the antenna beam pattern. We also summarize earlier work where the TGE was used to measure the C ? of the DGSE using 150 MHz GMRT data. Earlier measurements of C ? are restricted to \(\ell \le \ell _{\max } \sim 10^{3}\) for the DGSE, the signal at the larger ? values is dominated by the residual point sources after source subtraction. The higher sensitivity of the upcoming SKA1 Low will allow the point sources to be subtracted to a fainter level than possible with existing telescopes. We predict that it will be possible to measure the C ? of the DGSE to larger values of \(\ell _{\max }\) with SKA1 Low. Our results show that it should be possible to achieve \(\ell _{\max }\sim 10^{4}\) and ~105 with 2 minutes and 10 hours of observations respectively.  相似文献   

18.
Based on the new catalogue of solar proton events (SPEs) for the period of 1997?–?2009 (Solar Cycle 23) we revisit the long-studied problem of the event-size distributions in the context of those constructed for other solar-flare parameters. Recent results on the problem of size distributions of solar flares and proton events are briefly reviewed. Even a cursory acquaintance with this research field reveals a rather mixed and controversial picture. We concentrate on three main issues: i) SPE size distribution for \({>}\,10~\mbox{MeV}\) protons in Solar Cycle 23; ii) size distribution of \({>}\,1~\mbox{GV}\) proton events in 1942?–?2014; iii) variations of annual numbers for \({>}\,10~\mbox{MeV}\) proton events on long time scales (1955?–?2015). Different results are critically compared; most of the studies in this field are shown to suffer from vastly different input datasets as well as from insufficient knowledge of underlying physical processes in the SPEs under consideration. New studies in this field should be made on more distinct physical and methodological bases. It is important to note the evident similarity in size distributions of solar flares and superflares in Sun-like stars.  相似文献   

19.
Previous analysis of magnetohydrodynamic-scale currents in high-speed solar wind near 1 AU suggests that the most intense current-carrying structures occur at electron scales and are characterized by average current densities on the order of \(1~\mbox{pA}/\mbox{cm}^{2}\). Here, this prediction is verified by examining the effects of the measurement bandwidth and/or measurement resolution on the analysis of synthetic solar wind signals. Assuming Taylor’s hypothesis holds for the energetically dominant fluctuations at kinetic scales, the results show that when \(\nu_{c}\gg \nu_{b}\), where \(\nu_{c}\) is the measurement bandwidth and \(\nu_{b} \approx 1/3~\mbox{Hz}\) is the break frequency, the average scale of the most intense fluctuations in the current density proxy is approximately \(1/\nu_{c}\), and the average peak current density is a weakly increasing function that scales approximately like \(\nu_{c}^{0.1}\).  相似文献   

20.
We give arguments for a basically unified formation mechanism of slow (Lynden-Bell) and fast (common) galactic bars. This mechanism is based on an instability that is akin to the well-known instability of radial orbits and is produced by the mutual attraction and alignment of precessing stellar orbits (so far, only the formation of slow bars has been explained in this way). We present a general theory of the low-frequency modes in a disk that consists of orbits precessing at different angular velocities. The problem of determining these modes is reduced to integral equations of moderately complex structure. The characteristic pattern angular velocities Ωp of the low-frequency modes are of the order of the mean orbital precession angular velocity \(\bar \Omega _{pr}\). Bar modes are also among the low-frequency modes; while \(\Omega _p \approx \bar \Omega _{pr}\) for slow bars, Ωp for fast bars can appreciably exceed even the maximum orbital precession angular velocity in the disk Ω pr max (however, it remains of the order of these precession angular velocities). The possibility of such an excess of Ωp over Ω pr max is associated with the effect of “repelling” orbits. The latter tend to move in a direction opposite to the direction in which they are pushed. We analyze the pattern of orbital precession in potentials typical of galactic disks. We note that the maximum radius of an “attracting” circular orbit rc can serve as a reasonable estimate of the bar length lb. Such an estimate is in good agreement with the available results of N-body simulations.  相似文献   

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

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