CN107292047A - 一种基于最大位移的复杂板壳厚度快速优化设计方法 - Google Patents
一种基于最大位移的复杂板壳厚度快速优化设计方法 Download PDFInfo
- Publication number
- CN107292047A CN107292047A CN201710536304.3A CN201710536304A CN107292047A CN 107292047 A CN107292047 A CN 107292047A CN 201710536304 A CN201710536304 A CN 201710536304A CN 107292047 A CN107292047 A CN 107292047A
- Authority
- CN
- China
- Prior art keywords
- mrow
- thickness
- msub
- shell
- plate
- 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.)
- Withdrawn
Links
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明涉及一种基于最大位移的复杂板壳厚度快速优化设计方法，属于机械结构快速优化设计领域。本发明根据有限元理论和板壳理论中应力、载荷、位移以及刚度矩阵之间的理论关系，推导出基于最大位移的最优板壳厚度计算公式。针对复杂板壳结构，由于不能忽略复杂结构特征对结构性能的影响，复杂板壳结构的厚度难以直接用板壳理论进行计算，该方法能快速计算出复杂板壳结构在满足相应临界最大位移要求下的最优解，具有快捷、计算量小、求解容易等优点。
Description
技术领域
本发明属于机械结构快速优化设计领域，涉及一种基于最大位移的复杂板壳厚度的快速优化设计方法。
背景技术
板壳是典型的工程构件，其厚度尺寸远小于其它特征尺寸。随着复杂板壳结构广泛应用于航空航天、仪表元件和海洋工程中，其安全性越来越引起人们的注意，而板壳厚度的设计直接关系着其机械性能。因此对设计人员来说，如何快速准确地确定板厚显得尤为重要。
由于不能忽略复杂结构特征对结构性能的影响，复杂板壳结构的厚度难以直接用板壳理论进行计算。目前，国内外对复杂板壳厚度设计主要采用经验法和类比法，这两种设计方法对经验的依赖度很高且设计多为冗余设计。
发明内容
本发明的目的是提供一种基于最大位移的复杂板壳厚度快速优化设计方法，其可快速计算出复杂板壳结构在满足相应临界位移条件下的最优厚度，以解决板壳结构复杂而难以直接利用板壳理论进行计算的问题。
本发明按以下技术方案实现：根据有限元理论和板壳理论中应力、载荷、位移以及刚度矩阵之间的理论关系，推导出基于最大位移的最优板壳厚度计算公式：首先，假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构，其中厚度h1某一任意已知的初始值，厚度hcr未知，hcr为满足最大位移μcr时的最优板壳厚度，即hcr为所求最优解；然后建立有限元模型，通过有限元分析，计算出厚度为h1时的最大位移μ1；进一步地，将μ1、h1、μcr带入公式便可以快速计算出满足临界最大等效应力条件下的最优厚度。
获取最优厚度的计算公式的具体方法如下：
首先，假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构。根据有限元理论，板壳位移和载荷关系可表示为：
Kμ＝F (1)
式中，K为总体刚度矩阵，μ为位移，F为载荷。
根据有限元理论，对完全约束的结构，K为正定矩阵，因此式(1)可表示为：
μ＝K-1F (2)
根据有限元理论和板壳理论，刚度矩阵K可表示为：
K＝h3Ku (3)
式中，Ku表示当壳厚度为单位1的刚度矩阵。
将式(3)代入式(2)，可得：
μ＝h-3(Ku)-1F (4)
式(4)中，对满足设计要求的位移μ为已知，Ku和F已知的情况下，可计算板壳厚度h。但由于Ku的获取比较困难，所以对于设计人员来说很难直接通过上式来计算。
假设某一任意已知的板壳厚度h＝h1时，位移向量μ1和厚度h1的关系为：
当板壳厚度h＝hcr，位移向量μcr和厚度hcr的关系式可表示为：
根据假设，对于相同材料模型，相同网格划分和相同单元类型。即，式(5)和式(6)中将两式相除可得：
由式(7)可知，当厚度分别为h1和hcr时对应的节点位移向量μ1和μcr之间的关系式。设计时，临界位移μcr往往为一维的标量(如：最大位移或某点的位移等)。
根据线性理论和有限单元理论，具有相同物理意义的位移标量μ1和设计要求的临界位移μcr，同样可满足式(7)中的关系式。因此，优化设计厚度hcr可表示为：
通过有限元分析，计算出厚度为初始值h1时的最大位移μ1，根据相应的已知约束条件μcr，将μ1、h1、μcr带入公式便可以快速计算出满足临界最大位移μcr时的最优板壳厚度。
本发明的有益效果是：针对复杂板壳结构，由于不能忽略复杂结构特征对结构性能的影响，复杂板壳的厚度难以直接用板壳理论进行计算，该方法能快速计算出满足相应临界最大位移条件下的最优解，具有计算量小、求解容易等优点。
附图说明
附图1为本发明的方法步骤图。
具体实施方式
下面结合附图和实施例，对本发明作进一步说明，但本发明的内容并不限于所述范围。
本发明的实施例：附图1为本发明的方法步骤图，根据有限元理论和板壳理论中应力、载荷、位移以及刚度矩阵之间的理论关系，推导出基于最大位移的最优板壳厚度计算公式，本发明方法的具体实现步骤为：步骤一，假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构，其中厚度h1某一任意已知的初始值，厚度hcr未知，hcr为满足最大位移μcr时的最优板壳厚度，即hcr为所求最优解；步骤二，通过有限元分析，计算出厚度为h1时的最大位移μ1；步骤三，将μ1、h1、μcr带入公式便可以快速计算出满足临界最大等效应力条件下的最优厚度。
获取最优厚度的计算公式的具体方法如下：
首先，假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构。根据有限元理论，板壳位移和载荷关系可表示为：
Kμ＝F (1)
式中，K为总体刚度矩阵，μ为位移，F为载荷。
根据有限元理论，对完全约束的结构，K为正定矩阵，因此式(1)可表示为：
μ＝K-1F (2)
根据有限元理论和板壳理论，刚度矩阵K可表示为：
K＝h3Ku (3)
式中，Ku表示当壳厚度为单位1的刚度矩阵。
将式(3)代入式(2)，可得：
μ＝h-3(Ku)-1F (4)
式(4)中，对满足设计要求的位移μ为已知，Ku和F已知的情况下，可计算板壳厚度h。但由于Ku的获取比较困难，所以对于设计人员来说很难直接通过上式来计算。
假设某一任意已知的板壳厚度h＝h1时，位移向量μ1和厚度h1的关系为：
当板壳厚度h＝hcr，位移向量μcr和厚度hcr的关系式可表示为：
根据假设，对于相同材料模型，相同网格划分和相同单元类型。即，式(5)和式(6)中将两式相除可得：
由式(7)可知，当厚度分别为h1和hcr时对应的节点位移向量μ1和μcr之间的关系式。设计时，临界位移μcr往往为一维的标量(如：最大位移或某点的位移等)。
根据线性理论和有限单元理论，具有相同物理意义的位移标量μ1和设计要求的临界位移μcr，同样可满足式(7)中的关系式。因此，优化设计厚度hcr可表示为：
通过有限元分析，计算出厚度为初始值h1时的最大位移μ1，根据相应的已知约束条件μcr，将μ1、h1、μcr带入公式便可以快速计算出满足临界最大位移μcr时的最优板壳厚度。
Claims (1)
1.一种基于最大位移的复杂板壳厚度快速优化设计方法，其特征在于：根据有限元理论和板壳理论中应力、载荷、位移以及刚度矩阵之间的理论关系，推导出基于最大位移的最优板壳厚度计算公式；假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构，其中厚度h1为某一任意已知的初始值，厚度hcr未知，hcr为满足最大位移μcr时的最优板壳厚度，即hcr为所求最优解；建立有限元模型，通过有限元分析，求解厚度为h1时的最大位移μ1；将μ1、h1、μcr带入公式便可以快速计算出满足临界最大位移条件下的最优厚度；
获取最优厚度的计算公式的具体方法如下：
假设厚度为h1和hcr的板壳结构为网络类型及网络划分、边界条件、载荷完全相同，材料各向异性，仅厚度不一样的板壳结构；根据有限元理论，板壳位移和载荷关系可表示为：
Kμ＝F (1)
式中，K为总体刚度矩阵，μ为位移，F为载荷；
根据有限元理论，对完全约束的结构，K为正定矩阵，因此式(1)可表示为：
μ＝K-1F (2)
根据有限元理论和板壳理论，刚度矩阵K可表示为：
K＝h3Ku (3)
式中，Ku表示当壳厚度为单位1的刚度矩阵；
将式(3)代入式(2)，可得：
μ＝h-3(Ku)-1F (4)
式(4)中，对满足设计要求的位移μ为已知，Ku和F已知的情况下，可计算板壳厚度h；
假设某一任意已知的板壳厚度h＝h1时，位移向量μ1和厚度h1的关系为：
<mrow>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>h</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>K</mi>
<mn>1</mn>
<mi>u</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>F</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
当板壳厚度h＝hcr，位移向量μcr和厚度hcr的关系式可表示为：
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
<mi>u</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>F</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
对于相同材料模型，相同网格划分和相同单元类型；即，式(5)和式(6)中将两式相除可得：
<mrow>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
由式(7)可知，当厚度分别为h1和hcr时对应的节点位移向量μ1和μcr之间的关系式；临界位移μcr为一维的标量；
根据线性理论和有限单元理论，具有相同物理意义的位移标量μ1和设计要求的临界位移μcr，满足式(7)中的关系式；优化设计厚度hcr可表示为：
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
1
通过有限元分析，计算出厚度为初始值h1时的最大位移μ1，根据约束条件μcr，将μ1、h1、μcr带入公式快速计算出满足临界最大位移μcr时的最优板壳厚度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710536304.3A CN107292047A (zh) | 2017-06-21 | 2017-06-21 | 一种基于最大位移的复杂板壳厚度快速优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710536304.3A CN107292047A (zh) | 2017-06-21 | 2017-06-21 | 一种基于最大位移的复杂板壳厚度快速优化设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107292047A true CN107292047A (zh) | 2017-10-24 |
Family
ID=60098622
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710536304.3A Withdrawn CN107292047A (zh) | 2017-06-21 | 2017-06-21 | 一种基于最大位移的复杂板壳厚度快速优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107292047A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108446419A (zh) * | 2018-01-17 | 2018-08-24 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
| CN108614922A (zh) * | 2018-04-04 | 2018-10-02 | 南京航空航天大学 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
-
2017
- 2017-06-21 CN CN201710536304.3A patent/CN107292047A/zh not_active Withdrawn
Non-Patent Citations (1)
| Title |
|---|
| 吴盛彬: "挖掘机动臂力学性能分析及板厚优化方法研究", 《中国优秀硕士学位论文全文数据库-工程科技II辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108446419A (zh) * | 2018-01-17 | 2018-08-24 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
| CN108446419B (zh) * | 2018-01-17 | 2022-02-22 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
| CN108614922A (zh) * | 2018-04-04 | 2018-10-02 | 南京航空航天大学 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
| CN108614922B (zh) * | 2018-04-04 | 2021-08-06 | 南京航空航天大学 | 一种基于应力杂交化后处理的板壳边缘效应的计算方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103279627B (zh) | 一种基于有限元的热-机械-磨损耦合分析数值模拟方法 | |
| Alkebsi et al. | Design of graded lattice structures in turbine blades using topology optimization | |
| Kumar et al. | Finite element analysis of laminated composite and sandwich shells using higher order zigzag theory | |
| He et al. | Partitioned subiterative coupling schemes for aeroelasticity using combined interface boundary condition method | |
| Xing et al. | An eigenelement method of periodical composite structures | |
| Liu et al. | Error and its upper bound estimation between the solutions of GM (1, 1) grey forecasting models | |
| CN107292047A (zh) | 一种基于最大位移的复杂板壳厚度快速优化设计方法 | |
| Dai et al. | A multiple scale time domain collocation method for solving non-linear dynamical system | |
| Hamed Meraji et al. | An efficient algorithm based on the differential quadrature method for solving Navier–Stokes equations | |
| CN202372981U (zh) | 一种适用于空间飞行器的动力学模型修正装置 | |
| CN107357986A (zh) | 一种基于有限元分析的复杂板壳厚度快速优化设计方法 | |
| Xiangyang et al. | NS-FEM/ES-FEM for contact problems in metal forming analysis | |
| CN107346358A (zh) | 一种基于最大等效应力的复杂板壳厚度快速优化设计方法 | |
| Kanber et al. | Investigation of RPIM shape parameter effects on the solution accuracy of 2D elastoplastic problems | |
| Fish et al. | Constitutive modeling based on atomistics | |
| CN107341307A (zh) | 一种基于有限元分析的结构材料快速选择方法 | |
| Wang et al. | Accurate bending strength analysis of the asymmetric gear using the novel es-pim with triangular mesh | |
| Zhang et al. | An implicit algorithm for high-order DG/FV schemes for compressible flows on 2D arbitrary grids | |
| Zhang et al. | A high-order flux reconstruction/correction procedure via reconstruction method for shock capturing with space-time extension time stepping and adaptive mesh refinement | |
| Vanli et al. | A Bayesian approach for integration of physical and computer experiments for quality improvement in nano‐composite manufacturing | |
| Duan et al. | A modified liu-storey conjugate gradient method and its global convergence for unconstrained optimization | |
| CN107315882A (zh) | 一种基于静态分析的结构材料快速选择方法 | |
| Zhu et al. | A QUASI-NEWTON TRUST REGION METHOD BASED ON A NEW FRACTIONAL MODEL. | |
| CN104268320A (zh) | 适用于卫星敏感载荷的新型隔振器刚度矩阵快速预估方法 | |
| Miles et al. | Numerical techniques to model fractional-order nonlinear viscoelasticity in soft elastomers |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20171024 |