宇航机构与结构大作业2使用matlab计算太阳冀展开的蜗卷弹簧

已知：一太阳翼由连接架、三块基板组成，如图1所示。图1各板间采用2个平面涡卷弹簧展开铰链连接，连接架及基板质量、尺寸参数如表1所示。表1质量参数连接架质量基板1质量基板2质量基板3质量=8.615kg=9.342kg=9.900kg=9.257kg尺寸参数连接架长度基板长度连接架偏心=0.704m=1.129m=0.172m惯性半径连接架惯性半径基板1惯性半径基板2惯性半径基板3惯性半径=0.3462m=0.6518m=0.7007m=0.6518m根部铰链、板间铰链展开结束时刻的阻力矩和各铰链阻力矩变化常数如表2所示。表2铰链位置根部铰链板间铰链（B轴线）板间铰链（C轴线）板间铰链（D轴线）展开结束时刻阻力矩-0.3817Nm-0.1946Nm-0.20598Nm-0.15685Nm阻力矩变化常数0.1611Nm/rad0.0428Nm/rad0.0573Nm/rad0.0312Nm/rad（1）设计根部铰链及各板间铰链的涡卷弹簧参数，并使展开结束时刻产生的冲击尽可能小；（2）用matlab进行编程求解，绘制太阳翼展开过程中，展开铰链随时间变化的角速度、角加速度曲线。注：提交作业中包括完整的计算过程和带有注释的程序源代码，及对计算结果的相关说明。平面涡卷弹簧设计所设计的太阳冀展开机构有三块板，有一条根部铰链线，三条板间铰链线。其中，根部铰链线上有一个用涡卷弹簧驱动的铰链，三条板间铰链线上各有两个铰链。MACROBUTTONMTEditEquationSection2EquationChapter(Next)Section1SEQMTEqn\r\hSEQMTSec\r1\hSEQMTChap\h满足太阳冀展开状态下的静态扭矩裕度要求主动驱动力矩： 其中和是第i条铰链线上的两个铰链的驱动力矩。 其中为第i条铰链线上的总的驱动力矩。 其中为根部铰链的驱动力矩。 总的驱动力矩为： 阻力矩：根部铰链线的阻力矩： 其中，为根部铰链的阻力矩，为根部电缆的总阻力矩，为联动机构对根部铰链线作用产生的总阻力矩，为根部铰链的温度影响系数，为电缆的温度影响系数，为联动机构的温度影响系数。各板间铰链线上的阻力矩： 其中，为第i条铰链线上板间铰链的阻力矩，，分别为第i条铰链线上的上、下铰链的阻力矩。 静力矩裕度： 满足涡卷弹簧强度要求弹簧在工作中的最大主动力矩满足强度要求： 其中，为材料的强度极限。 展开速度要求展开结束时刻角速度是衡量展开冲击的一个重要标志，要求控制锁定冲击载荷不超过SADA许用值，但是在本作业中不将速度要求作为考虑范围。平面涡卷弹簧的计算平面涡卷弹簧的计算是一个首先初步估算参数值，再检验设计要求，然后修正计算，再检验设计要求，反复若干次后最终确定参数的过程。通过反复计算平面涡卷弹簧的宽度、厚度、圈间距等参数，最终得到同时满足前面所述要求的弹簧设计参数。其计算过程如REF_Ref353826873\h图1所示图SEQ图\*ARABIC1平面涡卷弹簧参数设计过程设计载荷的确定假设根部铰链可以提供的主动力矩是一条板间铰链线上可以提供力矩的两倍，则一个根部铰链可以提供的主动力矩是一个板间铰链可以提供的主动力矩的四倍。即MACROBUTTONMTEditEquationSection2EquationSection(Next)SEQMTEqn\r\hSEQMTSec\h 将GOTOBUTTONZEqnNum234822REFZEqnNum234822\*Charformat\!(2.1)式代入GOTOBUTTONZEqnNum477261REFZEqnNum477261\*Charformat\!(1.10)可以得到 其中，为初始计算时刻板间铰链平面涡卷弹簧在展开结束时所应提供的最小主动力矩，为初始计算时刻根部铰链平面涡卷弹簧在展开结束时所应提供的最小主动力矩。而和则是由展开结束时的SADA许用值确定，可以根据展开结束时的最大速度来确定，不在本次大作业的设计范围中。根据上述理论和已知的阻力矩条件，可以选择主动力矩分别为给出的阻力矩的两倍，由此可得、、、。要使设计的弹簧参数使展开结束时产生的冲击尽可能小，将上述值取等号，可得：、、、。圆整后取、。（单位为Nm）平面涡卷弹簧参数的计算依据型号设计原型的经验，可以初步确定平面涡卷弹簧的材料，弹簧的外径、内径和圈间距。根据已知的展开结束时刻阻力矩与阻力矩变化常数可知根部铰链与板间铰链大约有两倍的关系，为了简化平面涡卷弹簧的设计，只需设计根部铰链处的涡卷弹簧与板间铰链的涡卷弹簧。为了保证满足平面涡卷弹簧的强度要求，所选择的和的值应当满足如下要求 其中，和的值分别为和的值，和分别为根部弹簧与板间弹簧系数。由GOTOBUTTONZEqnNum286215REFZEqnNum286215\*Charformat\!(2.4)、GOTOBUTTONZEqnNum889593REFZEqnNum889593\*Charformat\!(2.5)和GOTOBUTTONZEqnNum561238REFZEqnNum561238\*Charformat\!(1.11)式可以得到平面涡卷弹簧的最小厚度 其中，和分别为根部铰链和板间铰链平面涡卷弹簧的最小厚度。根据已知的阻力矩变化常数可以取弹簧系数大于等于阻力矩变化常数，故取和，则其中为1.0827Nm，为0.4688Nm。选取弹簧宽度为，计算得，，根据JB/T7366-1994取根部铰链弹簧厚度为，板间铰链弹簧厚度为。平面涡卷弹簧的刚度 则平面涡卷弹簧的工作长度 弹簧有效作用圈数为 其中，为平面涡卷弹簧的长度，为材料的弹性模量。平面涡卷弹簧的节距： 其中，为弹簧外圈直径，为弹簧内圈直径。计算所得弹簧需满足式GOTOBUTTONZEqnNum561238REFZEqnNum561238\*Charformat\!(1.11)强度要求。 取弹簧外圈直径，弹簧内圈直径，弹簧材料为热处理弹簧钢带II级，根据JB/T7366-1994可知该材料抗拉强度为。可以得到,取，弹性模量，取用matlab编程clcclearbg=10;bb=10;hg=0.8;hb