基于非结构化T样条的薄壳等几何分析

doi:10.3901/JME.2018.15.132

机械工程学报 ›› 2018, Vol. 54 ›› Issue (15): 132-140.doi: 10.3901/JME.2018.15.132

基于非结构化T样条的薄壳等几何分析

袁沛, 刘振宇, 裘辿, 谭建荣

浙江大学CAD&CG国家重点实验室杭州 310027

收稿日期:2017-08-13 修回日期:2018-01-23 出版日期:2018-08-05 发布日期:2018-08-05
通讯作者: 裘辿(通信作者),男,1983年出生,博士,讲师。主要研究方向为数字化设计、数字样机仿真与虚拟现实。E-mail:qc@zju.edu.cn
作者简介:袁沛,男,1988年出生,博士研究生。主要研究方向为数字样机仿真、等几何分析。E-mail:bifeng@zju.edu.cn;刘振宇,男,1974年出生,博士,教授,博士研究生导师。主要研究方向为数字化设计、数字样机仿真与虚拟现实。E-mail:liuzy@zju.edu.cn;谭建荣,男,1954年出生,博士,教授,中国工程院院士。主要研究方向为数字化设计、机械设计及理论。E-mail:egi@zju.edu.cn
基金资助:
国家自然科学基金资助项目（51475417，51490663，U1608256）。

Unstructured T-spline-based Isogeometric Analysis for Thin Shells

YUAN Pei, LIU Zhenyu, QIU Chan, TAN Jianrong

State Key Lab of CAD & CG, Zhejiang University, Hangzhou 310027

Received:2017-08-13 Revised:2018-01-23 Online:2018-08-05 Published:2018-08-05

摘要/Abstract

摘要： 等几何分析方法采用几何的样条基函数来构建分析模型，从而避免传统有限元法的网格离散过程。采用非结构化T样条构建复杂模型，结合Kirchhoff-Love薄壳理论构造基于非结构化T样条的薄壳单元，并研究其在模态计算和弹性变形分析中的应用。由于样条几何基函数缺乏插值性，等几何分析中单元所受载荷无法像传统有限元法一样直接均分到单元各个节点上，针对这一问题，将任意载荷作用区域分为两种基本形式，通过一定的映射关系将这两种形式由规则形状表示，从而将规则形状的高斯积分点同样映射到不规则区域上。数值算例结果表明，和传统有限元相比，基于非结构化T样条的薄壳等几何分析能够以更少的系统自由度获得精确解，而任意区域的积分映射方法也有效地解决了等几何分析中载荷施加的问题。

关键词: Bézier提取, 薄壳, 等几何分析, 非结构化T样条, 积分映射

Abstract: Isogeometric analysis(IGA) is a method aimed at avoiding meshing procedure of the traditional finite element analysis by using the same basis for analysis as is used to describe the geometry, thus integrating CAE into a general CAD modeling system. The unstructured T-splines are adopted to construct complex geometric models. By integrating the Kirchhoff-Love theory, an IGA thin-shell element is established for the modal analysis and elastic analysis. Due to the non-interpolatory nature of the geometric basis functions, the loads of IGA elements cannot be evenly distributed to the control points as with the traditional finite element. The loading areas are then classified to two basic types and an integration mapping method is proposed to map the Gauss quadrature points of the regular domain to the basic types. Numerical examples illustrated that the unstructured T-splines-based thin-shell element could achieve the convergence with much less degrees of freedom than the traditional finite element method and the integration mapping method solves the imposition problem of external loading effectively.

Key words: Bézier extraction, integration mapping, isogeometric analysis, thin-shell, unstructured T-splines

中图分类号:

TP391

袁沛, 刘振宇, 裘辿, 谭建荣. 基于非结构化T样条的薄壳等几何分析[J]. 机械工程学报, 2018, 54(15): 132-140.

YUAN Pei, LIU Zhenyu, QIU Chan, TAN Jianrong. Unstructured T-spline-based Isogeometric Analysis for Thin Shells[J]. Journal of Mechanical Engineering, 2018, 54(15): 132-140.

参考文献

[1] HUGHES T J R,COTTRELL J A,BAZILEVS Y. Isogeometric analysis:CAD,finite elements,NURBS,exact geometry and mesh refinement[J]. Comput. Methods Appl. Mech. Eng.,2005,194(39-41):4135-4195.
[2] BAZILEVS Y,HSU M-C,SCOTT M A. Isogeometric fluid-structure interaction analysis with emphasis on non-matching discretizations,and with application to wind turbines[J]. Comput. Methods Appl. Mech. Engrg.,2012,249-252:28-41.
[3] 吴紫俊,黄正东,左兵权,等. 等几何分析研究概述[J]. 机械工程学报,2015,51(5):114-129. WU Zijun,HUANG Zhengdong,ZUO Bingquan,et al. Perspectives on isogeometric analysis[J]. Journal of Mechanical Engineering,2015,51(5):114-129.
[4] BUFFA A,SANGALLI G,VÁZQUEZ R. Isogeometric methods for computational electromagnetics:B-spline and T-spline discretizations[J]. Journal of Computational Physics,2014,257:1291-1320.
[5] KIENDL J,SCHMIDT R,WÜCHNER R,et al. Isogeometric shape optimization of shells using semi-analytical sensitivity analysis and sensitivity weighting[J]. Comput. Methods Appl. Mech. Engrg.,2014,274:148-167.
[6] 蔡守宇,张卫红,李杨. 基于面片删减的带孔结构等几何形状优化方法[J]. 机械工程学报,2013,49(13):150-157. CAI Shouyu,ZHANG Weihong,LI Yang. Isogeometric shape optimization method with patch removal for holed structures[J]. Journal of Mechanical Engineering,2013,49(13):150-157.
[7] TAVAKKOLI S M,HASSANI B. Isogeometric topology optimization by using optimality criteria and implicit function[J]. International Journal of Optimization in Civil Engineering,2014,4(2):151-163.
[8] HASSANI B,KHANZADI M,TAVAKKOLI S M. An isogeometrical approach to structural topology optimization by optimality criteria[J]. Struct. Multidisc. Optim.,2012,45:223-233.
[9] BENSON D J,BAZILEVS Y,HSU M C,et al. Isogeometric shell analysis:The Reissner-Mindlin shell[J]. Comput. Methods Appl. Mech. Engrg.,2010,199:276-289.
[10] DORNISCH W,KLINKEL S,SIMEON B. Isogeometric Reissner-Mindlin shell analysis with exactly calculated director vectors[J]. Comput. Methods Appl. Mech. Engrg.,2013,253:491-504.
[11] KIENDL J,BLETZINGER K U,LINHARD J,et al. Isogeometric shell analysis with Kirchhoff-Love elements[J]. Comput. Methods Appl. Mech. Engrg.,2009,198:3902-3914.
[12] SCOTT M A,SIMPSON R N,EVANS J A,et al. Isogeometric boundary element analysis using unstructured T-splines[J]. Comput. Methods Appl. Mech. Engrg.,2013,254:197-221.
[13] KIENDL J,BAZILEVS Y,HSU M C,et al. The bending strip method for isogeometric analysis of Kirchhoff-Love shell structures comprised of multiple patches[J]. Computer Methods in Applied Mechanics and Engineering,2010,199:2403-2416.
[14] SEDERBERG T W,ZHENG J,BAKENOV A,et al. T-Splines and T-NURCCS[J]. ACM Trans. Graph.,2003,22(3):477-484.
[15] BAZILEVS Y,CALO V M,COTTRELL J A,et al. Isogeometric analysis using T-splines[J]. Comput. Methods Appl. Mech. Engrg.,2010,199:229-263.
[16] SCOTT M A,BORDEN M J,VERHOOSEL C V,et al. Isogeometric finite element data structures based on Bézier extraction of T-splines[J]. Int. J. Numer. Meth. Engng,2011,88:126-156.
[17] SCOTT M A,HUGHES T J R,SEDERBERG T W,et al. An integrated approach to engineering design and analysis using the autodesk t-spline plugin for rhino3d[C]//ICES REPORT,The Institute for Computational Engineering and Sciences,The University of Texas,Austin,2014,14-33.
[18] DIMITRI R,LORENZIS L D,SCOTT M A,et al. Isogeometric large deformation frictionless contact using T-splines[J]. Comput. Methods Appl. Mech. Engrg.,2014,269:394-414.
[19] BUFFA A,SANGALLI G,VÁZQUEZ R. Isogeometric methods for computational electromagnetics:B-spline and T-spline discretizations[J]. Journal of Computational Physics,2014,257:1291-1320.
[20] SEDERBERG T W,CARDON D L,FINNIGAN GT,et al. T-spline simplification and local refinement[J]. ACM Trans. Graph.,2004,23(3):276-283.
[21] WANG A,ZHAO G,LI Y-D. An influence-knot set based new local refinement algorithm for T-spline surfaces[J]. Expert Systems with Applications,2014,41:3915-3921.
[22] GUO Y,RUESS M. Weak Dirichlet boundary conditions for trimmed thin isogeometric shells[J]. Computers and Mathematics with Applications,2015,70:1425-1440.
[23] RUESS M,SCHILLINGER D,BAZILEVS Y,et al. Weakly enforced essential boundary conditions for NURBS-embedded and trimmed NURBS geometries on the basis of the finite cell method[J]. Int. J. Numer. Meth. Engng.,2013,95:811-846.
[24] WANG D,XUAN J. An improved NURBS-based isogeometric analysis with enhanced treatment of essential boundary conditions[J]. Comput. Methods Appl. Mech. Engrg.,2010,199:2425-2436.
[25] KIM H J,SEO Y D,YOUN S K. Isogeometric analysis for trimmed CAD surfaces[J]. Comput. Methods Appl. Mech. Engrg.,2009,198:2982-2995.
[26] REDDY J. Theory and analysis of elastic plates and shells[D]. Boca Raton:CRC Press,Taylor and Francis,2007.
[27] BELYTSCHKO T,STOLARSKI H,LIU W K,et al. Stress projection for membrane and shear locking in shell finite elements[J]. Comput. Methods Appl. Mech. Engrg.,1985,51:221-258.

基于非结构化T样条的薄壳等几何分析

Unstructured T-spline-based Isogeometric Analysis for Thin Shells

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 4

编辑推荐

Metrics

本文评价

[1]	吴紫俊, 黄正东, 左兵权, 刘清华, 殷小亮. 等几何分析研究概述[J]. 机械工程学报, 2015, 51(5): 114-129.
[2]	蔡守宇;张卫红;李杨. 基于面片删减的带孔结构等几何形状优化方法[J]. , 2013, 49(13): 150-157.
[3]	陈涛;莫蓉;万能. 等几何分析中Dirichlet边界条件的配点施加方法[J]. , 2012, 48(5): 157-164.
[4]	丁晓红;林建中;山崎光悦. 利用植物根系形态形成机理的加筋薄壳结构拓扑优化设计[J]. , 2008, 44(4): 201-205.