01
结构动力优化设计简介
工程结构在其工作期间将经常性地受到风浪、地震、振动设备等各种外界激励的影响,本质上讲,工程结构所受到的环境载荷都是动载荷,其正常服役性能会受到动力破坏或损伤的影响。船舶在正常服役中会受到主机、螺旋桨的振动和波浪的冲击等外界载荷的作用,海洋平台在其使用期间也要受到复杂的环境荷载和各种振动设备的影响,船舶、海洋平台结构在这些外界载荷的作用下,将会发生不同程度的振动。如若振动过大,则可导致结构产生疲劳破坏,相关人员的工作效率也会受到影响,甚至还会损害到人员的身体健康,仪器和设备的正常使用也很难以得到保证。欲使结构和设备在动力载荷作用下的相应力学特性得到改善,使其振动得到有效控制和在正常工作时能一直保持较好的性能,对其进行结构动力优化设计是非常有效的方法。
结构的动力优化设计相对于结构的动力分析来说,是其反问题或逆问题,相比正问题而言,其理论求解要复杂和困难得多。在结构动力优化设计中,同时涉及到结构动力学基本理论、灵敏度分析、结构有限元建模、结构重分析、数学优化算法、软件编程和实际工程应用等诸多方面,是诸多学科相互交叉、有机结合的产物。它包括近代的结构动力学、数学优化方法、数值计算方法、结构有限元法以及程序设计等,是现代工程结构设计领域中的一个新兴的分支。结构动力优化设计的主要任务是:根据在实际中对结构要求的不同,构建相应的结构动力分析与优化设计理论模型,把相关算法语言、各种优化数学算法与结构在动载荷作用下的分析计算方法相融合,对结构的布局参数、尺寸参数、材料参数等进行优化配置,达到改善结构动力特性、降低结构振动的目的。
根据优化目标的不同,结构动力优化可以分成两个部分:
即针对结构动力特性的优化设计和针对结构动力响应的优化设计。其中,针对结构动力特性的优化设计研究主要是以结构的固有频率、振型、阻尼和刚度与质量分布等方面为目标的优化设计研究。其中最早涉及到的课题是以结构固有频率为目标或约束的优化设计。设计阶段若通过合理配置结构的相应参数,使结构的固有频率远离激励频率,能够提高结构在正常工作时的安全性能。另外,针对结构动力响应的优化设计也是结构优化设计的重要组成部分,它主要是以结构在动载荷作用下振动响应的物理量,如:速度、加速度、振动功率流、声辐射能量等为目标或约束的结构优化设计。结构动力响应优化设计因为同时涉及到结构动力特性和动力响应分析计算以及相应数学优化算法,与结构动力特性优化设计相比较,其求解难度要更高,计算量更大。现阶段,结构的动力优化主要是指结构动力特性的优化设计。
02
结构动力优化设计的基本原理
2.1 结构模态分析方法
模态计算是用来计算线性结构的动力学特性:
-结构的固有频率;
-结构的模态振型;
-振型参与系数
-有效质量
模态分析所有动力学计算的最基础和最重要的分析;通过模态分析可以使结构避免共振或让结构在指定的频率下振动。
模态计算的假设和限制条件:
-结构是线性的,即具有恒定的总体质量矩阵和总体刚度矩阵
-结构没有外载荷(力,温度,压力等),即结构是自由振动。
2.1.1 无阻尼模态分析理论:
无阻尼线性结构自由振动的控制方程:
![3bc712d4c0b28796e4143cb0c99c0292.png](https://i-blog.csdnimg.cn/blog_migrate/789fe601a8520db4896fcf6f5add1de3.jpeg)
假设结构的运动简谐运动:
![740bb4ff32c0457147f7948dbb9add3b.png](https://i-blog.csdnimg.cn/blog_migrate/4a26853a32d1af0d8cb27105972eddb6.jpeg)
将结构运动的位移和速度,代入到控制方程中,可得
![36c585b9675631eab08b74a88c0799ad.png](https://i-blog.csdnimg.cn/blog_migrate/175af6c7b232bd7b04fd75603ea08b49.jpeg)
以下两种情况可以满足上述方程
![12b8247847d58df6048a4ecd816305ba.png](https://i-blog.csdnimg.cn/blog_migrate/c22bca029dc429fd3ad8f928f0c2ab36.jpeg)
-表明结构没有振动,这个情况不考虑舍去
![41ba34a9ea4893fe52cce7ead0e30ac0.png](https://i-blog.csdnimg.cn/blog_migrate/9581b085571b7dabaf087b98848de332.jpeg)
这个是一个特征值问题,可以求解出n个方程的根
![9eb5556e6d70e7e5f40028870d4babcc.png](https://i-blog.csdnimg.cn/blog_migrate/e38c1d699ae409abf6c1902dca476d71.jpeg)
这些根是这个方程的特征值;
对于每一根(特征值),都对应着一个特征向量
![ce16b002c24e44612dc3fe23dc06c103.png](https://i-blog.csdnimg.cn/blog_migrate/26273ac6ba95c7bebb4043fffcb2d888.jpeg)
ANSYS采用下式输出计算的固有频率:
![0945994772b089973394a375dd5bac9c.png](https://i-blog.csdnimg.cn/blog_migrate/2ae2f1a3e418086442ce6961c89cde3b.jpeg)
其中:fi的单位为Hz,即转/秒。如果模型的约束不足导致产生刚体运动,则总体刚度矩阵[K]为半正定型,则会出现固有频率为0的情况。
![b71c1b31fe6187a95e0b7cbf1f5def99.png](https://i-blog.csdnimg.cn/blog_migrate/e667b7d2b966e4fc2fdee67f20da1ef4.jpeg)
![998ab717d4674143c9f85cf67fb88550.png](https://i-blog.csdnimg.cn/blog_migrate/5f1d4bbaf00b05e87539889a4db69a7b.jpeg)
![c84b2ad8b12bde2d46f585e1d478b618.png](https://i-blog.csdnimg.cn/blog_migrate/ff127d1d58660d214ba55607132fca32.jpeg)
2.1.2 有阻尼模态分析理论:
有阻尼模态分析中假设结构没有外力作用,则控制方程变为
![84fba7064983d0e0b35b2d674793cf6d.png](https://i-blog.csdnimg.cn/blog_migrate/c0c27a97a4919d447456b81a62c18cfa.jpeg)
设解为
![71ea1c3f32f69f6253b140995d9355f4.png](https://i-blog.csdnimg.cn/blog_migrate/c13bfe50c25f3e3a6aa57ed1469ce4a4.jpeg)
代入方程(1)得
![c24d19caaf0f8c466ad6f1d0b6b8258e.png](https://i-blog.csdnimg.cn/blog_migrate/382ea17b9a67744f6539cff08bc3c169.jpeg)
矩阵
![57bcf90be49740fd4b8ae047afdf2986.png](https://i-blog.csdnimg.cn/blog_migrate/8b15b98bd8fc1b0760ca7e3d5f9ed728.jpeg)
称为系统的特征矩阵。方程(3)是一个“二次特征值”问
题,要(3)式有非零解的充要条件为
![1d11f707db029ad2b2f0ca8160b24318.png](https://i-blog.csdnimg.cn/blog_migrate/4aef0f1e747a5c059f496c24e30168fe.jpeg)
对于包含陀螺效应的旋转软化结构或需考虑阻尼的结构,则使用QR的表达式:Damped法求解模态振型和复特征值。特征值
![a3364395bac5f6f5439729b1e2a968a1.png](https://i-blog.csdnimg.cn/blog_migrate/3bc908f03a3fe1386d8b8d347ad07d59.jpeg)
的表达式
![eb16338dff620c5f16a9835f2b050e17.png](https://i-blog.csdnimg.cn/blog_migrate/3f2ca1545205f100d9488e15813d66f0.jpeg)
![45a0ced306e707fe401d399ae800b7f8.png](https://i-blog.csdnimg.cn/blog_migrate/eb6d2bd505b9820f93aef2557cbb2740.jpeg)
复数特征值的实部;
![85160bebcf5c5fca0638c5b89e25d4fb.png](https://i-blog.csdnimg.cn/blog_migrate/109c257d296c5a14c6f0a1a1617991ca.jpeg)
复数特征值的虚部
特征值的虚部
![846143c8d8cf754dcac49fa37e60ae07.png](https://i-blog.csdnimg.cn/blog_migrate/ff10bcfcc651b89fed5eaf45489061c5.jpeg)
代表系统的稳态角频率。特征值的实部
![b3c358671e09c5ef93987b9418f88bb1.png](https://i-blog.csdnimg.cn/blog_migrate/cdfeb4f4ee4626d2bb8d4c7743d5fff2.jpeg)
代表系统的稳定性 。
模态阻尼比由下式给出:
![062867ae248f46e8c72dfc98133f60b4.png](https://i-blog.csdnimg.cn/blog_migrate/ea8694fea612ca360a2d240120972c81.jpeg)
它的物理意义是实际阻尼与临界阻尼之比对数衰减率表示任意连续位移峰值的比值。它的表达式如下:
![69300e96702a204b913a73f922bcf18d.png](https://i-blog.csdnimg.cn/blog_migrate/39a914e35d7125ee52ca11c6e74d7140.jpeg)
![fa65155c4a6354e06148351265e9696d.png](https://i-blog.csdnimg.cn/blog_migrate/f101601da41d612477e19da8aacf4f4b.jpeg)
考虑阻尼的模态计算中的完整阻尼矩阵:
![ee1e2a7fe715314374f0ea3ec34d530e.png](https://i-blog.csdnimg.cn/blog_migrate/4895cdced841283eefa75d3eafa99466.jpeg)
对于阻尼矩阵,程序支持材料阻尼和总体阻尼控制,当需要考虑模型中不同的材料阻尼时,用户可以在工程数据模块为不同的材料定义刚度阻尼系数和质量阻尼系数
![f1bdba1a3050375460d54f7c9c2f90f1.png](https://i-blog.csdnimg.cn/blog_migrate/9c29ee80cdeb0e7f293382ffd5c41b3a.jpeg)
2.模态分析理论和术语
2.2 有阻尼模态分析理论:
用户可以在进入Model中定义结构的总体阻尼特性:
![44a2983ffb338a3c4e3eb35af726f327.png](https://i-blog.csdnimg.cn/blog_migrate/262bd582beeeadb78f7a1f0ae964a757.jpeg)
2.1.2 有阻尼模态分析理论:
考虑阻尼的模态计算输出的特征值是复数:
![a4cb7a78007f3d0c40f660c8cab41217.png](https://i-blog.csdnimg.cn/blog_migrate/24014e12d6f3490d6063f3f435b14043.jpeg)
![4ad60b2cbc18b35f70f2e029048ecd72.png](https://i-blog.csdnimg.cn/blog_migrate/0accb27e0e4c7ddf9fe60f2e3d4edf96.jpeg)
模态计算的方程数=无约束的节点*节点的自由度,即能够提取的最大模态阶数。
2.1.3固有频率和模态振型
在模态计算中:
1.特征值等于结构自然频率(固有频率)值的平方;
2.特征向量,对应于结构的振型,即结构的每个固有频率都关联一个模态形状向量
![c61d2c932eb198fcb2caa841605b3e72.png](https://i-blog.csdnimg.cn/blog_migrate/8281af66514da12774f2cfaa2bb0fa59.jpeg)
3.模态振型可以关于质量矩阵进行归一化:
![356cf62e-6122-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/356cf62e-6122-eb11-8da9-e4434bdf6706.png)
模态振型也可以关于单位矩阵进行归一化,基于该归一化方法可以比较各个节点发生共振时的相对位移量,在关于单位矩阵进行归一化时,将模态的最大特征向量设置为1:
Workbench后处理显示的关于质量矩阵归一化的结果;
由于模态计算结果是一个归一化的相对值,因此只有每个节点自由度的振动形状具有真实意义。
特征值的平凡根等于结构的固有频率(rad/s)
![9eda826c6974e9ec7fa226a588027dbf.png](https://i-blog.csdnimg.cn/blog_migrate/5341ffd31ddfb04547dd7e1b8d66d377.jpeg)
ANSYS Workbench输入和输出的固有频率的单位为Hz,因为输入和输出时候已经除以了2π。
![2978fb5ba4055ae4f0468fad899835a1.png](https://i-blog.csdnimg.cn/blog_migrate/9c96e0300d1df3f11cf6884918114caf.jpeg)
模态计算中的特征向量表征了结构的模态振型,如图所示该形状即为假设结构按照频率249Hz振动时的形状。
![643721f5be4375c01b2f650c87768f9d.png](https://i-blog.csdnimg.cn/blog_migrate/d90b5a444f0c8c6cea1b0693886f19cb.jpeg)
2.1.4.参与系数,有效质量
模态计算后除了能够获取结构的固有频率和振型外,还有参与系数与有效质量,其中参与系数的计算公式:
![02b7d3a4be91c53a2d9e07c529e55675.png](https://i-blog.csdnimg.cn/blog_migrate/1c2d0723c44d00943bd7a74c3e6ad821.jpeg)
其中,{D}为总体笛卡尔坐标系中三个平动和三个转动的单位位移谱。
(1)参与系数是表征每阶模态每个方向质量运动的总和,即在特定方向上某个振型在多大程度上参与了振动)
(2)某一个方向的参与系数值越大,则表明如果在这个方向受到动载荷的作用下,则结构的振动响应越大。
(3)比值(Ratio),就是某个方向所有模态阶数的参与系数关于最大的参与系数进行归一化的比值。
参与系数的计算公式:
![2cdf772b1ebf5afd19660b63376fc1e0.png](https://i-blog.csdnimg.cn/blog_migrate/47eed7b629d477d6407c0b5c8623dbb0.jpeg)
![7b563a66d45a42e2427c0438d176c3dd.png](https://i-blog.csdnimg.cn/blog_migrate/85f5d83f6faa2bc0dc7971b0985fb7a0.jpeg)
![7b14ad53be4da25d437dbc6996a00315.png](https://i-blog.csdnimg.cn/blog_migrate/7cfa8d4db19d794832cb302ef76d3713.jpeg)
模态计算中的有效质量计算公式:
![7f854c754bfef417f7bbee6cb57e8207.png](https://i-blog.csdnimg.cn/blog_migrate/5c3ce3d097947592f35d7f4ec0b36642.jpeg)
由于程序模态计算时,各个振型关于质量矩阵进行归一化,即
![5765eb38c57eca03105c9f02cb0d7a0c.png](https://i-blog.csdnimg.cn/blog_migrate/85a606d7d94da3bd3ae07e7b5b2af590.jpeg)
则,有效质量可简化为
![89e90a88146b055b444a92102c21005e.png](https://i-blog.csdnimg.cn/blog_migrate/41c3f94e2ff6cd9b988457d1cc222448.jpeg)
-理想情况下,在每个方向的所有有效质量之和等于结构的总质量,但是这个取决于模态计算提出的模态阶数;
-有效质量与结构总质量的比值对于确定提取的模态数量是否足够,非常有帮助。
模态计算中的有效质量计算公式:
![696c111192bcd7fa5019c9697635ca7a.png](https://i-blog.csdnimg.cn/blog_migrate/8067a53c4c8c1985ec0d21b0f738baa8.jpeg)
![4744756223e9cbceb0cb20f0a1daae4a.png](https://i-blog.csdnimg.cn/blog_migrate/d8f3611623a5993d638de0ae6db1be21.jpeg)
![f1ebcba8974be5934aefb58bffd84bb3.png](https://i-blog.csdnimg.cn/blog_migrate/02da7a46066968dd84a2843828030ab6.jpeg)
2.1.5 模态的提取方法
在大多数情况下,建议用户选用 Program Controlled选项,程序会自 动优化进行选择算法。
![83e05e32e0983dcb2b3831441588b142.png](https://i-blog.csdnimg.cn/blog_migrate/dc350ca14a7957e125a2282afb0a24a2.jpeg)
(1)Direct-Block Lanczos
-能够处理对称矩阵;
-是一种功能强大的方法,当提取中型到大型模型(50000 ~ 100000 个自由度)的大量振型时(40+),这种方法很有效;
-经常应用在具有实体单元或壳单元的模型中;
-可以很好地处理刚体振型;
-需要较高的内存。
(2)Iterative-PCG Lanczos
-能够处理对称矩阵,但是不用于求解屈曲模态;
-适合求解中等到大规模的模态计算问题,提取的模态阶数高于100阶;
-适合于网格划分形状较好的三维实体单元;
(3)Unsymmetric
-能够处理非对称矩阵;
-模态计算中使用完整的刚度和质量矩阵;
-适合求解K和M为非对称矩阵的问题,如流-固耦合的振动,声学振动;
-计算以复数表示的特征值和特征向量:
--实数部分就是自然频率;
--虚数部分表示稳定性,负值表示稳定,正值表示不确定。
![5f27cfe1b38ae3c6e888da27e85d0d0e.png](https://i-blog.csdnimg.cn/blog_migrate/94955f11e4c160d1d534bcf1ed825fcb.jpeg)
(4)Supernode
-能够处理对称矩阵,但是不用于求解屈曲模态;
-适合求解大规模的模态计算问题,提取的模态阶数高于100000阶;
-主要应用于二维平面,壳体/梁结构(提取模态阶数高于100)和三维实体
结构(提取模态阶数高于250);
如果结构中存在阻尼,则将阻尼选项 设置为yes,然后选择
相应的方法进 行求解。
(5)Full Damped
![13e4ff73289b3b8e27bdee0e8545d7fd.png](https://i-blog.csdnimg.cn/blog_migrate/792af9a9f29b3aaf10a478627dd3d72c.jpeg)
(6) Reduced Damped
QR阻尼法能够很好地求解大阻尼系统模态解,阻尼可以是任意阻尼类型,即 无论是比例阻尼或非比例阻尼。由于该方法的计算精度取决于提取的模态数目,所以建议提取足够多的基频模态,特别是阻尼较大的系统更应当如此,这样才能保证得到好的计算结果。该方法不建议用于提取临界阻尼或过阻尼 系统的模态。该方法输出实部和虚部特征值(频率),但仅仅输出实特征向量(模态振型)。
![0c7d9a774d63b3816c1825703717d223.png](https://i-blog.csdnimg.cn/blog_migrate/4de2183044dd84c72808dd7e767512db.jpeg)
2.2 结构动力优化模型
以结构重量为目标函数,以某阶频率或频率禁区
![ecb20c3a376d0c7ee205f68b583bd831.png](https://i-blog.csdnimg.cn/blog_migrate/d9c44cc91a1e5ee292ec73a6c052947f.jpeg)
为约束,要求重量最小。
![fd3665686a58f43b5b1a4143f5c4f156.png](https://i-blog.csdnimg.cn/blog_migrate/5899cea8f0b20fadbbebc6f0acc278f1.jpeg)
![e9083d6add3e16d18070f6530f41e22e.png](https://i-blog.csdnimg.cn/blog_migrate/5b2d125992193d853878b790070e8114.jpeg)
数。其中设计变量的选择根据各种优化问题的具体情况,如杆件的截面面积、板的厚度设计变量梁的截面惯性矩等等,如果考虑形状的优化,节点坐标也可以作为在动力问题中,非结构质量的位置和大小也可以作为设计变量加以考虑。实际工程优化问题还应同时考虑应力、位移等其他的约束条件。
2.3 频率优化中的灵敏度分析
灵敏度分析是一种度量。它是一种评价因设计变量或参数的改变而引起结构响应特性变化率的方法。结构系统灵敏度的研究是一个很特别的领域,它是当前计算力学和结构工程领域的主要研究方向之一。灵敏度分析是结构优化设计方法的基础。在对离散化的结构模型进行优化时,除非应用少数几种比较初等的方法,否则,不论是用准则设计法还是数学规划法,在修改结构的尺寸使之向更优方向变化时,都必须作灵敏度分析。在结构动力优化设计中,寻求有效的结构自振频率及振型关于设计变量灵敏度的方法十分重要。关于频率优化中特征值和特征向量的灵敏度的求解,特征根问题可以描述如下:
![cfdc812e4fd24b38870a70f845e1c3ad.png](https://i-blog.csdnimg.cn/blog_migrate/151bec82463fb3edbfb01ff81129b277.jpeg)
关于频率优化中特征值和特征向量的灵敏度的求解,特征根问题可以描述如下:
![4f9ba7387fa3711cdd398a5762cc71d0.png](https://i-blog.csdnimg.cn/blog_migrate/ffc3fafc28da4a5e2a62c47d486d7cbb.jpeg)
其中,K为结构的刚度矩阵,M为结构的质量矩阵,
![48acac9e39e31034b48f2f5424c4c0b8.png](https://i-blog.csdnimg.cn/blog_migrate/a129d488906db1eaa73d4fb05d2f039a.jpeg)
为特征根,
![9f1df50cd6aee96cc1774a924259366e.png](https://i-blog.csdnimg.cn/blog_migrate/c865eae31b27dca80f39d82a2557970e.jpeg)
为特征向量,K’、M’分别为刚度矩阵和质量矩阵关于设计变量的导数。
03
结构动力优化设计的ANSYS实现
![59cc8c35d6289f5f7969204efd07ab0e.png](https://i-blog.csdnimg.cn/blog_migrate/4f65bbb9ee315f802efb0808e1e40274.jpeg)