Effect of boundary conditions on shakedown analysis of heterogeneous materials

The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.

Least Square Finite Element Model for Analysis of Multilayered Composite Plates under Arbitrary Boundary Conditions

Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, anisotropic rectangular plate subjected to arbitrary boundary conditions on one edge and simply supported on other edge. The rectangular laminate consists of anisotropic and homogeneous laminae of arbitrary thicknesses. This study presents the elastic analysis of laminated composite plates subjected to sinusoidal mechanical loading under arbitrary boundary conditions. Least square finite element solutions for displacements and stresses are investigated using a mathematical model, called a state-space model, which allows us to simultaneously solve for these field variables in the composite structure’s domain and ensure that continuity conditions are satisfied at layer interfaces. The governing equations are derived from this model using a numerical technique called the least-squares finite element method (LSFEM). These LSFEMs seek to minimize the squares of the governing equations and the associated side conditions residuals over the computational domain. The model is comprised of layerwise variables such as displacements, out-of-plane stresses, and in- plane strains, treated as independent variables. Numerical results are presented to demonstrate the response of the laminated composite plates under various arbitrary boundary conditions using LSFEM and compared with the 3D elasticity solution available in the literature.

Boundary Conditions for Momentum and Vorticity at an Interface between Two Fluids

Boundary conditions for momentum and vorticity have been precisely derived, paying attention to the physical meaning of each mathematical expression of terms rigorously obtained from the basic equations: Navier-Stokes equation and the equation of vorticity transport. It has been shown first that a contribution of fluid molecules crossing over a conceptual surface moving with fluid velocity due to their fluctuating motion is essentially important to understanding transport phenomena of momentum and vorticity. A notion of surface layers, which are thin layers at both sides of an interface, has been introduced next to elucidate the transporting mechanism of momentum and vorticity from one phase to the other at an interface through which no fluid molecules are crossing over. A fact that a size of δV, in which reliable values of density, momentum, and velocity of fluid are respectively defined as a volume-averaged mass of fluid molecules, a volume-averaged momentum of fluid molecules and a mass-averaged velocity of fluid molecules, is not infinitesimal but finite has been one of the key factors leading to the boundary conditions for vorticity at an interface between two fluids. The most distinguished characteristics of the boundary conditions derived here are the zero-value conditions for a normal component of momentum flux and tangential components of vorticity flux, at an interface.

Dynamic characterization of tailing dam using fully coupled dynamic analysis with different boundary conditions-a case study

The seismic response analysis of a tailing dam is studied using a fully coupled effective stress approach in conjunction with an advanced multi yield surface plastic constitutive model for tailing material.Strain controlled static and cyclic triaxial tests were carried out to obtain the constitutive model for the tailing material.The tailing materials were collected from the Rampura Agucha tailing dam(Rajasthan State,India).A 2D nonlinear finite element(FE)model was then developed using different boundary conditions from the tailing embankment constructed using the downstream and upstream method of rising using OpenSees software.In first case,the model boundary was fixed in both the X and Y directions,and in the second case,viscous dashpots were introduced for both side and horizontal boundaries.The model was validated with experimental results on tailing material.Analyses were carried out considering five different earthquake motions,which were applied at the base.Comparisons of the different boundary conditions in terms of displacement flow vectors,pore pressure and stress-strain curves during shaking are presented.From the analysis,it was observed that the viscous boundary condition replicates the actual field conditions more accurately than the fixed boundary condition.In addition,it was found that the tailing embankment constructed by the downstream and upstream method of rising is not susceptible to liquefaction and lateral spreading for earthquake motions,even for a magnitude>5.5.

Effect of spin on the instability of THz plasma waves in field-effect transistors under non-ideal boundary conditions

Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.

Free Vibration Analysis of Rectangular Plate with Cutouts under Elastic Boundary Conditions in Independent Coordinate Coupling Method

Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.

One-and Multi-dimensional CWENOZ Reconstructions for Implementing Boundary Conditions Without Ghost Cells

We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.

Existence of Positive Solutions for Systems of Second-order Nonlinear Singular Differential Equations with Integral Boundary Conditions on Infinite Interval

By using cone theory and the Monch fixed theorem combined with a monotone iterative technique,we investigate the existence of positive solutions for systems of secondorder nonlinear singular differential equations with integral boundary conditions on infinite interval and establish the existence theorem of positive solutions and iterative sequence for approximating the positive solutions.The results in this paper improve some known results.

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS WITH NON-SEPARATED TYPE INTEGRAL BOUNDARY CONDITIONS

In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.

Application of inverse method to estimation of boundary conditions during investment casting simulation

Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient （IHTC） between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.

A novel high resolution model without open boundary conditions applied to the China Seas: first investigation on tides

We developed a Global Ocean Circulation and Tide Model （GOCTM） with coarse grids in the open deep ocean degrading ‘smoothly’ into the highly resolved China Seas （CS） of refined grids to study the tides and circulation there.GOCTM is based on the framework of the Finite Volume approach for better mass conservation through improved transports across the discrete individual control volume.It also takes a full advantage of the geometric flexibility of unstructured mesh using a realistic global topography including the Arctic Ocean.The CS are given a special focus by refining the unstructured grids,but they are embedded into global domain naturally.Furthermore,GOCTM not only successfully avoids the treatment of the open boundaries,but also optimizes the trade-off between computational cost and model accuracy.Meanwhile,GOCTM is driven by the astronomical tide-generating potential and the secondary tide-generating potential directly,together with the wind stress and heat flux.GOCTM succeeds in reproducing the global eight principal tidal harmonic constants.Particularly,the simulated tidal results in the CS are improved compared to some other regional models with the discrepancy of 3.9 cm for M 2 tide.This idea of GOCTM can also be referred for other regional ocean study.

Boundary conditions of the dynamicnumerical model for the Chinese mainland lithosphere

The numerical dynamic model of Chinese mainland lithosphere's stress and strain field was constructed withelasto-viscous creep constitive relation. The most recent data of the stress field of Chinese mainland and thehorizontal movement velocity of the crust blocks were used as constraint conditions. The values of the boundalsforce were computed by trial-and-error method. The effect of the Qinghai-Xizang Plateau's excess potential energyto the movement of Chinese mainland was studied also in this model. The results of the numerical computing showthat, recent rapid raising of the Qinghai-Xizang Plateau and the generation of normal faults in the southem part ofthe plateau resulted from the convergence of Indian Plate to the Eurasian Plate, and also from the rapid convectivethinning of the lower lithosphere. Horizontal extension of the Qinghai-Xizang Plateau is the main dynamic factorto torm the present tectonic framework of the Chinese mainland. The compressive loads on the eastern boundaryof the model were mainly applied by the subduction of the Pacific Plate. The compression from the PhilippinePlate maybe slight.

Transient Hydroelastic Response of VLFS by FEM with Impedance Boundary Conditions in Time Domain

A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastic responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.

GENERAL DECAY FOR A VISCOELASTIC EQUATION OF VARIABLE COEFFICIENTS WITH A TIME-VARYING DELAY IN THE BOUNDARY FEEDBACK AND ACOUSTIC BOUNDARY CONDITIONS

A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay results of the energy are established via suitable Lyapunov functionals and some properties of the convex functions. Our result is obtained without imposing any restrictive growth assumption on the damping term and the elements of the matrix A and the kernel function g.

MHD flow of nanofluids over an exponentially stretching sheet in a porous medium with convective boundary conditions

This article concentrates on the steady magnetohydrodynamic （MHD） flow of viscous nanofluid. The flow is caused by a permeable exponentially stretching surface. An incompressible fluid fills the porous space. A comparative study is made for the nanoparticles namely Copper （Cu）, Silver （Ag）, Alumina （A1203） and Titanium Oxide （TiO2）. Water is treated as a base fluid. Convective type boundary conditions are employed in modeling the heat transfer process. The non-linear partial differential equations governing the flow are reduced to an ordinary differential equation by similarity transformations. The obtained equations are then solved for the development of series solutions. Convergence of the obtained series solutions is explicitly discussed. The effects of different parameters on the velocity and temperature profiles are shown and analyzed through graphs.

Frequency Analysis of Multiple Layered Cylindrical Shells under Lateral Pressure with Asymmetric Boundary Conditions

Natural frequency characteristics of a thin-walled multiple layered cylindrical shell under lateral pressure are studied. The multiple layered cylindrical shell configuration is formed by three layers of isotropic material where the inner and outer layers are stainless steel and the middle layer is aluminum. The multiple layered shell equations with lateral pressure are established based on Love＇s shell theory. The governing equations of motion with lateral pressure are employed by using energy functional and applying the Ritz method. The boundary conditions represented by end conditions of the multiple layered cylindrical shell are simply supported-clamped（SS-C）, free-clamped（F-C） and simply supported-free（SS-F）. The influence of different lateral pressures, different thickness to radius ratios, different length to radius ratios and effect of the asymmetric boundary conditions on natural frequency characteristics are studied. It is shown that the lateral pressure has effect on the natural frequency of multiple layered cylindrical shell and causes the natural frequency to increase. The natural frequency of the developed multilayered cylindrical shell is validated by comparing with those in the literature. The proposed research provides an effective approach for vibration analysis shell structures subjected to lateral pressure with an energy method.

THE EXISTENCE OF NONTRIVIAL SOLUTIONS OF HAMILTONIAN SYSTEMS WITH LAGRANGIAN BOUNDARY CONDITIONS

Some theorems are obtained for the existence of nontrivial solutions of Hamiltonian systems with Lagrangian boundary conditions by the minimax methods.

On a Discrete Fractional Boundary Value Problem with Nonlocal Fractional Boundary Conditions

In this paper, we investigate the nonlinear fractional difference equation with nonlocal fractional boundary conditions. We derive the Green's function for this problem and show that it satisfies certain properties. Some existence results are obtained by means of nonlinear alternative of Leray-Schauder type theorem and Krasnosel-skii's fixed point theorem.

A Comparison Study of Deep Galerkin Method and Deep Ritz Method for Elliptic Problems with Different Boundary Conditions

Recent years have witnessed growing interests in solving partial differential equations by deep neural networks,especially in the high-dimensional case.Unlike classical numerical methods,such as finite difference method and finite element method,the enforcement of boundary conditions in deep neural networks is highly nontrivial.One general strategy is to use the penalty method.In the work,we conduct a comparison study for elliptic problems with four different boundary conditions,i.e.,Dirichlet,Neumann,Robin,and periodic boundary conditions,using two representative methods:deep Galerkin method and deep Ritz method.In the former,the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter.Therefore,it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions.However,by a number of examples,we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides,in some cases,when the boundary condition can be implemented in an exact manner,we find that such a strategy not only provides a better approximate solution but also facilitates the training process.

SPECTRAL PROPERTIES OF DISCRETE STURM-LIOUVILLE PROBLEMS WITH TWO SQUARED EIGENPARAMETER-DEPENDENT BOUNDARY CONDITIONS

In this article,we consider a discrete right-definite Sturm-Liouville problems with two squared eigenparameter-dependent boundary conditions.By constructing some new Lagrange-type identities and two fundamental functions,we obtain not only the existence,the simplicity,and the interlacing properties of the real eigenvalues,but also the oscillation properties,orthogonality of the eigenfunctions,and the expansion theorem.Finally,we also give a computation scheme for computing eigenvalues and eigenfunctions of specific eigenvalue problems.