ansys耦合和约束方程教程.ppt

第16章耦合和约束方程,正如自由度约束能约束模型中确定的节点一样，耦合和约束方程可以建立节点间的位移关系。本章主要讨论何时需要建立、怎样建立节点间的耦合和约束方程。主要内容:A.耦合B.约束方程,16.1耦合,耦合是使一组节点具有相同的自由度值除了自由度值是由求解器计算而非用户指定外，与约束相类似例如:如果节点1和节点2在UX方向上耦合，求解器将计算节点1的UX值并简单地把该值赋值给节点2的UX一个耦合设置是一组被约束在一起，有着同一方向的节点(即一个自由度)一个模型中可以定义多个耦合，但一个耦合中只能包含一个方向的自由度,16.1.1耦合设置的特点,只有一个自由度卷标如：ux,uy或temp可含有任意节点数任意实际的自由度方向ux在不同的节点上可能是不同的主、从自由度的概念加在主自由度上的载荷,16.1.2一般应用,施加对称条件无摩擦界面铰接,如：用耦合施加循环对称性，在循环对称切面上的对应位置实施自由度耦合,16.1.3施加对称条件,耦合自由度常被用来施加移动或循环对称性条件。这可以保证平面截面依然是平面。例如：-对圆盘扇区模型(循环对称)，应使两个对称边界上的对应节点在各个自由度上耦合。对锯齿形模型的半齿模型(平移对称)，应使一个边上的节点在各自由度上耦合,16.1.3施加对称条件(续),考虑在均匀轴向压力下的空心长圆柱体，此3D结构可用下面右图所示的2D轴对称模型表示,16.1.4无摩擦界面,如果满足下列条件，则可用耦合自由度模拟接触面表面保持接触几何线性分析(小变形)忽略摩擦在两个面上，节点是一一对应的通过耦合垂直于接触面的重合节点来模拟接触面分析仍然是线性的无间隙收敛性问题,在竖向耦合每对节点,16.1.5铰接,耦合可用来模拟铰接，如：万向节、铰链借助力矩释放可模拟铰接：只耦合连接节点间的位移自由度，不耦合旋转自由度例如：下图中，若A处重合两节点在UX、UY方向上耦合，旋转不耦合，则A连接可模拟成铰接,节点1和节点2重合，为了看清分开显示,16.1.6创建耦合设置,根据使用不同，可用多种方法进行耦合设置将节点进行同方向耦合：选择所需要的设置接着使用CP命令或orPreprocessorCoupling/CeqnCoupleDOFs例如,cp,ux,all是把所有选择节点在UX方向上耦合,输入耦合设置参考号，选择自由度卷标,16.1.6创建耦合设置(续),在零偏移量的一组节点之间生成附加耦合关系:MainMenu:PreprocessorCoupling/CeqnGenw/SameNodes,3.单击OK,1.输入现存耦合设置的参考号,2.对每个设置指定新的自由度卷标,16.1.6创建耦合设置(续),同一位置节点间的耦合：首先确保所有要耦合的节点都被选择接着使用命令CPINTF或PreprocessorCoupling/CeqnCoincidentNodes例如：cpintf,uy同一位置的所有节点在UY上耦合(包含0.0001的缺省误差),16.1.6创建耦合设置(续),不在同一位置节点间的耦合，如循环对称:首先确保所有要耦合的节点都被选择。然后使用命令CPCYC或PreprocessorCoupling/CeqnOffsetNodes例如：cpcyc,all,1,0,30,0把圆心角相差30的对应节点的各自由度进行耦合(注：当前KCN选项是总体柱坐标系),16.1.7关于耦合的说明,记忆要点:耦合中的自由度方向(UX,UY,等)是节点坐标系中的方向求解器只保留耦合中的第一个自由度，并把它作为主自由度，而不保留其余自由度施加在耦合节点上的载荷(在耦合自由度方向)求和后作用在主节点上耦合自由度上的约束只能施加在主节点上,16.1.8练习：耦合循环对称边界,在此练习中，由生成耦合DOF设置来模拟有循环对称性的模型的接触问题1.建模并在图形窗口中画单元2.在总体柱坐标系下，生成具有Y的增量为30的节点复制件a.将当前坐标系变为总体柱坐标系b.在当前坐标系中，以Y=30的增量拷贝所有的结点,16.1.8练习：耦合循环对称边界(续),3.在同一位置的节点上生成适当的耦合关系a.Choosecouplecoincidentnodesb.ChooseAllAppropriate4.不选择附在单元上的节点a.选择entity,nodeattachedtob.选择unselect，并单击apply5.将新节点拷贝回原始位置(DY=-30,INC=0)a.以Y30的增量拷贝所有节点b.对节点号增量输入06.选择everything,16.1.8练习：耦合循环对称边界(续),7.对所有处于同一位置的节点进行merge操作a.NumberingcontrolsMergeitemsb.关掉警告信息8.将所有的节点坐标系转到总体柱坐标系a.MainMenu:Preprocessor-Modeling-Move/ModifyRotatenodeCStoactiveCSb.拾取all9.求解并进行后处理,16.2约束方程,约束方程定义了节点自由度间的线性关系若两个自由度耦合,它们的简单关系是UX1=UX2约束方程是耦合的更一般形式，允许写诸如UX1+3.5*UX2=10.0的约束方程在一个模型中可以定义任意多个约束方程另外，一个约束方程可以包含任意数量的节点和自由度的集合。约束方程的一般形式是:Coef1*DOF1+Coef2*DOF2+Coef3*DOF3+.=Constant,16.2.1约束方程的特点,约束方程的特点自由度卷标的任意组合任意节点号任意实际的自由度方向在不同的节点上ux可能不同,16.2.2一般应用,连接不同的网格实体与实体的界面2D或3D相同或相似的单元类型单元面在同一表面上，但结点位置不重合连接不同类型的单元壳与实体垂直于壳或实体的梁建立刚性区过盈装配,16.2.3连接不同的网格,两个已划分网格的实体部分在某个面相连接，若它们的节点不相同，可以通过建立约束方程来建立连接处理此类情况最容易的方法是使用CEINTF命令(PreprocessorCoupling/CeqnAdjacentRegions),首先选择网格划分较细的对象的节点和另一方的单元自动计算所有必要的系数和常数适合于实体单元对实体单元，2-D或3-D,16.2.4连接不同类型的单元,如果需要连接自由度集不同的单元类型，则要求写出约束方程以便于从一类单元向另一类单元传递载荷:梁与实体或垂直于壳的梁壳与实体命令：CE命令(PreprocessorCoupling/CeqnConstraintEqn),建立转动自由度和移动自由度之间的关系,16.2.5建立刚性区,约束方程通常被用来模拟刚性区作用在节点(主节点)上的载荷将被恰当地分配到刚性区的其它节点上使用CERIG命令(或PreprocessorCoupling/CeqnRigidRegion),在某些特殊情况下，全刚性区给出了约束方程的另一种应用全刚性区和部分刚性区的约束方程都可由程序自动生成,16.2.6过盈装配,同接触耦合相类似，但在两界面间允许有过盈量或间隙典型方程:0.01=UX(node51)-UX(node251),16.2.7建立约束方程的过程,人工建立约束方程的菜单路径：MainMenu:PreprocessorCoupling/CeqnConstraintEqn,2.单击OK,1.输入常数项、节点号、自由度卷标和方程系数,16.2.7建立约束方程的过程(续),以现有的约束方程为基础生成约束方程：1.生成第一个约束方程:MainMenu:PreprocessorCoupling/CeqnConstraintEqn2.生成其余的约束方程:MainMenu:PreprocessorCoupling/CeqnGenw/SameDOF,生成的约束方程数.,现存约束方程中的节点增量,3.选择OK,生成的约束方程的起始序号，终止序号和增量,16.2.7建立约束方程的过程(续),通过“刚性区”来建立约束方程：MainMenu:PreprocessorCoupling/CeqnRigidRegion拾取将要连在一起的结点，然后单击OK,1.选择将要使用的刚性区的类型(自由度设置),2.单击OK,16.2.7建立约束方程的过程(续),在相邻的区域生成约束方程：1.从网格较密的区域中选择节点2.从网格较稀的区域中选择单元MainMenu:PreprocessorCoupling/CeqnAdjacentRegions,指定容差，此容差作为单元区域中最小单元长度的比率在约束方程中将要使用的自由度,5.单击OK,16.2.8练习：在蜗轮叶片上建立约束方程,在此练习中，将使用约束方程将具有不同单元类型和不同网格的两部分连接起来。这两部分分别是涡轮叶片段及叶片连接的基座,16.2.8练习：在蜗轮叶片上建立约束方程(续),1.建模并划分单元2.选择基座上的单元(mat2)3.选择叶片底面上的节点a.首先，unselect附在底座单元上的节点(接第2步)b.然后，在位置Z0处reselect节点4.在所选的相邻区域生成约束方程.5.选择everything6.求解并进行后处理,16.2.9耦合练习-叶轮叶片,说明：对叶轮的30扇区使用耦合。确定叶片在绕Z轴1000弧度/秒角速度载荷下的vonMises应力分布。,16.2.9耦合练习-叶轮叶片(续),载荷和材料特性,16.2.9耦合练习-叶轮叶片(续),1.按教师指定的工作目录，用“cp-blade”作为作业名，进入ANSYS。2.恢复“cp-blade.db1”数据库文件:UtilityMenuFileResumefrom或使用命令:RESUME,cp-blade,db1,16.2.9耦合练习-叶轮叶片(续),3.进入前处理器，分别定义单元类型1为SOLID95,单元类型2为MESH200。对MESH200单元设置KEYOPY(1)=5：MainMenuPreprocessorElementTypeAdd/Edit/DeleteAdd.选择“StructuralSolid”和“Brick20node95”,然后按Apply选择“NotSolved”and“MeshFacet200”,然后按OK选择Options.SetK1=“TRIA6-NODE”,然后按OKClose或使用命令:/PREP7ET,1,SOLID95ET,2,MESH200KEYOPT,2,1,5,16.2.9耦合练习-叶轮叶片(续),4.使用VSWEEP对体volume2进行网格剖分:MainMenuPreprocessorMeshTool选择“Hex”（六面体）和“Sweep”（扫掠）,然后选择Sweep或使用命令:VSWEEP,2,16.2.9耦合练习-叶轮叶片(续),5.选择“智能尺寸”等级4并用MESH200单元对1号面剖分网格（扇区底侧边界）:MainMenuPreprocessorMeshTool智能尺寸”等级置为4Mesh置为Areas选择“Tri”和“Free”,然后按Mesh或使用命令:SMRT,4AMESH,1,16.2.9耦合练习-叶轮叶片(续),6.拷贝1号面的网格到11号面(扇区高段一侧边界):MainMenuPreprocessor-Modeling-CopyAreaMesh+拾取1号面(或者在ANSYS输入窗口键入“1”后按Enter键)OK拾取11号面(或者在ANSYS输入窗口键入“11”后按Enter键)在拾取对话框中选择OK设置KCN=1设置DY=30按OK或使用命令:MSHCOPY,AREA,1,11,1,0,30,16.2.9耦合练习-叶轮叶片(续),7.使用SOLID95对1号体剖分网格:MainMenuPreprocessorMeshTool或使用命令:VMESH,1,16.2.9耦合练习-叶轮叶片(续),8.将SOLID95退化为SOLID92单元:MainMenuPreprocessor-Meshing-ModifyMeshChangeTets.或使用命令:TCHG,95,92,3,16.2.9耦合练习-叶轮叶片(续),9.在柱坐标系(CSYS,1)中，把边界低侧节点自由度耦合到边界高侧节点的自由度:MainMenuPreprocessorCoupling/CeqnOffsetNodes设置KCN=1设置DY=30选择OK或使用命令:CPCYC,ALL,0.0001,1,0,30,16.2.9耦合练习-叶轮叶片(续),10.在X=0(或者在柱坐标系中r=0)的节点上约束UX和UY:10a.在X=0选择节点:UtilityMenuSelectEntities.或使用命令:NSEL,S,LOC,X,010b.约束节点UX、UY自由度:MainMenuPreprocessorLoads-Loads-Apply-Structural-DisplacementOnNodes+或使用命令:D,ALL,UX,UY,16.2.9耦合练习-叶轮叶片(续),11.为了防止在Z方向的刚体运动（轴向），约束坐标原点处节点(2426号节点)UZ自由度:11a.再选择Z=0处节点子集:UtilityMenuSelectEntities.或使用命令:NSEL,R,LOC,Z,011b.约束所选节点UX,UY,和UZ自由度:MainMenuPreprocessorLoads-Loads-Apply-Structural-DisplacementOnNodes+或使用命令:D,ALL,UZ,16.2.9耦合练习-叶轮叶片(续),12.选择所有节点并把节点的坐标系改变到总体柱坐标系:12a.选择所有节点:UtilityMenuSelectEverything或使用命令:NSEL,ALL,16.2.9耦合练习-叶轮叶片(续),12b.把激活坐标系设置为总体柱坐标系:UtilityMenuWorkPlaneChangeActiveCStoGlobalCylindrical或使用命令:CSYS,112c.改变节点坐标系到总体柱坐标系:MainMenuPreprocessor-Modeling-Move/Modify-RotateNodeCS-ToActiveCS+PickAll或使用命令:NROTAT,ALL,16.2.9耦合练习-叶轮叶片(续),13.关闭节点耦合符号:UtilityMenuPlotCtrlsSymbols在对话框中选择“ForIndividual:”和“Miscellaneous”,然后选择Ok设置CP=“Off”选择OK或使用命令:/PBC,CP,0,16.2.9耦合练习-叶轮叶片(续),14.检查单元:MainMenuPreprocessor-Meshing-CheckMeshSelBadElems在对话框中选择OK选择CloseUtilityMenuPlotElements选择CloseUtilityMenuSelectEverythingUtilityMenuPlotElements或使用命令:CHECK,ESEL,WARNEPLOTESEL,ALLEPLOT,16.2.9耦合练习-叶轮叶片(续),15.存储数据库并获取解答:Pickthe“SAVE_DB”buttonintheToolbar(orselect:UtilityMenuFileSaveasJobname.db)MainMenuSolution-Solve-CurrentLS查看“/STATUSCommand”然后关闭对话框选择OK选择Close-关闭黄色信息框完成求解或使用命令:SAVE/SOLUSOLVE,16.2.9耦合练习-叶轮叶片(续),16.进入后处理器，画出vonMises应力:MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:/POST1PLNSOL,S,EQV,16.2.9耦合练习-叶轮叶片(续),17.画出2号体（叶片）的vonMises应力:UtilityMenuSelectEntitiesSeleBelow(toselecteverythingbelowselectedvolumes)MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:VSEL,S,2ALLSEL,BELOW,VOLUPLNSOL,S,EQV,16.2.9耦合练习-叶轮叶片(续),18.画出1号体（基座）的vonMises应力:UtilityMenuSelectEntitiesSeleBelow(toselecteverythingbelowselectedvolumes)Replot或使用命令:VSEL,S,1ALLSEL,BELOW,VOLUPLNSOL,S,EQV,16.2.9耦合练习-叶轮叶片(续),19.选择全部实体并把结果扩展360度:UtilityMenuSelectEverythingUtilityMenuPlotCtrlsStyleSymmetryExpansionUser-SpecifiedExpansion设置NREPEAT=12设置TYPE=Polar设置DY=30,然后选择OKUtilityMenuPlotCtrlsPan,Zoom,Rotate.选择ISO选择Fit或使用命令:/EXPAND,12,PLOAR,FULL,30/VIEW,1,1,1,1/AUTO,1/REPLOT,16.2.9耦合练习-叶轮叶片(续),20.存储并退出ANSYS:Pickthe“QUIT”buttonintheToolbar或使用命令:SAVEFINISH/EXIT,16.2.10约束方程练习-涡轮叶片,说明使用约束方程连接涡轮叶片和叶片支座，然后施加X方向1000in/sec2的加速度进行应力分析。,16.2.10约束方程练习-涡轮叶片(续),载荷与材料特性,16.2.10约束方程练习-涡轮叶片(续),1.用“ce-blade”作为作业名，进入ANSYS。2.恢复数据库文件“ce-blade.db1”:UtilityMenuFileResumefrom或使用命令:RESUME,ce-blade,db1,16.2.10约束方程练习-涡轮叶片(续),3.选择基座单元(材料号为2)然后画出单元:UtilityMenuSelectEntitiesUtilityMenuPlotElements或使用命令:ESEL,S,MAT,2EPLOT,16.2.10约束方程练习-涡轮叶片(续),4.选择叶片表面在基座上的节点:4a.选择Z=0位置的节点然后画出它们:UtilityMenuSelectEntitiesUtilityMenuPlotNodes或使用命令:NSEL,S,LOC,Z,0NPLOT,16.2.10约束方程练习-涡轮叶片(续),4b.不选择与单元相连的节点:UtilityMenuSelectEntitiesUtilityMenuPlotNodes或使用命令:NSLE,UNPLOT,16.2.10约束方程练习-涡轮叶片(续),5.进入前处理器，生成“邻近区域”的约束方程:MainMenuPreprocessorCoupling/CeqnAdjacentRegions.认可缺省的单元容差(TOLER=0.25)选择OKUtilityMenuPlotElements或使用命令:/PREP7CEINTF,0.25EPLOT,16.2.10约束方程练习-涡轮叶片(续),6.选择全部实体，并约束全部基座前面（25号面）上的全部自由度:UtilityMenuSelectEverythingUtilityMenuPlotCtrlsNumbering设置AREA=“On”设置Elem/Attribnumbering=“Nonumbering”设置/NUM=Colors&numbers选择OKUtilityMenuPlotAreasMainMenuPreprocessorLoads-Loads-Apply-Structural-DisplacementOnAreas+,16.2.10约束方程练习-涡轮叶片(续),或使用命令:ALLSEL,ALL/PNUM,AREA,1/NUMBER,0APLOTDA,25,ALL,16.2.10约束方程练习-涡轮叶片(续),7.画出单元，然后存储数据库:UtilityMenuPlotElements在工具条上拾取“SAVE_DB”或使用命令:EPLOTSAVE8.求解:MainMenuSolution-Solve-CurrentLS或使用命令:/SOLUSOLVE,16.2.10约束方程练习-涡轮叶片(续),9.进入通用后处理器，画出vonMises应力:MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:/POST1PLNSOL,S,EQV,16.2.10约束方程练习-涡轮叶片(续),10.存储并退出ANSYS:Pickthe“QUIT”buttonintheToolbar选择“SaveEverything”选择OK或使用命令:SAVEFINISH/EXIT,ALL,16.2.11约束方程练习-摆杆刚性区,说明用“刚性区”(约束方程)替代摆杆右边2/3模型.刚性摆杆(EX=2.11e6kgf/cm2),右端施有500kgf的竖向力作用(在X=33cm处).注释:1kgf=9.81N,1kgf/cm2=98066.5Pa,16.2.11约束方程练习-摆杆刚性区(续),载荷和材料特性,16.2.11约束方程练习-摆杆刚性区(续),1.用“swaybar-cerig”作为作业名，进入ANSYS。2.读入文件“swaybar1.inp”:UtilityMenuFileReadInputfrom或使用命令:/INP,swaybar1,inp,16.2.11约束方程练习-摆杆刚性区(续),3.在X=33,Y=0,Z=0处定义节点(节点号为8000):MainMenuPreprocessor-Modeling-CreateNodes或使用命令:N,8000,33,16.2.11约束方程练习-摆杆刚性区(续),4.设置单元类型为2并在8000号节点上建立MASS21单元:MainMenuPreprocessor-Attributes-DefineDefaultAttribs.MainMenuPreprocessor-Modeling-CreateElements-AutoNumbered-ThruNodes+或使用命令:TYPE,2E,8000,16.2.11约束方程练习-摆杆刚性区(续),5.在8000号节点上施加FY方向-500kgf的集中载荷:MainMenuPreprocessorLoads-Loads-Apply-Structural-Force/MomentOnNodes+或使用命令:F,8000,FY,-500,16.2.11约束方程练习-摆杆刚性区(续),6.选择X=10的全部节点和8000号节点:UtilityMenuSelectEntitiesUtilityMenuPlotNodes(OrpickPLOTonSelectionTool).或使用命令:NSEL,S,LOC,X,10NSEL,A,8000NPLOT,在这种情况下,不能在菜单中用“ALL”拾取。直接输入命令会更方便。,16.2.11约束方程练习-摆杆刚性区(续),7.在X=10的全部节点和8000号节点之间生成自由度UX、UY的刚性连线:MainMenuPreprocessorCoupling/CeqnRigidRegion+或使用命令:CERIG,8000,ALL,UX,UY,16.2.11约束方程练习-摆杆刚性区(续),8.选择所有节点,打开力和约束方程符号并画出单元:UtilityMenuSelectEntitiesUtilityMenuPlotCtrlsSymbolsUtilityMenuPlotElements或使用命令:NSEL,ALL/PBC,F,1/PBC,CE,1EPLOT,16.2.11约束方程练习-摆杆刚性区(续),9.存储数据库，然后用PCG优化求解器求解:MainMenuSolution-AnalysisType-SolnControl.UtilityMenuFileSaveasJobname.dbMainMenuSolution-Solve-CurrentLS或使用命令:/SOLUEQSLV,PCGSAVESOLVE,16.2.11约束方程练习-摆杆刚性区(续),10.不选MASS21单元，然后进入后处理器POST1画出SEQV（vonMises应力）的节点解:UtilityMenuSelectEntities.MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:ESEL,U,ENAME,21/POST1PLNSOL,S,EQV,16.2.11约束方程练习-摆杆刚性区(续),11.画出位移:MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:PLNSOL,U,SUM,16.2.11约束方程练习-摆杆刚性区(续),12.如果时间允许，求解整个模型，然后与使用刚性连接模型的结果进行比较:UtilityMenuFileReadInputfromSelectthe“swaybar2.inp”file,thenOK或使用命令:/INP,swaybar2,inp,16.2.11约束方程练习-摆杆刚性区(续),13.求解完成后,进入后处理器POST1并画出SEQV（vonMises应力）的节点解:MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:/POST1PLNSOL,S,EQV,16.2.11约束方程练习-摆杆刚性区(续),14.选择体1到体15以及与体1到体15有关的所有单元和节点，然后画出节点上的SEQV（vonMises应力）解:UtilityMenuSelectEntitiesMainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:VSEL,S,1,5ALLSEL,BELOW,VOLUPLNSOL,S,EQV/AUTO,1/REPLOT,16.2.11约束方程练习-摆杆刚性区(续),15.画出位移:MainMenuGeneralPostprocPlotResults-ContourPlot-NodalSolu.或使用命令:PLNSOL,U,SUM,16.2.11约束方程练习-摆杆刚性区(续),16.存储并退出ANSYS:在工具条上拾取“QUIT”。选择“SaveEverything”选择OK或使用命令:SAVEFINISH/EXIT,ALL,