Jianchen Wu¹, Yujie Guo^1,*, Fangli Wang^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.128, No.3, pp. 1007-1031, 2021, DOI:10.32604/cmes.2021.016475

Abstract Spacecraft flexible appendages may experience thermally induced vibrations (TIV) under sudden heating loads, which in consequence will be unable to complete their intended missions. Isogeometric analysis (IGA) utilizes, in an isoparametric concept, the same high order and high continuity non-uniform rational B-splines (NURBS) to represent both the geometry and the physical field of the structure. Compared to the traditional Lagrange polynomial based finite element method where only C0-continuity across elements can be achieved, IGA is geometrically exact and naturally fulfills the C1-continuity requirement of Euler–Bernoulli (EB) beam elements, therefore, does not need extra rotational degrees-of-freedom. In this paper, we present… More >

Weihua Fang¹, Tiantang Yu^2,*, Yin Yang³

CMES-Computer Modeling in Engineering & Sciences, Vol.126, No.3, pp. 1315-1332, 2021, DOI:10.32604/cmes.2021.014575

Abstract Steady-state heat transfer problems in heterogeneous solid are simulated by developing an adaptive extended isogeometric analysis (XIGA) method based on locally refined non-uniforms rational B-splines (LR NURBS). In the XIGA, the LR NURBS, which have a simple local refinement algorithm and good description ability for complex geometries, are employed to represent the geometry and discretize the field variables; and some special enrichment functions are introduced into the approximation of temperature field, thus the computational mesh is independent of the material interfaces, which are described with the level set method. Similar to the approximation of temperature field, a temperature gradient recovery… More >

Houlin Yang¹, Bingquan Zuo^1,2,*, Zhipeng Wei^1,2, Huixin Luo^1,2, Jianguo Fei^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.126, No.3, pp. 1033-1052, 2021, DOI:10.32604/cmes.2021.014493

Abstract 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 from the subdivision progress still leads to the solution efficiency problems. This paper presents a multigrid method based on geometric multigrid (GMG) to solve the matrix system of IGA. This method extracts the required computational data for multigrid method from the IGA process, which also can be used to improve the traditional… More >

Leilei Chen¹, Kunpeng Li¹, Xuan Peng², Haojie Lian^3,4,*, Xiao Lin⁵, Zhuojia Fu⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.126, No.1, pp. 125-146, 2021, DOI:10.32604/cmes.2021.012821

Abstract This paper presents an isogeometric boundary element method (IGABEM) for transient heat conduction analysis. The Non-Uniform Rational B-spline (NURBS) basis functions, which are used to construct the geometry of the structures, are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations. B´ezier extraction technique is employed to accelerate the evaluation of NURBS basis functions. We adopt a radial integration method to address the additional domain integrals. The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis. More >

Daoyuan Yu, Shouyu Cai^*, Wenya Fan, Lan Zhang^*

CMES-Computer Modeling in Engineering & Sciences, Vol.124, No.3, pp. 787-806, 2020, DOI:10.32604/cmes.2020.09936

Abstract In this study, a patch removing based Isogeometric analysis (PR-IGA) method is proposed to conduct the holed structural analysis with only one parametric domain, in which there are also no trimmed elements. The theoretical foundation of this novel patch removing approach is that any holed structure can be obtained by removing sub-patches (i.e., the holes) from an intact base patch. Since the parametric domains of these patches are all meshed by rectangular grids, the elements in the resulted holed structural parametric domain could all be untrimmed rectangles under certain mapping conditions. To achieve the special condition, a systematic technique consisting… More >

Yingjun Wang^1,*, Zhenpei Wang², Xiaowei Deng³, David J. Benson⁴, Damiano Pasini⁵, Shuting Wang⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.124, No.3, pp. 783-785, 2020, DOI:10.32604/cmes.2020.013234

Abstract This article has no abstract. More >

Hamid Ghasemi¹, Harold S. Park², Xiaoying Zhuang^{3, 4, *}, Timon Rabczuk^{5, 6}

CMC-Computers, Materials & Continua, Vol.65, No.2, pp. 1157-1179, 2020, DOI:10.32604/cmc.2020.08358

Abstract Flexoelectricity is a general electromechanical phenomenon where the electric polarization exhibits a linear dependency to the gradient of mechanical strain and vice versa. The truncated pyramid compression test is among the most common setups to estimate the flexoelectric effect. We present a three-dimensional isogeometric formulation of flexoelectricity with its MATLAB implementation for a truncated pyramid setup. Besides educational purposes, this paper presents a precise computational model to illustrate how the localization of strain gradients around pyramidal boundary shapes contributes in generation of electrical energy. The MATLAB code is supposed to help learners in the Isogeometric Analysis and Finite Elements Methods… More >

Jinlan Xu^*, Ningning Sun, Gang Xu

CMES-Computer Modeling in Engineering & Sciences, Vol.124, No.2, pp. 627-642, 2020, DOI:10.32604/cmes.2020.010341

Abstract This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications. In isogeometric analysis, multi-patch structure is not easy to achieve high continuity between neighboring patches which will reduce the advantage of isogeometric analysis in a sense. The proposed method can achieve high continuity at surface stitching region with low geometric error, and this technique exploits constructing the approximate surface with several control points are from original surfaces, which guarantees the local feature of the surface can be well-preserved with high precision. With the proposed approximating method, isogeometric analysis… More >

Yang Xia^*, Luting Deng, Jian Zhao

CMES-Computer Modeling in Engineering & Sciences, Vol.124, No.2, pp. 605-626, 2020, DOI:10.32604/cmes.2020.010204

Abstract Geometric fitting based on discrete points to establish curve structures is an important problem in numerical modeling. The purpose of this paper is to investigate the geometric fitting method for curved beam structure from points, and to get high-quality parametric model for isogeometric analysis. A Timoshenko beam element is established for an initially curved spacial beam with arbitrary curvature. The approximation and interpolation methods to get parametric models of curves from given points are examined, and three strategies of parameterization, meaning the equally spaced method, the chord length method and the centripetal method are considered. The influences of the different… More >

Shige Wang¹, Zhongwang Wang¹, Leilei Chen¹, Haojie Lian^2,3,*, Xuan Peng⁴, Haibo Chen⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.124, No.2, pp. 585-604, 2020, DOI:10.32604/cmes.2020.09904

Abstract The paper applied the isogeometric boundary element method (IGABEM) to thermoelastic problems. The Non-Uniform Rational B-splines (NURBS) used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation. Due to the existence of thermal stress, the domain integral term appears in the boundary integral equation. We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral. In this way, IGABEM can maintain its advantages in dimensionality reduction and more importantly, seamless integration of CAD and numerical analysis based on boundary representation. The algorithm is verified… More >