发明名称 一种基于复合人因势场的口盖尺寸优化设计方法
摘要 本发明公开了一种人因势场的口盖尺寸优化设计方法,包括以下步骤:一、建立复合人因势场模型;二、进行人因分析验证维修人员初始位置;三、计算满足应力强度的最小倒圆;四、基于人因势场理论优化求解满足维修性要求的口盖参数。本发明可以计算得到最小的口盖面积,在避免应力集中的同时,保证良好的可视性和可达性水平。该方法物理意义明确,且可以通过计算机自动完成,避免了通过人工经验设计所导致的口盖大小不合理问题,也避免了需反复修改所导致的效率低下问题,为飞机、船舶等口盖的维修性设计提供了科学的技术途径。
申请公布号 CN103761381B 申请公布日期 2015.02.18
申请号 CN201410020674.8 申请日期 2014.01.16
申请人 中国人民解放军国防科学技术大学 发明人 葛哲学;杨拥民;陈健;吴福章;官凤娇;董孟高;罗旭;胡政;陆智淼
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 李振
主权项 一种基于复合人因势场的口盖尺寸优化设计方法,其特征在于包括以下步骤:(一)、建立复合人因势场模型:复合人因势场为可视性和可达性两个因素所构成的人工势场;定义复合人因势能函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>U</mi><mi>F</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mo>[</mo><mn>1</mn><mo>-</mo><msub><mi>C</mi><mi>&theta;</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>P</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>C</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>P</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>DOF</mi></munderover><msub><mi>&omega;</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>q</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup></mrow><mrow><msubsup><mi>q</mi><mi>i</mi><mi>U</mi></msubsup><mo>-</mo><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mtd><mtd><msubsup><mi>q</mi><mi>i</mi><mi>U</mi></msubsup><mo>&GreaterEqual;</mo><msub><mi>q</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup></mtd></mtr><mtr><mtd><mo>[</mo><mn>1</mn><mo>-</mo><msub><mi>C</mi><mi>&theta;</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>P</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>C</mi><mi>l</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>P</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>DOF</mi></munderover><msub><mi>&omega;</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>q</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup></mrow><mrow><msubsup><mi>q</mi><mi>i</mi><mi>L</mi></msubsup><mo>-</mo><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mtd><mtd><msubsup><mi>q</mi><mi>i</mi><mi>N</mi></msubsup><mo>&GreaterEqual;</mo><msub><mi>q</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msubsup><mi>q</mi><mi>i</mi><mi>L</mi></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000011.GIF" wi="1556" he="333" /></maths>其中:q<sub>i</sub>为从人体腰部到眼部和手部的各关节角度,<img file="FDA0000600596620000015.GIF" wi="64" he="79" />为中心角度,q<sub>imax</sub>为从人体腰部到眼部和手部的各关节的最大活动范围,q<sub>imin</sub>为从人体腰部到眼部和手部的各关节的最小活动范围;ω<sub>i</sub>为各关节舒适度的权值,满足<img file="FDA0000600596620000012.GIF" wi="228" he="143" />DOF为各关节自由度之和;C<sub>l</sub>(X<sub>P</sub>)为可视度视距函数;C<sub>θ</sub>(X<sub>P</sub>)为视野影响函数,由垂直视野函数和水平视野函数构成;定义垂直视野函数<img file="FDA0000600596620000016.GIF" wi="71" he="79" />为:<img file="FDA0000600596620000013.GIF" wi="1525" he="665" />式(2)中,θ<sub>y</sub>为垂直视线角:<img file="FDA0000600596620000014.GIF" wi="461" he="149" />θ<sub>y1min</sub>,θ<sub>y1max</sub>为最佳视野的最小垂直视线角和最大垂直视线角,θ<sub>y0min</sub>,θ<sub>y0max</sub>为最大视野的最小垂直视线角和最大垂直视线角;定义<img file="FDA0000600596620000024.GIF" wi="67" he="65" />水平视野函数表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>C</mi><msub><mi>&theta;</mi><mi>x</mi></msub></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>&theta;</mi><mi>x</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>|</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>|</mo><mo>></mo><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>0</mn><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>-</mo><mn>0.1</mn><mo>&CenterDot;</mo><mfrac><mrow><mo>|</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>|</mo></mrow><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>1</mn><mi>max</mi></mrow></msub></mfrac><mo>+</mo><mn>1</mn></mtd><mtd><mn>0</mn><mo>&lt;</mo><mo>|</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>|</mo><mo>&le;</mo><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>1</mn><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0.9</mn><mo>&CenterDot;</mo><mfrac><mrow><mo>|</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>|</mo><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>0</mn><mi>max</mi></mrow></msub></mrow><mrow><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>1</mn><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>0</mn><mi>max</mi></mrow></msub></mrow></mfrac></mtd><mtd><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>1</mn><mi>max</mi></mrow></msub><mo>&lt;</mo><mo>|</mo><msub><mi>&theta;</mi><mi>x</mi></msub><mo>|</mo><mo>&le;</mo><msub><mi>&theta;</mi><mrow><mi>x</mi><mn>0</mn><mi>max</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000021.GIF" wi="1470" he="472" /></maths>式中,θ<sub>x</sub>为水平视线角:<img file="FDA0000600596620000022.GIF" wi="462" he="139" />θ<sub>x1max</sub>为最佳视野内最大水平视线角,θ<sub>x0max</sub>为最大视野内最大水平视线角;C<sub>l</sub>(X<sub>P</sub>)为可视度视距函数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>l</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>l</mi><mo>=</mo><msub><mi>l</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0.2</mn><mo>&CenterDot;</mo><mfrac><mrow><mi>l</mi><mo>-</mo><msub><mi>l</mi><mn>1</mn></msub></mrow><mrow><msub><mi>l</mi><mn>1</mn></msub><mo>-</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>min</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>+</mo><mn>1</mn></mtd><mtd><msub><mi>l</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub><mo>&le;</mo><mi>l</mi><mo>&le;</mo><msub><mi>l</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>-</mo><mn>0.2</mn><mo>&CenterDot;</mo><mfrac><mrow><mi>l</mi><mo>-</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub></mrow><mrow><msub><mi>l</mi><mn>1</mn></msub><mo>-</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub></mrow></mfrac><mo>+</mo><mn>1</mn></mtd><mtd><msub><mi>l</mi><mn>1</mn></msub><mo>&le;</mo><mi>l</mi><mo>&le;</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0.8</mn><mo>&CenterDot;</mo><mfrac><mi>l</mi><msub><mi>l</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub></mfrac></mtd><mtd><mn>0</mn><mo>&lt;</mo><mi>l</mi><mo>&le;</mo><msub><mi>l</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0.8</mn><mo>&CenterDot;</mo><mfrac><mrow><mi>l</mi><mo>-</mo><msub><mi>l</mi><mrow><mn>0</mn><mi>max</mi></mrow></msub></mrow><mrow><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>l</mi><mrow><mn>0</mn><mi>max</mi></mrow></msub></mrow></mfrac></mtd><mtd><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub><mo>&lt;</mo><mi>l</mi><mo>&le;</mo><msub><mi>l</mi><mrow><mn>0</mn><mi>max</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>l</mi><mo>></mo><msub><mi>l</mi><mrow><mn>0</mn><mi>max</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000023.GIF" wi="1446" he="749" /></maths>式中,l<sub>1</sub>为最优视距,l<sub>1min</sub>和l<sub>1max</sub>为正常作业视距的最小值和最大值,l<sub>0max</sub>表示最大作业视距;二、通过人因分析验证人员位置:依据国家标准“中国成年人人体尺寸”(GB10000‑1988),取第90百分位的人体手臂尺寸作为设计参考,确定人体前端面与设备面板的距离,保证所有内部待修设备满足可视性和可达性要求;三、计算满足应力强度的最小倒圆:开口会产生应力集中,需要对存在直角过渡的口盖进行倒圆,使得最大应力满足强度约束;分析口盖应力情况,按照1mm为步长确定倒圆角尺寸,建立口盖应力分析的有限元模型,进行有限元分析;在保留尺寸余量Δ的情况下确定满足应力集中要求的最小倒圆半径η;四、基于人因势场理论优化求解面积最小情况时的口盖尺寸参数:(1)基于人因势场的口盖尺寸优化模型的建立建立基于人因势场的口盖尺寸优化模型,包括变量的确定与描述、目标函数的建立和约束条件:①变量的确定与描述考虑设备的尺寸和位置参数;设计变量为设备的位置参数X,即:X={X<sub>e</sub>,X<sub>s</sub>,X<sub>h</sub>}={(x<sub>e</sub>,y<sub>e</sub>,z<sub>e</sub>),(x<sub>1</sub>,…,x<sub>i</sub>,…x<sub>n</sub>),(x<sub>h</sub>,y<sub>h</sub>,z<sub>h</sub>)}    (5)其中,X<sub>e</sub>=(x<sub>e</sub>,y<sub>e</sub>,z<sub>e</sub>)为描述人眼睛的位置参数,X<sub>h</sub>=(x<sub>h</sub>,y<sub>h</sub>,z<sub>h</sub>)为口盖的中心坐标,X<sub>s</sub>=(x<sub>1</sub>,…,x<sub>i</sub>,…x<sub>n</sub>)(i=1,2,…n)为描述口盖的形状特征参数,n为形状特征参数个数;②约束条件的处理口盖的形状设计过程中的约束条件有:a.人的视线能够通过口盖观察到设备关键人因点;以Ω表示口盖经过人的视线投影在设备底面内部的各点集合(区域);P<sub>i</sub>(x<sub>i</sub>,y<sub>i</sub>,z<sub>i</sub>)表示各关键人因点的坐标;则有:P<sub>i</sub>(x<sub>i</sub>,y<sub>i</sub>,z<sub>i</sub>)∈Ω,i=1,2,…m                (6)其中m表示设备关键人因点的个数;b.观察点的人因势能值应不高于要求的极限值;较高的人因势能值意味着较差的维修性;形状优化的人因势能值约束为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>U</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><msub><mi>U</mi><mi>Fj</mi></msub><mo>&le;</mo><msub><mi>U</mi><mrow><mi>F</mi><mi>max</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000041.GIF" wi="1152" he="142" /></maths>其中,U<sub>Fj</sub>、ω<sub>j</sub>为设备第j个关键人因点复合人因势能值及其权重值,m表示设备关键人因点的个数;U<sub>Fmax</sub>为要求的极限人因势能值;ω<sub>j</sub>之间满足归一化条件:<img file="FDA0000600596620000042.GIF" wi="268" he="136" />(0<ω<sub>j</sub><1);U<sub>F</sub>为复合人因势场值;③目标函数的确立定义设备的关键人因点,优化设计变量,使人的眼睛通过口盖能够观察到关键人因点,并且满足复合人因势场势能值不大于0.2的要求;按照上面第三步的口盖强度要求,保证最小倒圆半径为η;同时希望口盖愈小愈好,也即口盖的面积最小:minS=minS(x<sub>1</sub>,…,x<sub>i</sub>,…x<sub>n</sub>)              (8)(2)求解最优的口盖尺寸参数采用自适应粒子群算法进行求解;具体过程如下:①采用Matlab软件初始化模型参数,包括维修设备相关的尺寸和位置、人体站立的姿势及空间位置;②根据公式(1)‑(8)建立基于复合人因势场、以口盖面积为目标的开口尺寸优化模型;采用Matlab软件初始化随机产生符合约束条件的粒子群位置和速度;计算初始化粒子群中在尺寸优化模型中最优值,确定初始最优位置参数;③更新速度和粒子群位置并检验是否满足约束条件,如不是继续更新速度和粒子群位置直至满足条件;粒子的速度和位置根据如下方程进行更新:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>v</mi><mi>id</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>v</mi><mi>id</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msubsup><mi>rand</mi><mn>1</mn><mi>k</mi></msubsup><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mi>id</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>id</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>&CenterDot;</mo><msubsup><mi>rand</mi><mn>2</mn><mi>k</mi></msubsup><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mi>gd</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>id</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000043.GIF" wi="1567" he="100" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>id</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>id</mi><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mi>id</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000600596620000044.GIF" wi="1238" he="94" /></maths>式中,<img file="FDA0000600596620000051.GIF" wi="149" he="82" />和<img file="FDA0000600596620000052.GIF" wi="144" he="82" />为(0~1)内均匀分布的随机数;<img file="FDA0000600596620000053.GIF" wi="67" he="83" />是粒子i在k次迭代中第d维的当前位置;c<sub>1</sub>和c<sub>2</sub>是学习因子或加速系数;p<sub>id</sub>为粒子群经历的当前最优值;p<sub>gd</sub>则为全局最优值;④计算面积目标函数,以目标函数变化率作为停止准则,检验是否满足停止准则,循环更新粒子速度和位置直至满足停止准则,得到最优结果。
地址 410073 湖南省长沙市开福区德雅路109号