CN109241559B - 一种基于子结构的复合材料弹性参数识别方法 - Google Patents
一种基于子结构的复合材料弹性参数识别方法 Download PDFInfo
- Publication number
- CN109241559B CN109241559B CN201810860694.4A CN201810860694A CN109241559B CN 109241559 B CN109241559 B CN 109241559B CN 201810860694 A CN201810860694 A CN 201810860694A CN 109241559 B CN109241559 B CN 109241559B
- Authority
- CN
- China
- Prior art keywords
- substructure
- matrix
- composite material
- modal
- residual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 79
- 230000035945 sensitivity Effects 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000005457 optimization Methods 0.000 claims abstract description 8
- 238000006068 polycondensation reaction Methods 0.000 claims abstract description 7
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims description 44
- 238000004088 simulation Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 239000012792 core layer Substances 0.000 description 7
- 238000013461 design Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种基于子结构的复合材料弹性参数识别方法，建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚；缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息；提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度；将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数。本发明通过考虑了子结构的复合材料建模，将模型待识别部分定义为残余结构，通过全模型模态频率对残余结构弹性参数的相对灵敏度分析和模态振型匹配，采用优化迭代算法识别复合材料待识别参数，节省计算资源，提高计算效率，具有十分重要的工程意义。
Description
技术领域
本发明涉及一种复合材料，具体涉及一种复合材料的参数识别方法。
背景技术
复合材料通常是由两种甚至两种以上的不同材料复合而成，且复合材料一般具有各向异性，因其具有高比强、高比模量、可设计性强等优点，广泛应用于航空航天、土木工程、机械工程等领域。为满足工程结构设计需求，需要对复合材料的材料参数给予准确信息，但因制作工艺和制作环境的影响，复合材料成品性能和理论性能相差较大，需要对其弹性参数进行识别，为工程材料设计提供基础信息。
传统的复合材料参数识别过程中，只需对整体结构中材料参数未知部分进行分析，对结构整体的参数识别效率低，产生资源浪费。
发明内容
发明目的：本发明的目的在于针对现有技术的不足，提供一种基于子结构的复合材料弹性参数识别方法，采用子结构理论和基于模态频率的参数识别方法，在保证参数识别精度的同时极大的提高计算效率具有十分重要的工程意义。
技术方案：本发明提供了一种基于子结构的复合材料弹性参数识别方法，包括以下步骤：
(1)建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚；
(2)将缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息；
(3)提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度；
(4)将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数。
进一步，步骤(1)包括以下步骤：
(1.1)建立复合材料精细化有限元模型，施加分析边界条件与载荷；将复合材料结构中材料参数已知部分定义为子结构，将材料参数未知部分定义为残余结构，分别对子结构和残余结构赋予确定性材料参数p和初始参数并将两部分有限元模型分别输出,其中为残余结构待识别的材料参数pr的初始值；
(1.2)根据所给出的材料参数p和利用有限元方法分别求解子结构和残余结构特征矩阵，并得到子结构在物理坐标u下的动力学方程：
式中，M、C、K分别表示子结构的质量矩阵、阻尼矩阵和刚度矩阵,F为子结构在物理坐标下的外载荷；
子结构分为内部自由度o集和边界自由度b集，将子结构动力学方程写成分块矩阵形式：
式中，下角标oo和bb分别表示特征矩阵的内部分块和外部分块，ob和bo表示耦合分块；
根据子结构动力学方程的分块形式，求解子结构的约束模态矩阵ψ和固定界面模态矩阵Φoo：
其中，Ibb为与Φob维度相同的单位矩阵，ωo为固定界面特征方程对应的主模态特征频率；
将约束模态矩阵和固定界面模态矩阵组装为转换矩阵H：
(1.3)复合材料子结构模型动力学缩聚，得到模态坐标下的特征矩阵和动力学方程；由子结构在物理坐标u下运动方程转换到缩减的模态坐标p下的运动方程为：
其中，
式中，分别表示子结构在模态坐标下的质量矩阵，阻尼矩阵和刚度矩阵；p为模态坐标；分别为子结构在模态坐标下的外部载荷。
进一步，步骤(2)利用子结构间位移协调条件及力平衡条件，将所有子结构和残余结构进行模态综合，得到全模型在广义坐标q下的运动方程：
其中，
式中，Mr、Cr、Kr为残余结构特征矩阵，为在广义坐标q下的全模型特征矩阵，为广义坐标下的全模型外部载荷，T为转换矩阵；
根据广义坐标下的运动方程，得到全模型在广义坐标下的特征方程：
根特征方程求解复合材料全模型模态ω。
进一步，步骤(3)包括以下步骤：
(3.1)提取全模型模态频率ω和残余结构模态振型Φr；
(3.2)根据子结构理论和虚功原理确定整体刚度矩阵对残余结构材料弹性参数的偏导：
式中，e表示有限元模型中的单元，为残余结构单元刚度矩阵，pr为残余结构待识别的弹性参数；
其中，为残余结构单元刚度矩阵对材料参数的偏导的累加：
式中，B为单元应变矩阵，为残余结构单元的体积，Dr为弹性系数矩阵；
(3.3)计算全模型模态频率对残余结构材料弹性参数的偏导：
式中，S为模态频率对残余结构待识别参数的灵敏度矩阵；
(3.4)通过无量纲处理得到相对灵敏度矩阵：
S'＝ω-1Spr
式中，S'为相对灵敏度矩阵。
进一步，步骤(4)包括以下步骤：
(4.1)将试验和有限元模拟的残余结构模型模态振型进行MAC值匹配，调整模态频率阶次匹配；
(4.2)将试验和有限元模拟的模态频率残差的二范数作为目标函数：
F(pr)＝||ω*-ω||2
其中，F(pr)为目标函数，ω*表示实验测量的模态频率，ω为有限元求解的模态频率；
(4.3)选择最小二乘法为基于子结构的复合材料弹性参数识别的优化算法，最小化目标函数；当运算结果不满足收敛准则时，更新材料参数，重新进行模态频率对残余结构弹性参数的灵敏度分析迭代计算，当满足收敛准则时，认为目标函数达到最小，迭代终止；材料参数初始值为第k次迭代的材料参数通过求解下式得到：
式中，S'k为第k次迭代过程相对灵敏度，和分别为第k次和第k-1次迭代残余结构材料参数；
收敛准则定义为：
式中，ε为给定的精度要求。
有益效果：本发明通过考虑了子结构的复合材料建模，将模型待识别部分定义为残余结构，通过全模型模态频率对残余结构弹性参数的相对灵敏度分析和模态振型匹配，采用优化迭代算法识别复合材料待识别参数，节省计算资源，提高计算效率，具有十分重要的工程意义。
附图说明
图1为本发明方法流程图；
图2为复合材料有限元全模型；
图3为复合材料子结构模型；
图4为复合材料残余结构模型；
图5为模态频率误差图；
图6为参数识别过程示意图。
具体实施方式
下面对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
一种基于子结构的复合材料弹性参数识别方法，如图1所示，包括以下步骤：
(1)建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚，包括以下步骤：
(1.1)建立复合材料精细化有限元模型，施加分析边界条件与载荷；将复合材料结构中材料参数已知部分定义为子结构，将材料参数未知部分定义为残余结构，分别对子结构和残余结构赋予确定性材料参数和初始参数，并将两部分有限元模型分别输出；其中复合材料整体有限元模型、子结构模型、残余结构模型分别如图2、图3、图4所示，本实施例中不包含阻尼项，且结构处于自由边界条件，无外部载荷；
复合材料整体有限元模型包括上层板、芯层、下层板，均为复合材料；几何参数为90*45*11.5mm，上面板厚度为1mm，芯层厚度为10mm，下面板厚度为0.5mm；
复合材料子结构模型包括上层板、芯层、下层板，均为复合材料，材料参数为确定性正交各向异性材料；几何参数为90*45*11.5mm-30*15*11.5mm，上面板厚度为1mm，芯层厚度为10mm，下面板厚度为0.5mm；
复合材料残余结构模型包括上层板、芯层、下层板，均为复合材料，芯层材料剪切模量G23和G13为待识别参数；几何参数为30*15*11.5mm，上面板厚度为1mm，芯层厚度为10mm，下面板厚度为0.5mm。
(1.2)根据有限元方法分别求解子结构和残余结构特征矩阵，并得到自由状态下无阻尼子结构在物理坐标u下的动力学方程：
式中，M、K分别表示子结构的质量矩阵和刚度矩阵；
子结构分为内部自由度(o集)和边界自由度(b集)，将子结构动力学方程写成分块矩阵形式：
根据子结构动力学方程的分块形式，求解子结构的约束模态矩阵ψ和固定界面模态矩阵Φoo：
其中，Ibb为与Φob维度相同的单位矩阵，ωo为固定界面特征方程对应的主模态特征频率；
将约束模态矩阵和固定界面模态矩阵组装为转换矩阵H：
(1.3)复合材料子结构模型动力学缩聚，得到模态坐标下的特征矩阵和动力学方程；由子结构在物理坐标u下运动方程转换到缩减的模态坐标p下的运动方程为：
其中，
式中，分别表示子结构在模态坐标下的质量矩阵和刚度矩阵；p为模态坐标。
(2)将缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息，包括以下步骤：
(2.1)利用子结构间位移协调条件及力平衡条件，将所有子结构和残余结构进行模态综合，得到全模型在广义坐标q下的运动方程：
其中，
式中，Mr、Kr为残余结构特征矩阵，为在广义坐标q下的全模型特征矩阵，T为转换矩阵；
根据广义坐标下的运动方程，得到全模型在广义坐标下的特征方程：
根据特征方程求解复合材料全模型模态ω。
(3)提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度，包括以下步骤:
(3.1)提取全模型模态频率ω和残余结构模态振型Φr，提取前10阶模态参数。
(3.2)根据子结构理论和虚功原理确定整体刚度对残余结构材料弹性参数的偏导，其中残余结构单元刚度矩阵对材料参数的偏导：
其中，为残余结构单元刚度矩阵，B为单元应变矩阵，为残余结构单元的体积，D为弹性系数矩阵，pr为残余结构待识别的弹性参数；
整体刚度矩阵对残余结构材料弹性参数偏导：
式中e表示有限元模型中的单元。
(3.3)计算全模型模态频率对残余结构材料弹性参数的偏导：
其中，S为模态频率对残余结构待识别参数的灵敏度矩阵。
(3.4)通过无量纲处理得到相对灵敏度矩阵：
S'＝ω-1Spr
其中，S'为相对灵敏度矩阵。
(4)将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数,包括以下步骤：
(4.1)将实验和有限元模拟的残余结构模型模态振型进行MAC值匹配，调整模态频率阶次匹配；残余结构实验模态振型取前3阶。
(4.2)将试验和有限元模拟的模态频率残差的二范数作为目标函数：
F(pr)＝||ω*-ω||2
其中，F(pr)为目标函数，ω*表示实验测量的模态频率，取前3阶，ω为对应有限元求解的模态频率；复合材料前3阶模态频率迭代误差如图5所示。
(4.3)选择最小二乘法为基于子结构的复合材料弹性参数识别的优化算法，最小化目标函数。当运算结果不满足收敛准则时，更新材料参数，重新进行模态频率对残余结构弹性参数的灵敏度分析迭代计算，当满足收敛准则时，认为目标函数达到最小，迭代终止。第k次迭代的材料参数通过求解下式得到：
式中，S'k为第k次迭代过程相对灵敏度，和分别为第k次和第k-1次迭代残余结构材料参数；
收敛准则定义为：
式中，ε为给定的精度要求。
复合材料弹性参数识别结果如图6所示。
Claims (4)
1.一种基于子结构的复合材料弹性参数识别方法，其特征在于：包括以下步骤：
(1)建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚；
(2)将缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息；
(3)提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度；
(4)将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数；
步骤(1)包括以下步骤：
(1.1)建立复合材料精细化有限元模型，施加分析边界条件与载荷；将复合材料结构中材料参数已知部分定义为子结构，将材料参数未知部分定义为残余结构，分别对子结构和残余结构赋予确定性材料参数p和初始参数并将两部分有限元模型分别输出,其中为残余结构待识别的材料参数pr的初始值；
(1.2)根据所给出的材料参数p和利用有限元方法分别求解子结构和残余结构特征矩阵，并得到子结构在物理坐标u下的动力学方程：
式中，M、C、K分别表示子结构的质量矩阵、阻尼矩阵和刚度矩阵,F为子结构在物理坐标下的外载荷；
子结构分为内部自由度o集和边界自由度b集，将子结构动力学方程写成分块矩阵形式：
式中，下角标oo和bb分别表示特征矩阵的内部分块和外部分块，ob和bo表示耦合分块；
根据子结构动力学方程的分块形式，求解子结构的约束模态矩阵ψ和固定界面模态矩阵Φoo：
其中，Ibb为与Φob维度相同的单位矩阵，ωo为固定界面特征方程对应的主模态特征频率；
将约束模态矩阵和固定界面模态矩阵组装为转换矩阵H：
(1.3)复合材料子结构模型动力学缩聚，得到模态坐标下的特征矩阵和动力学方程；由子结构在物理坐标u下运动方程转换到缩减的模态坐标p下的运动方程为：
其中，
式中，分别表示子结构在模态坐标下的质量矩阵，阻尼矩阵和刚度矩阵；p为模态坐标；为子结构在模态坐标下的外部载荷。
2.根据权利要求1所述的基于子结构的复合材料弹性参数识别方法，其特征在于：步骤(2)利用子结构间位移协调条件及力平衡条件，将所有子结构和残余结构进行模态综合，得到全模型在广义坐标q下的运动方程：
其中，
式中，Mr、Cr、Kr为残余结构特征矩阵，为在广义坐标q下的全模型特征矩阵，为广义坐标下的全模型外部载荷，T为转换矩阵；
根据广义坐标下的运动方程，得到全模型在广义坐标下的特征方程：
根据特征方程求解复合材料全模型模态频率ω。
3.根据权利要求2所述的基于子结构的复合材料弹性参数识别方法，其特征在于：步骤(3)包括以下步骤：
(3.1)提取全模型模态频率ω和残余结构模态振型Φr；
(3.2)根据子结构理论和虚功原理确定整体刚度矩阵对残余结构材料弹性参数的偏导：
式中，e表示有限元模型中的单元，为残余结构单元刚度矩阵，pr为残余结构待识别的弹性参数；
其中，为残余结构单元刚度矩阵对材料参数的偏导的累加：
式中，B为单元应变矩阵，为残余结构单元的体积，Dr为弹性系数矩阵；
(3.3)计算全模型模态频率对残余结构材料弹性参数的偏导：
式中，S为模态频率对残余结构待识别参数的灵敏度矩阵；
(3.4)通过无量纲处理得到相对灵敏度矩阵：
S'＝ω-1Spr
式中，S'为相对灵敏度矩阵。
4.根据权利要求3所述的基于子结构的复合材料弹性参数识别方法，其特征在于：步骤(4)包括以下步骤：
(4.1)将试验和有限元模拟的残余结构模型模态振型进行MAC值匹配，调整模态频率阶次匹配；
(4.2)将试验和有限元模拟的模态频率残差的二范数作为目标函数：
F(pr)＝||ω*-ω||2
其中，F(pr)为目标函数，ω*表示实验测量的模态频率，ω为全模型模态频率；
(4.3)选择最小二乘法为基于子结构的复合材料弹性参数识别的优化算法，最小化目标函数；当运算结果不满足收敛准则时，更新材料参数，重新进行模态频率对残余结构弹性参数的灵敏度分析迭代计算，当满足收敛准则时，认为目标函数达到最小，迭代终止；材料参数初始值为第k次迭代的材料参数通过求解下式得到：
式中，S'k为第k次迭代过程相对灵敏度，和分别为第k次和第k-1次迭代残余结构材料参数；
收敛准则定义为：
式中，ε为给定的精度要求。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810860694.4A CN109241559B (zh) | 2018-08-01 | 2018-08-01 | 一种基于子结构的复合材料弹性参数识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810860694.4A CN109241559B (zh) | 2018-08-01 | 2018-08-01 | 一种基于子结构的复合材料弹性参数识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109241559A CN109241559A (zh) | 2019-01-18 |
| CN109241559B true CN109241559B (zh) | 2019-06-18 |
Family
ID=65073332
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810860694.4A Active CN109241559B (zh) | 2018-08-01 | 2018-08-01 | 一种基于子结构的复合材料弹性参数识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109241559B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110008530B (zh) * | 2019-03-15 | 2020-04-24 | 东南大学 | 一种空间柔性复合材料分布式概率建模方法 |
| CN110188417B (zh) * | 2019-05-14 | 2023-03-24 | 东南大学 | 一种基于多级超单元的复杂结构有限元模型修正方法 |
| CN110837714A (zh) * | 2019-11-14 | 2020-02-25 | 宁波财经学院 | 用于智能机器人的压电材料全矩阵参数识别方法 |
| CN111274704B (zh) * | 2020-01-20 | 2022-04-15 | 湖北工业大学 | 一种添加子结构后结构动力学特性重分析方法 |
| TWI758953B (zh) * | 2020-11-16 | 2022-03-21 | 中國鋼鐵股份有限公司 | 耐火泥彈性模數的測量方法 |
| CN112861404B (zh) * | 2021-02-07 | 2023-08-15 | 中国汽车工程研究院股份有限公司 | 一种基于最小二乘法的材料弹性模量参数的求解方法 |
| CN113111547B (zh) * | 2021-03-15 | 2024-07-05 | 南京航空航天大学 | 基于缩减基的频域有限元模型修正方法 |
| CN114970292A (zh) * | 2022-07-29 | 2022-08-30 | 吉林大学 | 一种数字孪生驱动增材制造的仿生产品建模系统及方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
| CN107657073A (zh) * | 2017-08-18 | 2018-02-02 | 东南大学 | 一种三明治复合材料结构的参数识别方法 |
| CN107862170A (zh) * | 2017-12-25 | 2018-03-30 | 华中科技大学 | 一种基于动态缩聚的有限元模型修正方法 |
-
2018
- 2018-08-01 CN CN201810860694.4A patent/CN109241559B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
| CN107657073A (zh) * | 2017-08-18 | 2018-02-02 | 东南大学 | 一种三明治复合材料结构的参数识别方法 |
| CN107862170A (zh) * | 2017-12-25 | 2018-03-30 | 华中科技大学 | 一种基于动态缩聚的有限元模型修正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于灵敏度分析的复合材料组分参数识别方法;秦福溶 等;《复合材料学报》;20180212;全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109241559A (zh) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109241559B (zh) | 一种基于子结构的复合材料弹性参数识别方法 | |
| Gravenkamp et al. | A numerical approach for the computation of dispersion relations for plate structures using the scaled boundary finite element method | |
| WO2019011026A1 (zh) | 基于聚类分析的复合材料结构有限元模型修正方法 | |
| Dey et al. | Free vibration analyses of multiple delaminated angle-ply composite conical shells–A finite element approach | |
| CN104149340B (zh) | 基于3d打印技术获得金属零件性能的方法 | |
| CN115630558B (zh) | 一种复合材料构件装配变形预测方法 | |
| CN104484502A (zh) | 一种基于正向子结构的有限元模型修正方法 | |
| Grover et al. | An efficient C0 finite element modeling of an inverse hyperbolic shear deformation theory for the flexural and stability analysis of laminated composite and sandwich plates | |
| CN111159934B (zh) | 一种铝基复合材料动力学模拟方法 | |
| CN116011301B (zh) | B样条等几何状态空间有限元方法 | |
| Jiang et al. | An edge-based smoothed XFEM for fracture in composite materials | |
| Huang et al. | Three-dimensional vibrations of cracked rectangular parallelepipeds of functionally graded material | |
| CN117556568B (zh) | 一种考虑装配偏差的复合材料螺接结构建模方法 | |
| Yuan et al. | An equivalent modeling method for honeycomb sandwich structure based on orthogonal anisotropic solid element | |
| Sahoo | Free vibration behavior of laminated composite stiffened elliptic parabolic shell panel with cutout | |
| Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
| CN111027254A (zh) | 一种用于ecc双轴受压分析的本构模型构建方法 | |
| Jaworski et al. | Implementation features of composite materials effective mechanical characteristics finding method based on microlevel cellular structural models | |
| CN106650073B (zh) | 一种平衡力系边界条件子模型分析方法 | |
| Haldar et al. | Free vibration of skew laminates–A brief review and some benchmark results | |
| Perez et al. | Prediction of displacement and stress fields of a notched panel with geometric nonlinearity by reduced order modeling | |
| Michopoulos et al. | Advances on the constitutive characterization of composites via multiaxial robotic testing and design optimization | |
| Liu et al. | Base force element method (BFEM) on complementary energy principle for linear elasticity problem | |
| CN102930161A (zh) | 基于多舱段航天器混合模态预示方法 | |
| Shah et al. | Comparative Study of Moving Least Square and Point Interpolation Meshless Technique for Layered Composite Beam Subjected to Transverse Loading. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |