A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of ...A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of revolutions. All of the interface integral items participating in forming the potential and complementary energy functionals of the composite materials are expressed in terms of intrinsic quantities of the ellipsoidal of revolutions. Based on this, the upper and the lower bound for the effective elastic moduli of the composite materials with inclusions described above have been derived. Under three limiting conditions of sphere, disk and needle shaped inclusions, the results of this paper will return to the bounds obtained by Hashin([6]) (1992).展开更多
This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method (EFGM). The EFGM is employed as an alternative scheme to numerically so...This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method (EFGM). The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization (SIMP) scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes (MMA) belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.展开更多
This paper describes geometrical essentials of some iteration methods (e.g. Newton iteration, secant line method, etc.) for solving nonlinear equations and advances some geometrical methods of iteration that are fle...This paper describes geometrical essentials of some iteration methods (e.g. Newton iteration, secant line method, etc.) for solving nonlinear equations and advances some geometrical methods of iteration that are flexible and efficient.展开更多
In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical featu...In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.展开更多
A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) n...A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) not only the stress fields of a circular disk can be displayed, but the u -and v -displacement fields can also be produced.展开更多
In this paper, a new statistical averaging technique is proposed for finding an optimal solution to a multi-objective linear fractional programming problem (MOLFPP) and multi-objective linear programming problem (MOLP...In this paper, a new statistical averaging technique is proposed for finding an optimal solution to a multi-objective linear fractional programming problem (MOLFPP) and multi-objective linear programming problem (MOLPP) by using new arithmetic averaging method and new geometric averaging method. It is significantly noticeable same characteristics among all the technique while taking maximum or minimum among all optimized values for multi-objective functions using simplex algorithm. The characteristics provided from the problems are verified by the numerical examples.展开更多
The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consis...The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.展开更多
To balance inventory cost with diverse demand,an optimal investment decision on necessary process improvement for delayed product differentiation is studied. A two-stage flexible manufacturing system is modeled as a c...To balance inventory cost with diverse demand,an optimal investment decision on necessary process improvement for delayed product differentiation is studied. A two-stage flexible manufacturing system is modeled as a continuous time Markov chain. The first production stage manufactures semifinished products based on a make-to-stock policy. The second production stage customizes semi-finished products from the first production stage on a make-to-order policy. Various performance measures for this flexible manufacturing system are evaluated by using matrix geometric methods. An optimization model to determine the level of investment on process improvement that minimizes the manufacturer ’s total cost is established. The results show that,a higher investment level can reduce both the expected customer order fulfillment delay and the expected semi-finished products inventory. When the initial order penetration point is 0. 4,the manufacturer ’s total cost is reduced by 15. 89% through process investment. In addition, the optimal investment level increases with the increase in the unit time cost of customer order fulfillment delay,and decreases with the increase in the product value and the initial order penetration point.展开更多
The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rat...The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.展开更多
A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilater...A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method (QAC) based membrane element AGQ6- II, and a Timoshenko beam function (TBF) method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory (FSDT), for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.展开更多
The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite el...The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite element model for two-dimensional dynamic analysis of armored cable is presented. This model accounts for the geometric nonlinearities of large displacement of the armored cable, and effects of axial load and bending stiffness. The governing equations are derived by consistent linearization and finite element discretization of the total weak form of the armored cable system, and solved by the Newmark time integration method. To make the solution procedure avoid falling into the local extreme points, a simple adaptive stepping strategy is proposed. The presented model is validated via actual measured data. Results for dynamic configurations, motion and tension of both ends of the armored cable, and resonance-zone are presented for two numerical cases, including the dynamic analysis under the case of only ship heave motion and the case of joint action of ship heave motion and ocean current. The dynamics analysis can provide important reference for the design or product selection of the armored cable in a deep-sea ROV system so as to improve the safety of its marine operation under the sea state of 4 or above.展开更多
Multi-pass slab vertical-horizontal (V-H) rolling process with variable edging roll shape have been simulated with explicit dynamic finite element method and updating geometric method. The distributions of plastic s...Multi-pass slab vertical-horizontal (V-H) rolling process with variable edging roll shape have been simulated with explicit dynamic finite element method and updating geometric method. The distributions of plastic strain contour in slab daring rolling process with different edging roll and under different rolling stage have been obtained. The results show that there exist two thin strain assembling zones in slab when the flat edging roll is used, and there just exist one strain assembling zone in slab when the edging roll with groove is used. And compared the deformation equality between flat edging roll and edging roll with groove, the lateris better than the former, which supplies the theory prove to the slab deformation distribution during V-H rolling process and is helpful for predicting the slab texture.展开更多
We extend the theory of global geometrical optics method, proposed originally for the linear scalar high-frequency wave-like equations in [Commun. Math. Sci., 2013, 11(1): 105-140], to the more general vector- valu...We extend the theory of global geometrical optics method, proposed originally for the linear scalar high-frequency wave-like equations in [Commun. Math. Sci., 2013, 11(1): 105-140], to the more general vector- valued Schrodinger problems in the semi-classical regime. The key ingredient in the global geometrical optics method is a moving frame technique in the phase space. The governing equation is transformed into a new equation but of the same type when expressed in any moving frame induced by the underlying Hamiltonian flow. The classical Wentzel-Kramers-Brillouin (WKB) analysis benefits from this treatment as it maintains valid for arbitrary but fixed evolutionary time. It turns out that a WKB-type function defined merely on the underlying Lagrangian submanifold can be obtained with the help of this moving frame technique, and from which a uniform first-order approximation of the wave field can be derived, even around caustics. The general theory is exemplified by two specific instances. One is the two-level SchrSdinger system and the other is the periodic SchrSdinger equation. Numerical tests validate the theoretical results.展开更多
PL homotopy metheds are effective methods to locate zerces(or fixed points) of highly nonlinearmappirgs. Due to the Jexicographical system, the methods are feasible without exceptions Thispaper presents a geemetrical ...PL homotopy metheds are effective methods to locate zerces(or fixed points) of highly nonlinearmappirgs. Due to the Jexicographical system, the methods are feasible without exceptions Thispaper presents a geemetrical interpretation of the without-exception feasibility.展开更多
The objective of this article is to introduce a generalized algorithm to produce the m-point n-ary approximating subdivision schemes(for any integer m, n ≥ 2). The proposed algorithm has been derived from uniform B-s...The objective of this article is to introduce a generalized algorithm to produce the m-point n-ary approximating subdivision schemes(for any integer m, n ≥ 2). The proposed algorithm has been derived from uniform B-spline blending functions. In particular, we study statistical and geometrical/traditional methods for the model selection and assessment for selecting a subdivision curve from the proposed family of schemes to model noisy and noisy free data. Moreover, we also discuss the deviation of subdivision curves generated by proposed family of schemes from convex polygonal curve. Furthermore, visual performances of the schemes have been presented to compare numerically the Gibbs oscillations with the existing family of schemes.展开更多
Capacitance extraction is one of the key issues in integrated circuits and also a typical electrostatic prob- lem. The dual discrete geometric method (DGM) is investigated to provide relative solutions in two-dimens...Capacitance extraction is one of the key issues in integrated circuits and also a typical electrostatic prob- lem. The dual discrete geometric method (DGM) is investigated to provide relative solutions in two-dimensional unstructured mesh space. The energy complementary characteristic and quick field energy computation thereof based on it are emphasized. Contrastive analysis between the dual finite element methods and the dual DGMs are presented both from theoretical derivation and through case studies. The DGM, taking the scalar potential as unknown on dual interlocked meshes, with simple form and good accuracy, is expected to be one of the mainstreaming methods in associated areas.展开更多
Sensitivity analysis methods help to deal with the challenges of process variation in extraction of para- sitic capacitances in an integrated circuit. The dual discrete geometric methods (DGMs), which have been rece...Sensitivity analysis methods help to deal with the challenges of process variation in extraction of para- sitic capacitances in an integrated circuit. The dual discrete geometric methods (DGMs), which have been recently utilized to extract parasitic capacitances, are reviewed. The computation method based on the dual DGMs for sen- sitivities of capacitances with respect to the given process parameters is presented. As the dual DGMs utilize scalar electric potential is unknown, the capacitances are obtained effectively, and then the sensitivities are calculated conveniently.展开更多
In this paper,we consider a GI/M/1 queue operating in a multi-phase service environment with working vacations and Bernoulli vacation interruption.Whenever the queue becomes empty,the server begins a working vacation ...In this paper,we consider a GI/M/1 queue operating in a multi-phase service environment with working vacations and Bernoulli vacation interruption.Whenever the queue becomes empty,the server begins a working vacation of random length,causing the system to move to vacation phase 0.During phase 0,the server takes service for the customers at a lower rate rather than stopping completely.When a vacation ends,if the queue is non-empty,the system switches from the phase 0 to some normal service phase i with probability qi,i=1,2,⋯,N.Moreover,we assume Bernoulli vacation interruption can happen.At a service completion instant,if there are customers in a working vacation period,vacation interruption happens with probability p,then the system switches from the phase 0 to some normal service phase i with probability qi,i=1,2,⋯,N,or the server continues the vacation with probability 1−p.Using the matrix geometric solution method,we obtain the stationary distributions for queue length at both arrival epochs and arbitrary epochs.The waiting time of an arbitrary customer is also derived.Finally,several numerical examples are presented.展开更多
The unforeseen mobile data explosion as well as the scarce of spectrum resource pose a major challenge to the performance of today's cellular networks which are in urgent need of novel solutions to handle such volumi...The unforeseen mobile data explosion as well as the scarce of spectrum resource pose a major challenge to the performance of today's cellular networks which are in urgent need of novel solutions to handle such voluminous mobile data. Long term evolution-unlicensed (LTE-U), which extends the LTE standard operating on the unlicensed band, has been proposed to improve system throughput. In LTE-U system, arriving users will contend the unlicensed spectrum resource with wireless fidelity (WiFi) users to transmit data information. Nevertheless, there is no clear consensus as to the benefits of transmission using unlicensed bands for LTE users. To this end, in this paper an analytical model is presented based on a queue system to understand the performance achieved by unlicensed based LTE system taking quality of services (QoS) and LTE-U users' behaviors into account. To obtain the stead-state solutions of the queue system, a matrix geometric method is used to solve it. Then, the average delay and utilization of unlicensed band for the LTE-U users is derived by using the queuing model. The performance of LTE-U coexistence is evaluated with WiFi using the proposed model and provide some initial insights as to the advantage of LTE-U in practice.展开更多
In this paper,we study the Cauchy problem of an integrable evolution system,i.e.,the n-dimensional generalization of third-order symmetry of the well-known Landau-Lifshitz equation.By rewriting this equation in a geom...In this paper,we study the Cauchy problem of an integrable evolution system,i.e.,the n-dimensional generalization of third-order symmetry of the well-known Landau-Lifshitz equation.By rewriting this equation in a geometric form and applying the geometric energy method with a forth-order perturbation,we show the global well-posedness of the Cauchy problem in suitable Sobolev spaces.展开更多
文摘A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of revolutions. All of the interface integral items participating in forming the potential and complementary energy functionals of the composite materials are expressed in terms of intrinsic quantities of the ellipsoidal of revolutions. Based on this, the upper and the lower bound for the effective elastic moduli of the composite materials with inclusions described above have been derived. Under three limiting conditions of sphere, disk and needle shaped inclusions, the results of this paper will return to the bounds obtained by Hashin([6]) (1992).
基金the National '973' Key Fundamental Research Projects of China(No.2003CB716207)the National '863' High-Tech Development Projects of China(No.2006AA04Z162)also the Australian Research Council(No.ARC-DP0666683).
文摘This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method (EFGM). The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization (SIMP) scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes (MMA) belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.
文摘This paper describes geometrical essentials of some iteration methods (e.g. Newton iteration, secant line method, etc.) for solving nonlinear equations and advances some geometrical methods of iteration that are flexible and efficient.
基金This project is supported by Provincial Project Foundation of Science and Technology of Guangdong, China(No.2002104040101).
文摘In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.
文摘A geometrical moir method for displaying directly isochromatics and isoclinics of a diametrically pressed circular disk is presented. It is demonstrated that by using two identical or different Wulff nets (or grids) not only the stress fields of a circular disk can be displayed, but the u -and v -displacement fields can also be produced.
文摘In this paper, a new statistical averaging technique is proposed for finding an optimal solution to a multi-objective linear fractional programming problem (MOLFPP) and multi-objective linear programming problem (MOLPP) by using new arithmetic averaging method and new geometric averaging method. It is significantly noticeable same characteristics among all the technique while taking maximum or minimum among all optimized values for multi-objective functions using simplex algorithm. The characteristics provided from the problems are verified by the numerical examples.
基金supported by the Natural Science Foundation of Hubei Province(CN)(Grant No.2019CFB693)the Research Foundation of the Education Department of Hubei Province(CN)(Grant No.B2019003)the open Foundation of the Key Laboratory of Metallurgical Equipment and Control of Education Ministry(CN)(Grant No.2015B14).
文摘The isogeometric analysis method(IGA)is a new type of numerical method solving partial differential equations.Compared with the traditional finite element method,IGA based on geometric spline can keep the model consistency between geometry and analysis,and provide higher precision with less freedom.However,huge stiffness matrix fromthe subdivision progress still leads to the solution efficiency problems.This paper presents amultigrid method based on geometric multigrid(GMG)to solve the matrix system of IGA.This method extracts the required computational data for multigrid method fromthe IGA process,which also can be used to improve the traditional algebraic multigrid method(AGM).Based on this,a full multigrid method(FMG)based on GMG is proposed.In order to verify the validity and reliability of these methods,this paper did some test on Poisson’s equation and Reynolds’equation and compared the methods on different subdivision methods,different grid degrees of freedom,different cyclic structure degrees,and studied the convergence rate under different subdivision strategies.The results show that the proposed method is superior to the conventional algebraic multigrid method,and for the standard relaxed V-cycle iteration,the method still has a convergence speed independent of the grid size at the same degrees.
基金The National Natural Science Foundation of China(No.71661147004)
文摘To balance inventory cost with diverse demand,an optimal investment decision on necessary process improvement for delayed product differentiation is studied. A two-stage flexible manufacturing system is modeled as a continuous time Markov chain. The first production stage manufactures semifinished products based on a make-to-stock policy. The second production stage customizes semi-finished products from the first production stage on a make-to-order policy. Various performance measures for this flexible manufacturing system are evaluated by using matrix geometric methods. An optimization model to determine the level of investment on process improvement that minimizes the manufacturer ’s total cost is established. The results show that,a higher investment level can reduce both the expected customer order fulfillment delay and the expected semi-finished products inventory. When the initial order penetration point is 0. 4,the manufacturer ’s total cost is reduced by 15. 89% through process investment. In addition, the optimal investment level increases with the increase in the unit time cost of customer order fulfillment delay,and decreases with the increase in the product value and the initial order penetration point.
基金the National Natural Science Foundation of China (Nos. 11472041,11532002,11772049,and 11802320)。
文摘The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.
文摘A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method (QAC) based membrane element AGQ6- II, and a Timoshenko beam function (TBF) method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory (FSDT), for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.
基金Project(2008AA09Z201)supported by the National High Technology Research and Development Program of China
文摘The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite element model for two-dimensional dynamic analysis of armored cable is presented. This model accounts for the geometric nonlinearities of large displacement of the armored cable, and effects of axial load and bending stiffness. The governing equations are derived by consistent linearization and finite element discretization of the total weak form of the armored cable system, and solved by the Newmark time integration method. To make the solution procedure avoid falling into the local extreme points, a simple adaptive stepping strategy is proposed. The presented model is validated via actual measured data. Results for dynamic configurations, motion and tension of both ends of the armored cable, and resonance-zone are presented for two numerical cases, including the dynamic analysis under the case of only ship heave motion and the case of joint action of ship heave motion and ocean current. The dynamics analysis can provide important reference for the design or product selection of the armored cable in a deep-sea ROV system so as to improve the safety of its marine operation under the sea state of 4 or above.
文摘Multi-pass slab vertical-horizontal (V-H) rolling process with variable edging roll shape have been simulated with explicit dynamic finite element method and updating geometric method. The distributions of plastic strain contour in slab daring rolling process with different edging roll and under different rolling stage have been obtained. The results show that there exist two thin strain assembling zones in slab when the flat edging roll is used, and there just exist one strain assembling zone in slab when the edging roll with groove is used. And compared the deformation equality between flat edging roll and edging roll with groove, the lateris better than the former, which supplies the theory prove to the slab deformation distribution during V-H rolling process and is helpful for predicting the slab texture.
基金Acknowledgements This work was supported by the National Natural Science Foundation of China (Grant Nos. 11371218, 91630205).
文摘We extend the theory of global geometrical optics method, proposed originally for the linear scalar high-frequency wave-like equations in [Commun. Math. Sci., 2013, 11(1): 105-140], to the more general vector- valued Schrodinger problems in the semi-classical regime. The key ingredient in the global geometrical optics method is a moving frame technique in the phase space. The governing equation is transformed into a new equation but of the same type when expressed in any moving frame induced by the underlying Hamiltonian flow. The classical Wentzel-Kramers-Brillouin (WKB) analysis benefits from this treatment as it maintains valid for arbitrary but fixed evolutionary time. It turns out that a WKB-type function defined merely on the underlying Lagrangian submanifold can be obtained with the help of this moving frame technique, and from which a uniform first-order approximation of the wave field can be derived, even around caustics. The general theory is exemplified by two specific instances. One is the two-level SchrSdinger system and the other is the periodic SchrSdinger equation. Numerical tests validate the theoretical results.
基金The work is supported in part by the Foundation of Zhongshan University, Advanced Research Centre and in part by the National Natural Science Foundation of China
文摘PL homotopy metheds are effective methods to locate zerces(or fixed points) of highly nonlinearmappirgs. Due to the Jexicographical system, the methods are feasible without exceptions Thispaper presents a geemetrical interpretation of the without-exception feasibility.
基金supported by the National Research Program for Universities(No.3183)
文摘The objective of this article is to introduce a generalized algorithm to produce the m-point n-ary approximating subdivision schemes(for any integer m, n ≥ 2). The proposed algorithm has been derived from uniform B-spline blending functions. In particular, we study statistical and geometrical/traditional methods for the model selection and assessment for selecting a subdivision curve from the proposed family of schemes to model noisy and noisy free data. Moreover, we also discuss the deviation of subdivision curves generated by proposed family of schemes from convex polygonal curve. Furthermore, visual performances of the schemes have been presented to compare numerically the Gibbs oscillations with the existing family of schemes.
基金supported by the National Science Foundation of China(No.51207150)the Director Foundation of Institute of Microelectronics of Chinese Academy of Sciences(Nos.Y2SF017001,Y3SZ0701)
文摘Capacitance extraction is one of the key issues in integrated circuits and also a typical electrostatic prob- lem. The dual discrete geometric method (DGM) is investigated to provide relative solutions in two-dimensional unstructured mesh space. The energy complementary characteristic and quick field energy computation thereof based on it are emphasized. Contrastive analysis between the dual finite element methods and the dual DGMs are presented both from theoretical derivation and through case studies. The DGM, taking the scalar potential as unknown on dual interlocked meshes, with simple form and good accuracy, is expected to be one of the mainstreaming methods in associated areas.
基金supported by the National Natural Science Foundation of China(Nos.61574167,51407181)
文摘Sensitivity analysis methods help to deal with the challenges of process variation in extraction of para- sitic capacitances in an integrated circuit. The dual discrete geometric methods (DGMs), which have been recently utilized to extract parasitic capacitances, are reviewed. The computation method based on the dual DGMs for sen- sitivities of capacitances with respect to the given process parameters is presented. As the dual DGMs utilize scalar electric potential is unknown, the capacitances are obtained effectively, and then the sensitivities are calculated conveniently.
基金the National Natural Science Foundation of China(No.61773014)。
文摘In this paper,we consider a GI/M/1 queue operating in a multi-phase service environment with working vacations and Bernoulli vacation interruption.Whenever the queue becomes empty,the server begins a working vacation of random length,causing the system to move to vacation phase 0.During phase 0,the server takes service for the customers at a lower rate rather than stopping completely.When a vacation ends,if the queue is non-empty,the system switches from the phase 0 to some normal service phase i with probability qi,i=1,2,⋯,N.Moreover,we assume Bernoulli vacation interruption can happen.At a service completion instant,if there are customers in a working vacation period,vacation interruption happens with probability p,then the system switches from the phase 0 to some normal service phase i with probability qi,i=1,2,⋯,N,or the server continues the vacation with probability 1−p.Using the matrix geometric solution method,we obtain the stationary distributions for queue length at both arrival epochs and arbitrary epochs.The waiting time of an arbitrary customer is also derived.Finally,several numerical examples are presented.
基金supported by Beijing Municipal Commission of Education (201501001)Beijing Municipal Science and Technology Commission (Z16111000500000)+1 种基金the National Natural Science Foundation of China (61671073)supported by Beijing Laboratory of Advanced Information Network
文摘The unforeseen mobile data explosion as well as the scarce of spectrum resource pose a major challenge to the performance of today's cellular networks which are in urgent need of novel solutions to handle such voluminous mobile data. Long term evolution-unlicensed (LTE-U), which extends the LTE standard operating on the unlicensed band, has been proposed to improve system throughput. In LTE-U system, arriving users will contend the unlicensed spectrum resource with wireless fidelity (WiFi) users to transmit data information. Nevertheless, there is no clear consensus as to the benefits of transmission using unlicensed bands for LTE users. To this end, in this paper an analytical model is presented based on a queue system to understand the performance achieved by unlicensed based LTE system taking quality of services (QoS) and LTE-U users' behaviors into account. To obtain the stead-state solutions of the queue system, a matrix geometric method is used to solve it. Then, the average delay and utilization of unlicensed band for the LTE-U users is derived by using the queuing model. The performance of LTE-U coexistence is evaluated with WiFi using the proposed model and provide some initial insights as to the advantage of LTE-U in practice.
基金supported by National Basic Research Program of China(Grant No.2006CB805902)
文摘In this paper,we study the Cauchy problem of an integrable evolution system,i.e.,the n-dimensional generalization of third-order symmetry of the well-known Landau-Lifshitz equation.By rewriting this equation in a geometric form and applying the geometric energy method with a forth-order perturbation,we show the global well-posedness of the Cauchy problem in suitable Sobolev spaces.