ansys使用函数加载实例1

ANSYS 中使用函数加载的一个简单例子

2013-10-28 07:34:10| 分类：ANSYS 实例 - GUI | 标签：ansys 函数定义函数加载示例|字号订阅本文将通过一个具体实例说明在 ANSYS 中如何使用函数加载，后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。

实例的具体说明：

一个 1/4 圆柱，内半径 30 mm，外半径 42 mm，长度 100mm，如图 1 所示：

所用材料为双线性弹塑性材料，其机械性能为：

弹性模量 E = 201000 Mpa；泊松比μ=0.3

屈服应力σ= 200 Mpa；切线模量 Et = 2010

使用单元类型 solid185 (8 节点六面体单元)。

取整体单元边长 4 mm，然后可以直接对该几何模型划分 MAP 网格，划分网格结果如图 2：

约束条件为：

轴向两个截面为对称边界条件；一个端面约束轴向位移 Uz。

载荷条件为：

在外表面施加变化的压力载荷，载荷函数为：

P (y) = 8e7 + 7E7 * (X/42)

即：

X = 0 (最高点) 时，P = 15E7；

X = 42 （最低点）时，P = 8E7。

我们采用函数方式来施加这一压力载荷，首先定义函数：

在 Solution 模块中，点击菜单路径：

Solution > Define Loads > Apply > Functions > Define/Edit 将会弹出一个函数编辑器，可以在其中定义所需的函数。

在函数编辑器中，函数类型选择为Single equation，即单值函数；计算函数值时使用的插值坐标系 ( (x,y,z) interpreted in CSYS) 选择 0，即总体直角坐标系，如图 3 所示：

然后，在函数编辑器中间位置的“Result = “ 小窗口中输入要定义

的函数表达式，如果表达式中有 x, y, z, time 等变量 (供定义函数时使

用的“自变量”)，可以用 {X},{Y},{Z},{TIME} 等的形式输入；或者点击下面一个小窗口右边的小箭头，会出现一个下拉列表，列出可以选择的

变量，然后从该列表中选择某个自变量，则该自变量会按照上述格式

写入函数中，如图 5 所示：

接下来最好检查一下函数定义是否正确。做法是：点击函数编辑器中的 Graph 按钮，弹出一个“Plot information” 对话框，其中显示了所定义的函数方程，列出了方程中所有的变量。然后可以激活要变化的变量 (对于本函数，只有一个变量– X，会自动激活)，并在下面的输入窗口中输入该变量的最小和最大值，以及在该变量范围中的计算点数，再点击 Graph 按钮，则会绘制函数的曲线；或者点击 List 按钮，则列出各计算点上的函数值。

对于本函数，点击 Graph 后所显示的函数曲线如图 7 所示，可以看到它是一个线性函数：

点击 List 后列出的函数值：

如果我们不希望外表面的压力函数是一个线性函数，而是一个三角函数，使得压力从下向上以余弦方式变化，可以将当前坐标系转换为总体圆柱坐标系 (csys,1)，然后压力函数定义为：

P = 8E7 + 7E7 * cos(Y)

在总体圆柱坐标系中，Y 为角度，这里需要确定 Y 的单位是角度(Degree) 还是弧度 (Radians)。为了方便，在函数编辑器中，将角度单位设置为 Degree。

下面使用函数编辑器来定义这一函数：

然后同样对所定义的函数进行检查。在函数编辑器中点击 GRAPH 按钮，然后在弹出的对话框中输入 X-Axis Range （X 轴范围）为 0 到 90，再点击 Graph 按钮以显示压力载荷曲线，点击 List 按钮列表压力载荷值，结果如图 10 所示：

可以看到，压力不再是线性变化的。

如果我们接受了这一函数定义，需要将它保存起来，以便后续定义载荷时使用。

回到函数编辑器，点击上部菜单中的 File > Save：

图 11 保存函数

然后弹出一个“另存为” 对话框，在其中，输入文件名 (这里输入文件名为 press_test)，(如有需要还可以修改保存文件的目录) 再点击

“保存” 按钮，就可以把所定义的函数保存为以 .func 为扩展名的函数文件。这里将函数保持为文件 press.func：

将ansys作为子程序调用

将ANSYS作为子程序调用 对于优化或参数化设计，可以在VC或FORTRAN中将ANSYS作为子程序调用。具体调用方法如下： 1.在VC中调用ANSYS ::WinExec("d:/ANSYS57/BIN/INTEL/ANSYS57 -b -p ansys_product_feature -i input_file -o output_file",SW_SHOWNORMAL); 2.在FORTRAN中调用ANSYS LOGICAL(4) result RESULT=SYSTEMQQ('d:\ANSYS57\BIN\INTEL\ANSYS57 -b -p ansys_product_feature -i input_file -o output_file') 3.说明 1和2中，input_file为用APDL语言编写的ANSYS输入文件。 ansys_product_feature为你的ANSYS产品特征代码。 需要注意的是，在VC中调用ANSYS时，需要加一条判断语句，以确定ANSYS 已经执行完毕。 在FORTRAN中不需要判断，FORTRAN会等ANSYS执行完毕才继续执行下一条语句。 在VC中，我没有找到与FORTRAN类似的函数，只好加一条循环判断语句。 如果谁能找着这样的函数，请告诉我，谢谢！ 判断方法很简单，只需判断错误文件file.err是否可写就可以了。 因为当ANSYS在运行时，file.err是不可写的，只有当它运行完毕，此文件才可写。 数据文件(假设输出的数据文件名为opt.out): *dim,out1,,2,1 out1(1)=dmax !目标函数 out1(2)=1-eymax !约束条件1 *cfopen,opt,out *vwrite,out1(1),out1(2) (2f10.6) *cfclos 在VC中相应的显示数据文件命令为： result=system("notepad opt.out"); 图形文件（假设ANSYS工作文件名为test,输出jpg图形文件,具体信息请参考命令/show）：/SHOW,JPEG JPEG,QUAL,75, JPEG,ORIENT,HORIZ JPEG,COLOR,2 JPEG,TMOD,1

ANSYS命令流实例

/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容

ansys旋转经典命令流

1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0

ansys使用函数加载实例1

ANSYS 中使用函数加载的一个简单例子 2013-10-28 07:34:10| 分类：ANSYS 实例 - GUI | 标签：ansys 函数定义函数加载示例|字号订阅本文将通过一个具体实例说明在 ANSYS 中如何使用函数加载，后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。 实例的具体说明： 一个 1/4 圆柱，内半径 30 mm，外半径 42 mm，长度 100mm，如图 1 所示： 所用材料为双线性弹塑性材料，其机械性能为： 弹性模量 E = 201000 Mpa；泊松比μ=0.3 屈服应力σ= 200 Mpa；切线模量 Et = 2010 使用单元类型 solid185 (8 节点六面体单元)。

取整体单元边长 4 mm，然后可以直接对该几何模型划分 MAP 网格，划分网格结果如图 2： 约束条件为： 轴向两个截面为对称边界条件；一个端面约束轴向位移 Uz。 载荷条件为： 在外表面施加变化的压力载荷，载荷函数为： P (y) = 8e7 + 7E7 * (X/42) 即： X = 0 (最高点) 时，P = 15E7； X = 42 （最低点）时，P = 8E7。 我们采用函数方式来施加这一压力载荷，首先定义函数： 在 Solution 模块中，点击菜单路径：

Solution > Define Loads > Apply > Functions > Define/Edit 将会弹出一个函数编辑器，可以在其中定义所需的函数。 在函数编辑器中，函数类型选择为Single equation，即单值函数；计算函数值时使用的插值坐标系 ( (x,y,z) interpreted in CSYS) 选择 0，即总体直角坐标系，如图 3 所示： 然后，在函数编辑器中间位置的“Result = “ 小窗口中输入要定义 的函数表达式，如果表达式中有 x, y, z, time 等变量 (供定义函数时使 用的“自变量”)，可以用 {X},{Y},{Z},{TIME} 等的形式输入；或者点击下面一个小窗口右边的小箭头，会出现一个下拉列表，列出可以选择的 变量，然后从该列表中选择某个自变量，则该自变量会按照上述格式 写入函数中，如图 5 所示：

Ansys常见命令流

Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,

K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放

ansys实例命令流-弹塑性分析命令流

/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL

!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,

2 ANSYS加载时间函数的方法

2 ANSYS加载时间函数的方法 Apply/Functions/Define/Edit打开函数编辑器 ●Functions Type：选择函数类型。选择单个方程或多值函数。如果选择后者，必须键入状态变量名，也就是管理函数中方程的变量。当选择一个多值函数时，六状态表格将被激活。 ●Degrees/Radians：选择度或弧度，这一选择仅决定方程如何被运算，而不会影响*AFUN 设置。 ●使用初始变量方程和键区定义结果方程（单个方程）或描述状态变量的方程（多值函数），出如果定义单方程函数，保存方程。如果是定义多值函数，则继续下面的步骤。 ●单击Regime1，键入在函数表格下定义的状态变量的相应的最大最小值限制。 ●定义这个状态的方程。 ●单击Regime2，注意状态变量的最小值限制已被定义并且不可更改，这一特征确保状态保

持连续而无间隙。定义这个状态的最高值限制。 ●定义这个状态的方程。 ●在六个状态中连续如上操作。在每个状态里，不必储存或保存单个方程，除非想在另一状态中重用某个方程。 ●输入一个注释描述函数（可选）选择File/Comments。 ●计算器区域 使用计算器，你可以在输入表达式时，加入标准的数学操作符和函数调用，你只需点击序列数字，运算符或者函数等按钮，就可把函数加入表达式中，点击INV按钮，可轮流改变部分按钮的函数功能。 ?按钮“（”与“）”按钮，成对使用圆括号强制改变表达式中的运算顺序。 ?MAX/MIN按钮：查找变量中最大值/变量中最小值。 ?COMPLEX/CONJUGATE按钮：形成一个复变量/对一个复数变量执行共轭运算，利用INV 按钮进行函数功能切换。 ?LN/e^X按钮：求一个变量的自然对数/求变量的e次幂，利用INV按钮进行函数功能切换。?STO/RCL按钮：将表达式区域信息存储在内存中/从内存中恢复重复使用的表达式，利用

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT

ansys命令流

第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：

ansys实例命令流-谱分析命令流

/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。

MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF

几个ansys经典实例(长见识)

平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end

ANSYS建模apdl命令流实例应用

大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。

ansys加载函数载荷func

Apply/Functions/Define/Edit打开函数编辑器 ●Functions Type：选择函数类型。选择单个方程或多值函数。如果选择后者，必须键入状态变量名，也就是管理函数中方程的变量。当选择一个多值函数时，六状态表格将被激活。 ●Degrees/Radians：选择度或弧度，这一选择仅决定方程如何被运算，而不会影响*AFUN设置。 ●使用初始变量方程和键区定义结果方程（单个方程）或描述状态变量的方程（多值函数），出如果定义单方程函数，保存方程。如果是定义多值函数，则继续下面的步骤。 ●单击Regime1，键入在函数表格下定义的状态变量的相应的最大最小值限制。 ●定义这个状态的方程。 ●单击Regime2，注意状态变量的最小值限制已被定义并且不可更改，这一特征确保状态保持连续而无间隙。定义这个状态的最高值限制。 ●定义这个状态的方程。 ●在六个状态中连续如上操作。在每个状态里，不必储存或保存单个方程，除非想在另一状态中重用某个方程。 ●输入一个注释描述函数（可选）选择File/Comments。

●计算器区域 使用计算器，你可以在输入表达式时，加入标准的数学操作符和函数调用，你只需点击序列数字，运算符或者函数等按钮，就可把函数加入表达式中，点击INV按钮，可轮流改变部分按钮的函数功能。?按钮“（”与“）”按钮，成对使用圆括号强制改变表达式中的运算顺序。 ?MAX/MIN按钮：查找变量中最大值/变量中最小值。 ?COMPLEX/CONJUGATE按钮：形成一个复变量/对一个复数变量执行共轭运算，利用INV按钮进行函数功能切换。 ?LN/e^X按钮：求一个变量的自然对数/求变量的e次幂，利用INV按钮进行函数功能切换。 ?STO/RCL按钮：将表达式区域信息存储在内存中/从内存中恢复重复使用的表达式，利用INV 按钮进行功能切换。 ?CVAR按钮：计算两个变量之间的协方差（covariance），只适用于PSD求解。 ?RPSD按钮：计算响应PSD，只适用于PSD求解。 ?RESP按钮：计算响应功率谱密度，只适用于PSD求解。 ?LOG按钮：求一个变量的普通对数。 ?ABS/INS MEM按钮：求实变量的绝对值或复变量的模/将内存区域的内容插入到表达式中，利用I

ansys实例命令流-实体梁分析命令流

/FILNAME,SolidBeam ,1 !定义工作文件名。/TITLE,SolidBeam Analysis !定义工作标题。/PREP7 ET,1,SOLID95 !定义材料属性。 MPTEMP,,,,,,,, !定义材料属性。MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 !建立几何模型 K,1,,,, K,2,450,,, K,3,450,-55,, K,4,,-100,, FLST,2,4,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 A,P51X VOFFST,1,45, , !网格划分。 FLST,5,4,4,ORDE,2 FITEM,5,9 FITEM,5,-12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,4,4,ORDE,4

FITEM,5,2 FITEM,5,4 FITEM,5,6 FITEM,5,8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,1 FITEM,5,3 FITEM,5,5 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,30, , , , ,1 !* CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加载。

ANSYS_结构稳态(静力)分析之经典实例-命令流格式

/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !

命令流ansys经典实例

ansys钢筋混凝土建模实例 finish /clear /prep7 et,1,solid65,,,,,,,1 et,2,link8 et,3,185 et,4,solid45 !************************定义材料属性***************************** !混凝土材料属性 mp,ex,1,1.596e10 mp,prxy,1,0.2 mp,dens,1,2400 fc=1.68e7 !c30混凝土轴心抗压强度设计值 ft=1.86e6 tb,concr,1 tbdata,,0.5,0.95,ft,-1 tb,miso,1,,11 tbpt,,0.0002,fc*0.19 tbpt,,0.0004,fc*0.36 tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64 tbpt,,0.0010,fc*0.75 tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91 tbpt,,0.0016,fc*0.96 tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc tbpt,,0.0033,fc !57到68 v,6,7,56,54,10,11,60,58 vgen,4,57,57,1,0,0.38,0 vgen,2,57,60,1,0.25,0,0 vgen,2,57,60,1,-0.25,0,0 !69到80 v,189,190,234,233,191,192,238,237 vgen,4,69,69,1,0,0.38,0

vgen,3,69,72,1,0,0,0.25 !121到136 vgen,2,105,120,1,0,0,0.107 !193到208 k,1178,0.81,0,0.107 k,1179,0.823,0,0.107 k,1180,0.81,0,0.417 k,1181,0.81,0.16,0.107 k,1182,0.823,0.16,0.107 k,1183,0.81,0.16,0.417 k,1184,0.81,0.38,0.107 k,1185,0.823,0.38,0.107 k,1186,0.81,0.38,0.417 v,1178,1179,187,1180,1181,1182,189,1183 v,1181,1182,189,1183,1184,1185,191,1186 vgen,4,193,194,1,0,0.38,0 vgen,2,193,200,1,0.107,0,0 !209到232 k,1271,-0.02,0,-0.013 k,1272,-0.02,0,0 k,1273,-0.02,0.16,-0.013 k,1274,-0.02,0.16,0 k,1275,-0.02,0.38,-0.013 k,1276,-0.02,0.38,0 v,677,1,1272,1271,679,5,1274,1273 v,679,5,1274,1273,681,9,1276,1275 vgen,4,209,210,1,0,0.38,0 vgen,3,209,216,1,0.25,0,0 !233到256 vgen,2,209,232,1,0,0,0.107 vgen,2,257,280,1,0.107,0,0 k,5001,-0.1,1.3,0 k,5002,-0.02,1.3,0 k,5003,-0.02,1.3,0.094

ANSYS中加载与求解中的难点和陌生点

ANSYS中加载与求解中的难点和陌生点 1：对于加载的流程可以再有限元模型上加载也可以在几何模型上，但是是有区别的，几何模型的荷载独立于即将成型的有限元模型，可以自由修改网格，但是有限元模型上的荷载不行，必须先删除荷载，再修改网格，然后再重新输入荷载，原因很简单，如果你在有限元模型上建立荷载，那么已经网格划分完毕的节点被荷载依附，这样如果重新修改网格，有限元网格节点发生变化，那么荷载也就失去了整体性依附的特质，所以显然无法直接修改网格。 2：自由度约束如果是施加在线，面上，注意不要误以为施加在几何模型上，而是施加在有限元模型上，其实是施加在相应图素的节点上，所以必须先划分网格形成有限元模型以后才能施加对应的线面约束，注意是不是说约束就只能在有限元模型上施加呢，不是，因为还有关键点施加自由度约束，所以两种模型施加还是成立的，由此可见，施加自由度约束的对象是点是没有变的，关键看是关键点还是节点，节点按是不是一个图素的所有节点。 3：对于施加荷载或者删除荷载或者建立生死单元等等一系列针对对象的操作必须在选择之后完成，求解器进入前或者进入后处理前必须完成选择和FINISH命令。 4：施加集中荷载注意是针对节点坐标系的，只有集中荷载时，NTROTA，一般默认的是总体坐标系，这里注意整体坐标系和局部坐标系，柱坐标系，直角坐标系，球坐标系的区别。前面两者是两类，两类中又可以分别包含后三类，后面一类还多加环坐标系一类。 5：施加面荷载是难点，因为ANSYS中的面荷载时包括线荷载的。所以针对不同的单元类型，面荷载的单位和加载命令式不同的，比如对于2D单元不能使用SFA命令，对于3D单元不能使用SFL命令，对于梁单元施加单元荷载只能使用唯一的加载命令SFBEAM，另外对于SFE命令非常容易出错，因为要在对应的面号上进行施加，所以一般我们采用SFL，SFA命令，单位也不同，比如2D单元的面荷载是除以面积，壳单元的SFL是除以长度。 6：荷载步，子步，时间，平衡迭代是难点。在求解控制中最难的就是分清楚何种分析类型，注意几个范畴的区别，比如静力分析，动力分析，非线性分析等区别和命令的差异。其实分析的类型包括静力分析，瞬态分析，模态分析，谐分析，谱分析，屈曲分析，子结构分析等等。而非线性和线性分析师隶属于其各分析类型下的，只能说有的分析类型只能进行线性的分析。比如谐分析。对于静力线性分析比较简单，如果有非线性的因素有些命令必须掌握，