基于ANSYS的四杆机构分析详细
直列四缸发动机曲轴ansys分析
直列四缸发动机曲轴a n s y s分析本页仅作为文档封面，使用时可以删除This document is for reference only-rar21year.March有限元分析课程报告直列四缸发动机曲轴有限元分析姓名：学号：分数：年月日目录1.引言 .......................................................................................................错误!未定义书签。
2.曲轴有限元模型的建立........................................................................错误!未定义书签。
3.曲轴网格划分 .......................................................................................错误!未定义书签。
确定物理场和网格划分法.............................................................错误!未定义书签。
确定全局网格参数设置.................................................................错误!未定义书签。
确定局部网格参数设置.................................................................错误!未定义书签。
网格质量检查 ................................................................................错误!未定义书签。
调整网格划分 ................................................................................错误!未定义书签。
四杆机构实验报告
四杆机构实验报告四杆机构实验报告引言四杆机构是一种常见的机械结构，由四个连杆组成，通过铰链相连。
在工程学中，四杆机构被广泛应用于各种机械设备中，如发动机、机械手臂和运动机构等。
本实验旨在通过对四杆机构的研究和实验，深入了解其运动特性和应用。
一、实验目的本实验的主要目的是通过实验观察和测量四杆机构的运动规律，深入了解四杆机构的结构和工作原理。
具体目标包括：1. 理解四杆机构的基本构造和运动方式；2. 掌握四杆机构的运动规律和特性；3. 分析四杆机构的应用领域和优缺点。
二、实验装置和方法1. 实验装置：本实验使用由四个连杆和铰链组成的四杆机构模型，以及相应的测量仪器和工具。
2. 实验方法：a. 调整四个连杆的长度和角度，使其满足四杆机构的运动要求；b. 使用测量仪器测量和记录四个连杆的长度和角度的变化；c. 分析测量数据，绘制运动曲线和角度变化图表；d. 观察和记录四杆机构的运动规律和特点。
三、实验过程和结果1. 实验过程a. 调整四杆机构的连杆长度和角度，使其形成一个闭合的四边形；b. 使用测量工具测量和记录四个连杆的长度和角度的变化；c. 通过改变连杆长度和角度，观察和记录四杆机构的运动规律；d. 根据测量数据和观察结果，分析四杆机构的运动特性和应用。
2. 实验结果a. 通过测量仪器记录的数据，绘制了四个连杆的长度和角度的变化图表；b. 观察和记录了四杆机构的运动规律，如连杆的旋转和平移运动；c. 分析了四杆机构的运动特性和应用，如转动和传动功能。
四、实验分析和讨论1. 实验分析通过实验观察和测量，我们可以清楚地看到四杆机构的运动规律和特性。
连杆的长度和角度的变化直接影响着四杆机构的运动方式和速度。
通过改变连杆的长度和角度，我们可以实现不同的运动效果，如旋转、平移和传动等。
2. 实验讨论四杆机构具有广泛的应用领域和优点。
它可以用于实现复杂的运动和动力传递，如发动机和机械手臂。
同时，四杆机构的结构简单，易于制造和维护。
四连杆机运动学分析报告
栏杆机四杆机构运动学分析1 四杆机构运动学分析1.1 机构运动分析的任务、目的和方法曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。
对四杆机构进行运动分析的意义是：在机构尺寸参数已知的情况下，假定主动件（曲柄）做匀速转动，撇开力的作用，仅从运动几何关系上分析从动件（连杆、摇杆）的角位移、角速度、角加速度等运动参数的变化情况。
还可以根据机构闭环矢量方程计算从动件的位移偏差。
上述这些内容，无论是设计新的机械，还是为了了解现有机械的运动性能，都是十分必要的，而且它还是研究机械运动性能和动力性能提供必要的依据。
机构运动分析的方法很多，主要有图解法和解析法。
当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。
而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来，以便于机构的优化设计。
1.2 机构的工作原理在平面四杆机构中，其具有曲柄的条件为：a.各杆的长度应满足杆长条件，即：最短杆长度+最长杆长度≤其余两杆长度之和。
b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时，四杆机构为曲柄摇杆机构；当最短杆为机架时，则为双曲柄机构）。
三台设备测绘数据分别如下：第一组(2代一套)四杆机构L1=125.36mm，L2=73.4mm， L3=103.4mm，L4=103.52mm最短杆长度+最长杆长度(125.36+73.4) <其余两杆长度之和(103.4+103.52)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-1 II-1型栏杆机机构测绘及其运动位置图第二组(2代二套)四杆机构L1=125.36mm，L2=50.1mm，L3=109.8mm，L4=72.85mm最短杆长度+最长杆长度(125.36+50.1) <其余两杆长度之和(109.8+72.85)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-2 II-2型栏杆机机构测绘及其运动位置图第三组(3代)四杆机构L1=163.2mm，L2=64.25mm，L3=150mm，L4=90.1mm最短杆长度+最长杆长度(163.2+64.25) <其余两杆长度之和(150+90.1)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-3 III型栏杆机机构测绘及其运动位置图在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。
基于ANSYS的平面四杆机构仿真的结构优化设计
基于ANSYS的平面四杆机构仿真的结构优化设计作者：焦晨航来源：《硅谷》2014年第24期摘要随着计算机技术的迅猛发展，结构优化设计方法也随之变更。
本文简要介绍了有限元分析的典型步骤和流程，并运用ANSYS软件，以曲柄摇杆机构为例，对平面四杆机构仿真的结构优化设计进行探讨，以期为进一步的动力学分析打下基础。
关键词 ANSYS；平面四杆机构；曲柄摇杆机构；优化设计中图分类号：TH112 文献标识码：A 文章编号：1671-7597（2014）24-0005-01平面四杆机构是连杆机构中最常见的机构组成，广泛应用于工程机械和农业机械中，对四杆机构仿真的结构优化设计进行研究具有重要意义。
平面四杆机构仿真的结构优化设计就是根据实际情况在计算机编程的帮助下建立起机构的数学模型，并通过运用一定的优化算法寻找既能很好解决约束条件又能使目标函数最优的设计方案，最终达到优化设计的目的。
有限元方法是结构优化设计的重要方法之一，实际应用中，我们可通过自行编制有限元程序或采用通用的有限元分析软件来进行。
基于ANSYS的平面四杆机构仿真的结构优化设计时有效的、可行的。
1 有限元分析的典型步骤和流程有限元法是一种高效能、常用的、离散化的数值分析方法，其典型步骤为：①将连续体离散成有限个单元（杆系的单元是每一个杆件，连续体的单元是各种形状）；②选择单元类型或位移模型；③利用变分原理推导单元刚度矩阵；④单元场函数的集合；⑤建立有限个待定参量的代数方程组；⑥求解方程组，得到位移矢量；⑦由节点位移计算出单元的应变和应力。
进行有限元分析，首先要决定分析项目、分析的几何结构、外界条件和外力，获取材料性质；其次，建立有限元模型，包括单元类型、材料性质，直接或间接生成有限元网格；再次，加载并求解，输出分析结果，若结果不合理，重新建立有限元模型，若合理，则进行改进处理，解决问题，得到最佳设计。
2 基于ANSYS平面四杆机构仿真的结构优化设计1）结构优化设计方法。
【ANSYS算例】3.2.5(4)-四杆桁架结构的有限元分析(GUI)及命令流
四杆桁架结构的有限元分析下面针对【典型例题】(1)的问题，在ANSYS 平台上，完成相应的力学分析。
即如图3-8所示的结构，各杆的弹性模量和横截面积都为4229.510N/mm E,E=29.5X10 2100mm A ，基于ANSYS 平台，求解该结构的节点位移、单元应力以及支反力。
图3-8 四杆桁架结构解答 对该问题进行有限元分析的过程如下。
以下为基于ANSYS 图形界面( graphic user interface ，GUI)的菜单操作流程；注意：符号“→”表示针对菜单中选项的鼠标点击操作。
关于ANSYS 的操作方式见附录B 。
1． 基于图形界面的交互式操作(step by step)(1) 进入ANSYS(设定工作目录和工作文件)程序 →ANSYS → ANSYS Interactive →Working directory (设置工作目录) →Initial jobname (设置工作文件名): planetruss →Run → OK(2) 设置计算类型ANSYS Main Menu : Preferences… → Structural → OK(3) 选择单元类型ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete… →Add… →Link ：2D spar 1 →OK (返回到Element Types 窗口) →Close(4) 定义材料参数ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic → Isotropic: EX:2.95e11 (弹性模量)，PRXY: 0 (泊松比) → OK → 鼠标点击该窗口右上角的“ ”来关闭该窗口(5) 定义实常数以确定单元的截面积ANSYS Main Menu: Preprocessor →Real Constant s… →Add/Edit/Delete →Add →Type 1→ OK →Real Constant Set No: 1 (第1号实常数), AREA: 1e-4 (单元的截面积) →OK →Close(6) 生成单元 ANSYS Main Menu: Preprocessor →Modeling →Creat →Nodes → In Active CS →Node number 1 → X:0,Y:0,Z:0 →Apply →Node number 2 → X:0.4,Y:0,Z:0 →Apply →Node number 3 → X:0.4,Y:0.3,Z:0→Apply →Node number 4 → X:0,Y:0.3,Z:0→OKANSYS Main Menu: Preprocessor → Modeling → Create → Elements →Elem Attributes （接受默认值）→Usernumbered→Thru nodes→OK→选择节点1，2→Apply→选择节点2，3→Apply→选择节点1，3→Apply→选择节点3，4→Apply→OK(7)模型施加约束和外载添加位移的约束，分别将节点1 X和Y方向、节点2 Y方向、节点4的X和Y方向位移约束。
机械原理课件之四杆机构受力分析
通过解方程，求解出各个连杆的受力大小和方向。
四杆机构受力分析的案例研究
案例1
案例2
分析一台工业机械中的四杆机构， 确定各个连杆的受力情况。
在一个机器人手臂中应用四杆机 构，研究其受力和应力分析。
案例3
通过受力分析，优化四杆机构的 设计，提高其工作效率。
结论和总结
四杆机构受力分析是机械工程领域的重要研究方向之一。它不仅可以帮助我 们了解四杆机构的工作原理，还可以指导我们设计更优秀的机械系统。
四杆机构的组成和基本结构
连杆
四杆机构由四根连杆组成，包括两个边连杆和两个角连杆。
铰链
连杆通过铰链连接，使得四杆机构能够实现运动。
驱动装置
驱动装置为四杆机构提供动力，使其能够完成特定任务。
四杆机构的运动分析
1
自由度
四杆机构的自由度取决于连杆的个数和铰链的类型。
2
运动类型
四杆机构可以实现旋转、平动和复杂的运动。
3
工作轨迹
通过对四杆机构的运动分析，可以得到工作轨迹的方程。
四杆机构受力分析的基本原理
四杆机构受力分析的基本原理是根据静力学的原理，通过分析力的平衡条件 来确定各个连杆的受力情况。
四杆机构受力分析的方法和步骤1 建立坐标系确定来自适的坐标系，便于受力分析的计算。
2 列写平衡方程
根据力的平衡条件，列写各个连杆的受力方程。
机械原理课件之四杆机构 受力分析
这篇课件将详细介绍四杆机构的受力分析。从概述四杆机构的基本原理开始， 到运动分析和受力分析的具体方法，最后通过案例研究加深理解。让我们一 起来探索吧！
四杆机构的概述
四杆机构是一种常见的机械连杆机构，由四根连杆组成。它具有简单的结构 和广泛的应用领域，是研究机械原理的重要组成部分。
机械毕业设计(论文)基于ansys的连杆机构的有限元分析【全套设计】
湘潭大学兴湘学院毕业设计论文题目：连杆机构的有限元分析全套设计，加153893706专业：机械设计制造及其自动化学号： 2010963028 姓名：指导教师：完成日期： 2014 年 5 月 25 日湘潭大学兴湘学院毕业论文（设计）任务书论文（设计）题目：连杆机构的有限元分析学号： 2010963028姓名：专业：机械设计制造及其自动化指导教师：系主任：一、主要内容及基本要求1、总结连杆机构设计方法研究和连杆机构研究的发展状况和发展趋势，在总结前人研究成果的基础上，结合当前的技术发展趋势，采用有限元方法来进行开展研究。
2、阐述学习理论基础，即瞬态动力学分析,简要论述瞬态参数,识别原理。
3、简要论述有限元方法和动力学分析的基本求解过程，建立连杆机构中的曲柄滑块机构的有限元模型，合理的确定曲柄长度及转速、连杆长度和转速，偏距，选定和创建单元类型，指点单元属性，创建铰链单元，采用瞬态动力学分析瞬态分析类型对其进行瞬态分析，与图解法进行比较，验证有限元瞬态求解功能。
4、联系工程实际，对受力连杆进行结构静力学学习。
二、重点研究的问题1、 ANSYS的线性静力分析2 、构建几何模型3、在三维铰链单元COMBIN7的创建4、单元类型选择和网络划分5、 ANSYS瞬态动力学分析和静力学分析三、进度安排四、应收集的资料及主要参考文献[1]高耀东，刘学杰.ANSYS机械工程应用精华50例（第三版）.- 北京：电子工业出版社，2011.[2]孙波.毕业设计宝典.-西安：西安电子科技大学出版社，2008.[3]温正，张文电.ANSYS14.0有限元分析权威指南.-北京：机械工业出版社，2013.[4]欧阳周，汪振华，刘道德.毕业论文和毕业设计说明书写作指南.-长沙：中南工业大学出版社，1996.[5]华大年，华志宏.连杆机构设计与应用创新.-北京：机械工业出版社，2008.[6]胡仁喜，康士廷.机械与结构有限元分析从入门到精通.-北京：机械工业出版社，2012.[7]李红云，赵社戌，孙雁.ANSYS10.0基础及工程应用.北京：机械工业出版社，2008.[8]唐家玮，马喜川.平面连杆机构运动综合.-哈尔滨：哈尔滨工业大学出版社，1995.[9]潘存云，唐进元.机械原理.-长沙：中南大学出版社，2011.[10]李皓月，周田朋，刘相新.ANSYS工程计算应用教程.-北京:中国铁道出版社，2003湘潭大学兴湘学院毕业论文（设计）评阅表学号2010963028 姓名谭磁安专机械设计制造及其自动化毕业论文（设计）题目：连杆机构的有限元分析湘潭大学兴湘学院毕业论文（设计）鉴定意见学号2010963028 姓名谭磁安专业机械设计制造及其自动化毕业论文77 页图表30 张目录摘要............................................................................................ 错误！未定义书签。
ANSYS分析
机械臂是面向工业领域的多关节机械手，是自动执行工作的机器装置，是靠自身动力和控制能力来实现各种功能的一种机器。
它可以接受人类指挥，也可以按照预先编排的程序运行，现代的机械臂还可以根据人工智能技术制定的原则纲领行动。
机械臂由主体、驱动系统和控制系统三个基本部分组成。
主体即机座和执行机构，包括臂部、腕部和手部。
大多数机械臂有3～6个运动自由度，其中腕部通常有1～3个运动自由度；驱动系统包括动力装置和传动机构，用以使执行机构产生相应的动作；控制系统是按照输入的程序对驱动系统和执行机构发出指令信号，并进行控制。
图1 六自由度机械臂如图1所示为六自由度机械臂。
而机械臂在工作中需要承受一定的载荷，这会引起杆件的弹性变形，从而导致机械臂工作时产生一定的误差。
为了保证机械臂在运动中的定位误差，机械臂杆件结构需要有较高的刚性。
下面就机械臂杆件刚性的提高，从材料的选择、结构设计等方面进行有限元分析二、模型建立并导入到ansys使用solidworks进行机械臂三维模型的建立。
根据实验室相关尺寸建立模型如下图2、3。
其中图2形象的展现了其三维外观，而图3的前视图方便说明起尺寸大小。
图2 机械臂的三维模型图图3 机械臂的三维模型前视图由于是为了研究机械臂末端的应力、应变，所以为了在ansys中分析方便，在保留主体结构设计的前提下，可将该机械臂结构进行简化，得到如图4的第一种结构模型图。
且图4中为简化后结构的主体，其后在对不同结构设计的讨论中，需对该模型进行修改。
其中模型的体积为V A=18528979.98mm3=0.01853m3。
进行简化，得到如图4的第一种结构模型图。
且图4中为简化后结构的主体，其后在对不同结构设计的讨论中，需对该模型进行修改。
其中模型的体积为V A=18528979.98mm3=0.01853m3。
图4 简化后的结构模型1的三维图和基本尺寸将图4所示的模型导入到ansys中：(1)先在solidworks中把模型另存为Parasolid(*.x_t)格式（注意模型名字必须是英文，Ansys不接受中外，可保存为jxb.x_t）；(2)然后打开Ansys，在File下选择import-PARA，找到之前保存的Parasolid(*.x_t)格式的模型，将其导入；(3)在顶上菜单栏对话框中选择PoltCtrls—Style—Solid Model Facets，如图6所示。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ANSYS作业一．问题描述图1所示平行四边形机构，曲柄长200mm，连杆长l根据各自学号后2位乘以10，各杆截面为40x5，其中宽度为40，厚度为5。
现在连杆上表面加1MPa的三角分布的压力（铰接结构自习设计，加载面去除2端铰接结构），求各杆的强度和变形。
要求按报告格式写，写出主要步骤、注意事项、关键程序、结果及其评价（材料按Q235A），铰接处结构是否合理。
二．有限元分析本文采用ANSYS编程语言APDL，编制参数化程序。
简介方便，便于重复分析，节省大量的工作量。
1.定义材料、单元finish/clear/prep7 ！进入前处理器et,1,solid185 ！定义8节点实体单元solid185 mp,ex,1,2.08e5 ! 弹性模量mp,prxy,1,0.277 ! 泊松比mp,dens,1,7.86e-6 ! 密度et,2,conta173 ！定义接触单元conta173et,3,targe170 ！定义目标单元targe170 keyopt,2,5,3 ！设置接触单元选项keyopt,2,7,1keyopt,2,9,0keyopt,2,12,02.四杆机构建模参数p0=1 ！载荷a=40 ！连杆截面宽度b=5 ！连杆截面厚度l1=200 ！曲柄长度l2=90 ！连杆长度3.四杆机构参数化建模block,0,b,8,l1-8,0,a ！生成长方体块block,0,b,0,8,0,10block,0,b,0,8,a-10,ablock,0,b,l1-8,l1,0,10block,0,b,l1-8,l1,a-10,avadd,all ！将以上长方体布尔相加得到曲柄block,8,l2-8,0,b,0,a ！生成长方体块block,0,8,0,b,10,30block,l2-8,l2,0,b,10,30vadd,1,2,3 ！将以上长方体布尔相加得到机架vgen,2,6,6,0,l2-5 ！复制曲柄得到第四杆vgen,2,4,4,0,,l1-5 ！复制机架得到连杆wpoffs,2.5,2.5,0 ！工作平面沿x、y正向各平移2.5mm cylind,0,2,0,40,0,360 ！在铰接处画圆柱半径2mm、长度40mm vsbv,6,3,sepo ！布尔运算、曲柄减去圆柱生成曲柄铰链孔cylind,0,2,0,40,0,360 ！生成曲柄铰链的销vsbv,4,3,sepo ！布尔运算、机架减去圆柱生成机架铰链孔cylind,0,2,0,40,0,360vgen,2,3,3,0,l2-5 ！将圆柱销复制得到另外三个铰接出圆柱vsbv,6,4,sepo ！布尔运算得到铰链孔vgen,2,3,3,0,l2-5 ！生成圆柱vsbv,1,4,sepo ！布尔运算得到铰链孔vgen,2,3,3,0,l2-5 ！生成销vgen,2,3,3,0,,l1-5vsbv,5,4,sepo ！布尔运算得到铰链孔vgen,2,3,3,0,,l1-5 ！生成圆柱vsbv,2,4,sepo ！布尔运算得到铰链孔vgen,2,3,3,0,,l1-5 ！生成销vgen,2,3,3,0,l2-5,l1-5vsbv,5,4,sepo ！布尔运算得到铰链孔 vgen,2,3,3,0,l2-5,l1-5 ！生成圆柱vsbv,6,4,sepo ！布尔运算得到铰链孔 vgen,2,3,3,0,l2-5,l1-5 ！生成销四杆机构实体模型如图1所示：图1（a ） 图1（b ） 图1（c ）图1. 四杆机构、铰链处销实体模型4. 有限元网格划分wpoffs,-2.5,-2.5,0 ！工作平面移动到总体坐标系原点处 wpoffs,8,195,0 ！移动工作平面wprota,,,90 ！工作平面绕y 轴旋转90°vsbw,all ！用工作平面切割连杆与机架一端 wpoffs,,,l2-16 ！移动工作平面vsbw,all ！用工作平面切割连杆与机架另一端 wpoffs,,,-(l2-16) ！移动工作平面 wpoffs,,-3,0 ！移动工作平面wprota,,90 ！工作平面绕x 轴旋转90° vsbw,all ！用工作平面切割体 wpoffs,,,184 ！移动工作平面vsbw,all ！用工作平面切割体 esize,1,0, ！设置网格大小 mshape,1,3D ！设置单元形状 mshkey,0 ！网格划分方式vsweep,all ！扫略生成网格四杆机构有限元模型如图2所示：图2（a ） 图2（b ）图2（c）图2. 四杆机构有限元网格模型5.添加接触对asel,s,area,,3,4,1 ！选择第一个铰接处接触面（凸面）nsla,s,1 ！选择接触面上所有节点cm,jiechu1,node ！做成名为jiechu1的节点集合allsel,all ！选择所有asel,s,area,,69,70,1 ！选择第一个铰接处目标面（凹面）asel,a,area,,135,136,1asel,a,area,,21,23,2nsla,s,1 ！选择目标面上所有节点cm,mubiao1,node ！做成名为mubiao1的节点集合allsel,allasel,s,area,,11,12,1 ！选择第二个铰接处接触面（凸面）nsla,s,1 ！选择接触面上所有节点cm,jiechu2,node ！做成名为jiechu2的节点集合allsel,allasel,s,area,,18,26,8 ！选择第二个铰接处目标面（凹面）asel,a,area,,51,52,1asel,a,area,,77,78,1nsla,s,1 ！选择目标面上所有节点cm,mubiao2,node ！做成名为mubiao2的节点集合allsel,allasel,s,area,,20,22,2 ！选择第三个铰接处接触面（凸面）nsla,s,1 ！选择接触面上所有节点cm,jiechu3,node ！做成名为jiechu3的节点集合allsel,allasel,s,area,,111,112,1 ！选择第三个铰接处目标面asel,a,area,,121,122,1asel,a,area,,75,76,1nsla,s,1 ！选择目标面上所有节点cm,mubiao3,node ！做成名为mubiao3的节点集合allsel,allasel,s,area,,7,8,1 ！选择第四个铰接处接触面nsla,s,1 ！选择接触面上所有节点cm,jiechu4,node ！做成名为jiechu4的节点集合allsel,allasel,s,area,,97,100,1 ！选择第四个铰接处目标面asel,a,area,,14,16,2nsla,s,1 ！选择目标面上所有节点cm,mubiao4,node ！做成名为mubiao4的节点集合allsel,allr,1,,,1,0.1 ！定义接触属性：刚度渗透量等r,2,,,1,0.1r,3,,,1,0.1r,4,,,1,0.1type,2 ！生成第一个铰接处的接触对mat,1real,1cmsel,s,jiechu1,nodeesurf,toptype,3mat,1real,1cmsel,s,mubiao1,nodeesurf,toptype,2 ！生成第二个铰接处的接触对mat,1real,2cmsel,s,jiechu2,nodeesurf,toptype,3mat,1real,2cmsel,s,mubiao2,nodeesurf,toptype,2 ！生成第三个铰接处的接触对mat,1real,3cmsel,s,jiechu3,nodeesurf,toptype,3mat,1real,3cmsel,s,mubiao3,nodeesurf,toptype,2 ！生成第四个铰接处的接触对mat,1real,4cmsel,s,jiechu4,nodeesurf,toptype,3mat,1real,4cmsel,s,mubiao4,nodeesurf,top接触对模型如图3所示：图3（a）接触单元图3（b ）目标单元图3（c ）接触单元放大 图3（d ）目标单元放大图3. 四个铰链处接触对6. 加载与求解nsel,s,loc,x,5,l2-5 ！选择机架上所有节点，施加全约束 nsel,r,loc,y,0,5 d,all,allnsel,s,loc,x,8,l2/2,1！按坐标选择连杆上表面从左端点至中间位置的所有节点（铰接处除外） nsel,r,loc,y,199.9,200,0.1*get,nmax,node,,num,max, ！提取当前激活的最大节点数目编号nmax *get,nmin,node,,num,min, ！提取当前激活的最小节点数目编号nmin *dim,t1,array,nmax,1,1, ！定义名为t1的数组*do,j,nmin,nmax ！循环，j 从nmin 到nmax *if,nsel(j),eq,1,then ！if 判断 t1(j)=p0*abs(nx(j)-8)*1/l2/2！连杆从左端点到中间位置的载荷位置函数（一次函数）*else t1(j)=0*endif*enddo ！结束循环sffun,pres,t1(1)sf,all,pres,0 ！添加连杆上表面从左端点至中间位置的载荷local,12,0,l2,200,0 ！定义局部坐标系编号为12，远点位于连杆右端点上表面处csys,12allsel,allnsel,s,loc,x,(-1)*l2/2,-8,1nsel,r,loc,y,-0.1,0,0.1！按坐标选择连杆上表面从右端点至中间位置的所有节点（铰接处除外）*get,nmax,node,,num,max,*get,nmin,node,,num,min,*dim,t2,array,nmax,1,1,*do,i,nmin,nmax*if,nsel(i),eq,1,thent2(i)=p0*(abs(nx(i)+8)/l2/2)！连杆从右端点到中间位置的载荷位置函数（一次函数）*elset2(i)=0*endif*enddosffun,pres,t2(1)sf,all,pres,0allsel,allnummrg,node ！节点融合/solu ！进入求解器antype,static ！静态分析autots,on ! 使用自动时间步长neqit,200 ! 最大迭代次数200pred ! 跨越荷载步时不作预测nropt,full,,off ! 完全牛顿拉夫逊法，不使用自适应下降因子LNSRCH,onNLGEOM,on ! 考虑集合非线性EQSLV,PCG ！采用预条件共轭梯度迭代方程求解器nsubst,20,100,20 ！载荷步的子步数为20allsel,alloutpr,basic,all ！输出选项solve ！求解四杆机构加载以及约束如图4、5所示：图4. 机架全约束图5. 连杆上表面渐变的三角形分布载荷三．查看结果1.查看应力，应力图如图6、7所示。