In this paper, a new variational formulation for a reaction-diffusion problem in broken Sobolev space is proposed. And the new formulation in the broken Sobolev space will be proved that it is well-posed and equivalen...In this paper, a new variational formulation for a reaction-diffusion problem in broken Sobolev space is proposed. And the new formulation in the broken Sobolev space will be proved that it is well-posed and equivalent to the standard Galerkin variational formulation. The method will be helpful to easily solve the original partial differential equation numerically. And the method is novel and interesting, which can be used to deal with some complicated problem, such as the low regularity problem, the differential-integral problem and so on.展开更多
A mixed vuriational formulation for large deformation analysis of plates is introduced. In this formulation the equilibrium ami compatibility equations are satisfied identically by means of stress functions and displa...A mixed vuriational formulation for large deformation analysis of plates is introduced. In this formulation the equilibrium ami compatibility equations are satisfied identically by means of stress functions and displacement components, respectively, and the constitu,lye equations are satisfied in a least square sense. An example is solved and the results are compared with those available in the literature.Further, the functional is particularized for buckling analysis of plates and a simple example is solved to illustrate the theory.展开更多
In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave func...In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave functions together with plane wave functions to approximate the local properties of the field. We analyze the global convergence and give an error estimation of the method. Numerical examples are also presented to demonstrate the effectiveness of the strategy.展开更多
By introducing an image plane, the inverse heat conduction problem with free boundary is transformed into one with completely known boundaryt which is much simpler to handle.As a by-product, the classical Kirchhoff’s...By introducing an image plane, the inverse heat conduction problem with free boundary is transformed into one with completely known boundaryt which is much simpler to handle.As a by-product, the classical Kirchhoff’s transformation for accounting for variable conductivity is rederived and an invariance property of the inverse problem solution with respect to variable conductivity is indicated. Then a pair of complementary extremum principles are established on the image plane, providing a sound theoretical foundation for the Ritz’s method and finite element method (FEM).An example solved by FEM is also given.展开更多
The guided and leaky modes of a planar dielectric waveguide are eigensolutions of a singular Sturm-Liouville problem. The modes are the roots of a characteristicfunction which can be found with several methods that ha...The guided and leaky modes of a planar dielectric waveguide are eigensolutions of a singular Sturm-Liouville problem. The modes are the roots of a characteristicfunction which can be found with several methods that have been introduced in thepast. However, the evaluation of the characteristic function suffers from numericalinstabilities, and hence it is often difficult to find all modes in a given range. Here anew variational formulation is introduced, which, after discretization, leads either to aquadratic or a quartic eigenvalue problem. The modes can be computed with standardsoftware for polynomial eigenproblems. Numerical examples show that the method isnumerically stable and guarantees a complete set of solutions.展开更多
We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use fi...We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use first kind Bessel functions.We compare the performance of the two bases.Moreover,we show that it is possible to use coupled plane wave and Bessel bases in the same mesh.As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.展开更多
Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fre...Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.展开更多
Equations governing the vibrations and buckling of multilayered orthotropic graphene sheets can be expressed as a system of n partial differential equations where n refers to the number of sheets. This description is ...Equations governing the vibrations and buckling of multilayered orthotropic graphene sheets can be expressed as a system of n partial differential equations where n refers to the number of sheets. This description is based on the continuum model of the graphene sheets which can also take the small scale effects into account by employing a nonlocal theory. In the present article a variational principle is derived for the nonlocal elastic theory of rectangular graphene sheets embedded in an elastic medium and undergo- ing transverse vibrations. Moreover the graphene sheets are subject to biaxial compression. Rayleigh quotients are obtained for the frequencies of freely vibrating graphene sheets and for the buckling load. The influence of small scale effects on the frequencies and the buckling load can be observed qualiatively from the expressions of the Rayleigh quotients. Elastic medium is modeled as a combination of Winkler and Pasternak foundations acting on the top and bottom layers of the mutilayered nano-structure. Natural boundary con- ditions of the problem are derived using the variational principle formulated in the study. It is observed that free boundaries lead to coupled boundary conditions due to nonlocal theory used in the continuum formulation while the local （classical） elasticity theory leads to uncoupled boundary conditions. The mathematical methods used in the study involve calculus of variations and the semi-inverse method for deriving the variational integrals.展开更多
Sloshing-induced force and moment may affect the dynamic property of the liquid-contained system.Analytically presented linear Stokes-Joukowski potentials of fluid are usually needed for analytical study of sloshing i...Sloshing-induced force and moment may affect the dynamic property of the liquid-contained system.Analytically presented linear Stokes-Joukowski potentials of fluid are usually needed for analytical study of sloshing in liquid-filled tank under rotational(e.g.,pitching)excitations.To obtain the analytically approximate linear Stokes-Joukowski potentials of fluid in the rigid baffled tanks,a variational domain-decomposition scheme is proposed.This scheme includes three steps:(i)dividing the hydrostatic baffled fluid domain into simple sub-domains based on the positions of the baffles(i.e.,using the baffle as part of the boundaries of the sub-domain)by introducing artificial interfaces and densities of fluids in the different sub-domains or auxiliary normal fluid velocity functions on the artificial interfaces;(ii)expressing the solution for linear Stokes-Joukowski potential of each sub-domain as a linear combination of a class of harmonic functions with undetermined coefficients,and expressing the auxiliary normal fluid velocity functions on the artificial in terfaces as Fourier-type series with undetermined coefficients;(iii)solving the undetermined coefficients by the Trefftz method and the proposed variational formulations.The obtained semi-analytical linear Stokes-Joukowski potential agrees well with that published in literature or given by finite element method(FEM),and its applicability to study nonlinear sloshing problem is verified by applying it to a two-dimensional partially fluid-filled rectangular tank with a T-shaped baffle under pitching excitation.The present semi-analytical result is compared with that given by computational fluid dynamics(CFD)software or literature.展开更多
Vibration of a circular membrane in contact with a fluid has extensive applications in industry. The natural vibration frequencies for the asymmetric free vibra- tion of a circular membrane in contact with a bounded i...Vibration of a circular membrane in contact with a fluid has extensive applications in industry. The natural vibration frequencies for the asymmetric free vibra- tion of a circular membrane in contact with a bounded incompressible fluid are derived in this paper. Considering small oscillations induced by the membrane vibration in an incompressible and inviscid fluid, the velocity potential function is used to describe the fluid field. Two approaches are used to derive the free vibration frequencies of the sys- tem, which include a variational formulation and an approximate solution employing the Rayleigh quotient method. A good correlation is found between free vibration frequencies evaluated by these methods. Finally, the effects of the fluid depth, the mass density, and the radial tension on the free vibration frequencies of the coupled system are investigated.展开更多
On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualiza...On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, ＂mass＂ preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.展开更多
A family of one-dimensional(1D) elliptic boundary-value problems with periodic and rapidly-oscillating piecewise-smooth coefficients is considered. The coefficients depend on the local or fast variables corresponding ...A family of one-dimensional(1D) elliptic boundary-value problems with periodic and rapidly-oscillating piecewise-smooth coefficients is considered. The coefficients depend on the local or fast variables corresponding to two different structural scales. A finite number of imperfect contact conditions are analyzed at each of the scales. The reiterated homogenization method(RHM) is used to construct a formal asymptotic solution. The homogenized problem, the local problems, and the corresponding effective coefficients are obtained. A variational formulation is derived to obtain an estimate to prove the proximity between the solutions of the original problem and the homogenized problem. Numerical computations are used to illustrate both the convergence of the solutions and the gain of the effective properties of a three-scale heterogeneous 1D laminate with respect to their two-scale counterparts. The theoretical and practical ideas exposed here could be used to mathematically model multidimensional problems involving multiscale composite materials with imperfect contact at the interfaces.展开更多
In this paper,we propose and study neural network-based methods for solutions of high-dimensional quadratic porous medium equation(QPME).Three variational formulations of this nonlinear PDE are presented:a strong form...In this paper,we propose and study neural network-based methods for solutions of high-dimensional quadratic porous medium equation(QPME).Three variational formulations of this nonlinear PDE are presented:a strong formulation and two weak formulations.For the strong formulation,the solution is directly parameterized with a neural network and optimized by minimizing the PDEresidual.It can be proved that the convergence of the optimization problem guarantees the convergence of the approximate solution in the L^(1)sense.Theweak formulations are derived following(Brenier in Examples of hidden convexity in nonlinear PDEs,2020)which characterizes the very weak solutions of QPME.Specifically speaking,the solutions are represented with intermediate functions who are parameterized with neural networks and are trained to optimize the weak formulations.Extensive numerical tests are further carried out to investigate the pros and cons of each formulation in low and high dimensions.This is an initial exploration made along the line of solving high-dimensional nonlinear PDEs with neural network-based methods,which we hope can provide some useful experience for future investigations.展开更多
A novel wrist-inspired soft actuator,which is driven by a magneto-pneumatic hybrid system and based on a Kresling origami unit,is proposed.The geometric model,kinematic analysis model,and quasistatic analysis model of...A novel wrist-inspired soft actuator,which is driven by a magneto-pneumatic hybrid system and based on a Kresling origami unit,is proposed.The geometric model,kinematic analysis model,and quasistatic analysis model of the Kresling origami unit are presented.A key focus is on the formulation and investigation of the variation in rotation angle using the kinematic analysis model.A wrist-inspired soft actuator is designed,and its quasistatic characteristics are validated through various experiments.The paper proposes an innovative magneto-pneumatic hybrid actuation method,capable of achieving bidirectional torsion.This actuation method is experimentally validated,demonstrating the actuator's ability to maintain 3 steady states and its capability for bidirectional torsion deformation.Furthermore,the paper highlights the potential of the Kresling origami unit in designing soft actuators capable of achieving large rotation angles.For instance,an actuator with 6 sides(n=6)is shown to achieve a rotation angle of 239.5°,and its rotation ratio exceeds 277°,about twice the largest one reported in other literature.The actuator offers a practical and effective solution for bidirectional torsion deformation in soft robotic applications.展开更多
This paper presents a critical evaluation of the physical aspects of lift generation to prove that no lift can be generated in a steady inviscid flow.Hence,the answer to the recurring question in the paper title is ne...This paper presents a critical evaluation of the physical aspects of lift generation to prove that no lift can be generated in a steady inviscid flow.Hence,the answer to the recurring question in the paper title is negative.In other words,the fluid viscosity is necessary in lift generation.The relevant topics include D’Alembert’s paradox of lift and drag,the Kutta condition,the force expression based on the boundary enstrophy flux(BEF),the vortex lift,and the generation of the vorticity and circulation.The physi-cal meanings of the variational formulations to determine the circulation and lift are discussed.In particular,in the variational formulation based on the continuity equation with the first-order Tikhonov regularization functional,an incompressible flow with the artificial viscosity(the Lagrange multiplier)is simulated,elucidating the role of the artifi-cial viscosity in lift generation.The presented contents are valuable for the pedagogical purposes in aerodynamics and fluid mechanics.展开更多
A finite element model for piezoelectric smart beam in extension mode based on First-order Shear Deformation Theory(FSDT)with an appropriate through-thickness distribution of electric potential is presented.Accuracy o...A finite element model for piezoelectric smart beam in extension mode based on First-order Shear Deformation Theory(FSDT)with an appropriate through-thickness distribution of electric potential is presented.Accuracy of piezoelectric finite element formulations depends on the selection of assumed mechanical and electrical fields.Most of the conventional FSDT-based piezoelectric beam formulations available in the literature use linear through-thickness distribution of electric potential which is actually nonlinear.Here,a novel quadratic profile of the through-thickness electric potential is proposed to include the nonlinear effects.The results obtained show that the accuracy of conventional formulations with linear through-thickness potential approximation is affected by the material configuration,especially when the piezoelectric material dominates the beam cross section.It is shown that the present formulation gives the same level of accuracy for all regimes of material configurations in the beam cross section.Also,a modified form of the FSDT displacements is employed,which utilizes the shear angle as a degree of freedom instead of section rotation.Such a FSDT displacement field shows improved performance compared to the conventional field.The present formulation is validated by comparing the results with ANSYS 2D simulation.The comparison of results proves the improved efficiency and accuracy of the present formulation over the conventional formulations.展开更多
A class of steady-state metal-forming problems,with rigid-plastic,incompressible,strain-rate dependent material model and nonlocal Coulomb’s friction,is considered.Primal,mixed and penalty variational formulations,co...A class of steady-state metal-forming problems,with rigid-plastic,incompressible,strain-rate dependent material model and nonlocal Coulomb’s friction,is considered.Primal,mixed and penalty variational formulations,containing variational inequalities with nonlinear and nondifferentiable terms,are derived and studied.Existence,uniqueness and convergence results are obtained and shortly presented.A priori finite element error estimates are derived and an algorithm,combining the finite element and secant-modulus methods,is utilized to solve an illustrative extrusion problem.展开更多
This paper addresses the modeling of fracture in quasi-brittle materials using a phase-field approach to the description of crack topol-ogy.Within the computational mechanics community,several studies have treated the...This paper addresses the modeling of fracture in quasi-brittle materials using a phase-field approach to the description of crack topol-ogy.Within the computational mechanics community,several studies have treated the issue of modeling fracture using phase fields.Most of these studies have used an approach that implies the lack of a damage threshold.We herein explore an alternative model that includes a damage threshold and study how it compares with the most popular approach.The formulation is systematically explained within a rigorous variational framework.Subsequently,we present the corresponding three-dimensional finite element discretization that leads to a straightforward numerical implementation.Benchmark simulations in two dimensions and three dimensions are then presented.The results show that while an elastic stage and a damage threshold are ensured by the present model,good agreement with the results reported in the literature can be obtained,where such features are generally absent.展开更多
In this paper,we study the numerical solution of the Stokes system in deformed axisymmetric geometries.In the azimuthal direction the discretization is carried out by using truncated Fourier series,thus reducing the d...In this paper,we study the numerical solution of the Stokes system in deformed axisymmetric geometries.In the azimuthal direction the discretization is carried out by using truncated Fourier series,thus reducing the dimension of the problem.The resulting two-dimensional problems are discretized using the spectral element method which is based on the variational formulation in primitive variables.The meridian domain is subdivided into elements,in each of which the solution is approximated by truncated polynomial series.The results of numerical experiments for several geometries are presented.展开更多
An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of...An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of the governing Helmholtz equation.In this paper we are concerned with fast solver for the system generated by the method in[1].We propose a new preconditioner for such system,which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in[13].In our numerical experiments,this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations,and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.展开更多
基金Supported by the Natural Science Foundation of Henan Province（162300410031） Supported by the Excellent Youth Program of the Basic Research Operating Expenses Program of Henan Province （yqpy20140039）
文摘In this paper, a new variational formulation for a reaction-diffusion problem in broken Sobolev space is proposed. And the new formulation in the broken Sobolev space will be proved that it is well-posed and equivalent to the standard Galerkin variational formulation. The method will be helpful to easily solve the original partial differential equation numerically. And the method is novel and interesting, which can be used to deal with some complicated problem, such as the low regularity problem, the differential-integral problem and so on.
文摘A mixed vuriational formulation for large deformation analysis of plates is introduced. In this formulation the equilibrium ami compatibility equations are satisfied identically by means of stress functions and displacement components, respectively, and the constitu,lye equations are satisfied in a least square sense. An example is solved and the results are compared with those available in the literature.Further, the functional is particularized for buckling analysis of plates and a simple example is solved to illustrate the theory.
基金The authors would like to thank the reviewers and Dr.Zheng Enxi for many valuable suggcstions. This work is supported by the National Natural Science Foundation of China （Grant No. 11371172, 51178001）, Science and technology research project of the education department of Jilin Province （Grant No. 2014213）.
文摘In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave functions together with plane wave functions to approximate the local properties of the field. We analyze the global convergence and give an error estimation of the method. Numerical examples are also presented to demonstrate the effectiveness of the strategy.
文摘By introducing an image plane, the inverse heat conduction problem with free boundary is transformed into one with completely known boundaryt which is much simpler to handle.As a by-product, the classical Kirchhoff’s transformation for accounting for variable conductivity is rederived and an invariance property of the inverse problem solution with respect to variable conductivity is indicated. Then a pair of complementary extremum principles are established on the image plane, providing a sound theoretical foundation for the Ritz’s method and finite element method (FEM).An example solved by FEM is also given.
文摘The guided and leaky modes of a planar dielectric waveguide are eigensolutions of a singular Sturm-Liouville problem. The modes are the roots of a characteristicfunction which can be found with several methods that have been introduced in thepast. However, the evaluation of the characteristic function suffers from numericalinstabilities, and hence it is often difficult to find all modes in a given range. Here anew variational formulation is introduced, which, after discretization, leads either to aquadratic or a quartic eigenvalue problem. The modes can be computed with standardsoftware for polynomial eigenproblems. Numerical examples show that the method isnumerically stable and guarantees a complete set of solutions.
文摘We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use first kind Bessel functions.We compare the performance of the two bases.Moreover,we show that it is possible to use coupled plane wave and Bessel bases in the same mesh.As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.
文摘Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.
基金supported by research grants from the University of KwaZulu-Natal (UKZN)National Research Foundation (NRF) of South Africa
文摘Equations governing the vibrations and buckling of multilayered orthotropic graphene sheets can be expressed as a system of n partial differential equations where n refers to the number of sheets. This description is based on the continuum model of the graphene sheets which can also take the small scale effects into account by employing a nonlocal theory. In the present article a variational principle is derived for the nonlocal elastic theory of rectangular graphene sheets embedded in an elastic medium and undergo- ing transverse vibrations. Moreover the graphene sheets are subject to biaxial compression. Rayleigh quotients are obtained for the frequencies of freely vibrating graphene sheets and for the buckling load. The influence of small scale effects on the frequencies and the buckling load can be observed qualiatively from the expressions of the Rayleigh quotients. Elastic medium is modeled as a combination of Winkler and Pasternak foundations acting on the top and bottom layers of the mutilayered nano-structure. Natural boundary con- ditions of the problem are derived using the variational principle formulated in the study. It is observed that free boundaries lead to coupled boundary conditions due to nonlocal theory used in the continuum formulation while the local （classical） elasticity theory leads to uncoupled boundary conditions. The mathematical methods used in the study involve calculus of variations and the semi-inverse method for deriving the variational integrals.
基金the National Natural Science Foundation of China(Grant Nos.11572018 and 11772020).
文摘Sloshing-induced force and moment may affect the dynamic property of the liquid-contained system.Analytically presented linear Stokes-Joukowski potentials of fluid are usually needed for analytical study of sloshing in liquid-filled tank under rotational(e.g.,pitching)excitations.To obtain the analytically approximate linear Stokes-Joukowski potentials of fluid in the rigid baffled tanks,a variational domain-decomposition scheme is proposed.This scheme includes three steps:(i)dividing the hydrostatic baffled fluid domain into simple sub-domains based on the positions of the baffles(i.e.,using the baffle as part of the boundaries of the sub-domain)by introducing artificial interfaces and densities of fluids in the different sub-domains or auxiliary normal fluid velocity functions on the artificial interfaces;(ii)expressing the solution for linear Stokes-Joukowski potential of each sub-domain as a linear combination of a class of harmonic functions with undetermined coefficients,and expressing the auxiliary normal fluid velocity functions on the artificial in terfaces as Fourier-type series with undetermined coefficients;(iii)solving the undetermined coefficients by the Trefftz method and the proposed variational formulations.The obtained semi-analytical linear Stokes-Joukowski potential agrees well with that published in literature or given by finite element method(FEM),and its applicability to study nonlinear sloshing problem is verified by applying it to a two-dimensional partially fluid-filled rectangular tank with a T-shaped baffle under pitching excitation.The present semi-analytical result is compared with that given by computational fluid dynamics(CFD)software or literature.
文摘Vibration of a circular membrane in contact with a fluid has extensive applications in industry. The natural vibration frequencies for the asymmetric free vibra- tion of a circular membrane in contact with a bounded incompressible fluid are derived in this paper. Considering small oscillations induced by the membrane vibration in an incompressible and inviscid fluid, the velocity potential function is used to describe the fluid field. Two approaches are used to derive the free vibration frequencies of the sys- tem, which include a variational formulation and an approximate solution employing the Rayleigh quotient method. A good correlation is found between free vibration frequencies evaluated by these methods. Finally, the effects of the fluid depth, the mass density, and the radial tension on the free vibration frequencies of the coupled system are investigated.
文摘On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, ＂mass＂ preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.
基金Project supported by the Desenvolvimento e Aplicaoes de Mtodos Matemticos de Homogeneizaao(CAPES)(No.88881.030424/2013-01)the Homogeneizao Reiterada Aplicada a Meios Dependentes de Múltiplas Escalas con Contato Imperfeito Entre as Fases(CNPq)(Nos.450892/2016-6and 303208/2014-7)the Caracterizacin de Propiedades Efectivas de Tejidos Biolgicos Sanos y Cancerosos(CONACYT)(No.2016–01–3212)
文摘A family of one-dimensional(1D) elliptic boundary-value problems with periodic and rapidly-oscillating piecewise-smooth coefficients is considered. The coefficients depend on the local or fast variables corresponding to two different structural scales. A finite number of imperfect contact conditions are analyzed at each of the scales. The reiterated homogenization method(RHM) is used to construct a formal asymptotic solution. The homogenized problem, the local problems, and the corresponding effective coefficients are obtained. A variational formulation is derived to obtain an estimate to prove the proximity between the solutions of the original problem and the homogenized problem. Numerical computations are used to illustrate both the convergence of the solutions and the gain of the effective properties of a three-scale heterogeneous 1D laminate with respect to their two-scale counterparts. The theoretical and practical ideas exposed here could be used to mathematically model multidimensional problems involving multiscale composite materials with imperfect contact at the interfaces.
基金supported in part by National Science Foundation via grant DMS-2012286by Department of Energy via grant DE-SC0019449.
文摘In this paper,we propose and study neural network-based methods for solutions of high-dimensional quadratic porous medium equation(QPME).Three variational formulations of this nonlinear PDE are presented:a strong formulation and two weak formulations.For the strong formulation,the solution is directly parameterized with a neural network and optimized by minimizing the PDEresidual.It can be proved that the convergence of the optimization problem guarantees the convergence of the approximate solution in the L^(1)sense.Theweak formulations are derived following(Brenier in Examples of hidden convexity in nonlinear PDEs,2020)which characterizes the very weak solutions of QPME.Specifically speaking,the solutions are represented with intermediate functions who are parameterized with neural networks and are trained to optimize the weak formulations.Extensive numerical tests are further carried out to investigate the pros and cons of each formulation in low and high dimensions.This is an initial exploration made along the line of solving high-dimensional nonlinear PDEs with neural network-based methods,which we hope can provide some useful experience for future investigations.
基金supported by the National Natural Science Foundation of China(Grant Nos.91748209 and 11402229)the specialized research projects of Huanjiang Laboratory.
文摘A novel wrist-inspired soft actuator,which is driven by a magneto-pneumatic hybrid system and based on a Kresling origami unit,is proposed.The geometric model,kinematic analysis model,and quasistatic analysis model of the Kresling origami unit are presented.A key focus is on the formulation and investigation of the variation in rotation angle using the kinematic analysis model.A wrist-inspired soft actuator is designed,and its quasistatic characteristics are validated through various experiments.The paper proposes an innovative magneto-pneumatic hybrid actuation method,capable of achieving bidirectional torsion.This actuation method is experimentally validated,demonstrating the actuator's ability to maintain 3 steady states and its capability for bidirectional torsion deformation.Furthermore,the paper highlights the potential of the Kresling origami unit in designing soft actuators capable of achieving large rotation angles.For instance,an actuator with 6 sides(n=6)is shown to achieve a rotation angle of 239.5°,and its rotation ratio exceeds 277°,about twice the largest one reported in other literature.The actuator offers a practical and effective solution for bidirectional torsion deformation in soft robotic applications.
文摘This paper presents a critical evaluation of the physical aspects of lift generation to prove that no lift can be generated in a steady inviscid flow.Hence,the answer to the recurring question in the paper title is negative.In other words,the fluid viscosity is necessary in lift generation.The relevant topics include D’Alembert’s paradox of lift and drag,the Kutta condition,the force expression based on the boundary enstrophy flux(BEF),the vortex lift,and the generation of the vorticity and circulation.The physi-cal meanings of the variational formulations to determine the circulation and lift are discussed.In particular,in the variational formulation based on the continuity equation with the first-order Tikhonov regularization functional,an incompressible flow with the artificial viscosity(the Lagrange multiplier)is simulated,elucidating the role of the artifi-cial viscosity in lift generation.The presented contents are valuable for the pedagogical purposes in aerodynamics and fluid mechanics.
文摘A finite element model for piezoelectric smart beam in extension mode based on First-order Shear Deformation Theory(FSDT)with an appropriate through-thickness distribution of electric potential is presented.Accuracy of piezoelectric finite element formulations depends on the selection of assumed mechanical and electrical fields.Most of the conventional FSDT-based piezoelectric beam formulations available in the literature use linear through-thickness distribution of electric potential which is actually nonlinear.Here,a novel quadratic profile of the through-thickness electric potential is proposed to include the nonlinear effects.The results obtained show that the accuracy of conventional formulations with linear through-thickness potential approximation is affected by the material configuration,especially when the piezoelectric material dominates the beam cross section.It is shown that the present formulation gives the same level of accuracy for all regimes of material configurations in the beam cross section.Also,a modified form of the FSDT displacements is employed,which utilizes the shear angle as a degree of freedom instead of section rotation.Such a FSDT displacement field shows improved performance compared to the conventional field.The present formulation is validated by comparing the results with ANSYS 2D simulation.The comparison of results proves the improved efficiency and accuracy of the present formulation over the conventional formulations.
文摘A class of steady-state metal-forming problems,with rigid-plastic,incompressible,strain-rate dependent material model and nonlocal Coulomb’s friction,is considered.Primal,mixed and penalty variational formulations,containing variational inequalities with nonlinear and nondifferentiable terms,are derived and studied.Existence,uniqueness and convergence results are obtained and shortly presented.A priori finite element error estimates are derived and an algorithm,combining the finite element and secant-modulus methods,is utilized to solve an illustrative extrusion problem.
文摘This paper addresses the modeling of fracture in quasi-brittle materials using a phase-field approach to the description of crack topol-ogy.Within the computational mechanics community,several studies have treated the issue of modeling fracture using phase fields.Most of these studies have used an approach that implies the lack of a damage threshold.We herein explore an alternative model that includes a damage threshold and study how it compares with the most popular approach.The formulation is systematically explained within a rigorous variational framework.Subsequently,we present the corresponding three-dimensional finite element discretization that leads to a straightforward numerical implementation.Benchmark simulations in two dimensions and three dimensions are then presented.The results show that while an elastic stage and a damage threshold are ensured by the present model,good agreement with the results reported in the literature can be obtained,where such features are generally absent.
文摘In this paper,we study the numerical solution of the Stokes system in deformed axisymmetric geometries.In the azimuthal direction the discretization is carried out by using truncated Fourier series,thus reducing the dimension of the problem.The resulting two-dimensional problems are discretized using the spectral element method which is based on the variational formulation in primitive variables.The meridian domain is subdivided into elements,in each of which the solution is approximated by truncated polynomial series.The results of numerical experiments for several geometries are presented.
基金The second author was supported by the Major Research Plan of Natural Science Foundation of China G91130015the Key Project of Natural Science Foundation of China G11031006National Basic Research Program of China G2011309702.
文摘An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of the governing Helmholtz equation.In this paper we are concerned with fast solver for the system generated by the method in[1].We propose a new preconditioner for such system,which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in[13].In our numerical experiments,this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations,and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.
