Ansys中节点力提取

合集下载

ANSYS计算结果无难事，APDL经典命令让你的模型舞起来

ANSYS计算结果无难事，APDL经典命令让你的模型舞起来1、让你的ANSYS模型'舞'起来ANSYS计算结果的动画可采用ANTIME、ANMODE、ANCNTR、ANHARM等自动生成动画(AVI格式)，使结果展示更加生动直观，相信使用ANSYS的都会制作。

然而，几何模型或有限元模型则无动画显示功能，有时为展示模型本身，会从多个角度截取图片。

那么，模型能否也可制作动画呢？答案是肯定的。

利用ANSYS的图形存储命令/SEG可以实现此功能，让你的模型动起来。

具体过程详见命令流中及其注释，动画上传总是失败，自己生成不要观看吧。

Finish$/clear$/prep7!简单的创建几何模型以减少篇幅blc4,0,0,4,2,5cyl4,2,4,1,,2,,4!关闭图例信息/plopts,info,off!以下开始制作模型动画!删除当前储存的图形/seg,dele/seg,multi,jhdh,1 !独立存储且不覆盖，文件名为jhdh/auto,1 !自动计算与图形区合适显示方式!正视/view,1,0,0,1$vplot!侧视/view,1,1$vplot!俯视/view,1,,1$vplot!D视图/view,1,1,1,1$vplot!循环36次，每次改变10度视角*do,i,1,36$/ang,1,10,ys,1$/replot$*enddo!关闭图形存储操作，保存为jhdh.avi文件/seg,off$/anfile,save,jhdh,avi其实比较简单，一旦进入模型动画制作过程，所有的xPLOT(x=KLAVNE)绘制的图形都将进入动画序列，按显示过程形成一部连续的动画。

2、用一个命令解决ANSYS数据列表分页早年初学ANSYS时，经常用到xLIST(如NLIST、ELIST、KLIST、LLIST、ALIST、VLIST等命令)和PRxSOL（如PRNSOL、PRESOL、PRRSOL、PRETAB、PRPATH）等列表命令，并希望将这些内容保存到TXT文件中，然后再导入EXCEL中处理。

ansys中beam单元应力结果

ansys中beam单元应力结果ANSYS中的Beam单元（BEAM3或BEAM4）的应力结果可以通过查看在节点上的S1、S2和S3应力值来获得。

这些应力值表示相对于单元坐标系的正常应力，其中S1是沿单元长度方向的应力，S2是沿单元的切应力，S3是沿垂直于单元断面的应力。

要查看Beam单元的应力结果，请按照以下步骤操作：1. 在ANSYS中加载您的模型，并选择要查看应力结果的Beam单元。

2. 打开"Solution"菜单，并选择"Analysis Type"下的"Static Structural"。

3. 选择"Insert"菜单下的"Results"，然后选择"Engineering Data"。

4. 在Engineering Data窗口中，找到"Beam Stresses"项，并展开它。

5. 选择你要查看的Beam单元，然后单击"OK"。

6. 在主窗口中，选择"General Postproc"界面。

7. 单击"List"，然后选择"Nodal Solution"。

8. 在"Nodal Solution"窗口中，选择"Beam Stress"项，并单击"OK"。

9. 现在，您可以按节点编号或坐标对应的应力结果进行查看。

S1、S2和S3的值将以相应的单位（例如MPa）显示。

请注意，Beam单元的应力结果是基于Beam单元的假设和约束，因此可能不是准确的在复杂几何形状和应力条件下。

在解决复杂问题时，可能需要使用其他类型的单元来更准确地评估应力分布。

ANSYS后处理(结果查看)

一、显示某个时间点的温度云图1、General Postproc →Read Result →By Time/Freq2、在跳出的窗口中输入时间点，点击OK按钮3、然后点Plot Results按下图操作3、然后点击plot →Replot即可显示该时刻的云图二、提取某个节点的数值1、首先通过下列命令，选择部分单元nsel,s,loc,x,0,0.025esln,all然后读取所需节点的编号。

2、点击时间历程后处理器TimeHist postproc弹出如箭头所指对话框。

点击图对话框左上角的绿色增加按钮弹出对话框点击ok按钮，在弹出的对话框中输入节点编号，或者鼠标点击选择节点即可将新的数据读入对话框中如下图所示然后即可通过窗口上的按钮对数据进行操作处理。

/POST1set,last !定义数据集从结果文件中读出，last表示读取最后的数据集plnsol,s,eqv !以连续的轮廓线形式显示结果，S表示应力，EQV表示等效应力查看某个截面的云图!-----------------选取节点结果/post1!seltol,1.0e-10set,,,,,2.5!nsel,s,loc,y,0.1,0.1nsel,s,loc,x,0.02/page,99999,132,99999,240!-------------------显示某个截面wprota,,,90wpoffs,,,0.02/CPLANE,1 !指定截面为WP/TYPE,1,5 !结果显示方式选项工作平面移回全局坐标原点WPCSYS,-1nsel,s,loc,x,0,0.025esln,,1,ACTIVE。

ansys技巧总结_将面载荷转化为等效节点力施加的方法

ansys技巧总结_将面载荷转化为等效节点力施加的方法将面载荷转化为等效节点力施加的方法在进行分析时，有时候需要将已知的面载荷按照节点力来施加，比如载荷方向及大小不变的情况（ANSYS将面力解释为追随力，而将节点力解释为恒定力），那么，在只知道面力的情况下，如何施加等效于该面力的等效节点力呢？可以通过如下步骤给有限元模型施加与已知面载荷完全等效的节点力：1、在模型上施加与已知面力位置、大小相同但方向相反的面力。

Main Menu->Solution->Apply->Pressure->。

（注意：所施加面力要与已知力反号）。

2：将模型的所有节点自由度全部约束。

Main Menu->Solution->Apply->Displacement->On Nodes3：求解模型。

Main Menu->Solution->Current LS（这一步会生成结果文件Jobname.rst）4：开始新的分析：Main Menu->Solution->New Analysis5：删除前两步施加的面力和约束。

Main Menu->Solution->Delete->Pressure->Main Menu->Solution->Delete-> Displacement->On Nodes 6：从Jobname.rst中保存的支反力结果施加与已知面力完全等效的节点力。

Main Menu->Solution->Apply->Force/Moment->From Reaction7：施加其它必要的载荷和约束，然后求解。

这个方法基于ANSYS计算的支反力。

由于第二步将模型的所有自由度全部约束，所以通过第三步的求解，模型内不会产生任何的应力，为了保持载荷作用面上节点的平衡，这些节点上将产生与第一步施加面力大小完全等效但方向相反的节点支反力。

Ansys中节点力提取

Ansys中节点力提取几个问题的说明对于ansys中节点力提取的命令，一般有如下命令可以用，*GET,Par, NODE, N, RF,FX(FY/FZ/MX/MY/MZ)这组命令是我们最开始用的，用来提取节点反力，但是有个缺陷，节点反力只在有约束位置才能提取，如果在结构中任何一个节点处提取此节点所受合力，界面操作有两种方法。

Main Menu>General Postproc>Nodal Calcs>Total Force SumMain Menu>General Postproc>Nodal Calcs>Sum @ Each Node但是执行上面两个操作有个前提，需要选出对应的单元和节点，下面举例说明：如下图：800臂架结构由于要对连接架+塔帽进行单独的详细分析，需要提取旋转架与塔帽连接处铰点对塔帽的作用力。

而且为了在详细模型中施加载荷的时候方便，提取结果的坐标系需要是X向沿着主臂的局部坐标系，见示图1。

运用Main Menu>General Postproc>Nodal Calcs>Total Force Sum 或者Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node这两个操作可以实现。

下面就这两个操作的的结果进行比较。

就如上模型，研究塔帽和旋转架连接左侧铰耳处受力提取，见示图2。

将塔帽上与此铰点相连接的两个单元选出来，选择此节点，见示图3。

读取结果文件，设置结果坐标系为要求的局部坐标系（文件中为局部坐标系11）。

执行Main Menu>General Postproc>Nodal Calcs>Total Force Sum，选取Active Rsys，结果界面如下：执行Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node，结果界面如下：从结果界面上可以看出，了；两个结果界面中显示的结果来看，“结果显示1”与“结果显示3”的结果数据相同，而“结果显示2”与其他两个在在力的显示数据上相同，而矩的显示数据上并不相同。

ANSYS单元特性之命令流算例

EX1.1 （LINK1‎‎）（1）进入后处理‎模块,显示节点位‎移和杆件内‎力MID_N‎O DE = NODE (A/2,-B,0 )! 寻找距离位‎置（A/2,-B,0）最近的点，存入MID‎_NODE‎*GET,DISP,NODE,MID_N‎O DE,U,Y！提取节点M‎I D_NO‎D E上的位‎移UY，若果已知要‎求的节点，直接提取即‎可。

LEFT_‎E L = ENEAR‎N (MID_N‎O DE)! 需找距离节‎点MID_‎N ODE最‎近的单元，存入LEF‎T_ELETABL‎E,STRS,L S,1! 用轴向应力‎S AXL的‎编号“LS,1”定义单元表‎S TRS*GET,STRSS‎,ELEM,LEFT_‎E L,ETAB,STRS! 从单元表S‎‎TRS中提‎取LEFT‎‎_EL单元‎的应力结果‎，存入变量S‎‎TRSS。

注意：提取的轴向‎应力结果具‎体到指定的‎单元。

（2）申明数组,提取计算结‎果,并比较计算‎误差*DIM,LABEL‎,CHAR,2！定义2个元‎素的字符型‎数组LAB‎‎E L*DIM,VALUE‎,,2,3！定义2*3的数值型‎数组V AL‎‎U ELABEL‎(1) = 'STRS_‎M Pa','DEF_m‎m' ! 给字符型数‎组的第1个‎元素赋值*VFILL‎,VALUE‎(1,1),DATA,1,-0.05498‎! 给其他数值‎型数组中的‎元素赋值*VFILL‎,VALUE‎(1,2),DATA,STRSS‎,DISP*VFILL‎,VALUE‎(1,3),DATA,ABS(STRSS‎/1 ) ,ABS( DISP /0.05498‎)/OUT,EX1_1‎,out ！将输出内容‎重定向到文‎件EX1_‎‎1.out/COM ! 以注释形式‎输出内容/COM,------------------- EX1.1 RESUL‎T S COMPA‎R ISON‎---------------------/COM,/COM, | TARGE‎T | ANSYS‎ | RATIO‎/COM,*VWRIT‎E,LABEL‎(1),VALUE‎(1,1),VALUE‎(1,2),VALUE‎(1,3)(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)/COM,----------------------------------------------------------------/OUT ! 结束数据重‎定向，关闭输出文‎件FINIS‎H*LIST,EX1_1‎,out ! 列表显示文‎件EX1_‎‎1.out的内‎容EX1.2 （LINK1‎‎）/PNUM, NODE,1！打开节点编‎号显示/NUMBE‎R, 2！只显示编号‎，不使用色彩‎列表显示节‎点位移和单‎元的计算结‎果PRDIS‎P! 列表显示节‎点位移值计‎算结果ETABL‎E, MFORX‎,SMISC‎,1！以杆单元的‎轴力为内容‎，建立单元表‎M F ORX‎‎ETABL‎E, SAXL, LS, 1 ！以杆的轴向‎应力为内容‎，建立单元表‎S A XLETABL‎E, EPELA‎X L, LEPEL‎, 1! 以杆单元的‎轴向应变为‎内容，建立单元表‎E PELA‎‎X LPRETA‎B! 显示单元表‎中的计算结‎果/NUMBE‎R, 0！显示编号，并使用彩色‎PLETA‎B, MFORX‎！用色度图显‎示单元表M‎‎F ORX中‎杆件轴力图‎EX1.3 （LINK1‎‎）NSEL,S,LOC,Y,1.0 ！选择所有位‎于Y=1.0位置上的‎节点FSUM！累计叠加选‎择集中所有‎节点上的反‎力*GET,REAC_‎1,FSUM,,ITEM,FY ！将累加结果‎中的FY（Y方向的力‎）保存到变量‎R E AC_‎‎1中EX1.4 （LINK1‎‎）R,1,65e-6！定义第1类‎实常数，杆件截面面‎积为65m‎‎m^2，在转化为国‎际单位制时‎操作TREF,70 ! 设定参考温‎度为70度‎BFUNI‎F,TEMP,80 ! 温度从原来‎的70度均‎匀上升到8‎‎0度（TREF+10）EX1.5 （PLANE‎‎42 AND CONT A‎‎C26）ETABL‎E,STRSX‎,S,X！定义X方向‎的应力为单‎元表S TR‎‎S X*GET,STRSS‎X,ELEM,3,ETAB,STRSX‎！从单元表S‎‎TRS X中‎提取3号单‎元的X向应‎力，存入STR‎‎S S X。

Ansys中ETABLE命令详解(1)

Ansys中ETABLE命令详解1、ETABLE简介原来ansys求解完问题之后，会把计算结果保存在一个表里面，表的行是单元的编号，表的列是单元的计算结果，如节点位移、节点力、应力、应变等等。

ETABLE这个命令就是把表中的某个列取出来，赋值给某个自定义的向量，再通过*GET命令可以指定某个具体的单元，就可以把该单元的对应计算结果提取出来了。

ETABLE, Lab, Item, Comp 将单元值形成一个表以便进一步的处理。

2、命令选项解释：1）Lab用户自定义的表名，用于后续命令或输出的标题，最多可使用8个字母，不可与预定义的表名称重复。

默认的表名是Item和Comp项的前四个字母组合而成的8个字母。

如果与用户之前定义的表名相同，本次结果将被包括在同一表中。

最多可定义200个不同的表名。

以下表名是ANSYS预定义的，不可用作用户自定表名：REFL, STAT, 和ERAS.Lab = REFL以ETABLE的最新选项重写所有ETABLE命令预定义的表，但保留字段将被忽略，这个命令在载荷步改变后重写表时很方便。

Lab = STAT将显示储存的表的值。

Lab = ERAS将删除整个表。

2）Item选项名称。

常用的选项名称见后表。

某些选项需要栏目名。

Item = ERAS将删除表中的某一栏。

3）Comp选项的栏目名（如果需要的话）。

常用的栏目名见后表。

4）说明定义单元值的表以便后续处理。

单元值表可以被认为是工作表，其行代表所有被选择的单元。

其列代表通过ETABLE命令输入表中的单元值。

每一列数据有一个用户定义的名称，用于列表和显示。

将数据输入单元表后，你不仅可以列出和显示你的数据，还可以对数据进行许多操作，例如列相加或列相乘[SADD, SMULT]，为安全计算定义允许的应力[SALLOW]，或者将一列数据和另一列相乘 [SMULT]。

更多的细节请看ANSYS Basic Analysis Guide。

节点应力单元应力

2ansys提取六面体单元切（x向）应力，提取得到是8个节点的切应力，求这个单元切应力应怎么处理。
检举 | 2012-6-24 11:04 提问者： kaiti123 | 悬赏分：100 | 浏览次数：65次
我是这么提取的General postproc ——list result——Element solution——stress——X component of stress，得到的单元应力是八个节点的应力，我要的是单元应力，单元应力应该如何让处理，单元应力与节点应力有什么关系，请教大家了。
做ansys静力分析时，提取的单元应力和节点应力不一样，单元应力值是296MPa,节点应力值是240MPa,想问问这个分析强度时，选哪个应力为准呢？单元应力和节点应力为什么会有差别呢？
问题补充：所取应力均为von mises stress,谢谢
1你稍微翻看下有限元教材就很容易搞懂了。
ansys的初始解是节点变形量。由节点变形量用形函数方程求出单元应力。在有限元模型中，同一个节点周围有若干个单元。这些单元的应力均值就是该节点的应力值。
追问
非常感谢你的回答！我要算单元的应力就是把八个节点的应力加起来再除以8，是这样做的吗？
所以，无论什么软件，有限元计算的初始解永远是节点变形量。接下来的导出解是单元应力值。接下来的导出解是节点应力值。
至于选择哪个不是问题的重点。因为网格划分总够密集，两者的值会逐步接近的。如果差距较大，很可能是你网格划分不够。建议逐步细化，应力值变化不大时，就可以了。如果不想这么麻烦，我建议还是参考大值吧，稳妥起见。
所以，无论什么软件，有限元计算的初始解永远是节点变形量。接下来的导出解是单元应力值。接下来的导出解是节点应力值。

ansys读取某点应力和应力集中系数的问题

ansys读取某点应力和应力集中系数的问题2009-10-13 10:20提问者：我浩然哦|浏览次数：3086次热应力计算结束后后处理我想查看模型某点的热应力，该怎么操作？还有应力集中部位有多个但我想查看其中某个位置的应力集中系数可以查看吗？怎么操作？我来帮他解答精彩回答2009-10-13 22:33GUI操作：在General Postproc——Query Results——Subgrid Solu，选择你想显示的节点。

命令流：1. 最简单的办法是使用NSORT，打印出结果，可以通过控制使其输出到文件2. 使用apdl能复杂一点，下面是以前经常用的一段命令流，参考着修改一下吧*CREATE,GET_node_inf,mac,*GET,Nnod,NODE,0,COUNT !获取所选择的节点总数*DIM,S_Xyz,ARRAY,NNOD,5 !定义1个数组存放数据*GET,Nd,NODE,0,NUM,MIN !获取最小的节点编号*DO,I,1,Nnod,1S_Xyz(I,1)=Nd !将节点列表放数组第1列S_Xyz(I,2)=NX(Nd) !节点的X坐标放数组第2列S_Xyz(I,3)=NY(Nd) !节点的Y坐标放数组第3列S_Xyz(I,4)=NZ(Nd) !节点的Z坐标放数组第4列!*GET,S_Xyz(I,5),NODE,ND,S,EQV !节点的von mises值放数组第5列*GET,S_Xyz(I,5),NODE,ND,U,SUM !节点的总变形值值放数组第5列Nd=NDNEXT(Nd) !读出下一个节点编号*ENDDO*END*CREATE,OUT_node_inf,mac,*CFOPEN,node_info,txt,,*VWRITE,S_Xyz(1,1),S_Xyz(1,2),S_Xyz(1,3),S_Xyz(1,4),S_Xyz(1,5)(F10.0,3F15.4,E15.5)*CFCLOS*ENDGET_node_infOUT_node_inf/delete,GET_node_inf,mac/delete,OUT_node_inf,mac另附1.先对节点的值进行SORT，在提取最大的值即可。

ansys常见问题解答

ANSYS常见问题及应用技巧本篇开始讲述ANSYS在使用过程中常见的问题和在使用时一些常用的使用技巧，对与初学者来说，理解和弄清楚这些问题的根源和掌握这些使用技巧，能够更好的理解ANSYS这个软件本身。

1.ANSYS中的等效应力是什么物理含义？1）ANSYS中等效应力最大应力s1有什么区别，平常讨论应力分布，应该用等效应力还是最大应力s1呢？2）计算等效应力时是否需要输入等效泊松比呢？3）在实际的应用中，例如在讨论平板上的圆孔应力集中的应力分布问题时，应该用等效应力来描述应力集中的现象，还是采用主应力s1来反应集中的程度呢？还是采用一个单方向的sx来说明问题呢？答：1）这个等效应力应该就是弹塑性力学里的VonMises应力，他主要考察的是材料在各个方向上的应力差值，因为在实验室里获得材料强度都是单向载荷作用下的强度（当然现在也有三轴应力实验仪），所以有时候材料所受的单向载荷可能很大，但并没有造成破坏，这是就是看他的等效应力，具体计算公式是: σ等效=sqrt{0.5[(σ1-σ2)^2+(σ2-σ3)^2+(σ3-σ1)^2]}2）等效应力是三项主应力的组合如s，int即为max(si-sj),si,sj为三项主向应力。

i，j=1，2，3 i≠j即tresca型s,eqv为sqrt(0.5*∑(si-sj)**2)，i，j=1，2，3 i≠j即mises型3）个人认为应该采用等小应力来描述应力集中的现象，因为在实际中很难找到真正的单轴拉压的情况，一般结构的受力都没有这么简单，所以在分析的时候需要用等效应力来将各主应力进行转化，因此应该用等效应力来描述应力集中的现象。

4）等效泊松比就是泊松比，等效应力计算时不会用到泊松比，不过在计算mises 等效应变时会用到。

对于泊松比的取值原则应遵循以下两条：a:对于elastic & thermal strains 泊松比取为材料的泊松比；b:对于plastic creep hyperelastic strains 泊松比取为0.5。

ANSYS积分法和节点位移法求解应力强度因子附命令流

K I = √2π K II = √2π 其中： G为剪切模量；
G ∆v 1 + κ √r G ∆u 1 + κ √r
κ为材料常数，对于平面应力问题，取
3−�� 1+��
；
∆u为裂纹面在某点处的水平相对位移； ∆v为裂纹面在某点处的垂直相对位移。
图 1-6
位移法图解
根据断裂力学对于三种裂纹的定义，当∆v>0 时，K I 为正，裂纹上下面相对位移为顺时针为正，即顺时针时，∆u>0，K II 为正；反之为负。理论上，当取上下裂纹面同一位置的点，当该点趋向于裂尖时，结果更精确，本算例取奇异单元上 1/4 处的节点的位移进行计算，计算模型同上。首先，先对有限元模型进行求解，然后进入到后处理层，求出在局部坐标系系下，所处裂纹上下面的奇异单元上 1/4 处节点的水平及竖直位移 ux，uy，然后求出裂纹面的相对位移∆u、∆v，最后代入上式即可。计算结果如图 1-7 所示：KI=223.84Mpa*(mm)1/2，KII=217.63Mpa*(mm)1/2。计算误差分别为：3.1%、0.25%。
FINISH /CLEAR /TITIE,INTERACTIVE INTEGRATION METHOD BY IDUTER-ANSYS /PREP7 /RGB,INDEX,100,100,100, 0 /RGB,INDEX, 80, 80, 80,13 /RGB,INDEX, 60, 60, 60,14 /RGB,INDEX, 0, 0, 0,15 /REPLOT !------------------!UNIFIED UNIT(N,MM) PI=ACOS(-1) *SET,H,80 *SET,W,50 *SET,A,0.12*W *SET,BETA,90-45 *SET,ALPH,(90-BETA)*PI/180 *SET,SIGMA,100 R1=1 R2=2 R3=3 !THE HEIGHT OF MODEL !THE WEIGHT OF MODEL !HALF LENGTH OF THE ANGLED CRACK !THE INCLINED ANGLE OF CRACK ! RADIAN SYSTEM !SIGMA !FIRST ROW OF ELEMENT RADIUS !THIRD ROW OF ELEMENT RADIUS !SIXTH ROW OF ELEMENT RADIUS

（完整版）ansysetable提取结果

（完整版）ansysetable提取结果ANSYS etable提取结果etable提取运算结果ansys求解完结束后，会把计算结果保存在结果表中，行是单元的编号，列是计算信息，如节点位移、节点力、应力、应变等。

ETABLE这个命令就是把表中的需要的信息取出来，赋值给某个自定义的向量，再通过*GET命令可以指定某个具体的单元，就可以把该单元的对应计算结果提取出来了。

ETABLE, Lab, Item, Comp命令选项解释：Lab自定义的表名，用于后续命令或输出的标题，最多可使用8个字母，不可与预定义的表名称重复。

默认的表名是Item和Comp项的前四个字母组合而成的8个字母。

如果与用户之前定义的表名相同，本次结果将被包括在同一表中。

最多可定义200个不同的表名。

以下表名是ANSYS预定义的，不可用作用户自定表名：REFL, STAT, 和ERAS. Lab = REFL以ETABLE的最新选项重写所有ETABLE命令预定义的表，但保留字段将被忽略，这个命令在载荷步改变后重写表时很方便。

Lab = STAT将显示储存的表的值。

Lab = ERAS将删除整个表。

Item选项名称。

常用的选项名称见后表。

某些选项需要栏目名。

Item = ERAS将删除表中的某一栏。

Comp选项的栏目名（如果需要的话）。

常用的栏目名见后表。

说明：定义单元值的表以便后续处理。

单元值表可以被认为是工作表，其行代表所有被选择的单元。

其列代表通过ETABLE命令输入表中的单元值。

每一列数据有一个用户定义的名称，用于列表和显示。

更多的细节请看ANSYS Basic Analysis Guide。

有很多不同类型的数据结果可以被存在单元表中。

ansys读取某点应力和应力集中系数的问题

一种读取Ansys rst文件中节点力信息的方法

一种读取Ansys rst文件中节点力信息的方法作者：马斌谷翠军徐先伟黄艳丽来源：《计算机辅助工程》2022年第03期摘要：为从Ansys rst文件中准确读取焊缝节点力信息，对研究对象进行建模、有限元加载和强度计算，得到包含焊缝单元及焊缝节点编号的rst文件。

用MATLAB软件读取rst文件，筛选出其中的焊缝单元节点力信息，展示焊缝节点力信息的获取方法。

将获得的节点力信息转化为等效结构应力，利用主S-N曲线计算应力损伤比，将计算结果与FE-WELD的计算结果进行对比，验证节点力的准确性，反映读取方法的有效性。

关键词：壳单元; 结果文件; 数据结构; 文件头; 有限元; 焊缝节点力; 坐标系; MATLAB中图分类号： TP391.9文献标志码： BA method for reading node force informationin Ansys rst filesMA Bin， GU Cuijun， XU Xianwei， HUANG Yanli（Chongqing R&D Center， CRRC Qingdao Sifang Locomotive and Rolling Stock Research Institute， Chongqing 401133， China）Abstract： In order to accurately read the weld node force information from the Ansys rst file，the modeling， finite element loading and strength calculation of the research object are carried out，and the rst file containing the weld element and weld node number is obtained. The rst file is read using MATLAB， the node force information of the weld element is screened out， and the method of obtaining the node force information of the weld is shown. The obtained node force information are converted into equivalent structural stress， and the stress damage ratio is calculated by using the principal S-N curve. The calculation results are compared with the calculation results of FE-WELD to verify the accuracy of the node force and to reflect the effectiveness of the reading method.Key words： shell element; result document; data structure; file header; finite element; weld joint force; coordinate system; MATLAB0 引言近年來，随着国民经济的发展，铁路车辆的速度越来越快、运量越来越多，这对列车结构提出更高要求。

如何正确理解ANSYS的节点坐标系

如何正确理解ANSYS的节点坐标系2006年10月01日星期日下午 08:15节点坐标系用以确定节点的每个自由度的方向，每个节点都有其自己的坐标系，在缺省状态下，不管用户在什么坐标系下建立的有限元模型，节点坐标系都是与总体笛卡尔坐标系平行。

节点力和节点边界条件（约束）指的是节点坐标系的方向。

时间历程后处理器 /POST26 中的结果数据是在节点坐标系下表达的。

而通用后处理器/POST1中的结果是按结果坐标系进行表达的。

例如: 模型中任意位置的一个圆，要施加径向约束。

首先需要在圆的中心创建一个柱坐标系并分配一个坐标系号码(例如CS,11)。

这个局部坐标系现在成为激活的坐标系。

然后选择圆上的所有节点。

通过使用 "Prep7>Move/Modify>Rotate Nodal CS to active CS", 选择节点的节点坐标系的朝向将沿着激活坐标系的方向。

未选择节点保持不变。

节点坐标系的显示通过菜单路径Pltctrls>Symbols>Nodal CS。

这些节点坐标系的X方向现在沿径向。

约束这些选择节点的X方向，就是施加的径向约束。

注意：节点坐标系总是笛卡尔坐标系。

可以将节点坐标系旋转到一个局部柱坐标下。

这种情况下，节点坐标系的X方向指向径向，Y方向是周向（theta)。

可是当施加theta方向非零位移时，ANSYS总是定义它为一个笛卡尔Y位移而不是一个转动（Y位移不是theta位移）。

有限元分析中的很多相关量都是在节点坐标系下解释的，这些量包括：输入数据：1 自由度常数2 力3 主自由度4 耦合节点5 约束方程等输出数据：1 节点自由度结果2 节点载荷3 反作用载荷等但实际情况是，在很多分析中，自由度的方向并不总是与总体笛卡尔坐标系平行，比如有时需要用柱坐标系、有时需要用球坐标系等等，这些情况下，可以利用ANSYS 的“旋转节点坐标系”的功能来实现节点坐标系的变化，使其变换到我们需要的坐标系下。

Ansys中节点力提取

Ansys中节点力提取几个问题的说明对于ansys中节点力提取的命令,一般有如下命令可以用,GET,Par, NODE, N, RF,FXFY/FZ/MX/MY/MZ这组命令是我们最开始用的,用来提取节点反力,但是有个缺陷,节点反力只在有约束位置才能提取,如果在结构中任何一个节点处提取此节点所受合力,界面操作有两种方法;Main Menu>General Postproc>Nodal Calcs>Total Force SumMain Menu>General Postproc>Nodal Calcs>Sum Each Node但是执行上面两个操作有个前提,需要选出对应的单元和节点,下面举例说明：如下图：800臂架结构由于要对连接架+塔帽进行单独的详细分析,需要提取旋转架与塔帽连接处铰点对塔帽的作用力;而且为了在详细模型中施加载荷的时候方便,提取结果的坐标系需要是X向沿着主臂的局部坐标系,见示图1;运用Main Menu>General Postproc>Nodal Calcs>Total Force Sum 或者Main Menu>General Postproc>Nodal Calcs>Sum Each Node这两个操作可以实现;下面就这两个操作的的结果进行比较;就如上模型,研究塔帽和旋转架连接左侧铰耳处受力提取,见示图2;将塔帽上与此铰点相连接的两个单元选出来,选择此节点,见示图3;读取结果文件,设置结果坐标系为要求的局部坐标系文件中为局部坐标系11;执行Main Menu>General Postproc>Nodal Calcs>Total Force Sum,选取Active Rsys,结果界面如下：执行Main Menu>General Postproc>Nodal Calcs>Sum Each Node,结果界面如下：从结果界面上可以看出,了；两个结果界面中显示的结果来看,“结果显示1”与“结果显示3”的结果数据相同,而“结果显示2”与其他两个在在力的显示数据上相同,而矩的显示数据上并不相同;之所以出现这种差别,原因在于在求矩的作用数据时,矩的作用中心不相同;“结果显示2”中所显示的数据其矩的作用中心为提取载荷点1600078,而“结果数据1”和“结果数据3”中矩的作用中心为结果坐标系的原点;要想使上述3个结果数据显示值相同,只需要将结果显示中矩的作用中心设置到提取点上就可以了;通过Main Menu>General Postproc>Nodal Calcs>Summation Pt>At Node将矩的作用点设置到1600078,再执行命令Main Menu>General Postproc>Nodal Calcs>Total Force Sum,结果界面如下执行Main Menu>General Postproc>Nodal Calcs>Sum Each Node,结果界面如下可以看出,设置对应的矩的作用点后,3组结果显示数据很好的对应,而这个数据,应该就是我们提取的在模型中其他结构传到对应节点上的作用载荷;注意：在本模型提取位置是一个绕Y轴可以转动的铰接位置,理论上MY应该为0,但是在ansys大变形计算中提取出来的结果此处MY并不为0,这个值是ansys内部算法所得,所幸这个值很小,所以可以忽略不再深究;对于节点力提取的程序如下：RSYS,11 指定结果坐标系ESEL,S,,,1600057 选择对应单元ESEL,A,,,1600092 选择对应单元NSEL,S,,,1600078 选择对应节点SPOINT,1600078 矩的作用中心节点设置FSUM,rsys 显示“Total Force Sum”GET,TM_XZJZX, FSUM,0, Item1,FX 提取FXGET,TM_XZJZY, FSUM,0, Item1,FY 提取FYGET,TM_XZJZZ, FSUM,0, Item1,FZ 提取FZGET,TM_XZJZX, FSUM,0, Item1,MX 提取MXGET,TM_XZJZY, FSUM,0, Item1,MY 提取MYGET,TM_XZJZZ, FSUM,0, Item1,MZ 提取MZ。

APDL获取节点和单元的结果

APDL获取节点和单元的结果APDL获取ANSYS结果⽬录1. 获取节点结果1.1 获取节点应⼒结果序号结果APDL注释1X⽅向应⼒*GET,SX,NODE,1,S,X*GET,变量名,类型,编号,应⼒,X⽅向2Y⽅向应⼒*GET,SY,NODE,1,S,Y*GET,变量名,类型,编号,应⼒,Y⽅向3Z⽅向应⼒*GET,SZ,NODE,1,S,Z*GET,变量名,类型,编号,应⼒,Z⽅向4XY剪应⼒*GET,SXY,NODE,1,S,XY*GET,变量名,类型,编号,应⼒,XY⽅向5YZ剪应⼒*GET,SYZ,NODE,1,S,YZ*GET,变量名,类型,编号,应⼒,YZ⽅向6XZ剪应⼒*GET,SXZ ,NODE,1,S,XZ*GET,变量名,类型,编号,应⼒,XZ⽅向7第⼀主应⼒*GET,S1,NODE,1,S,1*GET,变量名,类型,编号,应⼒,第⼀主应⼒8第⼆主应⼒*GET,S2,NODE,1,S,2*GET,变量名,类型,编号,应⼒,第⼆主应⼒9第三主应⼒*GET,S3,NODE,1,S,3*GET,变量名,类型,编号,应⼒,第三主应⼒10von Mises应⼒*GET,SVON,NODE,1,S,EQV*GET,变量名,类型,编号,应⼒,von Mises应⼒1.2 获取节点位移结果序号结果APDL注释1X⽅向位移*GET,UX,NODE,1,U,X*GET,变量名,类型,编号,位移,X⽅向2Y⽅向位移*GET,UY,NODE,1,U,Y*GET,变量名,类型,编号,位移,Y⽅向3Z⽅向位移*GET,UZ,NODE,1,U,Z*GET,变量名,类型,编号,位移,Z⽅向4总位移*GET,USUM,NODE,1,U,SUM*GET,变量名,类型,编号,位移,总位移1.3 获取节点应变结果--总应变总机械应变(Total Mechanical Strain)，即为弹性应变+塑性应变+蠕变应变的和，不包括热应变。

Ansys中节点力提取

ANSYS计算结果无难事，APDL经典命令让你的模型舞起来

ansys中beam单元应力结果

ANSYS后处理(结果查看)

ansys技巧总结_将面载荷转化为等效节点力施加的方法

Ansys中节点力提取

ANSYS单元特性之命令流算例

Ansys中ETABLE命令详解(1)

节点应力单元应力

ansys读取 某点应力和应力集中系数的问题

ansys常见问题解答

ANSYS积分法和节点位移法求解应力强度因子附命令流

（完整版）ansysetable提取结果

ansys读取 某点应力和应力集中系数的问题

一种读取Ansys rst文件中节点力信息的方法

如何正确理解ANSYS的节点坐标系

Ansys中节点力提取

APDL获取节点和单元的结果

ansys读取某点应力和应力集中系数的问题

ansys读取某点应力和应力集中系数的问题