数学建模网络挑战赛
承 诺 书
我们仔细阅读了第五届“认证杯”数学中国数学建模网络挑战赛的竞赛规则。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。
我们允许数学中国网站(www.madio.net ) 公布论文，以供网友之间学习交流，数学中国网站以非商业目的的论文交流不需要提前取得我们的同意。
我们的参赛队号为： 参赛队员 (签名) ：
队员2：
队员3：
参赛队教练员 (签名) ：
参赛队伍组别：
队员1：
数学建模网络挑战赛 编 号 专 用 页
参赛队伍的参赛队号：（请各个参赛队提前填写好）：
竞赛统一编号（由竞赛组委会送至评委团前编号）：
竞赛评阅编号（由竞赛评委团评阅前进行编号）：
2012年第五届“认证杯”数学中国
数学建模网络挑战赛
题 目 对数螺线型蜘蛛网状的结构分析
关 键 词 蜘蛛网 对数螺线 蒙特卡洛方法 ANSYS分析法
摘要
本文针对蜘蛛网合适结构的问题，考虑吐丝量一定，外界环境较理想条件下，建立以对数螺线为核心的数学模型，追求蜘蛛网结构最优。运用蒙特卡洛方法，模拟昆虫触网的过程，考虑了在蜘蛛丝长度一定的条件下，对数螺旋比圆围成的面积大，但疏而不漏，应用随机过程近似昆虫触网的过程，得出了对数螺线更利于捕食的结论。另一方面，也对对数螺线型面联接理论和联接界面强度进行了分析与计算，利用ANSYS 进行接触分析，得出了对数螺线型面联接的接触应力和接触强度条件的表达式。采用随机数产生算法，利用MATLAB 7.0.1和C++编程，分别对模型进行求解，并对所得结果进行分析比较，以此来帮助设计最有蜘蛛网结构。
参赛队号
所选题目 A
Abstract
Our article aims to study the question about the best structure of the spider webs ,it is on the condition of certain output of the spinning the and quite ideal conditions ,establish mathematical model in the core of the logarithmic spiral to find the best way of the spider webs .We also analyze Logarithm of solenoid type surface connection theory, Interface connection strength and ANSYS to get the expression. we apply Monte Carlo method to simulate the process about Net insert and adopt the Random number produce algorithm ,we also use the software of Matlab 7.0.1 、Mathematica and Microsoft Visual C++ 6.0 to give the answer to the question about the model and analyze about the result from model ,so we establish the best structure of the spider webs by means of these datas.
一、问题重述
世界上生存着许多种类的蜘蛛，而其中的大部分种类都会通过结网来进行捕食。通过对蜘蛛网所形成的结构的分析，通过建立模型，设计一种更为合适的蜘蛛网结构。
二、问题分析
题目中主要研究的是：蜘蛛网织成怎样的结构才是最合适的。因此我们通过查阅资料，了解蜘蛛网的结构等方面内容，根据结构形状的不同, 蛛网可以分为片网、不规则网和圆网等几种类型。由于圆网在蛛网进化上的地位特殊，且结构简单、规则。因此，到目前为止对蛛网的研究大都集中在圆网上。圆网并不是标准的圆，而是近似于数学上的螺旋线。所谓合适的蜘蛛网结构就是利于蜘蛛捕食、防御、繁殖。根据题目的要求，我们提出以下几个问题：
1、为什么蜘蛛网是螺旋线状，而不是标准的同心圆； 2、对数螺线型面联接理论和联接界面强度分析与计算。
三、符号说明
ρ：模拟对数螺线的极径
θ：模拟对数螺线的极角
Φ：对数螺线型型面轴旋转角
σp ：面轴上接触应力
l ： 轴孔之间的轴向配合长度
γ：P 点接触应力σp 与极径ρ之间的夹角 f ：轴孔之间的摩擦系数 T ： 扭矩
S 1：面积就是最大圆的面积 L 1：四圈的长度为
四、模型假设：
1、假设蜘蛛网是规则的对数螺线；
2、不考虑蜘蛛网受到风雨等天气情况的影响;
3、假设昆虫飞向蜘蛛网时，落在网内每点的概率相同；
五、模型建立与求解：
蜘蛛网的中心和圆周之间呈辐射状的半径线，自外向里是螺旋线，愈近中心，每圈间的距离也愈小，直到不可辨认的地步，这正符合数学上的对数螺线的情况。因此，我们建立对数螺线的模型，近似代替蜘蛛网，研究其性质。
图1
㈠ 对数螺线的定义和性质
数学上对数螺线定义如下：动点的运动方向始终与极径保持定角θ的动点轨迹，称为对数螺线。如图1所示，其极坐标方程为：
ρ=ae m θ （1）
式中：a , m 为常数（λ=arctan (1/m )））；θ为极角，ρ为极径。
图2 对数螺线
对数螺线在渐屈、渐伸、垂迹、回光线等各种变换下的不变性质, 体现出自身的高度和谐、对称和统一性。 ㈡对数螺旋线与圆形蜘蛛网的比较
将四个标准圆形与对数螺旋线放入同一坐标系中，如下图3
图3
1、 四圈圆形蜘蛛网
面积就是最大圆的面积：S 1=82π=64⨯3.14=201.056 四圈的长度为： L 1=2(1+3.6+5.8+8.2)π=111.784 2、四圈螺旋线蜘蛛网
该对数螺旋线的方程为：ρ=0.02e 3θ，0≤θ≤8π
面积：由于对数螺旋线是一条不封闭的曲线，所以用下图中最外面的曲线和一条线段组成的封闭图形表示该螺旋线所包围的面积。通过数该封闭图形内的方格数估计面积，不足一格按半格记。
图4
共有188个正方形，45个不足一格的，所以面积为198+22.5=220.5 长度：
用Mathematica 计算该曲线长，
输入 Integrate[0.02e^3x,{ x,0,8Pi}] ， 得出结果L= 126.871
结论：
在蜘蛛丝长度一定的情况下，螺旋线所围成蜘蛛网的面积大。这样更利于蜘蛛捕食。
蒙特卡罗方法：
用蒙特卡罗方法模拟昆虫飞向蜘蛛网上的过程：
假设昆虫飞向蜘蛛网时是一个随机过程，此过程中不考虑环境因素（风向、风速等）的影响
编写C 语言程序，生成二维随机数。 程序及运行结果见附录。
将这些随机数在下图5中描点，为处理简单，以第一象限为例，其他象限相同。
图5
图5.1 图5.2 图5.3 图5.4 图5.5
由此可见，昆虫更可能碰到对数螺线。也就是说，对数螺线形的蜘蛛网更有利于捕食
㈢、对数螺线型面联接理论和联接界面强度分析与计算 3.1对数螺线曲线
以型面轴截面曲线为例, 选取三段对数螺线进行分析。θ从0°到110°的一段曲线组成, 其中各段曲线之间用直线圆滑联接, 减弱了应力集中现象。
图6三段对数螺线型面轴截面 2.2对数螺线曲线参数的确定
在图1所示对数螺线的方程ρ=ae m θ中, m的大小取决于型面联接轴与孔之间受力时的压力角α, 因α在整个曲线上是常数, 当α选定时, m 为常数。因此选择型面联接轴、孔截面形状时, 可选定a 值的大小来定轴、孔的尺寸, 再定出压力角, 即可确定式中m 值, 对数螺线形状也就随之确定, 同时型面轴、孔截形也就相应得到确定。
三段对数螺线组成的型面轴、孔有两个基本参数a 和α, 其中m 主要影响曲线的形状, a 主要影响曲线的大小。三段对数螺线型面联接如图7所示, 这是一种有间隙的配合(图7-a), 工作时通过一定量的相对旋转, 间隙补偿, 由于楔面的作用在接触面之间便产生正压力, 并摩擦闭锁, 形成可靠的联接(见图7-b), 其摩擦受力方向与运动方向所成角度a 的大小与参变量θ无关, 即对数螺线的压力角a 在任何位置都是相同的。
( a ) ( b) 图7对数螺线型面联接
3.2 型面轴孔接触初始位置确定
设型面轴孔的截面廓形曲线方程分别为ρ1=a 1e m θ、ρ2=a 2e m θ。假定固定轮毂，顺时针旋转型面轴Φ角后两者之间有初始接触，任取一接触点P ，则在P 点处有ρ1=ρ2，
即 a 1e m (θ+Φ)=a 2e m θ
可得 Φ=（2） （3） 1⎛a 2⎫ln ⎪ m ⎝a 1⎭
上式表明型面轴旋转角Φ仅与曲线的a , m 常数值有关，而与θ值无关，即与初始接触点P 的位置无关，即Φ为定值，说明轴孔工作表面之间同时发生接触，轴孔之间接触为面接触。
3.3利用ANSYS 进行接触分析
由于ANSYS 对复杂曲面建模具有一定的局限性，为了分析的准确性，利用Pr o /E 强大的三维建模功能，在Pr o /E 中建成模型后，再利用Pr o /E ,ANSYS 之间的接口程序，将模型导入, ANSYS 中。为了方便加载扭矩，在不影响分析结果性质的前提下，于轴中心建一半径为r 的圆形孔，在圆形孔边界节点处加载等效切向力F ，使得扭矩T=Fr，如图8所示。按照ANSYS 分析步骤设置属性、划分网格、加载、求解。在后处理器POST1中查看有限元模型在纯扭矩T 作用下的节点应力云图，如图9所示。
图8有限元模型
图9节点的应力云图
图9表明在轴孔各段工作表面上除了两端（刚进入接触与刚脱离接触的很小的一个区域）应力比较大之外，其它的区域所受到的应力分布都非常均匀。若在模型应力比较集中区域进行修磨处理，将各段曲线间联接用更圆滑的曲线过渡，可减少或消除应力集中现象。由此，认为对数螺线型面联接接触应力均匀分布是符合实际情况的。反复加载不同的扭矩对有限元模型进行分析求解，发现作用扭矩与接触应力之间成线性比例关系。
3.4接触应力计算
假定在无间隙无过盈的理想配合状态下，轴上作用纯扭矩，如图5所示。在对数螺线型面轴上工作表面共有三段，其上接触应力均匀分布，记为σp 。在AB
表面上任取一点P ，在包含P 点的微小弧段ds 上，接触应力的推动分力（即与该点的线速度平行的分力）及摩擦力对轴线o 之矩dt 应为
dT =ρ(σp sin γlds +σp f cos γlds )
=(sin γ+f cos γ)σp ρlds （4） 式中:ρ为P 点的极径; l 为轴孔之间的轴向配合长度（即工作长度）；γ为P 点接触应力σp 与极径ρ之间的夹角；f 为轴孔之间的摩擦系数（一般取0.1~0.15）。
由对数螺线性质知γ=α=π/2-λ, ρ=ae m θ, ds =ρd θ, 代入式（4）得
dT =(sin γ+f cos α)ρ2σp ld θ
=(sin α+f cos α)ρ2σp le 2m θd θ (5) 在传动轴匀速转动时，由型面轴廓形曲线的对称性知，求出一段工作表面上力对轴线o , 产生的力矩，再乘以3所得的总力矩应与作用扭矩T 平衡，即有 T =3⎰
211π180dT =3⎰11π180a σp l (sin α+f cos α)e 22m θd θ（6） m π⎛119⎫3a σp l (sin α+f cos α) e -1⎪⎝⎭=2m
即σp =2mT
⎛3a 2l (sin α+f cos α) e ⎝11m π9⎫-1⎪⎭ （7）
上式表明作用扭矩T 与σp 之间成线性比例关系，当作用扭矩一定时，接触
应力的大小取决于轴孔之间的轴向配合长度l 以及a , α等参数，增大a , l 或α均可减小对数螺线型面联接的接触应力σp 。
由于在实际中蜘蛛网面接触应力并不均匀。对于型面联接，要求在轴孔配合表面存在一定间隙，所以前面无间隙无过盈的假定与实际情况有一定出入。由于间隙的存在，各段工作表面并不能实现在其理论区域上全部接触。这样按式（7）计算的σp (小于其实际值。为使计算的σp 更加接近最大接触应力，可引入修正系数加以修正。因此对数螺线型面联接的接触强度条件可写成
σp =2kmT
⎛3a 2l (sinα+f cos α) e ⎝11m π9⎫-1⎪⎭≤[σ]p （8）
六、模型结果分析
通过建立对数螺线型近似蜘蛛网结构，运用蒙特卡洛方法模拟猎物触网过程，我们给出了在蜘蛛丝长度一定的情况下，螺旋线所围成蜘蛛网的面积大，更有利于捕食。对对数螺线型面轴廓形曲线的分析，在初步确定对数螺线型面轴孔联接初始接触位置及接触性质的基础上，应用有限元分析软件ANSYS 对对数螺线型面联接进行接触应力分布规律的研究不失为一种行之有效的方法。本文分析求解的接触应力的理论计算公式，由于各种综合因素考虑的不足，从本质上看，其计算结果与实际情况会有一定出入。然而该公式作为一种有效的实用算法，结构简单，使用方便，为此种联接在工程中的应用打下了良好的设计基础。
七、模型优缺点分析
本文中针对研究蜘蛛网的合适结构给出了对数螺线型模型，本模型较好的近似了蜘蛛网的结构，对于对数螺线型网状结构的优越性进行了较详细的说明，并给出了对数螺线型面联接理论和联接界面强度分析与计算。较好地说明了对数螺线型蜘蛛网是较为适合的蛛网结构。但是由于在模型建立时的诸多理想假设、主观原因和参数敏感度的影响，致使计算结果与实际有一定的偏差，在模型的最后给出了部分修正。
参考文献
[1] 郑友益. 对数螺线在机械设计与制造工程中的应用. 机械设计，
2003(1):69~70
[2] 金福. 关于对数螺线不变性之证明. 沈阳师范学院学报（自然科学
版）,1999(3):4~7
[3] 陈晓霞主编. ANSYS7.0高级分析. 北京：机械工业出版社，2004.
[4] 武良臣，郑友益. 型面联接制造理论与技术. 北京：煤炭工业出版社，2001.
[5] 濮良贵. 机械设计（第七版）. 北京：高等教育出版社，2001.
[6] 金福. 关于对数螺线不变性之证明. 沈阳师范学院学报.1999.( 3)
[7] 王忠. 对数螺线偏心夹紧机构在石材加工设备中的庙用[ J] .石材.2003(9).
[8] 程昆明. 对数螺线细粉分离器在新乡火电厂的应用. 节能, 2001(2).
[9] 韩旭东. 对数螺线叶片结构设计及应用. 煤矿机械, 1998(12).
[10] 王丕增. 修正对数螺线回零凸轮轮廓的设计. 哈尔滨工业大学学报, 1992(4).
附录
#include
using namespace std;
int main()
{
cout
for(int i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
}
数学建模网络挑战赛
承 诺 书
我们仔细阅读了第五届“认证杯”数学中国数学建模网络挑战赛的竞赛规则。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。
我们允许数学中国网站(www.madio.net ) 公布论文，以供网友之间学习交流，数学中国网站以非商业目的的论文交流不需要提前取得我们的同意。
我们的参赛队号为： 参赛队员 (签名) ：
队员2：
队员3：
参赛队教练员 (签名) ：
参赛队伍组别：
队员1：
数学建模网络挑战赛 编 号 专 用 页
参赛队伍的参赛队号：（请各个参赛队提前填写好）：
竞赛统一编号（由竞赛组委会送至评委团前编号）：
竞赛评阅编号（由竞赛评委团评阅前进行编号）：
2012年第五届“认证杯”数学中国
数学建模网络挑战赛
题 目 对数螺线型蜘蛛网状的结构分析
关 键 词 蜘蛛网 对数螺线 蒙特卡洛方法 ANSYS分析法
摘要
本文针对蜘蛛网合适结构的问题，考虑吐丝量一定，外界环境较理想条件下，建立以对数螺线为核心的数学模型，追求蜘蛛网结构最优。运用蒙特卡洛方法，模拟昆虫触网的过程，考虑了在蜘蛛丝长度一定的条件下，对数螺旋比圆围成的面积大，但疏而不漏，应用随机过程近似昆虫触网的过程，得出了对数螺线更利于捕食的结论。另一方面，也对对数螺线型面联接理论和联接界面强度进行了分析与计算，利用ANSYS 进行接触分析，得出了对数螺线型面联接的接触应力和接触强度条件的表达式。采用随机数产生算法，利用MATLAB 7.0.1和C++编程，分别对模型进行求解，并对所得结果进行分析比较，以此来帮助设计最有蜘蛛网结构。
参赛队号
所选题目 A
Abstract
Our article aims to study the question about the best structure of the spider webs ,it is on the condition of certain output of the spinning the and quite ideal conditions ,establish mathematical model in the core of the logarithmic spiral to find the best way of the spider webs .We also analyze Logarithm of solenoid type surface connection theory, Interface connection strength and ANSYS to get the expression. we apply Monte Carlo method to simulate the process about Net insert and adopt the Random number produce algorithm ,we also use the software of Matlab 7.0.1 、Mathematica and Microsoft Visual C++ 6.0 to give the answer to the question about the model and analyze about the result from model ,so we establish the best structure of the spider webs by means of these datas.
一、问题重述
世界上生存着许多种类的蜘蛛，而其中的大部分种类都会通过结网来进行捕食。通过对蜘蛛网所形成的结构的分析，通过建立模型，设计一种更为合适的蜘蛛网结构。
二、问题分析
题目中主要研究的是：蜘蛛网织成怎样的结构才是最合适的。因此我们通过查阅资料，了解蜘蛛网的结构等方面内容，根据结构形状的不同, 蛛网可以分为片网、不规则网和圆网等几种类型。由于圆网在蛛网进化上的地位特殊，且结构简单、规则。因此，到目前为止对蛛网的研究大都集中在圆网上。圆网并不是标准的圆，而是近似于数学上的螺旋线。所谓合适的蜘蛛网结构就是利于蜘蛛捕食、防御、繁殖。根据题目的要求，我们提出以下几个问题：
1、为什么蜘蛛网是螺旋线状，而不是标准的同心圆； 2、对数螺线型面联接理论和联接界面强度分析与计算。
三、符号说明
ρ：模拟对数螺线的极径
θ：模拟对数螺线的极角
Φ：对数螺线型型面轴旋转角
σp ：面轴上接触应力
l ： 轴孔之间的轴向配合长度
γ：P 点接触应力σp 与极径ρ之间的夹角 f ：轴孔之间的摩擦系数 T ： 扭矩
S 1：面积就是最大圆的面积 L 1：四圈的长度为
四、模型假设：
1、假设蜘蛛网是规则的对数螺线；
2、不考虑蜘蛛网受到风雨等天气情况的影响;
3、假设昆虫飞向蜘蛛网时，落在网内每点的概率相同；
五、模型建立与求解：
蜘蛛网的中心和圆周之间呈辐射状的半径线，自外向里是螺旋线，愈近中心，每圈间的距离也愈小，直到不可辨认的地步，这正符合数学上的对数螺线的情况。因此，我们建立对数螺线的模型，近似代替蜘蛛网，研究其性质。
图1
㈠ 对数螺线的定义和性质
数学上对数螺线定义如下：动点的运动方向始终与极径保持定角θ的动点轨迹，称为对数螺线。如图1所示，其极坐标方程为：
ρ=ae m θ （1）
式中：a , m 为常数（λ=arctan (1/m )））；θ为极角，ρ为极径。
图2 对数螺线
对数螺线在渐屈、渐伸、垂迹、回光线等各种变换下的不变性质, 体现出自身的高度和谐、对称和统一性。 ㈡对数螺旋线与圆形蜘蛛网的比较
将四个标准圆形与对数螺旋线放入同一坐标系中，如下图3
图3
1、 四圈圆形蜘蛛网
面积就是最大圆的面积：S 1=82π=64⨯3.14=201.056 四圈的长度为： L 1=2(1+3.6+5.8+8.2)π=111.784 2、四圈螺旋线蜘蛛网
该对数螺旋线的方程为：ρ=0.02e 3θ，0≤θ≤8π
面积：由于对数螺旋线是一条不封闭的曲线，所以用下图中最外面的曲线和一条线段组成的封闭图形表示该螺旋线所包围的面积。通过数该封闭图形内的方格数估计面积，不足一格按半格记。
图4
共有188个正方形，45个不足一格的，所以面积为198+22.5=220.5 长度：
用Mathematica 计算该曲线长，
输入 Integrate[0.02e^3x,{ x,0,8Pi}] ， 得出结果L= 126.871
结论：
在蜘蛛丝长度一定的情况下，螺旋线所围成蜘蛛网的面积大。这样更利于蜘蛛捕食。
蒙特卡罗方法：
用蒙特卡罗方法模拟昆虫飞向蜘蛛网上的过程：
假设昆虫飞向蜘蛛网时是一个随机过程，此过程中不考虑环境因素（风向、风速等）的影响
编写C 语言程序，生成二维随机数。 程序及运行结果见附录。
将这些随机数在下图5中描点，为处理简单，以第一象限为例，其他象限相同。
图5
图5.1 图5.2 图5.3 图5.4 图5.5
由此可见，昆虫更可能碰到对数螺线。也就是说，对数螺线形的蜘蛛网更有利于捕食
㈢、对数螺线型面联接理论和联接界面强度分析与计算 3.1对数螺线曲线
以型面轴截面曲线为例, 选取三段对数螺线进行分析。θ从0°到110°的一段曲线组成, 其中各段曲线之间用直线圆滑联接, 减弱了应力集中现象。
图6三段对数螺线型面轴截面 2.2对数螺线曲线参数的确定
在图1所示对数螺线的方程ρ=ae m θ中, m的大小取决于型面联接轴与孔之间受力时的压力角α, 因α在整个曲线上是常数, 当α选定时, m 为常数。因此选择型面联接轴、孔截面形状时, 可选定a 值的大小来定轴、孔的尺寸, 再定出压力角, 即可确定式中m 值, 对数螺线形状也就随之确定, 同时型面轴、孔截形也就相应得到确定。
三段对数螺线组成的型面轴、孔有两个基本参数a 和α, 其中m 主要影响曲线的形状, a 主要影响曲线的大小。三段对数螺线型面联接如图7所示, 这是一种有间隙的配合(图7-a), 工作时通过一定量的相对旋转, 间隙补偿, 由于楔面的作用在接触面之间便产生正压力, 并摩擦闭锁, 形成可靠的联接(见图7-b), 其摩擦受力方向与运动方向所成角度a 的大小与参变量θ无关, 即对数螺线的压力角a 在任何位置都是相同的。
( a ) ( b) 图7对数螺线型面联接
3.2 型面轴孔接触初始位置确定
设型面轴孔的截面廓形曲线方程分别为ρ1=a 1e m θ、ρ2=a 2e m θ。假定固定轮毂，顺时针旋转型面轴Φ角后两者之间有初始接触，任取一接触点P ，则在P 点处有ρ1=ρ2，
即 a 1e m (θ+Φ)=a 2e m θ
可得 Φ=（2） （3） 1⎛a 2⎫ln ⎪ m ⎝a 1⎭
上式表明型面轴旋转角Φ仅与曲线的a , m 常数值有关，而与θ值无关，即与初始接触点P 的位置无关，即Φ为定值，说明轴孔工作表面之间同时发生接触，轴孔之间接触为面接触。
3.3利用ANSYS 进行接触分析
由于ANSYS 对复杂曲面建模具有一定的局限性，为了分析的准确性，利用Pr o /E 强大的三维建模功能，在Pr o /E 中建成模型后，再利用Pr o /E ,ANSYS 之间的接口程序，将模型导入, ANSYS 中。为了方便加载扭矩，在不影响分析结果性质的前提下，于轴中心建一半径为r 的圆形孔，在圆形孔边界节点处加载等效切向力F ，使得扭矩T=Fr，如图8所示。按照ANSYS 分析步骤设置属性、划分网格、加载、求解。在后处理器POST1中查看有限元模型在纯扭矩T 作用下的节点应力云图，如图9所示。
图8有限元模型
图9节点的应力云图
图9表明在轴孔各段工作表面上除了两端（刚进入接触与刚脱离接触的很小的一个区域）应力比较大之外，其它的区域所受到的应力分布都非常均匀。若在模型应力比较集中区域进行修磨处理，将各段曲线间联接用更圆滑的曲线过渡，可减少或消除应力集中现象。由此，认为对数螺线型面联接接触应力均匀分布是符合实际情况的。反复加载不同的扭矩对有限元模型进行分析求解，发现作用扭矩与接触应力之间成线性比例关系。
3.4接触应力计算
假定在无间隙无过盈的理想配合状态下，轴上作用纯扭矩，如图5所示。在对数螺线型面轴上工作表面共有三段，其上接触应力均匀分布，记为σp 。在AB
表面上任取一点P ，在包含P 点的微小弧段ds 上，接触应力的推动分力（即与该点的线速度平行的分力）及摩擦力对轴线o 之矩dt 应为
dT =ρ(σp sin γlds +σp f cos γlds )
=(sin γ+f cos γ)σp ρlds （4） 式中:ρ为P 点的极径; l 为轴孔之间的轴向配合长度（即工作长度）；γ为P 点接触应力σp 与极径ρ之间的夹角；f 为轴孔之间的摩擦系数（一般取0.1~0.15）。
由对数螺线性质知γ=α=π/2-λ, ρ=ae m θ, ds =ρd θ, 代入式（4）得
dT =(sin γ+f cos α)ρ2σp ld θ
=(sin α+f cos α)ρ2σp le 2m θd θ (5) 在传动轴匀速转动时，由型面轴廓形曲线的对称性知，求出一段工作表面上力对轴线o , 产生的力矩，再乘以3所得的总力矩应与作用扭矩T 平衡，即有 T =3⎰
211π180dT =3⎰11π180a σp l (sin α+f cos α)e 22m θd θ（6） m π⎛119⎫3a σp l (sin α+f cos α) e -1⎪⎝⎭=2m
即σp =2mT
⎛3a 2l (sin α+f cos α) e ⎝11m π9⎫-1⎪⎭ （7）
上式表明作用扭矩T 与σp 之间成线性比例关系，当作用扭矩一定时，接触
应力的大小取决于轴孔之间的轴向配合长度l 以及a , α等参数，增大a , l 或α均可减小对数螺线型面联接的接触应力σp 。
由于在实际中蜘蛛网面接触应力并不均匀。对于型面联接，要求在轴孔配合表面存在一定间隙，所以前面无间隙无过盈的假定与实际情况有一定出入。由于间隙的存在，各段工作表面并不能实现在其理论区域上全部接触。这样按式（7）计算的σp (小于其实际值。为使计算的σp 更加接近最大接触应力，可引入修正系数加以修正。因此对数螺线型面联接的接触强度条件可写成
σp =2kmT
⎛3a 2l (sinα+f cos α) e ⎝11m π9⎫-1⎪⎭≤[σ]p （8）
六、模型结果分析
通过建立对数螺线型近似蜘蛛网结构，运用蒙特卡洛方法模拟猎物触网过程，我们给出了在蜘蛛丝长度一定的情况下，螺旋线所围成蜘蛛网的面积大，更有利于捕食。对对数螺线型面轴廓形曲线的分析，在初步确定对数螺线型面轴孔联接初始接触位置及接触性质的基础上，应用有限元分析软件ANSYS 对对数螺线型面联接进行接触应力分布规律的研究不失为一种行之有效的方法。本文分析求解的接触应力的理论计算公式，由于各种综合因素考虑的不足，从本质上看，其计算结果与实际情况会有一定出入。然而该公式作为一种有效的实用算法，结构简单，使用方便，为此种联接在工程中的应用打下了良好的设计基础。
七、模型优缺点分析
本文中针对研究蜘蛛网的合适结构给出了对数螺线型模型，本模型较好的近似了蜘蛛网的结构，对于对数螺线型网状结构的优越性进行了较详细的说明，并给出了对数螺线型面联接理论和联接界面强度分析与计算。较好地说明了对数螺线型蜘蛛网是较为适合的蛛网结构。但是由于在模型建立时的诸多理想假设、主观原因和参数敏感度的影响，致使计算结果与实际有一定的偏差，在模型的最后给出了部分修正。
参考文献
[1] 郑友益. 对数螺线在机械设计与制造工程中的应用. 机械设计，
2003(1):69~70
[2] 金福. 关于对数螺线不变性之证明. 沈阳师范学院学报（自然科学
版）,1999(3):4~7
[3] 陈晓霞主编. ANSYS7.0高级分析. 北京：机械工业出版社，2004.
[4] 武良臣，郑友益. 型面联接制造理论与技术. 北京：煤炭工业出版社，2001.
[5] 濮良贵. 机械设计（第七版）. 北京：高等教育出版社，2001.
[6] 金福. 关于对数螺线不变性之证明. 沈阳师范学院学报.1999.( 3)
[7] 王忠. 对数螺线偏心夹紧机构在石材加工设备中的庙用[ J] .石材.2003(9).
[8] 程昆明. 对数螺线细粉分离器在新乡火电厂的应用. 节能, 2001(2).
[9] 韩旭东. 对数螺线叶片结构设计及应用. 煤矿机械, 1998(12).
[10] 王丕增. 修正对数螺线回零凸轮轮廓的设计. 哈尔滨工业大学学报, 1992(4).
附录
#include
using namespace std;
int main()
{
cout
for(int i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
cout
for(i=0;i
{
cout
}