ansys使用函数加载实例1
ansys例题

ansys例题!实例1-1-1FINI/CLE/FILNAME,EX1-1-1/TITLE,EXAMPLE 1-1-1/UNITS,SI !使用国际单位制*AFUN,DEG !指定角度单位制为度!参数设定*SET,LL,0.085 !长*SET,BL,0.022 !宽*SET,RF,0.008 !漏斗半径*SET,AF,0.008 !漏斗处宽*SET,EXX,9.71E10 !弹性模量*SET,vv,0.3 !弹性模量!进入前处理器/PREP7ET,1,PLANE182 !定义单元类型为182号单元MP,EX,1,EXX !定义弹性模量MP,NUXY,1,vv !定义泊松比KEYOPT,1,3,0 !定义平面应力问题!定义应力应变数据点TB,MISO,1,1,25 !定义本构关系MISO TBTEMP,0 !定义温度TBPT,,0.000402,3.90E7 !在应力应变曲线上定义点TBPT,,0.002118,1.96E8TBPT,,0.003137,2.60E8TBPT,,0.004109,3.03E8TBPT,,0.004631,3.18E8TBPT,,0.005293,3.34E8TBPT,,0.006131,3.48E8TBPT,,0.006926,3.58E8TBPT,,0.007748,3.66E8TBPT,,0.008534,3.72E8TBPT,,0.009548,3.78E8TBPT,,0.010442,3.82E8TBPT,,0.011548,3.86E8TBPT,,0.012236,3.88E8TBPT,,0.014749,3.95E8TBPT,,0.016139,3.97E8TBPT,,0.018424,4.02E8TBPT,,0.020001,4.04E8TBPLOT,MISO,1 !图形显示应力应变曲线!建立四分之一几何模型RECTNG,0,RF+AF/2,0,BL/2 !建立矩形CYL4,,RF+AF/2,RF !建立圆面ASBA,1,2 !面相减得到漏斗RECTNG,RF+AF/2,LL/2,0,BL/2 !建立矩形平板AGLUE,1,3 !面粘接!划分网格MSHKEY,1 !打开映射网格开关MSHAPE,0 !使用四边形或者六面体网格/PNUM,LINE,1 !显示线的编号/PSYMB,LDIR,1 !显示线的方向LPLOT !绘制线LESIZE,9,,,6 !设置9号线划分份数为6 LESIZE,8,,,6LESIZE,11,,,12LESIZE,6,,,6LESIZE,7,,,6,3AMAP,2,8,1,7,9 !通过8,1,7,9关键点映射2号面LESIZE,3,,,12LESIZE,5,,,12LESIZE,4,,,6AMESH,1 !划分1号面!通过对称生成整体模型NSYM,X,175,ALL !通过X对称生成节点,节点增量175 ESYM,,175,ALL !生成单元EPLOTNSYM,Y,350,ALLESYM,,350,ALLEPLOTNUMMRG,ALL !合并所有元素NUMCMP,ALL !压缩编号FINI!进入求解器/SOLUNSEL,S,LOC,X,-LL/2 !选择左边界节点D,ALL,UX,0 !施加位移约束NSEL,R,LOC,Y,0 !选择左边界中轴线上节点D,ALL,ALL !全约束OUTRES,ALL,ALL !设定输出所有的结果TIME,1 !指定求解时间NSUBST,10 !指定载荷子步NLGEOM,ON !打开大变形选项NSEL,S,LOC,X,LL/2D,ALL,UX,0.00025 !施加位移载荷ALLSELSOLVEFINI!进入后处理器/POST1 !进入后处理器SET,LIST !显示所有计算的节果SET,1,LAST,1 !读取第一个加载步最后子步的节果PLDISP,1 !同时显示变形前与变形后的形状。
(完整版)疲劳分析的数值计算方法及ANSYS疲劳分析实例

第十四章疲劳分析的数值计算方法及实例第一节引言零件或构件由于交变载荷的反复作用,在它所承受的交变应力尚未达到静强度设计的许用应力情况下就会在零件或构件的局部位置产生疲劳裂纹并扩展、最后突然断裂。
这种现象称为疲劳破坏。
疲劳裂纹的形成和扩展具有很大的隐蔽性而在疲劳断裂时又具有瞬发性,因此疲劳破坏往往会造成极大的经济损失和灾难性后果。
金属的疲劳破坏形式和机理不同与静载破坏,所以零件疲劳强度的设计计算不能为经典的静强度设计计算所替代,属于动强度设计。
随着机车车辆向高速、大功率和轻量化方向的迅速发展,其疲劳强度及其可靠性的要求也越来越高。
近几年随着我国铁路的不断提速,机车、车辆和道轨等铁路设施的疲劳断裂事故不断发生,越来越引起人们的重视。
疲劳强度设计及其研究正在成为我国高速机车车辆设计制造中的一项不可缺少的和重要的工作。
金属疲劳的研究已有近150年的历史,有相当多的学者和工程技术人员进行了大量的研究,得到了许多关于金属疲劳损伤和断裂的理论及有关经验技术。
但是由于疲劳破坏的影响因素多而复杂并且这些因素互相影响又与构件的实际情况密切相关,使得其应用性成果尚远远不能满足工程设计和生产应用的需要。
据统计,至今有约90%的机械零部件的断裂破坏仍然是由直接于疲劳或者间接疲劳而引起的。
因此,在21世纪的今天,尤其是在高速和大功率化的新产品的开发制造中,其疲劳强度或疲劳寿命的设计十分重要,并且往往需要同时进行相应的试验研究和试验验证。
疲劳断裂是因为在零件或构件表层上的高应力或强度比较低弱的部位区域产生疲劳裂纹,并进一步扩展而造成的。
这些危险部位小到几个毫米甚至几十个微米的范围,零件或构件的几何缺口根部、表面缺陷、切削刀痕、碰磕伤痕及材料的内部缺陷等往往是这种危险部位。
因此,提高构件疲劳强度的基本途径主要有两种。
一种是机械设计的方法,主要有优化或改善缺口形状,改进加工工艺工程和质量等手段将危险点的峰值应力降下来;另一种是材料冶金的方法,即用热处理手段将危险点局部区域的疲劳强度提高,或者是提高冶金质量来减少金属基体中的非金属夹杂等材料缺陷等局部薄弱区域。
ANSYS加载命令流

ANSYS加载命令流⼀,在ANSYS⾥施加地震惯性⼒的⽅法在ANSYS⾥做地震分析时,需要对结构施加地震惯性荷载,地震惯性⼒是通过加速度的⽅式输⼊进结构的,然后与结构的质量⼀起形成动⼒计算时的惯性荷载,下⾯说⼀下在ANSYS ⾥施加地震惯性⼒的⽅法。
⾸先,将三个⽅向的地震加速度放到⼀个⽂本⽂件⾥,如accexyz.txt,在这个数据⽂件⾥共放三列数据,每列为⼀个⽅向的地震加速度值,这⾥仅给出数据⽂件中前⼏⾏的数据:-0.227109E-02 -0.209046E+00 0.467072E+01-0.413893E-02 -0.168195E+00 0.261523E+01-0.574753E-02 -0.157890E+00 0.809014E-01-0.731227E-02 -0.152996E+00 0.119975E+01-0.876865E-02 -0.138102E+00 0.130902E+01-0.101067E-01 -0.131582E+00 0.143611E+00 .......................然后,再建⼀个⽂本⽂件⽤来存放三个⽅向的地震加速度时间点,如time.txt,在这个数据⽂件⾥仅⼀列数据,对应于加速度数据⽂件⾥每⼀⾏的时间点,这⾥给出数据⽂件中前⼏⾏数据:0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.600000E-01.......................编写如下的命令流⽂件,并命名为acce.inp*dim,ACCEXYZ,TABLE,2000,3 !01⾏*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02⾏(3e16.6) !03⾏*vread,ACCEXYZ(1,0),time,txt !04⾏(e16.6) !05⾏ACCEXYZ(0,1)=1 !06⾏ACCEXYZ(0,2)=2 !07⾏,同上ACCEXYZ(0,3)=3 !08⾏,同上finish/SOLUANTYPE,transbtime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itimeAUTOTS,0NSUBST,1, , ,1KBC,1acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个⽅向的地震加速度SOLVE*ENDDO最后,在命令窗⼝⾥输⼊/input,acce,inp即可对结构进⾏地震动⼒分析。
ansys施加面荷载

1、如果是线性变化的,可以采用水压方式定义;1、如果可以用函数表示,则可以用函数来定义,2、也可以使用表面效应单元来定义;在ANSYS中如果要在一个面上施加沿某个方向变化的面荷载,需要有两步来完成:这里以一个在圆筒内表面加内水压力的例子进行说明。
第一步,设置面荷载变化规律。
如果面荷载沿Z向变化,后面指定面荷载从Z=100开始变化,并按斜率为-9800进行变化,可用如下语句sfgrad,pres,,z,100,-9800 !也就是准备在高100米的圆柱加内水压力吧第二步,施加面荷载。
在指定的面上施加按第一步设置的面荷载变化规律的面荷载。
SFA,P51X,1,PRES,0这个语句相当于在指定面上施加法向荷载(选圆筒体内表面),在Z=100时荷载值为0,随Z坐标变化荷载值以变化率-9800进行变化,这样在Z=0时荷载值为-9800*100每次用sfgrad进行设置后仅对随后的sfa命令有效,直倒下次再用sfgrad进行设置。
在面上施加荷载后,对模型剖分后可以执行以下命令来查看加的面荷载是否正确/PSF,PRES,NORM,2,0,1 以箭头方式显示面荷载sftran 将面荷载转化到有限元模型上本文摘自《ANSYS工程分析进阶实例》---王呼佳、陈洪军主编,在此对本书作者表示感谢!一般可以通过两种方法施加面荷载,一是在表面上覆盖一层表面效应单元SURF153或SURF154;二是通过apdl语言编程施加。
基本思路如下:人为将面上压力荷载换算成集中力并施加到节点上。
施加集中力时,将合力分解为X,Y,Z方向的分力。
(1)选中所要施加压力的表面,在面上生成一层shell63单元。
(2)对生成的shell63单元,使用循环语句逐步进行以下操作。
(3)得到每个单元的面积及单元中心的X,Y,Z坐标值。
(4)将坐标值代入压力随坐标变化的函数式,得到单元中心点处的压力值,并乘以面积得到单元所受的合力。
(5)将合力平均后,施加到单元的各个节点。
ANSYS算例100个

1.考虑热应力的模态分析2.ansys中抓图的命令流3.生死单元的例子4.正在建设中5.热应变与结构应变6.表面单元的使用7.修改ANSYS的默认工具条8.表格加载函数加载9.特征值屈曲分析10.包辛格效应11.塑性理论12.使用*vwrite和*vread13.ansys数值变量转字符变量14.ansys中保存视图设置15.ansys中最值问题16.apdl命令分类17.ANSYS宏加密18.用路径得到任一点的应力值19.ansys中的函数20.塑性应变实例21.残余应力22.beam4与solid45的连接23.塑性分析算例24.模态缩减法25.离心应力刚化26.ansys中泊松比效应27.比较弹性模量算例28.获得单元积分点的结果29.等向强化与随动强化小算例30.网格密度与计算结果31.ansys工况组合32.谐响应算例33.安全系数云图34.常应变单元35.ANSYS网格装配36.ANSYS子模型技术37.ansys中螺栓联接的模拟38.ANSYS接触过盈分析39.ANSYS板金冲压算例40.ANSYS节点解单元解单元表41.ANSYS单元生死42.过盈装配算例43.ANSYS超单元子结构44.ANSYS点到面热辐射算例45.ANSYS面到面热辐射46.ANSYS蠕变算例47.ANSYS蠕变算例248.ANSYS转子坎贝尔图分析49.ANSYS P单元算例50.ansys路径算例51.循环对称结构的模态分析52.ansys复合材料算例53.ansys复合材料算例254.扭矩和转动惯量对模态的影响55.ANSYS接触反力提取56.ANSYS热分析后处理57.ANSYS热结构耦合算例58.ANSYS瞬态动力学小算例59.ANSYS谐响应动力学小算例60.ANSYS模态动力学小算例61.ANSYS旋转软化应力刚化算例62.link10-gap单元63.循环对称结构的模态分析64.Ansys疲劳算例65.ANSYS材料实验小算例66.ANSYS各向异性材料算例67.ANSYS焊接模拟68.港口起重机有限元实例69.ansys中由网格生成几何模型70.ansys相变分析算例71.ANSYS输出单元刚度矩阵72.ANSYS模拟单摆运动73.ansys点点接触算例74.ANSYS弹性地基梁分析75.ANSYS平面应力plane4276.平面应变plane4277.轴对称plane4278.ANSYS螺旋的建模79.ANSYS弹簧建模80.ANSYS两个DB合并81.Hill_VS_vonmises82.ANSYS将文本数据给数组83.ansys中定义变截面梁84.ansys裂纹模拟85.ANSYS粘塑性86.shell63和shell187.ANSYS超弹性88.ANSYS剪切锁定算例89.ANSYS体积锁定算例90.ANSYS材料多线性与材料非线性91.工程应力与真实应力92.导热系数对热分布的影响93.ANSYS双线性随动/等向强化应94.ANSYS双线性随动/等向强化应95.ANSYS双线性随动/等向强化应96.ANSYS双线性随动/等向强化应97.Chaboche应力控制非循环对98.循环强化: Chaboche+非99.ANSYS数组插值。
ANSYS经典实例汇集(共72个例子)

47.非线性分析考虑刚度退化 .............................................................................................237 48.一个圆形水池的静力分析 .............................................................................................238 49.ANSYS 中混凝土模式预应力模拟的算例......................................................................239 50.悬臂梁受重力作用发生大变形求其固有频率.............................................................241 51.循环对称结构模态分析.................................................................................................243 52.三角平台受谐波载荷作用的结构响应.........................................................................245 53.三角平台受一地震谱激励的应力分布和支反力.........................................................247 54.三角平台受时程载荷作用的应力分布和变形过程.....................................................249 55.经典层合板理论.............................................................................................................251 56.定易圆轨迹的例子.......................................................................................................258 57.模拟门式刚架施工-单元生死...................................................................................258 58.钢筋混凝土整体式模型例子.......................................................................................261 59.在荷载步之间改变材料属性例子...............................................................................263 60.含预应力的特征值屈曲计算.......................................................................................264 61.振型叠加计算及工况组合例子...................................................................................266 62.柱子稳定分析算(预应力,特征值屈曲,初始缺陷)...........................................269 63. module MConcrete !混凝土模板 ............................................................................272 64.混凝土开裂实例.............................................................................................................280 65.螺栓网格划分...............................................................................................................281 66.自由液面的土石坝平面渗流分析...............................................................................282 67.导出刚度矩阵...............................................................................................................286 68.某混凝土拱坝工程施工期及运行期温度场仿真分析...............................................287 69.移动温度荷载计算.......................................................................................................294 70.SHSD 用于壳-实体装配实例 An .................................................................................296 71.ansys 显示-隐式-回弹分析实例 .........................................................................300 72.工况组合的经典例子...................................................................................................315
Ansys_Ls-dyna压力载荷加载方法
ANSYS/LS-DYNA软件压力载荷加载总结隐式分析:1、进入Solution项,找到Pressure项进行施加压力载荷,选取作用面(Areas)上。
2、设置Plotctrls/Symbal/定义载荷显示形式,如Pressure(压力)以Arrow(箭头表示方向),这样设置好,可以帮助判断施加载荷位置及方向是否正确。
3、一旦施加载荷后,压力的方向将以箭头形式出现。
显式分析:1、只有Element形式加载压力载荷。
2、定义压力载荷载体Element组元。
3、定义时间-压力载荷曲线,没有正负之分。
4、设置Plotctrls/Symbal/定义载荷显示形式,如Pressure(压力)以Arrow(箭头表示方向)。
可以帮助判断施加载荷位置及方向是否正确。
一旦施加载荷后,压力的方向将以箭头形式出现。
5、一般要以映射网格划分有限元模型,正常状态下单元有6个面,并且分配了面号。
施加压力载荷:Specify Load/PRES,面号需要通过Arrow方向来以手工方式逐个试出来(因为每个每次模型的方向和顺序是不一样的),选时间和载荷。
6、加载完毕。
隐式-显式分析步骤:1. 求解分析的隐式部分,从而得到预载求解分析的隐式部分(预加载荷)2. 改变现在的文件名进行显式求解部分Utility Menu > File > Change Jobname > Jobname2 (由原来的Jobname1改为Jobname2)…为了防止显式求解的结果覆盖隐式求解的结果3. 将隐式单元改为相应的显式单元Preprocessor > Element Type > Switch Elem Type …> Implic to Explic如果使用了非对应的单元,使用ETCHG, ITE 命令不能自动地将它们转变,而是用EMODIF命令手动将它们转变。
Preprocessor > Move/Modify > -Elements- Modify Attrib >Select elements to be modified > Elem Type –TYPE (STLOC field) > TYPE参考号与显式单元关联若LINK160, BEAM161, and LINK167 单元都需要第三个节点(方向点), 所以如果相应的隐式单元只定义了端点,那么必须增加第三个节点。
ANSYS函数
ANSYS函数⼀、ANSYS 中可以使⽤的数学函数在ANSYS 帮助系统中关于*SET 命令的注释下列出了ANSYS 中可以使⽤的数学函数。所有这些数学函数均可以在ANSYS 环境中使⽤,这些数学函数包括,其中三⾓函数的单位默认是弧度:
ABS(X)——xACOS(X)——)arccos(xASIN(X)——)arcsin(xATAN(X)——)(x arctgATAN2(X,Y)——)(xy arc COS(X) ——)cos(xCOSH(X)——)(x cohEXP(X)——x eGDIS(X,Y) 求以X 为均值,Y 为标准差的⾼斯分布,在使⽤蒙地卡罗法研究随机荷载和随机材料参数时,可以⽤该函数处理计算结果
LOG(X)——)ln(xLOG10(X)——)lg(xMOD(X,Y) 求 X/Y 的余数. 如果 Y=0, 函数值为 0NINT(X) 求最近的整数RAND(X,Y) 取随机数,其中X 是下限, Y 是上限SIGN(X,Y) 取 X 的绝对值并赋予Y 的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,.SIN(X) 正弦SINH(X) 双曲正弦SQRT(X)——xTAN(X)——)(x tg TANH(X) 双曲正切⼆、ansys 内嵌函数1. distnd( i,j) — I,j 两点的距离2. node(x,y,z)—提取距离位置(x,y,z)最近的节点号3. kp(x,y,z)—提取距离位置(x,y,z)最近的关键点号4.基本函数GDIS(x,y) Random sample of a Gaussian (normal) distribution with mean x and standard deviation y.VALCHR(CPARM) Numerical value of CPARM (if CPARM is non-numeric, returns 0.0).CHRVAL(PARM) Character value of numerical parameter PARM. Number of decimal places depends on magnitude.UPCASE(CPARM) Upper case equivalent of CPARM.LWCASE(CPARM) Lower case equivalent of CPARM.5. kx(i) 表⽰关键点i 的x坐标值;同理 ky(i);kz(i)6. nx(i)表⽰节点i 的x坐标值;同理 ny(i);nz(i)7. nsel(k)是节点k在就是1,不在就是0.8. NDNEXT(N) Next higher node number above N in selected set (or zero if none found).9. NELEM(ENUM,NPOS) returns the node number in position NPOS for element ENUM. Node number at position 1,2,... or20 of elementN, where npos is 1,2,...20.
一线工程师总结AnsysWorkbench之Mechanical应用——分析设置
一线工程师总结AnsysWorkbench之Mechanical应用——分析设置对于结构静力学中的简单线性问题,不需要对其进行设置,但是对于复杂的分析需要设置一些控制选项。
分析设置是在Mechanical分析树的Static Structural下的Anslysis Settings细节设置中。
本文主要对载荷步控制、求解器控制、重启控制、非线性控制、输出控制、分析数据管理进行介绍。
1 载荷步控制载荷步控制用于指定求解步数和时间。
在非线性分析时,用于控制时间步长。
载荷步控制也用于创建多载荷步,如螺栓预紧载荷。
1.1 载荷步与子步载荷步、子步和平衡迭代是控制加载求解过程的三个载荷时间历程节点。
1.1.1 载荷步在线性静力学分析或稳态分析中,可以使用不同的载荷步施加不同的载荷组合。
在瞬态分析中,可以将多个载荷步加载到同一加载历程曲线的不同时间点。
注意:载荷可以分步,约束不能分步。
实例1,固定矩形条一端,在另一端分3步加载载荷,第一步只加载100N的力,第二步只加载10000Nm的逆时针扭矩,第三步推力与扭矩共同作用,求每一步的变形。
Step1,设置零件材料,接触关系,网格划分,过程略。
Step2,分析设置,将载荷步设置为3,其余默认。
Step3,设置边界条件,如下图。
载荷默认都是渐增(斜坡)加载的,用一个载荷步将载荷从0增加到设定值。
选中分析树中的Force,在信息窗口中出现了Tabular Data表格和Graph图表,代表了Force的加载历程,在第一步中,力从0渐变到100,并在第二三步中保持。
对于静力学分析,渐增加载与恒定加载计算无区别,本例将力与扭矩都改为恒定加载,在表格第一行将数字改为设定值。
要想Force在第二步不起作用,只需要点击图表的第二步区域或表格对应行,右击选择Activate/Deactive at this step!(在此步激活/取消),此载荷便在第二步中消失。
同样设置Moment载荷,使它在第一步中不起作用。
ansys apdl结构变形函数输出
ansys apdl结构变形函数输出ANSYS APDL是一款应用程序开发语言,它作为ANSYS的核心模块,是专门应用于结构力学分析的软件工具。
我们可以通过ANSYS APDL快速地进行结构变形函数的“输出”,这将帮助我们更好地了解结构的崩溃点、破坏类型和承载能力。
下面,我将从以下几个方面来详细讲解如何使用ANSYS APDL输出结构变形函数:1. 建立模型和加载首先,我们需要在ANSYS中建立您所要分析的结构。
这通常需要使用设计软件来绘制原始图纸,并将其转换为ANSYS中可识别的文件格式。
在建立好模型后,我们需要加载该结构。
在ANSYS中,我们可以通过选择“LOAD”按钮来执行此操作。
2. 定义节点、单元在加载结构后,我们需要定义节点和单元。
节点是我们在结构中使用的点,单元是我们在结构中使用的元素。
为了更好地输出结构变形函数,我们需要确保每个节点和单元都已正确定义,并正确地与加载模型的结构相对应。
这可通过在ANSYS 命令字符串中键入一系列命令来实现。
3. 运行结构分析在定义好节点和单元后,我们可以通过ANSYS菜单中的“求解器”按钮来运行结构分析。
这将在ANSYS中运行一个代码,根据我们输入的模型和参数,计算出结构的强度和变形情况。
4. 输出结构变形函数在分析完结构的强度和变形情况后,我们可以轻松地输出结构变形函数。
在ANSYS中,我们可以通过键入以下命令来执行此操作:*VWRITE,'NODE',NODE(2),'Ux',U(1)其中,“NODE”表示节点的编号,而“Ux”则表示该节点的横向变形量。
5. 结束分析在输出结构变形函数后,我们需要结束分析。
在ANSYS中,我们可以通过选择“结束”按钮或输入“EXIT”来执行此操作。
综上所述,使用ANSYS APDL输出结构变形函数是一项非常棒的技能。
它将使我们能够更好地了解结构强度和变形情况,以便进行更准确的设计和结构优化。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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: