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.

SUFFICIENT AND NECESSARY CONDITIONS ON THE EXISTENCE AND ESTIMATES OF BOUNDARY BLOW-UP SOLUTIONS FOR SINGULAR p-LAPLACIAN EQUATIONS

Let?denote a smooth,bounded domain in R^(N)(N≥2).Suppose that g is a nondecreasing C^(1)positive function and assume that b(x)is continuous and nonnegative inΩ,and that it may be singular on■Ω.In this paper,we provide sufficient and necessary conditions on the existence of boundary blow-up solutions to the p-Laplacian problem△_(p)u=b(x)g(u)for x∈Ω,u(x)→+∞as dist(x,■Ω)→0.The estimates of such solutions are also investigated.Moreover,when b has strong singularity,the nonexistence of boundary blow-up(radial)solutions and infinitely many radial solutions are also considered.

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.

Solution of the Third Kind Boundary Value Problem of Laplace’s Equation Based on Conformal Mapping

In order to overcome the difficulty in solving the boundary value problem of electrostatic field with complex boundary and to give a new method for solving the third boundary value problem of Laplace’s equation, in this paper, the third boundary value problem of Laplace’s equation is studied by combining conformal mapping with theoretical analysis, the several analytical solutions of third boundary value problems of Laplace’s equation are gives, the correctness of its solution is verified through computer numerical simulation, and a new idea and method for solving the third boundary value problem of Laplace’s equation is obtained. In this paper, the boundary condition of the solving domain is changed by the appropriate conformal mapping, so that the boundary value problem on the transformed domain is easy to be solved or be known, and then the third kind boundary value of the Laplace’s equation can be solved easily;its electric potential distribution is known. Furthermore, the electric field line and equipotential line are plotted by using the MATLAB software.

Integral Operator Solving Process of the Boundary Value Problem of Abstract Kinetic Equation with the First Kind of Critical Parameter and Generalized Periodic Boundary Conditions

在这份报纸，边界的概念与第一种批评参数珍视抽象运动方程的问题 0 并且概括周期的边界条件在由函数组成，向量在一个一般 Banach 空格珍视的一个 Lebesgue 空格被介绍，然后描述这些的解决方案由 Volterra 的抽象线性不可分的操作符的抽象边界值问题类型。我们把这个过程称为解决过程的不可分的操作员。

Effect of boundary slip on electroosmotic flow in a curved rectangular microchannel

The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.

Existence of Entropy Solution for Degenerate Parabolic-Hyperbolic Problem Involving p(x)-Laplacian with Neumann Boundary Condition

We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.

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.

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

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.

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.

The boundary conditions for simulations of a shake-table experiment on the seismic response of 3D slope

Boundary conditions can significantly affect a slope＇s behavior under strong earthquakes. To evaluate the importance of boundary conditions for finite element （FE） simulations of a shake-table experiment on the slope response, a validated three-dimensional （3D） nonlinear FE model is presented, and the numerical and experimental results are compared. For that purpose, the robust graphical user-interface ＂SlopeSAR＂, based on the open-source computational platform OpenSees, is employed, which simplifies the effort-intensive pre- and post-processing phases. The mesh resolution effect is also addressed. A parametric study is performed to evaluate the influence of boundary conditions on the FE model involving the boundary extent and three types of boundary conditions at the end faces. Generally, variations in the boundary extent produce inconsistent slope deformations. For the two end faces, fixing the y-direction displacement is not appropriate to simulate the shake-table experiment, in which the end walls are rigid and rough. In addition, the influence of the length of the 3D slope＇s top face and the width of the slope play an important role in the difference between two types of boundary conditions at the end faces （fixing the y-direction displacement and fixing the （y, z） direction displacement）. Overall, this study highlights that the assessment of a comparison between a simulation and an experimental result should be performed with due consideration to the effect of the boundary conditions.

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.

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 THE COEFFICIENTS OF DIFFERENTIATED EXPANSIONS AND DERIVATIVES OF CHEBYSHEV POLYNOMIALS OF THE THIRD AND FOURTH KINDS

Two new analytical formulae expressing explicitly the derivatives of Chebyshev polynomials of the third and fourth kinds of any degree and of any order in terms of Chebyshev polynomials of the third and fourth kinds themselves are proved. Two other explicit formulae which express the third and fourth kinds Chebyshev expansion coefficients of a general-order derivative of an infinitely differentiable function in terms of their original expansion coefficients are also given. Two new reduction formulae for summing some terminating hypergeometric functions of unit argument are deduced. As an application of how to use Chebyshev polynomials of the third and fourth kinds for solving high-order boundary value problems, two spectral Galerkin numerical solutions of a special linear twelfth-order boundary value problem are given.

Discussion on electromagnetic crack face boundary conditions for the fracture mechanics of magneto-electro-elastic materials

This paper discusses electromagnetic boundary conditions on crack faces in magneto- electroelastic materials, where piezoelectric, piezomagnetic and magnetoelectric effects are coupled. A notch of finite thickness in these materials is also addressed. Four idealized electromagnetic boundary conditions assumed for the crack-faces are separately investigated, i.e. （a） electrically and magnetically impermeable （crack-face）, （b） electrically impermeable and magnetically permeable, （c） electrically permeable and magnetically impermeable, and （d） electrically and magnetically permeable. The influence of the notch thickness on important parameters, such as the field intensity factors, the energy release rate at the notch tips and the electromagnetic fields inside the notch, are studied and the results are obtained in closed-form. Results under different idealized electromagnetic boundary conditions on the crack-face are compared, and the applicability of these idealized assumptions is discussed.

Investigation on the choice of boundary conditions and shape functions for flexible multi-body system

The objective of this investigation is to examine the correctness and efficiency of the choice of boundary conditions when using assumed mode approach to simulate flexible multi-body systems. The displacement field due to deformation is approximated by the Rayleigh-Ritz assumed modes in floating frame of reference （FFR） formulation. The deformations obtained by the absolute nodal coordinate （ANC） formulation which are transformed by two sets of reference coordinates are introduced as a criterion to verify the accuracy of the simulation results by using the FFR formulation. The relationship between the deformations obtained from different boundary conditions is revealed. Nu- merical simulation examples demonstrate that the assumed modes with cantilevered-free, simply-supported and free- free boundary conditions without inclusion of rigid body modes are suitable for simulation of flexible multi-body system with large over all motion, and the same physical deformation can be obtained using those mode functions, differ only by a coordinate transformation. It is also shown that when using mode shapes with statically indeterminate boundary conditions, significant error may occur. Furthermore, the slider crank mechanism with rigid crank is accurate enough for investigating boundary condition problem of flexible multi-body system, which cost significant less simulating time.

APPLICABILITY OF THE CRACK FACE ELECTRICAL BOUNDARY CONDITIONS IN PIEZOELECTRIC MECHANICS

The electrical boundary conditions on the crack faces and their applicability in piezoelectric materials are discussed. A slit crack and a notch of ?nite thickness in piezoelectric materials subjected to combined mechanical and electrical loads is considered. Here, a crack is de?ned as a notch without thickness, which is ?lled with air or vacuum. The crack or notch is perpendicular to the poling direction of the medium. The ideal crack face electrical boundary conditions, i.e., the electrically permeable crack and the electrically impermeable crack, are in- vestigated ?rst. Then dependence of the ?eld intensity factors on notch thickness at the notch tips is analyzed to obtain a closed-form. The results are compared with the ideal crack solutions. Some useful results are found.