Ansys材料参数的定义问题
- 格式:docx
- 大小:14.82 KB
- 文档页数:2
ansys警告和错误(持续更新)(2010-12-20 11:05:16)1、The value of UY at node 1195 is 449810067.It is greater than the current limit of 1000000.This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model si properly constrained.错误的可能:1).出现了刚体位移,要增加约束2).求解之前先merge或者压缩一下节点3).有没有接触,如果接触定义不当,也会出现这样类似的情况4)材料属性设置不对会出现这种情况,例如密度设置的太离谱;2、Large negative pivot value...May be because of a bad temperature-dependent material property used in the model.出现这个错误很可能的原因是约束不够!请仔细检查模型!3、开始求解后出现以下提示,Solid model data is contaminated后来终于找到原因了有限元网格里包含一些未被划分网格的线,一般来说出现在面于面之间有重合的线,导致虽然面被划分了网格,却包含未被划分网格的线。
解决办法,把模型存为.cdb格式(去掉几何信息),然后再读取,就可以求解了命令:cdwrite,db,模型名,cdb听起来不错,不过也没遇到过,一般在划分后用一下NUMMRG命令,合并元素,以避免这种情况出现4、*** WARNING ***There are 79 small equation solver pivot terms.几个可能:1) 约束不够,但警告有79 个方程出现小主元,这一条可能性较小,但也不妨检查一下。
ansys分析出现问题NO.0052some contact elements overlap with the other contact element which can cause over constraint. 这是由于在同一实体上,即有绑定接触(MPC)的定义,又有刚性区或远场载荷(MPC)的定义,操作中注意在定义刚性区或远场载荷时避免选择不必要的DOF自由度,以消除过约束NO.0053Shape testing revealed that 450 of the 1500 new or modified elements violate shape warning limits.是什么原因造成的呢?单元网格质量不够好尽量,用规则化网格,或者再较为细密一点NO.0054在用Area Fillet对两空间曲面进行倒角时出现以下错误:Area 6 offset could not fully converge to offset distance 10. Maximum error between the two surfaces is 1% of offset distance.请问这是什么错误?怎么解决?其中一个是圆柱接管表面,一个是碟形封头表面。
ansys的布尔操作能力比较弱。
如果一定要在ansys里面做的话,那么你试试看先对线进行倒角,然后由倒角后的线形成倒角的面。
建议最好用UG、PRO/E这类软件生成实体模型然后导入到ansysNO.0055There are 21 small equation solver pivot terms.;SOLID45 wedges are recommended only in regions of relatively lowstress gradients.第一个问题我自己觉得是在建立contact时出现的错误,但自己还没有改正过来;第二个也不知道是什么原因。
Ansys workbench常用材料属性
1. isotropic secant coefficient of expansion 各向同性的热胀系数
需要输入基准温度、热膨胀系数。
基准温度,默认22度热膨胀系数
2. orthotropic secant coefficient of expansion 各向异性的热胀系数
需要输入基准温度、三个方向的热膨胀系数。
3. isotropic instantaneous coefficient of expansion 各向同性的热胀系数(随温度变化)需要输入基准温度、热膨胀系数。
(随温度变化)
4. orthotropic instantaneous coefficient of expansion 各向异性的热胀系数(随温度变化)需要输入基准温度、三个方向的热膨胀系数。
(随温度变化)
5. 阻尼系数、质量阻尼、刚度阻尼
6.Isotropic elasticity 各项同性的线弹性材料
需要输入弹性模量与泊松比
7.orthotropic elasticity 各项异性的线弹性材料
需要输入各方向的弹性模量与泊松比
8 Bilinear isotropic/kinematic hardening 双线性材料(非线性材料)需要输入屈服强度及切向模量,需要配合isotropic elasticity使用。
9.multilinear isotropic/kinematic hardening 多线性材料(非线性材料,应力应变曲线)需要配合isotropic elasticity使用,输入应力应变曲线。
1、Real constant 1 referenced by at least element types 1 and 2.原因:单元1和2共用了实常数1。
2、STURM number= 15 should be 0 Probably initial shift greater than first mode or Final Mode(s) is in a cluster.原因:所建模型刚度太小,你可以看一下你所求的频率几乎全部接近于0。
应该修改模型。
3、The calculated reference FORCE CONVERGENCE VALUE = 1.507996765E-04 is less than 1E-2 a threshold zero, A value of 1E-2 or specified MINREF is used. Check results carefully.原因:收敛准则没有定义。
4、The program chosen initial timestep/load-factor is arbitrary. It is necessary for the user to supply a suitable initi timestep/load-factor through the NSUB or DELTIM command for convergence and overall efficiency.原因:没有定义荷载子步。
子空间(Subspace)模态提取法的输出如果采用子空间模态提取法,那么在输出内容中可能会看到这样的警告:“STURM number=n should be m”,n和m是整数,表示某阶模态被漏掉了,或者第m阶和第n阶模态的频率相同而要求输出的只有第m阶模态。
你可以用下面的两个方法来检查被漏掉的模态:使用更多的迭代向量和改变特征值提取法的漂移点数值。
ANSYS中非线性材料的定义ANSYS中非线性材料的定义ANSYS中定义材料非线性包括如下步骤:1.定义材料的弹性模量(MP或MPDATA命令);2.激活非线性材料属性表并定义(TB族命令,包括:TB + TBTEMP + TBDATA或TBPT等)即:ANSYS中材料非线性定义命令流:1.定义材料的弹性模量:①MP, ! (该命令中应含有材料号)或①MPTEMP,MPDATA, ! (该命令中应含有材料号, MPTEMP+MPDATA是连续的)2. 激活非线性材料属性表并定义②TB, ! (该命令中应含有材料号, 三个命令是连续的)TBTEMP,TBDATA,或②TB, ! (该命令中应含有材料号, 三个命令是连续的)TBTEMP,TBPT,详述如下:1. 利用MP或MPDATA命令定义材料的弹性模量MP,Lab,MAT,C0,C1,C2,C3,C4说明:定义材料的属性(Material Property),材料属性为固定值时,其值为C0,当随温度变化时,由后四个参数控制。
MAT:对应ET所定义的号码(ITYPE),表示该组属性属于ITYPE。
Lab:材料属性类别,任何元素具备何种属性在元素属性表中均有说明。
例如:杨氏系数(Lab=EX,EY,EZ),密度(Lab=DENS),泊松比(Lab=NUXY,NUXYZ,NUZX),剪切模数(Lab=GXY,GYZ,GXZ),热膨胀系数(Lab=ALPX,ALPY,ALPZ)等。
2. 利用TB命令激活非线性材料属性表,并利用TBTEMP及TBDATA或TBPT命令定义属性表中数据TB, Lab, mat, ntemp,npts,tbopt,eosopt 激活非线性材料特性表的定义TBTEMP,temp,kmod 为材料表定义温度值(每一个温度对应一个材料非线性公式或应力-应变曲线)TBDATA, stloc, c1,c2,c3,c4,c5,c6 给当前数据表定义数据或TBPT,oper, x,y 在应力-应变曲线上定义一个点上述两个命令要配合TB及TBTEMP使用。
ansys常见错误ansys分析出现问题NO.0052some contact elements overlap with the other contact element which can cause over constraint. 这是由于在同一实体上,即有绑定接触(MPC)的定义,又有刚性区或远场载荷(MPC)的定义,操作中注意在定义刚性区或远场载荷时避免选择不必要的DOF自由度,以消除过约束NO.0053Shape testing revealed that 450 of the 1500 new or modified elements violate shape warning limits.是什么原因造成的呢?单元网格质量不够好尽量,用规则化网格,或者再较为细密一点NO.0054在用Area Fillet对两空间曲面进行倒角时出现以下错误:Area 6 offset could not fully converge to offset distance 10. Maximum error between the two surfaces is 1% of offset distance.请问这是什么错误?怎么解决?其中一个是圆柱接管表面,一个是碟形封头表面。
ansys的布尔操作能力比较弱。
如果一定要在ansys里面做的话,那么你试试看先对线进行倒角,然后由倒角后的线形成倒角的面。
建议最好用UG、PRO/E这类软件生成实体模型然后导入到ansysNO.0055There are 21 small equation solver pivot terms.; SOLID45 wedges are recommended only in regions of relatively lowstress gradients.第一个问题我自己觉得是在建立contact时出现的错误,但自己还没有改正过来;第二个也不知道是什么原因。
ansys警告和错误(持续更新)(2010-12-20 11:05:16)1、The value of UY at node 1195 is 449810067.It is greater than the current limit of 1000000.This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model si properly constrained.错误的可能:1).出现了刚体位移,要增加约束2).求解之前先merge或者压缩一下节点3).有没有接触,如果接触定义不当,也会出现这样类似的情况4)材料属性设置不对会出现这种情况,例如密度设置的太离谱;2、Large negative pivot value...May be because of a bad temperature-dependent material property used in the model.出现这个错误很可能的原因是约束不够!请仔细检查模型!3、开始求解后出现以下提示,Solid model data is contaminated后来终于找到原因了有限元网格里包含一些未被划分网格的线,一般来说出现在面于面之间有重合的线,导致虽然面被划分了网格,却包含未被划分网格的线。
解决办法,把模型存为.cdb格式(去掉几何信息),然后再读取,就可以求解了命令:cdwrite,db,模型名,cdb听起来不错,不过也没遇到过,一般在划分后用一下NUMMRG命令,合并元素,以避免这种情况出现4、*** WARNING ***There are 79 small equation solver pivot terms.几个可能:1) 约束不够,但警告有79 个方程出现小主元,这一条可能性较小,但也不妨检查一下。
1怎么在ansys中查看内部的应力分布最佳答案可以通过工作平面切片法看在PlotCtrls 下,选Style,选Hidden Lione Option,在Type of plot后选Capped hidden ,在cutting plane is后选Normal to view 确定后转动模型就可以看到相应截面的应力了。
也可以选其他选项有不同的切面法可自己尝试2 ANSYS中如何查看应力对一个受力情况进行了静力学分析然后产生的偏移现在要查看他的应力和切应力情况,怎么看?GUI操作最佳答案GUI-general postproc-plot results-contour plot下面有四个选项,根据你的具体情况选择是要节点还要是单元3如何查看ANSYS模型受力之后的应力应变曲线图比如说我想知道某点的,就是用应力表示纵坐标,应变表示横坐标那种曲线图,最佳答案我想,你说的应该是应力沿着座标的分布关系吧?进入后处理程序以后,建立映射路径,把应力数据映射到路径上去,可以得到应力沿着路径的分布曲线。
至于你说的“应力应变曲线”,严格地说,它是材料的本构关系,不是ansys计算出来的,而是由材料力学性能试验机所测得的。
4 ansys,模态分析完成后层间位移如何提取最佳答案建议看看结构动力学方面的基础书籍,那个位移数值其实就是位移,只不过这个位移不是结构的真实位移,而是结构对应于各阶频率的模态位移。
它反映的是结构的振动特性,比如,结构的一阶频率1.4HZ,一阶振动形态是个竖向的正弦半波。
那么这就是简支梁的典型振动特征。
值得注意的是,这个模态位移对于各阶频率来说,也只是个相对大小。
ANSYS默认的是关于质量矩阵归一化的模态。
而如何由得到的模态位移求结构在动力荷载下的真实位移,要用到模态叠加的方法去求。
模态分析只是结构动力分析的第一步,接下来还要继续分析,求出所谓的各阶的振型坐标。
才能最终通过线形叠加的方法得到结构在荷载下的响应。
【转】ANSYS实常数和材料定义总结2010-07-08 20:46定义实常数实常数用于描述那些用单元几何形状不能完全确定的几何参数。
壳单元通过四边形和三角形定义了壳的表面,实常数用来定义其厚度;而梁单元的实常数相对复杂。
主要包括截面积、截面对zz轴、yy轴的惯性短、沿z轴、y 轴的厚度(最大应力发生在离轴最远点)等。
对于简单截面梁,其几何特性这里不再赘述。
但对于实体结构复杂的复合梁,其截面特性的定义具有技巧。
在有限元建模过程中,为简化结构,减少单元数量,通常将其简化为单根梁。
如下图所示结构,经过受力分析可知,主要承力构件为4根立柱,其余斜杆只是起辅助支撑作用,因此其截面应简化如右图所示。
但是,经过计算会发现,计算结果数据中位移和应力明显偏小,与实际情况有出入。
经过分析不难发现,造成这种情况的原因是截面的选择只考虑了截面积和惯性矩,忽视了梁单元的质量,从而造成重力变形减小。
解决这个问题,不能简单增大截面积,那样会使计算应力不可信。
我们可以采取2种方法:(1)沿梁轴线均匀加载一个沿重力方向的线性载荷;(2)将梁单元材料密度乘一个系数。
上述2种方法均切实可行,也得到了工程实践的验证。
单元的材料特性定义绝大多数单元类型都需要材料特性。
根据应用的不同,材料特性可以是线性或非线性。
与单元类型、实常数一样,ANSYS软件对每一组材料特性有一个材料参考号。
但值得注意的是,材料库中的特性值是为了方便而提供的,这些数值是材料的典型值,供用户进行基本分析及一般应用场合,特殊情况用户应自己输人数据。
线性材料特性可以是常数或温度相关的,各向同性或正交异性的,对各向同性材料只需指定其一个方向的特性。
非线性材料特性通常是表格数据,如塑性数据、磁场数据、蛹变数据、膨胀数据、超弹性材料数据等。
材料特性主要由材料本身物理特性决定,在此不再赞述。
Ansys材料参数的定义问题参数的定义正确与否ANSYS计算结果的精度,不仅与模型,⽹格,算法紧密相关,⽽且材料参数对结果的可靠性也有决定性的作⽤,为⽅便⼤家的学习,本⼈就⽤过的⼀些材料模型,作出⼀些总结,并给出相关的命令操作,希望对从事ANSYS应⽤的兄弟姐妹们有所帮助,⽔平有限,不对之处还望及时纠正.先给出线性材料的定义问题,线性材料分为三类:1.isotropic:各向同性材料2.orthotropic:正交各向异性材料3.anisotropic:各向异性材料1. isotropic各向同性材料的定义:这种材料⽐较普遍,⽽且定义也⾮常简单,只需定义两个常数:EX, NUXYNUXY默认为0.3,剪切模量GXY默认为EX/(2(1+NUXY)),如果你定义的是各向同性的弹参数⼀般不⽤定义.如果要定义,⼀定要和公式: EX/(2(1+NUXY))的值匹性材料的话,这个参数配,否则出错,另泊松⽐的定义⼀般推荐不要超过0.5.相关命令,例如:mp,ex,1,300e9mp,nuxy,1,0.252.orthotropic:正交各向异性材料:这种材料也是⽐较常见的,不过定义起来稍微⿇烦⼀点,需定义的常数有: EX, EY, EZ, NUXY, NUYZ, NUXZ, GXY, GYZ, GXZ参数不定义的话,程序会提⽰出错,⽐如:XY注意:在这⾥没有默认值,就是说,如果你某些参数平⾯的平⾯应⼒问题,如果你只定义了EX, EY,程序将提⽰你,这是正交各向异性材料, GXY, NUXY是必须的.相关命令,例如:mp,ex,1,300e9mp,ey,1,200e9mp,nuxy,1,0.25mp,gxy,1,170e9…3.anisotropic:各向异性材料:材料定义起来较为复杂,这⾥我只作些简单的说明,更详细的资料,⼤家可以去各向异性材料定义看帮助.对于各向异性弹性材料的定义,需要定义弹性系数矩阵,这个矩阵是⼀个对称正定阵,因⽽输⼊的值⼀定要为正值.弹性常数矩阵如下图所⽰,各向异性体只有21个独⽴的弹性常数,因⽽我们也就只需输⼊参数即可,⽽且对于⼆维问题,弹性常数缩减为10个.弹性系数矩阵可以⽤刚度或柔21个参数度两种形式来定义,⾃⼰根据情况选⽤,输⼊的时候,可以通过菜单或者TB命令的TBOPT选项来控制.相关的命令流,例如:tb,anel,1tbdata,1, 110e6, 120e6, 130e6, 140e6, 150e6, 160e6tbdata,7, 220e6, 230e6, 240e6, 250e6, 260e6tbdata,12, 330e6, 340e6, 350e6, 360e6tbdata,16, 440e6, 450e6, 460e6tbdata,19, 550e6, 560e6tbdata,21, 660e6另:需注意⼀下各个参数参数的编号顺序和起始位置,不要搞错了,输⼊的时候,是按照上三⾓阵来录⼊的,即:D11,D12,D13,D14,D15,D16,D22,D23…MPMP, Lab, MAT, C0, C1, C2, C3, C4Defines a linear material property as a constant or a function of temperature.PREP7: MaterialsMP ME ST PR PRN DS DSS FL EM EH DY PP <> EME MFSProduct RestrictionsLabValid material property label. Applicable labels are listed under "Material Properties" in the input table for each element type in the Element Reference. See Linear Material Properties of the Element Reference for more complete property label definitions:EXElastic moduli (also EY, EZ).—ALPXSecant coefficients of thermal expansion (also ALPY, ALPZ).—CTEXInstantaneous coefficients of thermal expansion (also CTEY, CTEZ).—THSXThermal strain (also THSY, THSZ).—REFTReference temperature. Must be defined as a constant; C1 through C4 are ignored.—PRXYMajor Poisson's ratios (also PRYZ, PRXZ).—NUXYMinor Poisson's ratios (also NUYZ, NUXZ).—GXYShear moduli (also GYZ, GXZ).—DAMPK matrix multiplier for damping.—Note:If used in an explicit dynamic analysis, the value corresponds to the percentage ofdamping in the high frequency domain. For example, 0.1 roughly corresponds to 10%damping in the high frequency domain.DMPRConstant material damping coefficient.—MUCoefficient of friction.—DENSMass density.—CSpecific heat.—ENTHEnthalpy.—KXXThermal conductivities (also KYY, KZZ).—HFConvection or film coefficient.—EMISEmissivity.—QRATEHeat generation rate.—VISCViscosity.—SONCSonic velocity.—RSVXElectrical resistivities (also RSVY, RSVZ).—Electric relative permittivities (also PERY, PERZ).Electric relative permittivities (also PERY, PERZ).PERX—Note:If you enter permittivity values less than 1 for SOLID5, PLANE13, or SOLID98, theprogram interprets the values as absolute permittivity. Values input for PLANE223,SOLID226, or SOLID227 are always interpreted as relative permittivity.MURXMagnetic relative permeabilities (also MURY, MURZ).—MGXXMagnetic coercive forces (also MGYY, MGZZ).—LSSMMagnetic loss tangent.—LSSTElectric loss tangent.—SBKXSeebeck coefficients (also SBKY, SBKZ).—MATMaterial reference number to be associated with the elements (defaults to the current MAT setting [MAT]).C0Material property value, or if a property-versus-temperature polynomial is being defined, the constant term in the polynomial. C0 can also be a table name (%tabname%); if C0 is a table name, C1 through C4 are ignored.C1, C2, C3, C4Coefficients of the linear, quadratic, cubic, and quartic terms, respectively, in the property-versus-temperature polynomial. Leave blank (or set to zero) for a constant material property.。
ANSYS 查询函数(Inquiry Function)在ANSYS操作过程或条件语句中,常常需要知道有关模型的许多参数值,如选择集中的单元数、节点数,最大节点号等。
此时,一般可通过*GET命令来获得这些参数。
现在,对于此类问题,我们有了一个更为方便的选择,那就是查询函数— Inquiry Function。
Inquiry Function类似于ANSYS的 *GET 命令,它访问ANSYS数据库并返回要查询的数值,方便后续使用。
ANSYS每执行一次查询函数,便查询一次数据库,并用查询值替代该查询函数。
假如你想获得当前所选择的单元数,并把它作为*DO循环的上界。
传统的方法是使用*GET命令来获得所选择的单元数并把它赋给一个变量,则此变量可以作为*DO循环的上界来确定循环的次数*get, ELMAX,elem,,count*do, I, 1, ELMAX……*enddo现在你可以使用查询函数来完成这件事,把查询函数直接放在*DO循环内,它就可以提供所选择的单元数*do, I, ELMIQR(0,13)……*enddo这里的ELMIQR并不是一个数组,而是一个查询函数,它返回的是现在所选择的单元数。
括弧内的数是用来确定查询函数的返回值的。
第一个数是用来标识你所想查询的特定实体(如单元、节点、线、面号等等),括弧内的第二个数是用来确定查询函数返回值的类型的(如选择状态、实体数量等)。
同本例一样,通常查询函数有两个变量,但也有一些查询函数只有一个变量,而有的却有三个变量。
查询函数的种类和数量很多,下面是一些常用、方便而快速快捷的查询函数1 AREA—arinqr(areaid,key)areaid—查询的面,对于key=12,13,14可取为0;key—标识关于areaidr的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数=-1,材料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数…arinqr(areaid,key)的返回值对于key=1=0, areaid未定义=-1,areaid未被选择=1, areaid被选择…2 KEYPOINTS—kpinqr(kpid,key)kpid—查询的关键点,对于key=12,13,14为0 key —标识关于kpid的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数目=-1,数料号=-2,单元类型=-3,实常数=-4,节点数,如果已分网=-7,单元数,如果已分网kpinqr(kpid,key)的返回值对于key=1=-1,未选择=0,未定义=1,选择3 LINE—lsinqr(lsid,key)lsid—查询的线段,对于key=12,13,14为0 key—标识关于lsid的返回信息=1,选择状态=2, 长度=12,定义的数目=13,选择的数目=14,定义的最大数=-1,材料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数…4 NODE—ndinqr(node,key)node—节点号,对于key=12,13,14为0 key—标识关于node的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数=-2,超单元标记=-3,主自由度=-4,激活的自由度=-5,附着的实体模型ndinqr(node,key)的返回值对于key=1=-1,未选择=0,未定义=1,选择5 VOLUMES—vlinqr(vnmi,key)vnmi—查询的体,对于key=12,13,14为0key—标识关于vnmi的返回信息=1,选择状态=12,定义的数目=13,选择的数目=14,定义的最大数目=-1,数料号=-2,单元类型=-3,实常数=-4,节点数=-6,单元数=-8,单元形状=-9,中节点单元=-10,单元坐标系vlinqr(vnmi,key)的返回值对于key=1=-1,未选择=0,未定义=1,选择ANSYS能实现直接流-固耦合分析吗?ANSYS流固耦合分析有三种形式,可以实现全直接或半直接耦合分析:一: ANSYS/Mechanical模块或含该模块的软件包中的流固耦合分析功能,但此处的流体是非流动的流体,而是静流体,它计算流体由于重力、惯性力、波动压力等引起的分布压力载荷与结构的相互作用。
ansys警告和错误(持续更新)(2010-12-20 11:05:16)1、The value of UY at node 1195 is 449810067.It is greater than the current limit of 1000000.This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model si properly constrained.错误的可能:1).出现了刚体位移,要增加约束2).求解之前先merge或者压缩一下节点3).有没有接触,如果接触定义不当,也会出现这样类似的情况4)材料属性设置不对会出现这种情况,例如密度设置的太离谱;2、Large negative pivot value...May be because of a bad temperature-dependent material property used in the model.出现这个错误很可能的原因是约束不够!请仔细检查模型!3、开始求解后出现以下提示,Solid model data is contaminated后来终于找到原因了有限元网格里包含一些未被划分网格的线,一般来说出现在面于面之间有重合的线,导致虽然面被划分了网格,却包含未被划分网格的线。
解决办法,把模型存为.cdb格式(去掉几何信息),然后再读取,就可以求解了命令:cdwrite,db,模型名,cdb听起来不错,不过也没遇到过,一般在划分后用一下NUMMRG命令,合并元素,以避免这种情况出现4、*** WARNING ***There are 79 small equation solver pivot terms.几个可能:1) 约束不够,但警告有79 个方程出现小主元,这一条可能性较小,但也不妨检查一下。
目录引言 (2)一杆件受拉压的内力、应力、变形 (2)1.1轴向拉压的内力、轴力图 (2)1.2 轴向拉压杆横截面上的应力 (5)1.3 轴向拉压杆横截面上的变形 (7)1.4 圣维南原理 (9)1.5 工程结构实例分析 (11)二圆轴扭转 (15)2.1、扭转的力学模型及ANSYS建模 (15)2.2、圆轴扭转时,横截面上的内力偶矩------扭矩 (15)2.3、圆轴扭转时,横截面上的应力、强度条件 (15)(1) 横截面上的切应力 (15)(2) 极惯性矩与抗扭截面系数 (15)三、梁弯曲的内力、变形、应力 (20)3.1 梁的弯曲内力、变形 (20)3.2 弯曲应力 (27)3.3 工程实例: (31)四、压杆稳定 (35)4.1、压杆稳定的概念 (35)4.2、临界压力 (35)4.3、三类压杆的临界载荷 (36)4.4、压杆稳定性计算 (36)4.5 工程实例4 (38)引 言《材料力学》是机械、土木类工科学生重要的技术基础课,其计算方法和思想在工程计算中应用非常广泛。
为了使学生对课内知识体系有一个比较清晰的感性认识,锻炼学生的求真精神和实践动手能力,进一步培养学生的综合创造力,兴趣小组的学生们在教师的指导下基于ANSYS 有限元分析软件对《材料力学》的某些知识点进行数值计算与模拟,得到相关的数据、云图或动画,从而对理论公式进行形象验证,更开阔了学生的视野,提高了学生的CAE 水平。
本研究内容包括三部分:(1)对《材料力学》课程中的基本内容,包括拉压、剪切、扭转、弯曲的内力、应力、变形、压杆稳定、动载荷、疲劳强度、圣维南原理等重要理论知识点情况通过ANSYS 进行分析,得到内力、变形、应力、应变相关的数据、云图或动画;(2)对重要知识点的典型例题通过ANSYS 进行计算,并与理论计算结果进行对比验证。
(3)对《材料力学》理论知识能够解决的典型工程实际问题进行建模、分析与计算。
一 杆件受拉压的内力、应力、变形1.1轴向拉压的内力、轴力图在工程结构和机械中,发生轴向拉伸或压缩的构件是很常见的。
材料参数的定义问题
我想用过ANSYS的人都知道:ANSYS计算结果的精度,不仅与模型,网格,算法紧密相关,而且材料参数的定义正确与否对结果的可靠性也有决定性的作用,为方便大家的学习,本人就用过的一些材料模型,作出一些总结,并给出相关的命令操作,希望对从事ANSYS应用的兄弟姐妹们有所帮助,水平有限,不对之处还望及时纠正.
先给出线性材料的定义问题,线性材料分为三类:
1.isotropic:各向同性材料
2.orthotropic:正交各向异性材料
3.anisotropic:各向异性材料
1. isotropic各向同性材料的定义:
这种材料比较普遍,而且定义也非常简单,只需定义两个常数:EX, NUXY
NUXY默认为0.3,剪切模量GXY默认为EX/(2(1+NUXY)),如果你定义的是各向同性的弹性材料的话,这个参数一般不用定义.如果要定义,一定要和公式: EX/(2(1+NUXY))的值匹配,否则出错,另泊松比的定义一般推荐不要超过0.5.
相关命令,例如:
mp,ex,1,300e9
mp,nuxy,1,0.25
2.orthotropic:正交各向异性材料:
这种材料也是比较常见的,不过定义起来稍微麻烦一点,需定义的常数
有: EX, EY, EZ, NUXY, NUYZ, NUXZ, GXY, GYZ, GXZ
注意:在这里没有默认值,就是说,如果你某些参数不定义的话,程序会提示出错,比如:XY平面的平面应力问题,如果你只定义了EX, EY,程序将提示你,这是正交各向异性材料, GXY, NUXY是必须的.
相关命令,例如:
mp,ex,1,300e9
mp,ey,1,200e9
mp,nuxy,1,0.25
mp,gxy,1,170e9
…
3.anisotropic:各向异性材料:
各向异性材料定义起来较为复杂,这里我只作些简单的说明,更详细的资料,大家可以去看帮助.对于各向异性弹性材料的定义,需要定义弹性系数矩阵,这个矩阵是一个对称正定阵,因而输入的值一定要为正值.
弹性常数矩阵如下图所示,各向异性体只有21个独立的弹性常数,因而我们也就只需输入21个参数即可,而且对于二维问题,弹性常数缩减为10个.弹性系数矩阵可以用刚度或柔度两种形式来定义,自己根据情况选用,输入的时候,可以通过菜单或者TB命令的TBOPT选项来控制.
相关的命令流,例如:
tb,anel,1
tbdata,1, 110e6, 120e6, 130e6, 140e6, 150e6, 160e6
tbdata,7, 220e6, 230e6, 240e6, 250e6, 260e6
tbdata,12, 330e6, 340e6, 350e6, 360e6
tbdata,16, 440e6, 450e6, 460e6
tbdata,19, 550e6, 560e6
tbdata,21, 660e6
另:需注意一下各个参数的编号顺序和起始位置,不要搞错了,输入的时候,是按照上三角阵来录入的,即:D11,D12,D13,D14,D15,D16,D22,D23…。