基于等几何分析的参数化曲梁结构非线性动力学降阶模型研究

郭玉杰; 吴晗浪; 李薇; 李迎港; 吴剑晨

doi:10.6052/j.issn.1000-4750.2021.04.0285

基于等几何分析的参数化曲梁结构非线性动力学降阶模型研究

doi: 10.6052/j.issn.1000-4750.2021.04.0285

南京航空航天大学飞行器先进设计技术国防重点学科实验室，南京 210016

基金项目: 国家自然科学基金面上项目(11972187)

详细信息

作者简介:
吴晗浪(1993−)，男，贵州人，硕士生，主要从事等几何分析研究(E-mail: 473807649@qq.com)

李　薇(1998−)，女，河南人，硕士生，主要从事等几何结构动力学研究(E-mail: cyona@nuaa.edu.cn)

李迎港(1997−)，男，安徽人，硕士生，主要从事结构疲劳寿命研究(E-mail: liyinggang@nuaa.edu.cn)

吴剑晨(1997−)，男，安徽人，硕士生，主要从事航天器传热研究(E-mail: jianchenwu@nuaa.edu.cn)

通讯作者:
郭玉杰(1986−)，男，江苏人，副教授，博士，主要从事结构计算力学研究(E-mail: yujieguo@nuaa.edu.cn)

中图分类号: O342；TB122
计量
- 文章访问数: 488
- HTML全文浏览量: 156
- PDF下载量: 118
- 被引次数: 0
出版历程
- 收稿日期: 2021-04-16
- 修回日期: 2021-07-05
- 网络出版日期: 2021-08-06
- 刊出日期: 2022-08-01

MODEL ORDER REDUCTION FOR NONLINEAR DYNAMIC ANALYSIS OF PARAMETERIZED CURVED BEAM STRUCTURES BASED ON ISOGEOMETRIC ANALYSIS

Key Laboratory of Fundamental Science for National Defense-Advanced Design Technology of Flight Vehicle, College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

摘要

摘要: 模型降阶方法通过构造全阶模型的低阶近似模型有效提升了求解效率，同时也保留了原阶模型的主要信息从而保证了较高的计算精度。对于结构非线性以及参数化的模型降阶问题，常需要重复计算刚度矩阵等非线性以及参数依赖项，求解效率较低。此外，当参数化模型的几何形状改变时，往往需要重复进行CAD与有限元(FEA)模型的转换，这对于复杂结构较为耗时。等几何分析采用描述几何形状的非均匀有理B样条(NURBS)插值物理场，实现了CAD与FEA模型的统一，消除了两者之间繁琐的模型转换过程，其具有几何精确、高阶连续等优点，并且几何形状在细化过程中保持不变，非常适合于薄壁类结构的分析以及参数化表达。该研究结合等几何分析、特征正交分解(POD)以及离散经验插值方法(DEIM)研究参数化的平面曲梁结构的非线性动力学模型降阶问题。数值结果表明，基于等几何分析的POD-DEIM降阶模型能够显著提升平面曲梁结构的非线性动力学计算效率，并且该模型对于参数化以及变载荷等情形显示出了良好的适应性。
- 等几何分析 /
- 模型降阶 /
- 非线性动力学 /
- 离散经验插值 /
- 参数化模型
Abstract: Model order reduction (MOR) approach generates lower dimensional approximations to the original system while preserving model's essential information and computational accuracy. For nonlinear and parameterized structural problems, where the stiffness matrix is configuration dependent, an iterative solution procedure is inevitable and a revisit to all the elements is essential for updating the stiffness matrix. Considering the difference between CAD and finite element models, a preprocessing step is necessary for each parameterized geometry, which is time consuming for complex geometries. Isogeometric analysis (IGA) utilizes non-uniform rational B-spline (NURBS) basis for both the geometry description and physical filed interpolation, thusly eliminates the time-consuming preprocessing step between CAD and finite element models. IGA is a perfect candidate for the analysis and geometry parameterization of thin-walled structures due to its geometrically exact and high order continuous properties. The nonlinear dynamics of the parameterized planar curved beams are studied based on the isogeometric analysis and their model order reductions are investigated based on the proper orthogonal decomposition and discrete empirical interpolation method (POD-DEIM). Numerical results show that IGA-based POD-DEIM method significantly improves the computational efficiency of the nonlinear dynamic analysis of planar curved beams. Additionally, the proposed method also applies to different external loadings and time step sizes, which demonstrate the adaptivity of the method.
- isogeometric analysis /
- model order reduction /
- nonlinear dynamic analysis /
- discrete empirical interpolation method /
- parameterized models

HTML全文

图 1 NURBS曲线及NURBS基函数

Figure 1. NURBS curve and NURBS basis functions

下载: 全尺寸图片幻灯片

图 2 初始构型与当前构型

Figure 2. Reference and current configurations

下载: 全尺寸图片幻灯片

图 3 HHT-α时间积分算法流程图

Figure 3. HHT-α time integration scheme for nonlinear dynamic analysis

下载: 全尺寸图片幻灯片

图 4 两端固支半圆弧曲梁模型

Figure 4. Half-circle arc with clamped boundary conditions at both ends

下载: 全尺寸图片幻灯片

图 5 两端固支半圆弧曲梁位移收敛曲线(几何非线性)

Figure 5. Convergence studies of the half-circle arc beam(geometrically nonlinear)

下载: 全尺寸图片幻灯片

图 6 两端固支半圆弧曲梁奇异值

Figure 6. Singular values of displacement snapshots and stiffness matrix snapshots for polynomial degrees 2, 3, 4 for clamped half-circle arc

下载: 全尺寸图片幻灯片

图 7 两端固支半圆弧曲梁载荷-位移曲线

Figure 7. Load-displacement response of half-circle arc beam

下载: 全尺寸图片幻灯片

图 8 两端固支半圆弧曲梁中点轨迹对比图

Figure 8. Comparison of the paths traced by the mid-point of the half-circle arc beam obtained using the FOM and ROM

下载: 全尺寸图片幻灯片

图 9 两端固支半圆弧曲梁降阶与全阶模型位移差值图

Figure 9. Difference of mid-point displacement between FOM and ROM for half-circle arc beam

下载: 全尺寸图片幻灯片

图 10 两端固支半圆弧曲梁降阶、全阶模型计算时间对比

Figure 10. Comparison of computational time for different reduced order models of half-circle arc beam

下载: 全尺寸图片幻灯片

图 11 一端固支阿基米德螺旋梁

Figure 11. Archimedes spiral shaped beam with one end clamped

下载: 全尺寸图片幻灯片

图 12 一端固支阿基米德螺旋曲梁奇异值

Figure 12. Singular values of displacement snapshots and stiffness matrix snapshots for clamped Archimedes spiral beam

下载: 全尺寸图片幻灯片

图 13 一端固支阿基米德螺旋曲梁载荷位移曲线

Figure 13. Load-displacement curve of clamped Archimedes spiral arc beam

下载: 全尺寸图片幻灯片

图 14 一端固支阿基米德螺旋曲梁端点的轨迹对比图

Figure 14. Comparison of the paths traced by the end-point of the Archimedes spiral beam obtained using the FOM and ROM

下载: 全尺寸图片幻灯片

图 15 阿基米德螺旋曲梁端点降阶与全阶模型位移差值图

Figure 15. Difference of end-point displacement between FOM and ROM for the Archimedes spiral beam

下载: 全尺寸图片幻灯片

图 16 两端简支曲梁模型

Figure 16. Arc beam with simply supported boundary conditions at both ends

下载: 全尺寸图片幻灯片

图 17 两端简支曲梁奇异值

Figure 17. Singular values of displacement snapshots and stiffness matrix snapshots for simply supported arc beam

下载: 全尺寸图片幻灯片

图 18 两端简支曲梁载荷-位移曲线

Figure 18. Load-displacement curve of simply supported arc beam

下载: 全尺寸图片幻灯片

图 19 两端简支曲梁中点的轨迹对比图

Figure 19. Comparison of the paths traced by the mid-point of the arc beam obtained using the FOM and ROM

下载: 全尺寸图片幻灯片

图 20 两端简支曲梁降阶与全阶模型位移差值图

Figure 20. Difference of mid-point displacement between FOM and ROM for the arc beam

下载: 全尺寸图片幻灯片

图 21 两端简支曲梁频率响应对比图

Figure 21. Comparison of frequency response obtained using FOM and ROM

下载: 全尺寸图片幻灯片

图 22 两端简支曲梁载荷位移曲线($F = - 800\sin (36t)$)

Figure 22. Load-displacement curve of simply supported arc beam ($F = - 800\sin (36t)$)

下载: 全尺寸图片幻灯片

图 23 两端简支曲梁中点的轨迹对比图($F = - 800\sin (36t)$)

Figure 23. Comparison of the paths traced by the mid-point of the arc beam obtained using the FOM and ROM ($F = - 800\sin (36t)$)

下载: 全尺寸图片幻灯片

图 24 两端简支曲梁降阶与全阶模型位移差值图($F = - 800\sin (36t)$)

Figure 24. Difference of mid-point displacement between FOM and ROM for the arc beam ($F = - 800\sin (36t)$)

下载: 全尺寸图片幻灯片

图 25 两端简支曲梁频率响应对比图($F = - 800\sin (36t)$)

Figure 25. Comparison of frequency response obtained using FOM and ROM ($F = - 800\sin (36t)$)

下载: 全尺寸图片幻灯片

图 26 两端简支曲梁的参数化模型

Figure 26. Parametric model of simply supported curved beams

下载: 全尺寸图片幻灯片

图 27 控制点P₂抽样分布图

Figure 27. Sampling positions of the control point P₂

下载: 全尺寸图片幻灯片

图 28 参数化曲梁奇异值

Figure 28. Singular values of displacement snapshots and stiffness matrix snapshots for parameterized arc beam

下载: 全尺寸图片幻灯片

图 29 不同抽样方法所对应的位移载荷响应

Figure 29. Load-displacement response of half-circle arc beam

下载: 全尺寸图片幻灯片

图 30 参数化曲梁$\xi = 0.5$处的轨迹对比图

Figure 30. Comparison of the paths traced at the position of $\xi = 0.5$ for the parameterized arc beam obtained using the FOM and ROM

下载: 全尺寸图片幻灯片

图 31 参数化曲梁降阶与全阶模型位移差值图$(\xi = 0.5)$

Figure 31. Difference of the displacement between FOM and ROM for the parameterized arc beam $(\xi = 0.5)$

下载: 全尺寸图片幻灯片

图 32 参数化曲梁频率响应对比图

Figure 32. Comparison of frequency response obtained using FOM and ROM for the parameterized arc beam

下载: 全尺寸图片幻灯片

图 33 不同抽样方法所对应的位移载荷响应$(\xi = 0.5)$ ($F = - 800\sin (48t)$)

Figure 33. Comparisons of load-displacement histories of the parameterized arc beam at $(\xi = 0.5)$ ($F = - 800\sin (48t)$)

下载: 全尺寸图片幻灯片

图 34 不同抽样方法所对应的位移载荷响应$(\xi = 0.5)$ ($F = - 800\sin (48t)$, 时间步长：0.005 s)

Figure 34. Comparisons of load-displacement histories of the parameterized arc beam at $(\xi = 0.5)$($F = - 800\sin (48t)$, time step: 0.005 s)

下载: 全尺寸图片幻灯片

图 35 参数化曲梁$(\xi = 0.5)$处的轨迹对比图($F = - 800\sin (48t)$, 时间步长：0.005 s)

Figure 35. Comparison of the paths traced at the position of $(\xi = 0.5)$ for the parameterized arc beam obtained using the FOM and ROM ($F = - 800\sin (48t)$, time step: 0.005 s)

下载: 全尺寸图片幻灯片

图 36 参数化曲梁降阶与全阶模型位移差值图$(\xi = 0.5)$($F = - 800\sin (48t)$, 时间步长：0.005 s)

Figure 36. Difference of the displacement between FOM and ROM for the parameterized arc beam $(\xi = 0.5)$ ($F = - 800\sin (48t)$, time step: 0.005 s)

下载: 全尺寸图片幻灯片

图 37 参数化曲梁频率响应对比图($F = - 800\sin (48t)$, 时间步长：0.005 s)

Figure 37. Comparison of frequency response obtained using FOM and ROM for the parameterized arc beam ($F = - 800\sin (48t)$, time step: 0.005 s)

下载: 全尺寸图片幻灯片

表 1 两端固支半圆弧曲梁降阶模型计算时间对比

Table 1. Comparisons of computational time for reduced order models of half-circle arc beam

计算模型	全阶模型	POD降阶模型				POD-MDEIM 降阶模型
计算模型	全阶模型	k = 2	k = 3	k = 5	k = 15	k = 20 m = 15
计算时间/s	302.9	208.9	212.2	234.3	242.4	36.9
时间比	1	1.4	1.4	1.3	1.2	8.2

下载: 导出CSV

表 2 两端简支曲梁刚度矩阵MDEIM插值位置

Table 2. MDEIM Positions for stiffness matrix of simply supported arc beam

位置$ {\phi }_{i} $	行(i)	列(j)	单元编号e	−	位置$ {\phi }_{i} $	行(i)	列(j)	单元编号e
223	1	3	1		7057	64	64	19,20,21,22
225	3	3	1		7615	67	69	20,21,22,23
2689	25	25	6,7,8,9		7953	72	72	33,34,35,36
3697	34	34	9,10,11,12		8961	81	81	24,25,26,27
3699	36	34	9,10,11,12		10081	91	91	28,29,30,31
4371	42	40	11,12,13,14		10753	97	97	30,31,32,33
5599	49	51	14,15,16,17		11091	102	100	31,32,33,34
6607	58	60	17,18,19,20		−	−	−	−

下载: 导出CSV

表 3 两端简支曲梁刚度矩阵MDEIM插值位置

Table 3. MDEIM Positions for stiffness matrix of simply supported arc beam

位置$ {\phi }_{i} $	行i	列j	单元编号e	位置$ {\phi }_{i} $	行i	列j	单元编号e
243	3	3	1	8953	73	75	22,23,24,25
3752	32	32	8,9,10,11	11 738	98	98	30,31,32,33
4237	37	36	10,11,12	12 827	107	107	33,34,35,36
4841	41	41	11,12,13,14	13 795	115	115	36,37
5082	42	43	12,13,14	13 797	117	115	36,37
5814	54	49	15,16,17	14 037	117	117	36,37
7135	55	60	17,18,19	14 158	118	118	37
8108	68	68	20,21,22,23	−	−	−	−

下载: 导出CSV

表 4 参数化曲梁降阶模型计算时间对比

Table 4. Comparisons of computational time for reduced order models of parameterized arc beam

控制点P₂坐标	计算时间/s			速度提升率
	全阶模型	POD-MDEIM 降阶模型		抽样1	抽样2
	全阶模型	抽样1	抽样2	抽样1	抽样2
P₂ = (7.5,9.5) m F = −800 N	237.0	41.1	43.5	5.8	5.4
P₂ = (5.7,7.2) m F = −800 N	242.1	39.0	37.6	6.2	6.4
P₂ = (6,5.5) m F = −800sin(48t) N	402.0	78.1	71.6	5.1	5.6

下载: 导出CSV

参考文献(46)

[1]	Tiso P, Jansen E, Abdalla M. Reduction method for finite element nonlinear dynamic analysis of shells [J]. AIAA Journal, 2011, 49: 2259 − 2304.
[2]	Nickell R E. Nonlinear dynamics by mode superposition [J]. Computer Methods in Applied Mechanics and Engineering, 1976, 7: 107 − 129. doi: 10.1016/0045-7825(76)90008-6
[3]	Idelsohn S R, Cardona A A. A reduction method for nonlinear structural dynamic analysis [J]. Computer Methods in Applied Mechanics and Engineering, 1985, 49: 253 − 279. doi: 10.1016/0045-7825(85)90125-2
[4]	Idelsohn S R, Cardona A A. A load-dependent basis for reduced nonlinear structural dynamics [J]. Computers and Structures, 1985, 20: 203 − 210. doi: 10.1016/0045-7949(85)90069-0
[5]	Krysl P L, Marsden J. Dimensional model reduction in nonlinear finite element dynamics of solids and structures [J]. International Journal for Numerical Methods in Engineering, 2001, 51: 479 − 504. doi: 10.1002/nme.167
[6]	郑保敬, 梁钰, 高效伟, 等. 功能梯度材料动力学问题的POD模型降阶分析[J]. 力学学报, 2018, 50(4): 787 − 797. doi: 10.6052/0459-1879-18-069 Zheng Baojing, Liang Yu, Gao Xiaowei, et al. Analysis for dynamic response of functionally graded materials using POD based reduced order model [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(4): 787 − 797. (in Chinese) doi: 10.6052/0459-1879-18-069
[7]	Chaturanabut S. Nonlinear model reduction via discrete empirical interpolation [D]. USA: Rice University, 2011.
[8]	Astrid P. Reduction of process simulation models: A proper orthogonal decomposition approach [D]. The Netherlands: Eindhoven University of Technology, 2004.
[9]	Rewienski M, White J. A trajectory piecewise-linear approach to model order reduction and fast simulation of nonlinear circuits and micromachined devices [J]. IEEE Transactions on Computer Aided Design of Integrated Circuits and Systems, 2003, 22: 155 − 170. doi: 10.1109/TCAD.2002.806601
[10]	Gaonkar A K, Kulkarni S S. Model order reduction for dynamic simulation of slender beams undergoing large rotations [J]. Computational Mechanics, 2017, 59: 809 − 829. doi: 10.1007/s00466-017-1374-7
[11]	Gaonkar A K, Kulkarni S S. Model order reduction for dynamic simulation of beams with forcing and geometric nonlinearities [J]. Finite Elements in Analysis and Design, 2013, 76: 50 − 62. doi: 10.1016/j.finel.2013.08.001
[12]	满兴博, 伍晓红, 孙清. 采用非线性Galerkin方法的柔性梁模型降阶研究[J]. 西安交通大学学报, 2015, 49(7): 113 − 119. doi: 10.7652/xjtuxb201507019 Man Xingbo, Wu Xiaohong, Sun Qing. Study on the model order reduction of flexible beam based on nonlinear Galerkin method [J]. Journal of Xi’an Jiaotong University, 2015, 49(7): 113 − 119. (in Chinese) doi: 10.7652/xjtuxb201507019
[13]	Barrault M, Maday Y, Nguyen N C, et al. An "empirical interpolation" method: Application to efficient reduced basis discretization of partial differential equations [J]. Comptes Rendus Mathematique, 2004, 339: 667 − 672. doi: 10.1016/j.crma.2004.08.006
[14]	Grepl M A, Maday Y, Nguyen N C, et al. Efficient reduced-basis treatment of nonaffine and nonlinear partial differential equations [J]. Mathematical Modelling and Numerical Analysis, 2007, 41: 575 − 605. doi: 10.1051/m2an:2007031
[15]	Galbally D, Fidkowski K, Wilcox K, et al. Non-linear model reduction for uncertainty quantification in large-scale inverse problems [J]. International Journal for Numerical Methods in Engineering, 2010, 81: 1581 − 1608.
[16]	Hesthaven J S, Rozza G, Stamm B. Certified reduced basis methods for parameterized partial differential equations [M]. Switzerland: Springer, 2016.
[17]	Prud’homme C, Rovas D V, Veroy K, et al. Reliable real-time solution of parameterized partial differential equations: Reduced-basis output bound methods [J]. Journal of Fluids Engineering, 2002, 124: 70 − 80. doi: 10.1115/1.1448332
[18]	Balmes E. Parametric families of reduced finite element models: Theory and applications [J]. Mechanical Systems and Signal Processing, 1996, 10: 381 − 394. doi: 10.1006/mssp.1996.0027
[19]	Huynh P, Rozza G, Sen S, et al. A successive constraint linear optimization method for lower bounds of parametric coercivity and inf-sup stability constants [J]. Comptes Rendus Mathematique, 2007, 345: 473 − 478. doi: 10.1016/j.crma.2007.09.019
[20]	Rinaldi M. Reduced basis method for isogeometric analysis: Application to structural problems [D]. Italy: Polytecnico Di Milano, 2014.
[21]	Dedden R J. Model order reduction using the discrete empirical interpolation method [D]. The Netherlands: Delft University of Technology, 2012.
[22]	Salmoiraghi F, Ballarin F, Heltai L, et al. Isogeometric analysis-based reduced order modelling for incompressible linear viscous flows in parameterized shapes [J]. Advanced Modelling and Simulation in Engineering Sciences, 2016, 3: 21. doi: 10.1186/s40323-016-0076-6
[23]	Hughes T J R, Cottrell J A, Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement [J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194: 4135 − 4195. doi: 10.1016/j.cma.2004.10.008
[24]	徐岗, 李新, 黄章进, 等. 面向等几何分析的几何计算[J]. 计算机辅助设计与图形学报, 2015, 27(4): 570 − 581. Xu Gang, Li Xin, Huang Zhangjin, et al. Geometric computing for isogeometric analysis [J]. Journal of Computer-Aided Design & Computer Graphics, 2015, 27(4): 570 − 581. (in Chinese)
[25]	何宜谦, 王霄腾, 祝雪峰, 等. 求解粘弹性问题的时域自适应等几何比例边界有限元法[J]. 工程力学, 2020, 37(2): 23 − 33. doi: 10.6052/j.issn.1000-4750.2019.01.0099 He Yiqian, Wang Xiaoteng, Zhu Xuefeng, et al. A temporally piecewise adaptive isogeometric SBFEM for viscoelastic problems [J]. Engineering Mechanics, 2020, 37(2): 23 − 33. (in Chinese) doi: 10.6052/j.issn.1000-4750.2019.01.0099
[26]	刘涛, 汪超, 刘庆运, 等. 基于等几何方法的压电功能梯度板动力学及主动振动控制分析[J]. 工程力学, 2020, 37(12): 228 − 242. doi: 10.6052/j.issn.1000-4750.2020.04.0266 Liu Tao, Wang Chao, Liu Qingyun, et al. Analysis for dynamic and active vibration control of piezoelectric functionally graded plates based on isogeometric method [J]. Engineering Mechanics, 2020, 37(12): 228 − 242. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.04.0266
[27]	汪超, 谢能刚, 黄璐璐. 基于扩展等几何分析和混沌离子运动算法的带孔结构形状优化设计[J]. 工程力学, 2019, 36(4): 248 − 256. doi: 10.6052/j.issn.1000-4750.2018.03.0119 Wang Chao, Xie Nenggang, Huang Lulu. Design and shape optimization of holed structures by extended isogeometric analysis and chaotic ion motion optimization [J]. Engineering Mechanics, 2019, 36(4): 248 − 256. (in Chinese) doi: 10.6052/j.issn.1000-4750.2018.03.0119
[28]	Kiendl J, Bletzinger K -U, Linhard J, et al. Isogeometric shell analysis with Kirchhoff-Love elements [J]. Computer Methods in Applied Mechanics and Engineering, 2009, 198: 3902 − 3914. doi: 10.1016/j.cma.2009.08.013
[29]	Benson D J, Bazilevs Y, Hsu M C, et al. Isogeometric shell analysis: The Reissner-Mindlin shell [J]. Computer Methods in Applied Mechanics and Engineering, 2010, 199: 276 − 289. doi: 10.1016/j.cma.2009.05.011
[30]	Hosseini S, Remmers J J C, Verhoosel C V, et al. An isogeometric solid-like shell element for nonlinear analysis [J]. International Journal for Numerical Methods in Engineering, 2013, 95: 238 − 256. doi: 10.1002/nme.4505
[31]	Echter R, Oesterle B, Bischoff M. A hierarchic family of isogeometric shell finite elements [J]. Computer Methods in Applied Mechanics and Engineering, 2013, 254: 170 − 180. doi: 10.1016/j.cma.2012.10.018
[32]	Nagy A P, Abdalla M M, Gurdal Z. Isogeometric sizing and shape optimization of beam structures [J]. Computer Methods in Applied Mechanics and Engineering, 2010, 199: 1216 − 1230. doi: 10.1016/j.cma.2009.12.010
[33]	Weeger O, Narayanan B, Dunn M L. Isogeometric shape optimization of nonlinear, curved 3D beams and beam structures [J]. Computer Methods in Applied Mechanics and Engineering, 2019, 345: 26 − 51. doi: 10.1016/j.cma.2018.10.038
[34]	Maurin F, Dede L, Spadoni A. Isogeometric rotation-free analysis of planar extensible-elastica for static and dynamic applications [J]. Nonlinear Dynamics, 2015, 81: 77 − 96. doi: 10.1007/s11071-015-1974-8
[35]	夏阳, 廖科. 复杂三维曲梁结构的无闭锁等几何分析算法研究[J]. 工程力学, 2018, 35(11): 17 − 25. doi: 10.6052/j.issn.1000-4750.2017.08.0602 Xia Yang, Liao Ke. Locking-free isogeometric analysis of complex three-dimensional beam structures [J]. Engineering Mechanics, 2018, 35(11): 17 − 25. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.08.0602
[36]	Guo Y J, Zou Z H, Ruess M. Isogeometric multi-patch analyses for mixed thin shells in the framework of non-linear elasticity [J]. Computer Methods in Applied Mechanics and Engineering, 2021, 380: 113771. doi: 10.1016/j.cma.2021.113771
[37]	Guo Y J, Heller J, Hughes T J R, et al. Variationally consistent isogeometric analysis of trimmed thin shells at finite deformations, based on the STEP exchange format [J]. Computer Methods in Applied Mechanics and Engineering, 2018, 336: 39 − 79. doi: 10.1016/j.cma.2018.02.027
[38]	Guo Y J, Do H, Ruess M. Isogeometric stability analysis of thin shells: From simple geometries to engineering models [J]. International Journal for Numerical Methods in Engineering, 2019, 118: 433 − 458. doi: 10.1002/nme.6020
[39]	宋敏. 带参数椭圆型方程的等几何分析及POD模型降阶研究[D]. 浙江: 浙江大学, 2019. Song Min. IGA and proper orthogonal decomposition for parameter dependent elliptic problems [D]. Zhejiang: Zhejiang University, 2019. (in Chinese)
[40]	Manzoni A, Salmoiraghi F, Heltai L. Reduced basis isogeometric methods (RB-IGA) for the real-time simulation of potential flows about parameterized NACA airfoils [J]. Computer Methods in Applied Mechanics and Engineering, 2015, 284: 1147 − 1180. doi: 10.1016/j.cma.2014.11.037
[41]	Ding C S, Deokar R R, Lian H J, et al. Resolving high frequency issues via proper orthogonal decomposition based dynamic isogeometric analysis for structures with dissimilar materials [J]. Computer Methods in Applied Mechanics and Engineering, 2020, 359: 112753. doi: 10.1016/j.cma.2019.112753
[42]	Li S, Trevelyan J, Wu Z, et al. An adaptive SVD-Krylov reduced order model for surrogate based structural shape optimization through isogeometric boundary element method [J]. Computer Methods in Applied Mechanics and Engineering, 2019, 349: 312 − 338. doi: 10.1016/j.cma.2019.02.023
[43]	Khatir S, Wahab M A. Fast simulation for solving fracture mechanics inverse problems using POD-RBF XIGA and Jaya algorithm [J]. Engineering Fracture Mechanics, 2019, 205: 285 − 300. doi: 10.1016/j.engfracmech.2018.09.032
[44]	Luo K, Liu C, Tian Q, et al. An efficient model reduction method for buckling analyses of thin shells based on IGA [J]. Computer Methods in Applied Mechanics and Engineering, 2016, 309: 243 − 268. doi: 10.1016/j.cma.2016.06.006
[45]	Piegl L, Tiller W. The NURBS book – Monographs in visual communication [M]. Heidelberg: Springer-Verlag Publishing Company, 1997.
[46]	De Borst R, Crisfield M A, Remmers J J C, et al. Non-linear finite element analysis of solids and structures [M]. United Kingdom: John Wiley & Sons, 2012.