共查询到20条相似文献，搜索用时 171 毫秒
1.
The problem of free vibration of non-linear structures is considered initially. It is shown that this problem can be represented as a non-linear eigenvalue problem. Variational principles for non-linear eigenvalue problems are defined. These variational principles are implemented with finite element models to define numerical approximations for the free vibration problem. The solution of these approximate equations provides a set of non-linear modal vectors and natural frequencies which vary with the amplitude of the solution. The non-linear eigenvalue parameters can be used in modal expansion approximations for the non-linear transient or steady state response of structural systems. To demonstrate the proposed techniques the free vibration and steady state vibration characteristics of a geometrically non-linear circular plate are determined. 相似文献
2.
We study the dynamics of soil moisture in the presence of precipitation, evapotranspiration and leakage. This study utilizes a probabilistic framework, and considers the evolution of the probability density function of soil moisture. Our main objective is to assess when the stationary solution is a physically realizable state of the system. Mathematically, this is equivalent to analyzing when the stationary solutions are stable. We show the following main results: For a general rainfall and evapotranspiration model, in a multiple plant framework, the stationary solution is a stable attractor in L1. However, the stationary solution is at best conditionally stable in L∞. This has important implications in terms of (1) the validity of moment analysis of the transient problem, (2) the physical importance of transient solutions, and (3) the numerical solution of the stationary problem as a late time solution of the transient problem. 相似文献
3.
P.W. France 《Journal of Hydrology》1974,21(4):381-398
A numerical method is presented for analysing either steady state or transient three-dimensional groundwater flow problems. The governing equation is formulated in terms of the finite element process using the Galerkin approach, and cubic isoparametric elements are used to simulate the flow domain as these permit accurate modelling of curved boundaries. Particular attention is paid to the time dependent movement of the phreatic surface where an iterative technique based on the replacement of the original transient problem by a discrete number of steady state problems is used to effect a solution. Furthermore, in tracing the movement of the surface use is made of the element formulation theory in order to compute the normal to the boundary.The validity of the technique is first established by analysing a radially symmetrical problem for which an alternative analytical solution is available. Finally, a general three-dimensional flow system is studied for which there is no known analytical solution. It is shown that relatively few elements are required to yield practical solutions. 相似文献
4.
In most previous studies on the dynamic response of a long cylindrical cavity subjected to internal transient dynamic loads, the porous medium was usually assumed to be completely saturated by ground water. In practice, however, the full saturation condition does not always exist. In this paper the surrounding soil and the lining of the cavity are respectively treated as a nearly saturated porous medium and an elastic material, and the governing equations for the dynamic problem are derived. A set of exact solutions are obtained in the Laplace transform domain for three types of transient loads, i.e. suddenly applied constant load, gradually applied step load and triangular pulse load. By utilizing a reliable numerical method of inverse Laplace transforms, the time-domain solutions are then presented. The influence of the degree of saturation of the surrounding soil on the dynamic response of the lined cavity is examined for numerical examples. 相似文献
5.
Matthew J. Simpson 《Ground water》2018,56(2):337-342
Predicting the amount of time required for a transient groundwater response to take place is a practical question that is of interest in many situations. This time scale is often called the response time. In the groundwater hydrology literature, there are two main methods used to calculate the response time: (1) both the transient and steady‐state groundwater flow equations are solved, and the response time is taken to be amount of time required for the transient solution to approach the steady solution within some tolerance; and (2) simple scaling arguments are adopted. Certain limitations restrict both of these approaches. In this study, we outline a third method, based on the theory of mean action time. We derive the governing boundary value problem for both the mean and variance of action time for confined flow in two‐dimensional heterogeneous porous media. Importantly, we show that these boundary value problems can be solved using widely available software. Applying these methods to a test case reveals the advantages of the theory of mean action time relative to standard methods. 相似文献
6.
Saltwater intrusion problems have been usually tackled through analytical models because of its simplicity, easy implementation and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, which neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. This paper provides insight into the validity of a sharp‐interface approximation defined from a steady state solution when applied to transient seawater intrusion problems. The validation tests have been performed on a 3D unconfined synthetic aquifer, which include spatial and temporal distribution of recharge and pumping wells. Using a change of variable, the governing equation of the steady state sharp‐interface problem can be written with the same structure of the steady confined groundwater flow equation as a function of a single potential variable (?). We propose to approach also the transient problem solving a single potential equation (using also the ? variable) with the same structure of the confined groundwater flow equation. It will allow solving the problem by using the classical MODFLOW code. We have used the parameter estimation model PEST to calibrate the parameters of the transient sharp‐interface equation. We show how after the calibration process, the sharp‐interface approach may provide accurate enough results when applied to transient problems and improve the steady state results, thus avoiding the need of implementing a density‐dependent model and reducing the computational cost. This has been proved by comparing results with those obtained using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. The comparison was performed in terms of piezometric heads, seawater penetration, transition zone width and critical pumping rates. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
7.
Scattering of elastic waves by two dimensional multilayered dipping sediments of arbitrary shape embedded in an elastic half-sapce is investigated by using a bondary method. The displancement field is evaluated throughout the elastic media for both steady state and transient incident SH waves. The unknown scattered field is expressed in terms of wave functions which satisfy the equation of motion, traction-free boundary condition and appropariate radiation conditions. The transient response is constructed from the steady state solution by using the fast Fourier transform technique. The numerical results presented demonstrate that scattering of waves by subsurface irregularities may cause locally very large amplification of surface ground motion. The motion can be affected greatly by the scattered surface waves in the sediments. The results clearly indicate that the surface ground motion depends upon a number of parameters present in the problem, such as frequency and the angle of incidence of the incoming wave, impedance contrast between the layers and location of the observation point. 相似文献
8.
J. R. Sturgul 《Pure and Applied Geophysics》1967,68(1):66-82
Summary The techniques of Complex Variables are used to obtain a solution for a general notch problem having a load distributed over any portion of its free surface. This solution is very general in that it can be used for many different shapes of notches, including protrusions, by simply varying a few parameters. The same techniques are then used to solve the problem of a notched plate acted on by loads at infinity. This solution is of a very general nature. One particular solution holds for the problem studied in an earlier publication. 相似文献
9.
This paper presents an approach for coupling MODFLOW and MT3DMS for the simulation of variable-density ground water flow. MODFLOW routines were modified to solve a variable-density form of the ground water flow equation in which the density terms are calculated using an equation of state and the simulated MT3DMS solute concentrations. Changes to the MODFLOW and MT3DMS input files were kept to a minimum, and thus existing data files and data files created with most pre- and postprocessors can be used directly with the SEAWAT code. The approach was tested by simulating the Henry problem and two of the saltpool laboratory experiments (low- and high-density cases). For the Henry problem, the simulated results compared well with the steady-state semianalytic solution and also the transient isochlor movement as simulated by a finite-element model. For the saltpool problem, the simulated breakthrough curves compared better with the laboratory measurements for the low-density case than for the high-density case but showed good agreement with the measured salinity isosurfaces for both cases. Results from the test cases presented here indicate that the MODFLOW/MT3DMS approach provides accurate solutions for problems involving variable-density ground water flow and solute transport. 相似文献
10.
A model for investigations of ground motions due to continuously moving loads with constant and time-varying amplitudes is presented. The vertical displacements excited by moving load areas are obtained at a fixed observation point at the surface of the three-dimensional halfspace in time domain. The load is moving along a straight line with constant speed. To solve this nonaxisymmetric, initial boundary value problem a semi-analytical, discretized model is developed. It is based on Green's functions for a suddenly applied, stationary surface point load with Heaviside time dependency. These functions, also called influence functions of the halfspace, are valid for any homogeneous, isotropic and linear-elastic medium. The principle of superposition is used. Results are shown for the transient and the steady-state ground motions, and they are compared with analytical solutions. The load speed is varied in the subcritical range up to the propagation velocity of Rayleigh-waves. 相似文献
11.
Vertical wells with radial extension at the well bottom can improve the rate of water production. No study has yet investigated the effects of the transient state and anisotropy in directional hydraulic conductivities on the wellbore flux rate for this type of well. This study derives a semianalytical transient drawdown solution for constant-head pumping at a fully penetrating well radially extended at the bottom of a confined, anisotropic aquifer by applying Laplace transform and separation of variables as well as conducting a Fourier analysis. The results of this new solution indicate that transient and steady-state wellbore flux rates can be increased by a factor of two for greater radial extension of the well. Compared with an isotropic aquifer (a ratio of vertical and horizontal hydraulic conductivities equal to one), an anisotropic aquifer with the ratio less than one may produce a higher transient wellbore flux rate and lower steady-state wellbore flux rate. Moreover, the time required to achieve the steady-state wellbore flux rate can be substantially affected by anisotropy of the aquifer. 相似文献
12.
Kalyan Kumar Bagchi Sukla Mukherjee Lokenath Debnath 《Pure and Applied Geophysics》1973,111(1):2255-2263
Summary With the aid of the generalized function method, a study is made of the linearized theory of transient development of capillary-gravity waves in an inviscid, incompressible and homogeneous liquid of finite and infinite depth due to an arbitrary oscillating source situated at a finite depth below the undisturbed free surface of the liquid. The initial value problem is solved by using Laplace-Fourier transforms combined with asymptotic methods. The asymptotic solution is found to consist of the steady state and the transient components which are independently modified by surface tension. The latter decays more rapidly as timet due to the presence of surface tension than in the case where surface tension is neglected. It is predicted that the principal effect of surface tension is to increase both the phase and group velocity of the waves and make the energy more readily available among the rapidly travelling progressive surface waves. In addition to the effects of surface tension on the physical properties of the wave motions, our method of solution provides an interesting illustration of the applicability of generalized functions in water wave phenomena. 相似文献
13.
A Fourier transform approach is applied to the transient analysis of dynamic soil–structure interaction under SH-motion. The governing equations are formulated in the frequency domain using a Finite Element–Boundary Element (FE–BE) coupling method. After solving the transformed problem, the transient solution is obtained using the discrete inverse Fourier transform with a fast Fourier transform algorithm. Two examples are presented in order to show the numerical performance of the proposed technique. 相似文献
14.
Emil O. Frind 《Advances in water resources》1982,5(2):73-88
A numerical model for the economical simulation of long-term transient response in density-dependent transport problems is introduced. Although a classical Galerkin finite element approach is used, emphasis on optimum efficiency throughout the development results in a scheme that is found to be significantly less costly than comparable existing schemes. This advantage in efficiency increases the scope of simulation problems that can be handled within the constraints of a limited research budget. Some distinctive aspects are the elimination of static quantities in the fluid continuity equation, achieved by the introduction of equivalent freshwater head, and the elimination of numerical integration, achieved by the deliberate choice of linear elements. As a result of this choice, fluid velocities are discontinuous across the element boundaries. It is shown, however, that the solution obtained with discontinuous velocities approaches that obtained with continuous velocities as the grid is refined, and that the two types of solutions give essentially the same results when the elements are in the same size range. The model is applied to simulate the complete transient response for a well-known problem of seawater intrusion in a confined aquifer. The simulation is performed with both the constant dispersion coefficient used by previous researchers, and a more physically realistic velocity-dependent dispersion coefficient. Responses are found to be substantially different for the two types of coefficients, with the velocity-dependent dispersion coefficient producing much slower convergence to a state of dynamic equilibrium, and a much more pointed saltwater toe, which at the bottom of the aquifer tends to a sharp interface at equilibrium. Finally, it is shown by means of large-scale applications that the model is capable of efficiently simulating the long-term transient response in systems of practical significance. 相似文献
15.
Barry Saltzman 《Pure and Applied Geophysics》1968,69(1):237-259
Summary Equations governing the axially-symmetric time-average state of the atmosphere and the transient departures from this mean state are set down. As a first step toward a solution of this system for seasonal average conditions, a model is formulated based on the thermodynamical energy equation for the vertical average of the mean state, and on the perturbation solutions of the linearized equations governing the baroclinic growth of transient eddies. All forms of non-adiabatic heating within the atmosphere and at the earth's surface are parameterized. The resulting differential equation governing the axially-symmetric mean potential temperature distribution takes the form of a steadystate diffusion equation in surface spherical coordinates, with a variable Austausch coefficient which is to be determined iteratively as a dependent variable.Global solutions, for winter and summer equilibrium conditions, are obtained for the thermal structure, the heat balance components, the transient eddy variances of temperature and meridional wind speed, and the covariance representing the meridional eddy heat transport. These solutions are for different types of surface conditions (ocean, land), and for a successively more complete variety of modes of heat transfer ranging from pure radiation to a combination of radiation, latent heat processes, and conduction and convection within the atmosphere and the subsurface layers. The results for this latter complete case seem to be a reasonable first order approximation to the observed distributions. Suggestions are made for improving and generalizing the study. 相似文献
16.
On the solution of transient free-surface flow problems in porous media by the finite element method
A numerical procedure is presented to deal with solution of transient free-surface flows in porous media. The governing boundary-value problem for the piezometric potential is solved by the finite element method. The initial-value problem which describes the transient motion of the free-surface is solved by the method of quasi-linearization. The numerical scheme has been applied to isotropic and anisotropic earth dam problem and also to a ditch drainage problem. Excellent agreements have been reached when compared with known solutions. This computational procedure is shown to be stable and suitable for this class of problems with the aid of a digital computer. 相似文献
17.
Simon A. Mathias 《Ground water》2010,48(3):438-441
When seeking to predict plume geometry resulting from fluid injection through partially penetrating wells, it is common to assume a steady-state spherically diverging flow field. In reality, the flow field is transient. The steady-flow assumption is likely to cause overestimation of injection plume radius since the accommodation of fluid by increases in porosity and fluid density is ignored. In this paper, a transient solution is developed, resulting in a nonlinear ordinary differential equation expressing plume radius as a function of time. It is shown that the problem can be fully described by one type curve. A critical time, tc, is identified at which the percentage error of the steady-state flow solution compared to the fully dynamic problem is less than 1%. Only for large injection rates and low permeabilities, does tc become greater than 1 h. Nevertheless, an improved approximate solution is obtained by a simple linearization procedure. The critical time, tc for the new approximate solution is 0.3% of that required for the steady-state flow solution. 相似文献
18.
M. J. KUHN 《Geophysical Prospecting》1985,33(8):1103-1137
This paper considers propagation of elastodynamic waves in an imperfectly elastic half-space. Two different excitation modes are investigated: a buried source of compressional waves and a vertically directed areal load applied to the surface. Numerical integration of the analytical solution of the wave equation allows study of the vertical and horizontal components of displacement and/or particle velocity anywhere in the half-space. One case of particular interest concerns the examination of particle displacement and velocity at the surface in a circular area above the source. In another application seismograms generated by an explosive buried source are contrasted with seismograms generated by the transient application of a vertically directed load to the free surface. Still another application of considerable practical interest concerns the study of the nongeometrical pS—wave, in particular its characteristics as functions of range and depth. Finally, in the last application the behavior of a rarely observed wave (denoted here by the letter U) is studied in both elastic and visco-elastic half-spaces. 相似文献
19.
Generalized linear models may be used as a systematic and flexible alternative to existing practices in fitting sediment or solute concentration and load rating curves. Through widely available statistical software they offer a one‐step solution to the common and awkward problem that a bias correction is needed when estimating total loads. In particular, using a logarithmic link function has an advantage over logarithmic transformation of concentration or load, as generalized linear models return predictions directly on the original measured scale; there is thus no need for back‐transformation or subsequent correction factors. Example analyses are given for suspended sediment concentration for an upland stream in the Northern Pennines and for phosphorus load for the Illinois River. General criteria for choosing models and good practice in the use of figures of model merit and of residual plots are also discussed in detail. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
20.
An exact theoretical formulation is presented for the analysis of a thin-walled pile embedded in an elastic half-space under vertically-incident P-wave excitation. In the framework of three-dimensional elastodynamics and a shell theory, the axisymmetrical wave-scattering problem is shown to be reducible to a set of Fredholm boundary integral equations. With the incorporation of the singular characteristics of the wave-induced contact load distributions into the solution scheme, a computational boundary element method is developed for a rigorous treatment of the seismic soil-structure interaction problem. Typical results for the dynamic contact load distributions, displacements, complex-valued foundation input motion functions, and resonant pile foundation response are included for direct engineering applications. 相似文献