一、Pro/MECHANICA环境下的优化设计
% z8 B. `% }9 j8 s0 X 1.Pro/MECHANICA Structure模块简介& J4 C7 K& ^) w2 T* r
与其他专业的有限元分析软件一样,Pro/MECHANICA Structure模块可以完成有限元分析和模型的优化设计,它的设计研究种类主要有三种:标准分析(Standard)、灵敏度分析(Sensitivity)和优化设计分析(Optimization)。图1为Pro/MECHANICA Structure模块的主菜单。概括的说,Pro/MECHANICA Structure模块的分析任务为两类,第一类为设计验证或设计校核,例如进行设计模型的应力应变检验,和其他有限元分析软件一样,须通过创建几何模型、简化模型、设定单位和材料属性、定义约束、定义载荷、定义分析任务、运行分析、显示评价计算结果这样的工作流程;第二类为模型的设计优化,这是Pro/MECHANICA区别其他有限元软件最显著的特征。它的工作流程如图2所示。现以管道支架为例简述Pro/MECHANICA独特的优化设计功能。
& m: k) Y2 P) ~
^7 q! k' A; {+ }3 K. m) u
" W' ` ^, ^- ]图1 Pro/MECHANICAStructure模块主菜单
2 X; a6 m7 W" u# ?+ R0 `5 U9 j5 R5 j" v* ]
图2 Pro/MECHANICA优化设计工作流程图 2.管道支架的优化设计
1 v1 W) k7 R3 m 管道支架是一种广泛应用于建筑、电力等行业的支撑类零件,对其进行有限元分析和优化设计对实际生产有着重要的指导意义。
+ [- ?' w& J4 b' [. `6 N: z% W (1)几何模型的建立
6 S K' J; B! R. C* e3 E0 { 图3为某管道支架在Pro/E下的三维造型图,主要由中间板、边缘板、中间切槽、左端固定板组成,其左端面固定。由材料力学可知,该支架左端固定板与中间板、边缘板连接面为危险截面,截面上下位置处的点所受应力最大,为危险点。因此分别建立测量点PNT0、PNT1以监测危险点处的应力值。
0 d0 Q4 n; f) f+ o 5 l! A" P( ?5 m9 L: i
) p1 Y3 c# t8 @6 J; T2 S
图3 管道支架三维效果图 (2)设置模型材料属性
1 _3 [0 p' o. M6 ~: I9 a* R* r$ d% C. x 该模型所选材料为45钢,其密度为 ,弹性模量E=200GPa,泊松比 。根据材料属性,我们选择Pro/MECHANICA Structure模块中材料库中的材料模型steel,然后整理编辑,图4为定义的材料属性表。
1 V- O& ^; O3 [5 e4 r9 s O% }7 J- k3 O5 N% r
2 k$ x7 ^5 F, X2 t* I* ?! E图4 管道支架材料属性定义框 (3)约束及载荷的定义; @* u+ z; ]/ M, E' p0 q( h: a9 _
由于管道支架的左端面固定,因此定义约束时将该端面6个自由度完全固定;支架上端面承受管道重量1000N,因此选择载荷分布形式为均布(Uniform)的总载荷(Total Load),Y轴的分量-1000N,施加在50×100的面域上,如图5所示。
! J" C, q1 m( i* z5 O 2 R4 ^7 f$ l! B) O o9 A. S
. A5 `# Q- w0 @/ t
图5 管道支架的约束与载荷 (4)灵敏度分析 在建造模型过程中,设计人员定义了很多设计参数,如目标设计参数、物性参数等,当这些设计参数变化时,必然会对模型性能产生影响,如果这些设计参数都用于优化设计,必然加大计算量。考虑到不同设计参数对模型性能影响的程度不一样,如果能定量地表示这种影响程度,就能确定哪些是重要设计参数,并排除那些影响程度很小的设计参数,从而提高优化设计的效率,这是优化设计时要着重考虑的。而Pro/M的灵敏度分析恰恰提供了这样的定量分析功能,为优化设计高效率提供了保证。
- d3 G! V; m" _. j) q! r" B7 h 1)局部灵敏度分析( U: F( \- y9 Y1 r0 y, ]& A& @
局部灵敏度分析就是定量分析每个设计参数对模型性能影响的程度,从而可以确定哪些参数对模型性能影响最大,即重要设计参数,设计优化时作重点考虑。/ l3 {1 H1 A Y' ]8 D, s
在上述的管道支架中,设计参数主要有中间板的长度(Mid_Length)、宽度(Mid_Width)、板厚(Mid_Thickness)、边缘板的宽度(Edge_Width)、厚度(Edge_Thickness)、中间切槽的长度(Cut_ Length)和宽度(Cut_ Width)等。因此分别对这些参数进行局部灵敏度分析,最后确定对模型最大应力影响最深的重要参数,运行局部灵敏度分析后,图6~12分别为中间板的长度、宽度、厚度、边缘板的宽度、厚度、中间切槽的长度和宽度等7个参数对模型应力影响程度的曲线。从曲线图可以得知,管道支架模型应力对设计参数中间板宽度(Mid_Width)和板厚(Mid_Thickness)的变化最敏感,即中间板宽度和厚度的变化对模型应力影响最大,因此这是优化设计时重点考虑的参数;其他设计参数尺寸的变化对模型最大应力影响不大。, O8 ?( Y" h1 }" s$ p' u% i
2)全局灵敏度分析4 H$ E O( [2 `* f/ E6 v
全局灵敏度分析就是在局部灵敏度分析出对模型性能影响较大的参数的基础上,为这些重要设计参数确定用于优化设计的变化范围,也就是说全局灵敏度分析能准确地描述模型性能对用于优化设计的重要设计参数的灵敏度,确定其合理的参数变化范围,从而在这些参数的变化范围中寻求最佳设计。
% g' ]! o' t" c4 V9 Y0 u) T4 \ 由于在管道支架的局部灵敏度分析中的得出了中间板宽度和厚度对模型应力影响最大的结果,因此分别建立全局灵敏度分析任务(定义参数Mid_Width、Mid_Thickness的变化范围分别为160~240mm和15~25mm),分析运算后准确确定了管道支架模型应力对中间板宽度和厚度尺寸变化的灵敏度如图13、14所示。! r. h ]/ k8 {9 C( N" y
- V; x6 X/ M+ I. ~
& D( Z; E- c: K* @9 j3 S图6 中间板长度对模型应力的影响曲线 ! i" a. H" \, a: c' x! N+ z
, K+ E6 J5 u1 r' T图7中间板宽度对模型应力的影响曲线
0 k0 E. x% J$ Z$ ^% ?1 ?1 d r* I
7 b* e! E4 E: e2 [, c2 [6 @图8 中间板厚度对模型应力的影响曲线 ) }( p1 B+ R: ~( s
5 @( t& R1 K$ t' H7 S/ S6 J图9边缘板宽度对模型应力的影响曲线
3 I" T& z3 W6 }% M: ]' W& }0 g- D9 m5 x2 W
图10 边缘板厚度对模型应力的影响曲线 & w: J. N/ _* d1 ^, r
3 H/ O! e/ c e: U图11 中间切槽长度对模型应力的影响曲线 & | S7 l/ F/ {6 |
6 E5 H5 B2 @1 }( {2 C! E3 L9 {
2 ~7 P+ k& c! Y2 M1 N图12 中间切槽宽度对模型应力的影响曲线
/ ~5 T( U5 |' K) s% T6 P, h# u2 O7 F$ C3 k
图13中间板宽度全局灵敏度分析曲线 ! ?2 H9 K6 ?$ w, `( Y+ @
% y, J: X2 m! X; X( F图14 中间板厚参数全局灵敏度分析曲线 (5)优化设计
* k/ H0 D2 v# L# Y6 U, ~9 T! [" { 优化设计研究的内容主要就是保证模型的约束条件前提下,尽可能达到质量轻、体积小、形状合理、成本最低以及力学方面最大限度地减缓过渡区应力集中等目标条件。其数学模型一般可统一写成如下格式:
% c! k& \1 ]" H: _6 |% G# ^' @
6 S% v2 v5 }4 L1 d 式中D为设计变量X在欧式空间Rn的取值范围,称为可行域;“s.t.”为约束条件。即优化设计的数学模型包括三因素:设计变量X、目标函数f(X)和约束条件 。; c" p4 {/ R% r0 E
在本文的管道支架实例中,分别确定三要素如下:
/ A5 p4 w8 ]3 s/ R* [ w; C 目标函数:使支架质量最小;2 k% d. d7 e( w1 T1 h8 K" M
设计参数:中间板宽度(Mid_Width)、厚度(Mid_Thickness)。由灵敏度分析可知管道支架模型的最大应力对中间板的宽度和厚度变化最敏感,对其他参数的变化不敏感,因此中间板的宽度和厚度尺寸是着重研究的参数。原始尺寸Mid_Width=200mm,Mid_Thickness=20mm。
) N2 Y7 p8 _( Y/ N* j, q6 X$ Y 约束条件:强度条件 =40MPa,前面建立的测量点PNT0、PNT1可反应危险截面处的应力情况,在满足材料强度的前提下,必须控制危险截面处的最大应力值,减缓应力现象,即 , 。3 o# M. R; \# T& C
根据上述数学模型建立优化设计任务并运行得到了一组使管道支架重量最轻的优化结果。现将原始值与最优值作一个比较,如表1所示。
# T8 E2 f2 q$ }. z表1 最优值与原始值的比较 & P" W( _) f4 L/ l* ^6 [5 T1 Z
从上表中可以看出:经优化之后;支架的重量减少了0.00526(Tonne),减重率达12.6%;最大应力和位移比原始( v" O. P+ F0 L" {- V
值略有增大,但前提是满足所有的约束条件,对模型影响不大,图15为优化过程的历史结果曲线,图16为最后确定的优化管道支架模型。% | t( ]7 l, o6 i
6 ^" G3 m$ N# y) y1 t+ V/ W S) r- O. s( U. x
图15 优化历史结果曲线
2 ], N7 b" y! R4 O0 s) J
( Y8 `) |8 Z' Y1 T% ?9 p* ]8 V图16 优化后的管道支架几何模型几尺寸 |