一种平面结构设计方法技术

本发明专利技术涉及一种平面结构设计方法，初始控制点输入以生成多片初始模型，再拼接为一片，并提高模型的连续性；通过分析模块，分析该初始模型的刚度矩阵和位移，计算出设计变量的值；将初始控制点导入到优化模块中，求解出目标函数的结果；判断目标函数结果是否为最优设计变量值，求解目标函数的结果；对模型进行细化处理，以获取最优化显示控制点；将最优化控制点进行应力云图显示。在建模中解决了分片拼接的连续性问题，使得模型无论是位移还是应力的变化，都呈现平稳过渡；在优化中运用求解等效应力和等效应力对设计变量的导数；将灵敏度方法作为一个桥梁，把等效应力对设计变量的导数和控制点的位移连接在了一起。使得几何分析优化更简洁合理。

A design method of plane structure

The invention relates to a plane structure design method, in which the initial control point is input to generate multiple initial models, and then spliced into one, and the continuity of the model is improved; through the analysis module, the stiffness matrix and displacement of the initial model are analyzed, and the value of the design variable is calculated; the initial control point is imported into the optimization module, and the result of the objective function is solved; the objective function knot is judged If it is the optimal design variable value, solve the result of objective function; refine the model to obtain the optimal display control point; display the optimal control point by stress nephogram. In the modeling, the continuity problem of piecewise splicing is solved, so that the model presents a smooth transition no matter the change of displacement or stress; in the optimization, the derivative of the equivalent stress and the equivalent stress to the design variable is solved; the sensitivity method is used as a bridge, and the derivative of the equivalent stress to the design variable is connected with the displacement of the control point. It makes the optimization of geometric analysis more concise and reasonable.

【技术实现步骤摘要】
一种平面结构设计方法
本专利技术涉及一种结构设计技术，特别涉及一种平面结构设计方法。
技术介绍
机械设计是机械工业中最重要的环节，就是在各种限定的条件下(如材料，加工，理论知识和计算手段等)下设计出最好的机械，即做出优化设计。过去，设计的优化主要依靠设计者的知识，经验和远见。随着机械工程理论发展和计算机设计软件的应用，使机械设计更加高效，精准。之后等几何理论的发展，使得机械设计优化更加的简单。等几何分析的优点在于，将一个复杂难以表达的模型轻松的表达出来。现在流行的等几何分析建模方法主要有两种：裁剪法和分片法。裁剪法的缺点在于，虽然可以根据参数空间对模型进行修改，再映射到物理空间，但是处理边界问题时，就显得非常的复杂，对于接下来的分析优化，有着一定的难度。分片法是最近比较流行的方法，优点在于，可以将复杂的模型分解成简单的几片，表达出来之后，再进行拼接计算。现在分片法的技术问题在于，如何提高连续性，另在等几何分析里，优化的目标函数研究也是至关重要的，模型的集中应力的优化也必不可少。
技术实现思路
本专利技术是针对现在优化设计存在的问题，提出了一种平面结构设计方法，几何分析优化更简洁合理，并解决现在分片法中的连续性问题。本专利技术的技术方案为：一种平面结构设计方法，具体包括如下步骤：1)观察需要设计的模型，分解模型，将模型分解成常规可以表达的单片，再通过需要考虑的应力分布输入单片初始控制点以生成单片的初始模型，再通过程序使单片的控制点一一编号，将相邻的两片拼接为一片，在拼接完成的两片里，插入连续性控制点，再通过改变边界和内部曲线的节点矢量，表现出在图形上显示等参线是圆滑的；2)利用环环相扣的超参单元对传统等参单元进行替换，在参数空间计算形函数R对一个控制点坐标的xy两个方向的导数，通过雅克比矩阵的转换，求解出单元应力应变矩阵B，再通过单元应力应变矩阵B求解出单元刚度矩阵K，最后根据公式u＝K/f求解出单元中各个控制点的位移，u为控制点坐标列阵；K为单元刚度矩阵；f为载荷列阵；3)将等效应力作为目标函数，优化等效应力：采用形状尺寸优化，建立目标函数的解和控制点位置坐标的关系，将根据需要优化的控制点作为初始控制点，进行优化，得到优化后的控制点坐标；4)将步骤3)求得的优化后的控制点坐标输入，计算出控制点对应模型的等效应力，判断此等效应力是否是等效应力函数的极小值，如果不是，则返回步骤3)对控制点进行重新优化计算，反复进行迭代计算，直到得到等效应力函数的极小值，也就是最优控制点坐标，退出循坏，进行下一步；5)对模型进行细化处理，模型划分的更多的单元；6)根据最优控制点坐标求解出每个单元的等效应力，按照数值对应相应的颜色，再将其填充到模型里对应的细分后的每个单元位置，进行模型应力的显示。所述步骤3)中优化目标为结构的等效应力最小，约束是结构的面积小于等于给定值，优化方程如下：minC(|αi|)＝fTus.t.V≤V*[αi]min≤[αi]≤[αi]maxK/u＝f其中，C为优化目标等效应力；[αi]为优化设计变量序列，在此为优化形状边界控制点坐标；V为优化迭代中模型体积；V*为模型体积上限值；[αi]min、[αi]max分别为设计变量的下限值和上限值,f＝K/u为离散控制方程；K为单元刚度矩阵；u为控制点坐标列阵；f为载荷列阵。本专利技术的有益效果在于：本专利技术平面结构设计方法，将整个等几何分析优化，分成了几个6步，使得程序更加的简洁有逻辑,为以后程序的丰富打好了基础；在建模步骤里，解决了分片拼接的连续性问题，使得模型无论是位移还是应力的变化，都呈现平稳过渡；在优化步骤里，运用了一种新的方法求解等效应力和等效应力对设计变量的导数；将灵敏度方法作为一个桥梁，把等效应力对设计变量的导数和控制点的位移连接在了一起。附图说明图1为本专利技术平面结构的设计方法流程图；图2为本专利技术方法实施例中需要进行设计的模型图；图3为本专利技术图2中模型的分片图；图4为本专利技术图3中两片之间的拼接和连续性处理图；图5为本专利技术模型全部拼接后的等参线图；图6为模型的等效应力图；图7为本专利技术模型优化后的等效应力图。具体实施方式如图1所示平面结构的设计方法流程图，具体步骤如下：步骤S1：观察需要设计的模型，分解模型，将模型分解成常规可以表达的单片，如图2所示需要进行设计的模型图，如图3所示模型的分片图，再通过需要考虑的应力分布输入单片初始控制点以生成单片的初始模型，再通过程序使单片的控制点一一编号，将相邻的两片拼接为一片，在拼接完成的两片里，插入连续性控制点，再通过改变边界和内部曲线的节点矢量，达到c1连续。c1连续是指模型内部的nurbs基函数可以进行二次求导，表现出在图形上显示等参线是圆滑的，不会出现尖点的情况。图4是图3中的patch1和patch2两片之间的拼接和连续性处理图，进行拼接的时候，先将他们内部的控制点进行清除，拼接完成后再在他们的内部进行插入控制点，这里插入要求，在内部边界上不能存在公共控制点，因为nurbs在公共控制点上会有尖点，c1不连续。步骤S2：利用环环相扣的超参单元对传统等参单元进行替换，在参数空间计算形函数R对一个控制点坐标的xy两个方向的导数，通过雅克比矩阵(多变数向量函数的最佳线性逼近)的转换，求解出单元应力应变矩阵B，再通过单元应力应变矩阵B求解出单元刚度矩阵K，K＝∫Ω(B)TDBdΩ，Ω表示单元；D是弹性矩阵，属于常数矩阵；然后通过标准方法组装求解出总刚度矩阵，最后根据公式u＝K/f(u为控制点坐标列阵；K为单元刚度矩阵；f为载荷列阵)求解出单元中各个控制点的位移；如图5所示为模型全部拼接后的等参线图。步骤S3：等效应力可以清晰描述出一种结果在整个模型中的变化，从而使分析人员可以快速的确定模型中的最危险区域。因此将等效应力作为目标函数，优化等效应力。将根据需要优化的控制点作为初始控制点，带入相应公式，求出等效应力函数的解。由于本专利技术采用的优化是形状尺寸优化，所以需要建立目标函数的解和控制点位置坐标的关系。也就是将等效应力的导数和控制点坐标建立联系(将(4)式里的设计变量，换成了控制点坐标)。将求得的参数(参数为等效应力对控制点坐标的导数)输入到MMA程序中进行优化计算，得到优化后的控制点坐标。步骤S4：将步骤S3求得的控制点坐标输入到判断语句中，判断语句根据计算，计算出控制点对应的模型的等效应力，判断这时的等效应力是否为等效应力函数的极小值，如果不是，则返回步骤S3对控制点进行重新优化计算，反复进行迭代计算，直到得到等效应力函数的极小值，也就是最优解。退出循坏，进行下一步。步骤S5：细化对模型的质量有着重要的作用，细化次数越多，模型划分的单元就越多，求得的解就越逼近理论解。但是相应的，计算量会大大增加。所以对模型进行细化处理时，要选择合适的细化数，这里h细化和p细化(等几何分析中两种处理有限元单元的方法)本文档来自技高网...

【技术保护点】
1.一种平面结构设计方法，其特征在于，具体包括如下步骤：/n1)观察需要设计的模型，分解模型，将模型分解成常规可以表达的单片，再通过需要考虑的应力分布输入单片初始控制点以生成单片的初始模型，再通过程序使单片的控制点一一编号，将相邻的两片拼接为一片，在拼接完成的两片里，插入连续性控制点，再通过改变边界和内部曲线的节点矢量，表现出在图形上显示等参线是圆滑的；/n2)利用环环相扣的超参单元对传统等参单元进行替换，在参数空间计算形函数R对一个控制点坐标的xy两个方向的导数，通过雅克比矩阵的转换，求解出单元应力应变矩阵B，再通过单元应力应变矩阵B求解出单元刚度矩阵K，最后根据公式u＝K/f求解出单元中各个控制点的位移，u为控制点坐标列阵；K为单元刚度矩阵；f为载荷列阵；/n3)将等效应力作为目标函数，优化等效应力：采用形状尺寸优化，建立目标函数的解和控制点位置坐标的关系，将根据需要优化的控制点作为初始控制点，进行优化，得到优化后的控制点坐标；/n4)将步骤3)求得的优化后的控制点坐标输入，计算出控制点对应模型的等效应力，判断此等效应力是否是等效应力函数的极小值，如果不是，则返回步骤3)对控制点进行重新优化计算，反复进行迭代计算，直到得到等效应力函数的极小值，也就是最优控制点坐标，退出循坏，进行下一步；/n5)对模型进行细化处理，模型划分的更多的单元；/n6)根据最优控制点坐标求解出每个单元的等效应力，按照数值对应相应的颜色，再将其填充到模型里对应的细分后的每个单元位置，进行模型应力的显示。/n...

【技术特征摘要】
1.一种平面结构设计方法，其特征在于，具体包括如下步骤：
1)观察需要设计的模型，分解模型，将模型分解成常规可以表达的单片，再通过需要考虑的应力分布输入单片初始控制点以生成单片的初始模型，再通过程序使单片的控制点一一编号，将相邻的两片拼接为一片，在拼接完成的两片里，插入连续性控制点，再通过改变边界和内部曲线的节点矢量，表现出在图形上显示等参线是圆滑的；
2)利用环环相扣的超参单元对传统等参单元进行替换，在参数空间计算形函数R对一个控制点坐标的xy两个方向的导数，通过雅克比矩阵的转换，求解出单元应力应变矩阵B，再通过单元应力应变矩阵B求解出单元刚度矩阵K，最后根据公式u＝K/f求解出单元中各个控制点的位移，u为控制点坐标列阵；K为单元刚度矩阵；f为载荷列阵；
3)将等效应力作为目标函数，优化等效应力：采用形状尺寸优化，建立目标函数的解和控制点位置坐标的关系，将根据需要优化的控制点作为初始控制点，进行优化，得到优化后的控制点坐标；
4)将步骤3)求得的优化后的控制点坐标输入，计算出控制点对应模型的等效应力，判断此等...

【专利技术属性】
技术研发人员：陈龙，徐力，王恺，
申请(专利权)人：上海理工大学，
类型：发明
国别省市：上海;31