CN107609221B - 一种基于遗传算法的含铰结构非线性参数识别方法 - Google Patents
一种基于遗传算法的含铰结构非线性参数识别方法 Download PDFInfo
- Publication number
- CN107609221B CN107609221B CN201710696203.2A CN201710696203A CN107609221B CN 107609221 B CN107609221 B CN 107609221B CN 201710696203 A CN201710696203 A CN 201710696203A CN 107609221 B CN107609221 B CN 107609221B
- Authority
- CN
- China
- Prior art keywords
- nonlinear
- hinged structure
- hinge
- degree
- freedom
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Genetics & Genomics (AREA)
- Evolutionary Biology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于遗传算法的含铰结构非线性参数识别方法,包括用一种非线性模型来描述含铰结构中铰链的非线性特性,建立铰链结构的动力学模型,得到含铰结构的M、K、C方程;将含铰结构的M、K、C方程转化为非线性的代数方程组;利用数值分析中的Newton迭代法对得到的代数方程组进行求解,得到含铰结构的频响数据;根据求解的频响数据与试验获得的频响数据构建目标函数,通过Matlab中的GA遗传算法进行全局化搜索,最终识别出最优的非线性参数值。本发明的含铰结构非线性参数的识别方法,可以识别含铰结构的非线性特性,克服了传统方法需要提供经验初值及只能保证局部最优的局限性。
Description
技术领域
本发明涉及含铰结构参数的识别方法,特别是涉及一种遗传算法的含铰结构中非线性参数识别方法。
背景技术
随着航天事业的迅速发展,可展结构在空间任务中得到了广泛应用,如太阳帆、太阳能电池阵和空间天线支撑机构等。
可展结构各部件之间通过铰链来连接,但铰链的构造复杂,由于间隙的存在,在外界荷载的作用下,部件之间会发生相互接触、摩擦和碰撞等非线性现象,使结构的动力学特性具有较强的非线性特征。因此,铰链对结构动力学特性的影响不可忽视,需要对铰链的力学特性展开研究,为设计提供可靠的评价和参考。
为了研究含铰结构的非线性特性,需要获得准确的非线性参数值。在工程中,许多非线性因素无法用理论或者试验的方法直接得到,需要根据试验数据反推结构的非线性参数值。动力学***的逆向问题——参数识别是现代控制和优化设计的关键环节,因此今年来受到国内外学者的广泛关注,并相继提出了各种参数识别方法,但是这些方法大多存在提供非线性参数的经验初值,且极易陷入局部极小的陷阱而出现“死循环”现象,使得迭代无法进行,以及只能保证参数值局部最优等局限性。
发明内容
发明目的:提供一种可有效识别含铰结构中铰链的非线性参数,克服了传统方法存在的需要给出经验初值及只能保证局部最优的局限性的基于遗传算法的含铰结构非线性参数识别方法。
技术方案:一种基于遗传算法的含铰结构非线性参数识别方法,包括以下步骤:
(1)用一种非线性模型来描述含铰结构中铰链的非线性特性,建立铰链结构的动力学模型,得到含铰结构的M、K、C方程;
(2)通过谐波平衡法对含铰结构中各自由度的位移、速度、加速度以及铰链非线性恢复力进行一次谐波展开,将含铰结构的M、K、C方程转化为非线性的代数方程组;
(3)利用数值分析中的Newton迭代法对步骤(2)中得到的代数方程组进行求解,得到含铰结构的频响数据;
(4)根据步骤(3)中求解的频响数据与试验获得的频响数据构建目标函数,通过Matlab中的GA遗传算法进行全局化搜索,最终识别出最优的非线性参数值。
进一步的,在步骤(1)中,包含以下步骤:
(11)根据实际情况选择合适的非线性模型来描述含铰结构中铰链的非线性特性,其中,根据铰链非线性恢复力与铰链部件相对位移xj的关系,用间隙非线性模型、立方非线性模型或分段线性模型来描述铰链的非线性特性,其中相对位移xj指的是铰链两侧自由度的位移差;
(12)对含铰结构进行动力学建模,得到其矩阵式动力学方程:
其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵,q为含铰结构的各自由度位移向量,为q的一阶导数,即含铰结构的各自由度速度向量,为q的二阶导数,即含铰结构的各自由度加速度向量,FNL为铰链引入的非线性向量,F为激励广义力向量。
更进一步的,所述步骤(11)中间隙非线性模型,由于铰链中存在间隙,在尚未发生接触时不存在刚度,非线性恢复力和位移表现为***特性,则铰链的非线性恢复力表达式为:
其中,f为非线性恢复力,xj为相对位移,δ为铰链的间隙,k为铰链的线性接触刚度。
更进一步的,所述步骤(11)中立方非线性模型,由非线性接触引起的铰链非线性恢复力具有奇函数特性,且有硬弹簧特性,因此采用立方非线性模型来描述铰链的力位移关系,则表达式为:
其中,FNL为非线性恢复力,xj为相对位移,kcs为刚度系数。
更进一步的,所述步骤(11)中分段线性模型,由于铰链中存在间隙,构件间的接触会导致刚度的变化,可以用分段线性模型来描述,则表达式为:
其中,f为非线性恢复力,ε为变刚度阀值,K1和K2为分段接触刚度,当位移小于等于ε时,接触刚度为K1,当位移大于ε,接触刚度为K2。
进一步的,所述步骤(2)中包括:
(21)基于谐波平衡理论,对***的响应和铰链的非线性恢复力进行一次谐波分解,将n自由度对应的位移、速度、加速度以及铰链的非线性恢复力表示为:
qn=an sinωt+bn cosωt
其中,an和bn分别为谐波展开后各自由度对应的正弦系数和余弦系数,jp和jq分别为铰链非线性恢复力谐波展开后的正弦系数和余弦系数,f为铰链非线性恢复力,qn、和分别表示n自由度对应的位移、速度、加速度,ω表示频率。
(22)将步骤(21)中的公式带入步骤(12)的动力学方程,并假设外部激励力为Fsinωt,根据方程两边正弦系数和余弦系数分别相等,将步骤(12)中的动力学方程转化为代数方程组:
其中,Cp为铰链非线性恢复力谐波展开后的正弦系数对应的矩阵,Cq为铰链非线性恢复力谐波展开后的余弦系数对应的矩阵,令铰链连接的自由度分别为n1和n2,有:
其中,A为各自由度对应的正弦系数向量,B为各自由度对应的余弦系数向量,有:
A=[a1 a2 … an],B=[b1 b2 … bn]
进一步的,所述步骤(3)包括:
(31)将步骤(22)中得到的代数方程组进行一个等式变换,得到下式:
(32)可以将步骤(31)中的等式视为f(x)=0,其中为由各自由度对应的正余弦向量组成的列向量,由此可以将步骤(31)中的等式构造成Newton迭代格式:
其中,it表示迭代次数,迭代求解时综合运算效率及计算精度选择合适的迭代次数,即it值。
进一步的,所述步骤(4)中包括:
(41)定义残差项为结构的位移响应的试验值与谐波平衡法计算值的差值:
R(p)=xe(ω)-xa(ω,p)
其中,p代表待识别的铰链非线性参数,xe和xa分别代表位移响应的试验值和计算值,试验中,并不能够获得幅频曲线的解析表达式,而是以分辨率表示的离散数据点,假设试验选取的频率点为ω1,ω2,...,ωn,并根据最小二乘法原则,若要计算所得的频响数据能够最佳拟合试验结果,定义目标函数为:
(42)利用GA遗传算法,对非线性参数进行全局搜索,直到识别出最优的非线性参数值pa,使得目标函数R(p)取极小值,则pa为非线性参数的精确值。
有益效果:与现有技术相比,本发明的含铰结构非线性参数的识别方法,可以准确的识别含铰结构的非线性参数值,由于遗传算法具有全局寻优的优势,克服了传统方法需要提供经验初值及只能保证局部最优的局限性。
附图说明
图1是基于遗传算法的含铰结构非线性参数识别流程图;
图2是间隙非线性模型示意图;
图3是立方非线性模型示意图;
图4是分段线性模型示意图;
图5是铰梁结构动力学模型示意图;
图6是非线性参数值n迭代情况示意图;
图7是识别后的模型与试验结构幅频曲线图。
具体实施方式
下面结合附图对本发明技术方案进行详细说明。
如图1所示,本发明的一种基于遗传算法的含铰结构非线性参数识别方法,本发明以铰梁结构为例进行说明,其中,梁截面为圆形,直径D为0.02m,每根梁长度le为1m,总长度L为2m,材料为铝,弹性模量E为70Gpa,密度ρ为2700kg/m3,令梁结构中的阻尼为比例阻尼,质量矩阵阻尼系数CM为0.0001,刚度矩阵阻尼系数CK为0.0001,铰链对应的广义位移为转动自由度,为方便研究不同的非线性参数和激振力幅值对结构幅频特性的影响,令铰链刚度系数kcs=nEI/le,激励力幅值Fc=0.5EI/L2。
该方法包括以下步骤:
步骤1:假定一种非线性模型来描述含铰结构中铰链的非线性特性,建立铰链结构的动力学模型,并最终得到含铰结构的M、K、C方程。其中,确定铰梁结构动力学方程包含以下步骤:
步骤1.1:根据实际情况选择合适的非线性模型来描述含铰结构中铰链的非线性特性,一般地,根据铰链非线性恢复力与铰链部件相对位移xj的关系,可以有三种非线性模型来描述铰链的非线性特性;间隙非线性模型、立方非线性模型或分段线性模型,其中相对位移xj指的是铰链两侧自由度的位移差。
(1)间隙非线性模型:由于铰链中存在间隙,在尚未发生接触时不存在刚度,非线性恢复力和位移表现为***特性,如图2所示,则铰链的非线性恢复力表达式为:
其中,f为非线性恢复力,xj为相对位移,δ为铰链的间隙,k为铰链的线性接触刚度。
(2)立方非线性模型:由非线性接触引起的铰链非线性恢复力具有奇函数特性,且有硬弹簧特性,因此采用立方非线性模型来描述铰链的力位移关系,如图3所示,则表达式为:
其中,f为非线性恢复力,xj为相对位移,kcs为刚度系数。
(3)分段线性模型:由于铰链中存在间隙,构件间的接触会导致刚度的变化,可以用分段线性模型来描述,如图4所示,则表达式为:
其中,f为非线性恢复力,ε为变刚度阀值,K1和K2为分段接触刚度,当位移小于等于ε时,接触刚度为K1,当位移大于ε,接触刚度为K2。
本实施例采用的铰链非线性模型为立方非线性模型。
步骤1.2:以图5的铰梁结构为对象,对含铰结构进行动力学建模,得到其矩阵形式的结构动力学方程为:
其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵,q为含铰结构的各自由度位移向量,为q的一阶导数,即含铰结构的各自由度速度向量,为q的二阶导数,即含铰结构的各自由度加速度向量,FNL为铰链引入的非线性向量,即,由铰链非线性恢复力f引入的非线性向量,其表达式为:F为激励广义力向量,代表的是由作用在每个自由度上面的激励力幅值构成的列向量。
梁单元的单元刚度矩阵和单元质量矩阵为:
其中EI为梁的抗弯刚度,le为单元长度,ρ为材料的密度,对其进行刚度组装,可以得到整体刚度矩阵和整体质量矩阵为:
阻尼矩阵可以由下面公式得到:
C=CM×M+CK×K (9)
如图5所示,假设激励力作用位置为6自由度,运用上述方法求解6自由度处的频域响应,铰链非线性弯矩转角关系中对应的转角应为4,5自由度的转动位移差,即位移幅值假设试验的频响数据为激振力幅值0.5Fc,铰链刚度为0.1kcs情况下获得。
步骤2:通过谐波平衡法对含铰结构中各自由度的位移、速度、加速度以及铰链非线性恢复力进行一次谐波展开,将含铰结构的M、K、C方程转化为非线性的代数方程组。其中,谐波平衡法求解动力学方程包含以下步骤:
步骤2.1基于谐波平衡理论,对***的响应和铰链的非线性恢复力进行一次谐波分解,***的响应指的即是含铰结构中各自由度的位移、速度和加速度。将n自由度对应的位移、速度、加速度以及铰链的非线性恢复力表示为:
其中,an和bn分别为谐波展开后各自由度对应的正弦系数和余弦系数,jp和jq分别为铰链非线性恢复力谐波展开后的正弦系数和余弦系数,f为铰链非线性恢复力,qn、和分别表示n自由度对应的位移、速度、加速度,ω表示频率。
铰链立方非线性模型的一次谐波展开系数为:
步骤2.2:将式(12)带入式(11)的动力学方程,并假设外部激励力为Fsinωt,带入式(4)根据方程两边正弦系数和余弦系数分别相等,将式(4)动力学方程转化为代数方程组:
其中,Cp为铰链非线性恢复力谐波展开后的正弦系数对应的矩阵,Cq为铰链非线性恢复力谐波展开后的余弦系数对应的矩阵,令铰链连接的自由度分别为n1和n2,有:
其中,A为各自由度对应的正弦系数向量,B为各自由度对应的余弦系数向量,有:
A=[a1 a2 … an],B=[b1 b2 … bn] (16)
至此,将结构动力学方程组转化为了非线性代数方程组。
步骤3:利用数值分析中的Newton迭代法对步骤(2)中得到的代数方程组进行求解,得到含铰结构的频响数据。Newton迭代法求解非线性代数方程组包含以下步骤:
步骤3.1:将步骤(22)中得到的代数方程组进行一个等式变换,得到下式:
将上式构造成Newton迭代格式并利用Newton迭代法求解,最终求得结构的频响数据。
步骤3.2:可以将式(17)中的等式视为f(x)=0,其中为由各自由度对应的正余弦向量组成的列向量,由此可以将式(17)的等式构造成Newton迭代格式:
其中,it表示迭代次数,迭代求解时综合运算效率及计算精度选择合适的迭代次数,即it值。
步骤4:根据数值方法求解的频响数据与试验获得的频响数据的比对结果,通过Matlab中的GA遗传算法进行全局化搜索,最终识别出最优的非线性参数值。
GA遗传算法进行全局最优搜索包含以下步骤:
步骤4.1:为了识别含铰结构的非线性参数值,即铰梁刚度值kcs=nEI/le中的n,此处的非线性参数n为本实施例中具体的非线性参数值,定义适应度函数,即目标函数,有:
步骤4.2:基于GA遗传算法,将非线性参数n作为GA的个体,为了减少搜索时间,定义个体搜索上下限[0.05,0.15],全局搜索最小适应度值,此时的个体即为识别出的非线性参数值。图6和图7分别为非线性参数值迭代曲线图及识别后的模型与试验结构幅频曲线,由图6可知,纵坐标为“非线性参数值”,横坐标为“代数”,表示迭代次数,因为遗传算法是模仿生物遗传得到的优化算法,因此习惯将迭代次数叫“代数”,类似于生物学里面父代和子代。可以看到,经过520代的振荡搜索后,非线性参数值逐渐趋于一个稳定值,即识别后的铰链非线性参数值,为0.1kcs,同时将识别后的结果与试验结果进行对比,如图7所示,可以看到两个结果吻合的很好,说明了识别后的参数值的准确性。
可见本发明公开的含铰结构非线性参数识别方法可以有效识别含铰结构中铰链的非线性参数值。
Claims (7)
1.一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,包括以下步骤:
(1)用一种非线性模型来描述含铰结构中铰链的非线性特性,建立铰链结构的动力学模型,得到含铰结构的M、K、C方程,其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵;
(2)通过谐波平衡法对含铰结构中各自由度的位移、速度、加速度以及铰链非线性恢复力进行一次谐波展开,将含铰结构的M、K、C方程转化为非线性的代数方程组;
(3)利用数值分析中的Newton迭代法对步骤(2)中得到的代数方程组进行求解,得到含铰结构的频响数据;
(4)根据步骤(3)中求解的频响数据与试验获得的频响数据构建目标函数,通过Matlab中的GA遗传算法进行全局化搜索,最终识别出最优的非线性参数值;具体为:
(41)定义残差项为结构的位移响应的试验值与谐波平衡法计算值的差值:
R(p)=xe(ω)-xa(ω,p)
其中,p代表待识别的铰链非线性参数,xe和xa分别代表位移响应的试验值和计算值,试验中,并不能够获得幅频曲线的解析表达式,而是以分辨率表示的离散数据点,假设试验选取的频率点为ω1,ω2,…,ωn,并根据最小二乘法原则,若要计算所得的频响数据能够最佳拟合试验结果,定义目标函数为:
(42)利用GA遗传算法,对非线性参数进行全局搜索,直到识别出最优的非线性参数值pa,使得目标函数R(p)取极小值,则pa为非线性参数的精确值。
2.根据权利要求1所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于:在步骤(1)中,包含以下步骤:
(11)根据实际情况选择非线性模型来描述含铰结构中铰链的非线性特性,其中,根据铰链非线性恢复力与铰链部件相对位移xj的关系,用间隙非线性模型、立方非线性模型或分段线性模型来描述铰链的非线性特性,其中相对位移xj指的是铰链两侧自由度的位移差;
(12)对含铰结构进行动力学建模,得到其矩阵式动力学方程:
其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵,q为含铰结构的各自由度位移向量,为q的一阶导数,即含铰结构的各自由度速度向量,为q的二阶导数,即含铰结构的各自由度加速度向量,FNL为铰链引入的非线性向量,F为激励广义力向量。
3.根据根据权利要求2所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,所述步骤(11)中间隙非线性模型,由于铰链中存在间隙,在尚未发生接触时不存在刚度,非线性恢复力和位移表现为***特性,则铰链的非线性恢复力表达式为:
其中,f为非线性恢复力,xj为相对位移,δ为铰链的间隙,k为铰链的线性接触刚度。
4.根据权利要求2所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,所述步骤(11)中立方非线性模型,由非线性接触引起的铰链非线性恢复力具有奇函数特性,且有硬弹簧特性,因此采用立方非线性模型来描述铰链的力位移关系,则表达式为:
其中,f为非线性恢复力,xj为相对位移,kcs为刚度系数。
5.根据权利要求2所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,所述步骤(11)中分段线性模型,由于铰链中存在间隙,构件间的接触会导致刚度的变化,可以用分段线性模型来描述,则表达式为:
其中,f为非线性恢复力,ε为变刚度阀值,K1和K2为分段接触刚度,当位移小于等于ε时,接触刚度为K1,当位移大于ε,接触刚度为K2,xj为相对位移。
6.根据权利要求1所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,所述步骤(2)中包括:
(21)基于谐波平衡理论,对***的响应和铰链的非线性恢复力进行一次谐波分解,将n自由度对应的位移、速度、加速度以及铰链的非线性恢复力表示为:
qn=α n sinωt+bn cosωt
其中,an和bn分别为谐波展开后各自由度对应的正弦系数和余弦系数,jp和jq分别为铰链非线性恢复力谐波展开后的正弦系数和余弦系数,f为铰链非线性恢复力,qn、和分别表示n自由度对应的位移、速度、加速度,ω表示频率;
(22)将步骤(21)中的公式带入动力学方程中,并假设外部激励力为Fsinωt,根据方程两边正弦系数和余弦系数分别相等,将动力学方程转化为代数方程组为:
其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵,F为激励广义力向量,Cp为铰链非线性恢复力谐波展开后的正弦系数对应的矩阵,Cq为铰链非线性恢复力谐波展开后的余弦系数对应的矩阵,q为含铰结构的各自由度位移向量,为q的一阶导数,即含铰结构的各自由度速度向量,为q的二阶导数,即含铰结构的各自由度加速度向量,FNL为铰链引入的非线性向量;令铰链连接的自由度分别为n1和n2,有:
其中,A为各自由度对应的正弦系数向量,B为各自由度对应的余弦系数向量,有:
A=[a1 a2 … an],B=[b1 b2 … bn]。
7.根据权利要求1所述的一种基于遗传算法的含铰结构非线性参数识别方法,其特征在于,所述步骤(3)包括:
(31)将代数方程组:
进行一个等式变换,得到下式:
其中,M为含铰结构的质量矩阵,C为含铰结构的阻尼矩阵,K为含铰结构的刚度矩阵,F为激励广义力向量,Cp为铰链非线性恢复力谐波展开后的正弦系数对应的矩阵,Cq为铰链非线性恢复力谐波展开后的余弦系数对应的矩阵,ω表示频率;
(32)将步骤(31)中的等式视为f(x)=0,其中为由各自由度对应的正余弦向量组成的列向量,由此可以将步骤(31)中的等式构造成Newton迭代格式:
其中,it表示迭代次数,迭代求解时综合运算效率及计算精度选择迭代次数,即it值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710696203.2A CN107609221B (zh) | 2017-08-15 | 2017-08-15 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
PCT/CN2018/083369 WO2019033771A1 (zh) | 2017-08-15 | 2018-04-17 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710696203.2A CN107609221B (zh) | 2017-08-15 | 2017-08-15 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107609221A CN107609221A (zh) | 2018-01-19 |
CN107609221B true CN107609221B (zh) | 2018-07-17 |
Family
ID=61065125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710696203.2A Active CN107609221B (zh) | 2017-08-15 | 2017-08-15 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107609221B (zh) |
WO (1) | WO2019033771A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107609221B (zh) * | 2017-08-15 | 2018-07-17 | 东南大学 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
CN108573111A (zh) * | 2018-04-27 | 2018-09-25 | 福州大学 | 一种外压容器设计中解析新算法参数的拟合方法 |
CN109190328B (zh) * | 2018-11-27 | 2020-03-31 | 东南大学 | 一种混合有限-复变差分的多参数结构动响应灵敏度分析方法 |
CN109614720B (zh) * | 2018-12-14 | 2020-07-28 | 厦门大学 | 在役隔震/减振装置的无模型非线性特性实时识别方法 |
CN109902357B (zh) * | 2019-01-31 | 2020-06-02 | 东南大学 | 一种复变差分的光滑非线性结构动响应灵敏度分析方法 |
CN110377965A (zh) * | 2019-06-26 | 2019-10-25 | 东南大学 | 一种含铰柔性结构非线性特征的辨识方法 |
CN110569540A (zh) * | 2019-07-31 | 2019-12-13 | 西北工业大学 | 一种舵机传动机构动力学分析方法 |
CN111144206B (zh) * | 2019-11-21 | 2021-02-09 | 东南大学 | 一种柔性结构中立方非线性***识别方法 |
CN111144205B (zh) * | 2019-11-21 | 2021-01-05 | 东南大学 | 一种航天器结构中间隙非线性***识别方法 |
CN112180762B (zh) * | 2020-09-29 | 2021-10-29 | 瑞声新能源发展(常州)有限公司科教城分公司 | 非线性信号***构建方法、装置、设备和介质 |
CN112487577B (zh) * | 2020-11-27 | 2023-06-02 | 重庆大学 | 一种结构稳态非线性动力学响应快速计算方法及***和存储介质 |
CN112784378A (zh) * | 2021-03-01 | 2021-05-11 | 西安交通大学 | 一种螺栓结合部非线性连接层建模方法 |
CN113297908B (zh) * | 2021-04-23 | 2024-05-03 | 东南大学 | 一种脉冲激励下基于模式搜索的非线性刚度识别方法 |
CN113297907B (zh) * | 2021-04-23 | 2024-04-30 | 东南大学 | 一种脉冲激励下基于数据驱动的非线性阻尼识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950317A (zh) * | 2010-09-03 | 2011-01-19 | 清华大学 | 基于模态分割和遗传算法的飞行器定阶次参数模型辨识方法 |
CN102540894A (zh) * | 2012-02-17 | 2012-07-04 | 南京电力设备质量性能检验中心 | 一种基于遗传算法的带有未知负载的机械手参数辨识方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004265069A (ja) * | 2003-02-28 | 2004-09-24 | Yaskawa Electric Corp | 仮想受動関節モデルのモデルパラメータ同定方法およびその制御方法 |
US9260147B2 (en) * | 2013-06-26 | 2016-02-16 | Wisconsin Alumni Research Foundation | Dynamic predictor for articulated mechanisms |
CN105243231B (zh) * | 2015-11-02 | 2019-06-25 | 广东工业大学 | 高速轻载机构非线性动态***结构拓扑参数混合优化方法 |
CN106021715B (zh) * | 2016-05-18 | 2019-05-24 | 中南大学 | 一种机构的转动副铰接间隙磨损的敏感性指标的确定方法 |
CN106407582B (zh) * | 2016-09-28 | 2019-10-01 | 中国工程物理研究院总体工程研究所 | 一种用于柔性机构静力学和动力学分析的高效计算方法 |
CN107609221B (zh) * | 2017-08-15 | 2018-07-17 | 东南大学 | 一种基于遗传算法的含铰结构非线性参数识别方法 |
-
2017
- 2017-08-15 CN CN201710696203.2A patent/CN107609221B/zh active Active
-
2018
- 2018-04-17 WO PCT/CN2018/083369 patent/WO2019033771A1/zh active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950317A (zh) * | 2010-09-03 | 2011-01-19 | 清华大学 | 基于模态分割和遗传算法的飞行器定阶次参数模型辨识方法 |
CN102540894A (zh) * | 2012-02-17 | 2012-07-04 | 南京电力设备质量性能检验中心 | 一种基于遗传算法的带有未知负载的机械手参数辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107609221A (zh) | 2018-01-19 |
WO2019033771A1 (zh) | 2019-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609221B (zh) | 一种基于遗传算法的含铰结构非线性参数识别方法 | |
Harne et al. | Harnessing bistable structural dynamics: for vibration control, energy harvesting and sensing | |
Kiani | Analysis of FG-CNT reinforced composite conical panel subjected to moving load using Ritz method | |
Martin et al. | Methodology for wind/wave basin testing of floating offshore wind turbines | |
Ghafoori et al. | Dynamic analysis of laminated composite plates traversed by a moving mass based on a first-order theory | |
CN106844914A (zh) | 一种空天飞行器机翼振动响应的快速仿真方法 | |
Raja et al. | Active control of tensegrity structures under random excitation | |
Saboori et al. | Static analysis of tapered FRP transmission poles using finite element method | |
CN109398020B (zh) | 一种基于非线性模型的车辆液电耦合式isd悬架的预测控制方法 | |
Shahgholi et al. | Nonlinear vibration analysis of a spinning shaft with multi-disks | |
CN105426641A (zh) | 基于大质量法的简谐加速度激励下结构拓扑优化设计方法 | |
Yang et al. | LPV control of a tensegrity-membrane system | |
Funes et al. | An efficient dynamic formulation for solving rigid and flexible multibody systems based on semirecursive method and implicit integration | |
Suresh et al. | Vibration and damping analysis of thin-walled box beams | |
Stanford et al. | Cost reduction techniques for the design of non‐linear flapping wing structures | |
Liu et al. | Application of artificial intelligence and numerical method to improve the natural oscillation of the sport mountain bike tires | |
Xu et al. | Magnetoelastic nonlinear free vibration of an annular plate under a nonuniform magnetic field of the long straight current‐carrying wire | |
CN113919104A (zh) | 旋转鼓筒的非线性动态响应获取方法及*** | |
Magri et al. | Example of a non-smooth Hopf bifurcation in an aero-elastic system | |
Sylla et al. | Dynamics of a rotating flexible and symmetric spacecraft using impedance matrix in terms of the flexible appendages cantilever modes | |
Suryanarayanan et al. | A procedure for the development of control-oriented linear models for horizontal-axis large wind turbines | |
Mottaghi et al. | Multiscale Analysis of Carbon Nanotube-Reinforced Curved Beams: A Finite Element Approach Coupled with Multilayer Perceptron Neural Network | |
Merz | STAS Aeroelastic 1.0-Theory Manual | |
Takarics et al. | Friction compensation of an aeroelastic wing—A TP model transformation based approach | |
Afanaseva et al. | Aeroelastic Vibrations of a Thin Ribbon in a Laminar Air Flow |
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 |