表面效应单元

【讨论】表面效应单元的使用方法

这是一个加表面效应单元的例子，你看一下有没有帮助。

FINISH

/CLEAR

!!!gauss加载

/title, the simulate of laser cladding process

*SET,P,2000， !the power of the laser

*SET,V,0.008 !scanning speed

!*SET, RADIUS,0.003 !radius of the laser beam

*SET, R,0.003

*set,k,0.6 !the coefficient of the absorb

*SET,L,0.05 !the length of the block

/prep7

!熔覆层高5mm=0.005m

!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m

BLOCK,0,0.05,0, -0.0105,-0.01,0,

Vsel,all

vplot

!

/USER, 1

/VIEW, 1, -0.198026943995E-01, 0.382677277552 , 0.923669829831

/ANG, 1, -0.107453019382

/REPLO

!

ET,1,SOLID70

ET,2,Surf152

DOF,TEMP

KEYOPT,2,3,0

KEYOPT,2,4,1

KEYOPT,2,5,0

KEYOPT,2,6,0

KEYOPT,2,7,0

KEYOPT,2,8,4 !根据固体表面温度计算对流系数KEYOPT,2,9,0

!用温度表方式定义材料性能

Mp,kxx,1,82.9

Mp,dens,1,8900

Mp,c,1,471

mptemp,1,20,300,600,750,900,1200,1500,1800,2000 mpdata,hf,1,1,6,50,120,180,200,250,378,700,850

!

!

Mshkey,0 !自由网格化

Mshape,1,3d

type,1

esize,0.002,

Vmesh,all

eplot

!

!

asel,s,loc,y,0 !对流边界条件要施加在5个面的周围asel,a,loc,z,-0.01

asel,a,loc,x,0

asel,a,loc,x,0.05

type,2 !定义属性

amesh,all

esel,s,type,,2 !选择表面效应单元

sf,all,conv,-1,25 !确定对流边界条件，表格加载

!sf,all,conv,200,25 !确定对流边界条件

!!底面也施加轻微的对流

allsel,all

asel,s,loc,y, -0.0105

type,2 !定义属性

amesh,all

esel,s,type,,2

sf,all,conv,3,25

allsel,all

!

nsel,all

/psf,hflux,,2,0 !以箭头方式显示载荷

!

finish

/SOLU

!!!!!

antype,transient,new

TRNOPT,FULL

Nropt,,,off !关闭自适应下降LUMPM,0

Solcontrol,on

Autots,on

Kbc,1

Outres,all,last,

Lnsrch,on

neqit,50 !每一子步中方程的迭代次数限值!load step1,initial conditions 25

time,0.001

deltim,0.001,0.001,0.001

allsel,all

timint,off

tunif,25 !均匀的节点初始温度25，只对分析的第一个子步有效  

solve

EPLOT

!

! load step2………,apply moving heat flux

allsel

nsel,all

t=0

dt=0.5

local,11,1, 0.001+0.005,0,0,0,-90,0

CSYS,11

esel,s,type,,1

nsel,r,loc,x, 0,0.003

nsel,r,loc,y,0,180

nsel,r,loc,z,0

sf,,hflux,7000000

!sf,,hflux,%flux2%

csys,0

t=t+dt

time,t !每一载荷步结束时的时间deltim,0.1,0.05,0.1,on

kbc,1 !设定为阶越载荷

!neqit,25 !每个子步的默认次数autots,on !打开自动时间步长timint,on,therm !时间积分效果TINTP,,,,1

Lnsrch,on

!outres,all,all

outres,all,last

!outres,all,10

allsel

solve

save

这是一个加表面效应单元的例子，你看一下有没有帮助。FINISH

/CLEAR

!!!gauss加载

/title, the simulate of laser cladding process

*SET,P,2000，!the power of the laser

*SET,V,0.008 !scanning speed

!*SET, RADIUS,0.003 !radius of the laser beam

*SET, R,0.003

*set,k,0.6 !the coefficient of the absorb

*SET,L,0.05 !the length of the block

/prep7

!熔覆层高5mm=0.005m

!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m BLOCK,0,0.05,0, -0.0105,-0.01,0,

Vsel,all

vplot

!

/USER, 1

/VIEW, 1, -0.198026943995E-01, 0.382677277552 , 0.923669829831 /ANG, 1, -0.107453019382

/REPLO

!

ET,1,SOLID70

ET,2,Surf152

DOF,TEMP

KEYOPT,2,3,0

KEYOPT,2,4,1

KEYOPT,2,5,0

KEYOPT,2,6,0

KEYOPT,2,7,0

KEYOPT,2,8,4 !根据固体表面温度计算对流系数

KEYOPT,2,9,0

!用温度表方式定义材料性能

Mp,kxx,1,82.9

Mp,dens,1,8900

Mp,c,1,471

mptemp,1,20,300,600,750,900,1200,1500,1800,2000

mpdata,hf,1,1,6,50,120,180,200,250,378,700,850

!

!

Mshkey,0 !自由网格化

Mshape,1,3d

type,1

esize,0.002,

Vmesh,all

eplot

!

!

asel,s,loc,y,0 !对流边界条件要施加在5个面的周围

asel,a,loc,z,-0.01

asel,a,loc,x,0

asel,a,loc,x,0.05

type,2 !定义属性

esel,s,type,,2 !选择表面效应单元

sf,all,conv,-1,25 !确定对流边界条件，表格加载

!sf,all,conv,200,25 !确定对流边界条件

!!底面也施加轻微的对流

allsel,all

asel,s,loc,y, -0.0105

type,2 !定义属性

amesh,all

esel,s,type,,2

sf,all,conv,3,25

allsel,all

!

nsel,all

/psf,hflux,,2,0 !以箭头方式显示载荷

!

finish

/SOLU

!!!!!

antype,transient,new

TRNOPT,FULL

Nropt,,,off !关闭自适应下降

LUMPM,0

Solcontrol,on

Autots,on

Kbc,1

Outres,all,last,

Lnsrch,on

neqit,50 !每一子步中方程的迭代次数限值

!load step1,initial conditions 25

time,0.001

deltim,0.001,0.001,0.001

allsel,all

timint,off

tunif,25 !均匀的节点初始温度25，只对分析的第一个子步有效solve

EPLOT

!

! load step2………,apply moving heat flux

allsel

nsel,all

t=0

dt=0.5

local,11,1, 0.001+0.005,0,0,0,-90,0

esel,s,type,,1

nsel,r,loc,x, 0,0.003

nsel,r,loc,y,0,180

nsel,r,loc,z,0

sf,,hflux,7000000

!sf,,hflux,%flux2%

csys,0

t=t+dt

time,t !每一载荷步结束时的时间

deltim,0.1,0.05,0.1,on

kbc,1 !设定为阶越载荷

!neqit,25 !每个子步的默认次数

autots,on !打开自动时间步长

timint,on,therm !时间积分效果

TINTP,,,,1

Lnsrch,on

!outres,all,all

outres,all,last

!outres,all,10

allsel

solve

save

问题：在热分析中，有的用户经常问，用表面效应单元和不用表面效应单元，计算得到的温度分布差异比较大，为什么？

解答：Ansys的表面效应单元是非常有特色的，主要用于表面载荷的处理。在进行热分析时

，模型的面载荷（如热流密度、热对流等）既可以直接施加在模型的表面，也可以在模型的表面建立表面热效应单元（SURF151、SURF152），将这些面载荷施加在表面效应单元上

。用表面效应单元和不用表面效应单元，计算得到的温度分布差异比较大，主要问题在于对表面效应单元使用不当。在应用表面效应单元，施加面载荷时，要注意单元选项。比如对于SURF152号单元，选项k8是指定施加在SURF152单元的载荷类型，K8=0,表示在这个表面单元上忽略任何面载荷；K8=1，表示在这个单元上施加热流密度；K8=2，表示在这个单元上施加对流载荷。我们要根据实际施加的表面载荷类型对选项K8进行设置。比如，有的用户已经将对流条件直接施加在模型的表面了，但是，为了计算表面的散热和吸热，又在这个表面建立了表面效应单元，并且将该表面效应单元选项K8=2，那么这就不对了，该SURF152单元就将这些边界条件再处理一次，所以计算结果就与实际情况不符合。对这种情况，正确的应该置K8=0或1。模型的某个表面和这个表面建立的表面效应单元可以同时施加表面载荷，程序在处理时将这些载荷按照一定的规则进行叠加，所以建议它们施加不同类型的面载荷，并且通过单元选项按实际情况进行设置。

Q为问A为回复

-------------------------------------------------------------------------------

Q:模态分析得到的结果是不是某个方向上的各阶频率啊我要得到各个方向的一阶频率能做到吗

A模态分析得到的结果是你所选择的自由度内的振型如果需要得到某一个方向内的振型用RUDECED法选择你所希望的主自由度如UX将所求的频率值设置为一即可

Q不过大型复杂结构的振型一般都不是一个方向的应该是一个方向为主其它方向为次如果把所要求的方向设为主自由度会不会跟真是实际结果有出入

Q对于复杂的大型实体其主自由度无法给出此时就无法应用reduce法可不可以改约束还是就认为一次得到前3阶频率就是x,y,z方向的一阶频率啊

A正是因为复杂结构的主自由度不好确定所以我还是倾向于用SUBSPACE和BLOCK LANCZOS法约束不要随便改它应该描述结构真实的状态不能说前3阶频率即为XYZ方向的第一阶频率你可以好好想想所用方程中[K]和[M]的意义

--------------------------------------------------------------------------------

Q轴承是用来支撑要研究的对象把轴承等价为弹簧将其刚度加在对象的节点上这样就要定义刚度的大小在加是否还要定义单元类型然后才能定义其刚度

A我曾经用LINK8单元模拟过根据K=EA/L确定刚度我当时取E=210e9,L=0.1,然后

A=KL/E经过检验是可行的

A1用link单元模拟其刚度k=EA/L通过假定E,A,L的值使之和k相等其密度要02用combin14模拟其刚度阻尼0

-------------------------------------------------------------------------------- QANSYS处理动力问题如土层地震响应或基础振动反应可否从边界施加加速度荷载若可以如何瞬态动力反应分析时间子步为>200?

A如果你所希望求的最高频率的周期为f则ANSYS取时间步长为1/20f为默认值

Q板单元的后处理的显示在单元坐标还是在总体坐标? 比如位移,应力

A应力位移在总体坐标弯矩不一定

A应力位移是显示在结果坐标系的顺便说明一下节点信息显示在显示坐标系

不过结果坐标系默认为全局坐标系而显示坐标系也是默认为全局坐标系

你可以这样修改结果坐标系

Main Menu:General Postprocessor>Options for Output...

你可以这样修改显示坐标系

Utility Menu:Workplane>Change Display CS to

A你可以用ETABLE命令先存储对应于各单元的位移等等再用ANSYS的柱状图或是别的画图软件PLOT出图形你认为如何

Q我现在需分析一个板梁结构板已用SHELL63单元划分好梁我是用板上的一条线划分单元并添加截面而生成的但现在运算时发现板和梁是分开的它们之间互不影响请教各位高手怎样将板和梁合并为一个整体

APreprocessor>Numbering Ctrls>Merge Items里element and node 合并

Q: 用SHELL63单元划分的板结构想加垂直于板面的匀布载荷?

A: 找到如下位置ANSYS HELP->ANSYS ELEMENTS REFERENCE->ELEMENT LIBRARY->SHELL63 ELASTIC SHELL

在INPUT SUMMARY 中有SURFACE LOADS PRESSURES

1

其中face 1 (I-J-K-L)(bottom,in +Z direction)的意思为当Load Key=1时,pressure 加在由节点IJKL确定的面上在单元坐标中以+Z方向为正也可以说垂直指向面IJKL时为正至于均布则将几个值都填成一样的就可以了或者只填主要参数如VALUE at I自己多看看帮助文件做一个例子验证一下

--------------------------------------------------------------------------------

Q: ANSYS可以自己定义输出格式吗?

A: 可以先打开一个文件*cfopen再用*vwrite,label

通过打开一个文件读出数据或以追加的方式写入数据并可控制输入和输出的格式

*cfopen,,

*vread,,,,,

(f12.6)

*vwrite,

(f12.6)

*cfclos,,

--------------------------------------------------------------------------------

Q:作模态分析时就是集中质量单元和刚度加法

A在选用求解方法时例如SUBSPACE和BLOCK LANCZOS会有一个选项[LUMPM]默认为选用一致质量矩阵你可以选为集中质量矩阵

A用mass21通过mesh keypoint来把点划分成mass21单元即可

Q我在原有的静态分析的基础上想继续作一下结构的模态分析应该采取那些步骤啊

我用solution>new analysis-modal,then current LS,出现如下的提示

1.the eigenvalue extraction method must be specified on MODOPT command. SUBSP or LANB are usually recommanded.

2.please expand modes if any downstream MODAL analysis based calculations need to be performed.

这是怎么一回事

A你最好先看看帮助文件中的ANSYS STRUCTURE REFERENCE中的BUCKLING ANALYSIS

如果你对理论不清楚我劝你不要盲目计算不然结果也许并不可信这是我的建议

A系统默认为0阶!

Q1.在ANSYS中如何设定小数的有效位数比如我输入的应力值为0.0004495568而且我的一组数据都是最后几位差别但是ANSYSY自做聪明把它四舍五入了结果弄的都一样了如何更改这个有效位数呢

2.ANSYS在输出时如何修改输出的数字显示模式比如用科学记数法或者将数字缩放多少倍那样显示在什么地方能够修改呢

3.ANSYS5.7.1中什么地方能够出现那种材料库表就是说不同ANSYS的什么material models对话框生成而是通过table输入材料属性诸如应力应变曲线之类的东西?

A用a.bEn来表示如2.45678e-9材料特性用tb

--------------------------------------------------------------------------------

Q我是个ANSYS初学者现在碰到下面几个问题请各位高手帮忙看一下

1是什么原因使得相邻ELEMENT的共用节点解得的力不一样应该是方法步骤的原因而不是

误差吧请简单介绍一下

2假如我想用ANSYS模拟陶瓷的烧结过程是不是有实际*作的可能性如果有其基本的

2

原理和模型怎样

3ANSYS可以用来模拟热循环下的应力分配不均造成的焊点63%Sn&37%Pb失效的过程我看到有一篇文献介绍用在Sn球里随机分布Pb相和其它共晶相的方法建模的那么这个随机过程在ANSYS中怎么实现的而且如果要生成足够多的晶相小颗粒是不是只能手动一个一个的建立

A可能是你提取数据的时候思路有问题比如说当你提取的是有多个连接单元的节点数据而实际上你提取的是各个单元的

A相邻ELEMENT的共用节点你MERGE了吗ANSYS是通过单元的节点传递数值的从而得出结果的

Q在进行模拟计算结构多阶段施工过程计算时,在ansys里面说可以用死活单元的方法但是死活单元每次只允许选择一个,我的模型比较大,如果采用这个方法不太现实,各位同志有什么比较好的方法没有请赐教

A涉及到两个方面的问题一以TABLE形式定义参数二DO-LOOP循环你可以参考ANSYS GENERAL GUIDE中SOLUTION中对多重载荷步的叙述很简单

Q比如在beam4和shell63中显示的输出结果有应力和变形以及节点等效荷载请问如何才能显示截面弯矩扭矩和剪力呢

A在ELEMENT TABLE 中首先定义变量注意一你要知道BEAM设置了几个积分点注意二要知道SEQUENCE NUMBER值所以你需要参考ELEMENT REFERENCE 中的单元介绍

Q最近在计算一个模型管状在地下分析中要施加地层给予的支撑我想简化为弹性支撑可是发现在ansys中找不到这样的支撑情况

A首先你要知道弹性支撑的弹性系数具体的计算方法你应该知道的然后用LINK8单元模拟要注意的是刚度k=EA/L由此确定LINK单元的面积A来体现单元的弹性刚度我对此作过验证结论是可行的

Q用弹簧可是问题很多因为是一个管体整个管壁都有支撑简化为弹簧出现的问题很多了比如用多少个合适隔多长合适线性支撑和非线性支撑能等同吗我最近看校验手册发现V135这个例子可以借鉴可是还是有许多麻烦哪位大侠如果有这方面经验请指点指点

--------------------------------------------------------------------------------

Q如何梁单元同时施加两个方向分布荷载竖向侧向

A同时施加两个方向的均布荷载相当与分别施加不同的均布荷载即分别加

主要是要选择好均布荷载的ID参考单元参考手册中的说明例如BEAM4 INPUT SUMMARY 中

Surface Loads Pressure: 有face2(i-j)(-Y normal direction) 这里2就是LOAD ID 以单元坐标的-Y轴为正你试一试看

Q请问如何打印图形结果

A用ANSYA中的功能就可以拉GUI菜单命令PLOTCTRLSCAPTURE IMAGES或RESTORE IMAGES

A我是用plotctrl>hard copy>to file,然后在剪贴板中就可以看见了

A我都是用抓图工具抓下来然后再打印的

A/ui,copy,save是一个截取屏幕图形的命令,现在总算又实现了gui肯定可以通过相应的命令方式来实现的观点

AAnsys 5.6已经做到了把图形窗口COPY到文件如JPG/BMP...等格式PLOTCTRL/HARD COPY/TO FILE很好用

--------------------------------------------------------------------------------

3

Q其实APDL语言是ANSYS的二次开发语言是在FORTRAN编程语言的基础上形成的你只要会FORTRAN语言这就很容易理解APDL语言了也有循环选择等编程结构把枯燥的*作用计算机自动实现然后形成一个宏*作拉

Q是的APDL语言确实很方便而且和FORTRAN语言类似很容易读懂在计算电机电磁场中电机的反电势更加方便前几天刚算了一个需要循环运算288次的程序使用APDL语言使的我可以让电脑自己运行程序保存结果数据而我的任务就是睡一觉醒来对结果数据进行处理

Q: 我看见好多人在讨论用命令流处理程序很不明白命令流和gui的形式有很大区别吗在什么情况下更适用命令呢

A很容易控制Ansys且结果文件小功能很强大可以进行循环控制分支等*作运算速度快充分利用APDL语言对ansys进一步开发

Q用ansys5.6.2计算一个压力容器,划分网格9万余个.运行环境win2000,cpu:1G,内存640M,硬盘40G,模型建立,网格划分均无问题.但是运行到占用硬盘大约9.2G的时候,就出错.提示内存无法read.请问该如何解决?是内存的问题?还是ansys本身问题?还是系统问题?硬件问题?模型建立实在无法简化.

A我的硬盘是30G的内存1GOS windowsNT,Total Workingspace是1400M,再高就控制不了Database 是56,我能解的有效节点数是140000,只是因为磁盘空间不够工作的分区是0G,当ansys探测到磁盘空间已满就自动退出尽量大的增大节点数并不是好办法因为运行速度太慢了利用对称性减小求解区域利用Lesize,esize等控制网格都是好的办法

A选用对称性的时候要注意边界条件但的确是个好办法

A以前在使用ALGOR时也遇到过类似情况当时有10万个结点对WIN2000来说已不存在WIN9X中文件最大只能2G的瓶颈但在处理超大虚拟内存需求时仍然会出问题处理方法有二

1使用PQMAGIC划分出一个20G的逻辑分区一定要跟WIN2000系统在不同的分区再在控制面板-系统中将虚拟内存文件指定到这个大分区上以尽量避免发生系统堆栈错误虚拟文

件大小仍交WINDOWS控制

2将至少两个分区的硬盘升级至NTFS由于NTFS文件可以跨分区存储且要小10%左右比FAT32要好用得多

Q我在求解是出现如下的提示

1. real constant 1 referenced by at least element types 1 and

2.

2. 2.coefficient ratio exceeds 1.0e8.

我全部忽略上述提示得出了结果结果能可信吗我有错在哪里呢模型我都检查了数遍了还是和初始文件名有关我是在一个初始名下面建立若干相关的模型

由于我的单元只是板梁的结合所以我重新检查了一下删除了一个影响不大的扁梁很扁尺寸不大系统可能认为是板了运行就通过了上述的提示可能是因为系统对某些单元有默认的尺寸限制

A有些单元类型在默认值状态下不需要输入REAL CONSTANT但是当你改变了某项选项你就需要重新输入REAL CONSTANT你可以用UTILITY>SELECT>ENTITIES命令中的SELECT ELEMENT BY ATTRIBUTE来查看单元材料REAL CONSTANT等等是否定义明确然后用APPLY REPLOT命令来显示每个不同材料的单元

A我试了一下按照你所说的步骤出现了如下提示a default command range is not

provided for the ESEL

4

command.the ESEL command is ignored.所以没有任何结果可见是不是需要对某些default的命令做些改动呢

A.real constant 1 referenced by at least element types 1 and 2.此句意思好象是说你为1号单元和2号单元或是别的单元定义的real constant 均为1你需要检查你的模型的单元类型和实常数部分

A第二个问题经常出现.第一个问题有点问题,你看看单元定义对没有如果单元都没有划分对,结果当然不可靠啦

-------------------------------------------------------------------------------- Qslab on girder highway bridge 桥板为SOLID单元由于girder上部flange嵌贴在板的下表面所以建模时我将桥板定义为SOLID单元桥板下表面AREA定义为SHELL63单元但SOLID和SHELL单元的自由度不同这样建模肯定有不对的地方运算结果与实验结果也相差了40%有什么好的解决办法

A你可以根据需要为不同的单元增加自由度DOF顺便问一句你的模型的两个部分的节点是否完全重合提高计算精度要用更细的网格划分另外要保证模型中每个单元的长宽比和理是否有使用大应变的需要在ANSYS非线形手册中有影响精度的说明

A本来SOLID体的下表面被我定义为SHELL这时面单元的节点同时也是体单元下表面的节点但由于两种不同的单元自由度不同所以我放弃了这个模型另建模型的SHELL单元和SOLID单元相贴面节点坐标完全相同节点号不同这样才可能把它们耦合但是结果仍不理想可有什么方法?

A比较麻烦,可以通过编程输入命令流!

QA我已经把SHELL和SOLID对应点的UXUYUZ全部耦合起来让这些NODE的XYZ三个方向的位移相同但是结果仍然不理想我想约束方程的功能也进乎如此有没有更进一步的建议十分感谢

-------------------------------------------------------------------------------- Qansys有没有评价一个网格划分好坏的参数或命令请告知菜单*作APREPROCESSOR>MESHING>CHECK MESH>两个命令中一个用来CHECK网格划分另一个用来CHECK连贯性只是CHECK完了别哭就好

APLOT BAD ELEMENT!

--------------------------------------------------------------------------------

Q(1)在gui下运行宏时,重新定义原来已有的数组,ansys就会有个提示,想不要这个提示,怎么办?

(/uis,msgpop,3这个行不通,我试过了)

(2)ansys中数组的定义为a*b*c a<10040 b<255 c<7 而我实际上要定义更大的有什么好办法,难道非要写到文件里再读进来?

Aq1:有两种办法1在ultility menu / menuctrls/ message controls下设(2)将

/uis,msgpop,3 加入到start.ans文件中去以后每一次运行ansys时都会默认这个设置--------------------------------------------------------------------------------

Q我现在做一个板壳的屈曲分析想要得到该结构在荷载作用下的极限承载力利用非线性屈曲分析请问极限承载力的荷载比例因子在什么命令下可以查到

A非线性屈曲分析与特征屈曲的不同之处是它一开始就挠屈变形随着荷载的增大变形日益增大我们认为结构最危险点的变形时间历程曲线中变形若从某点开始发生较大的转折与先前完

全不一样就发生了屈曲此时的荷载就是屈曲荷载故不存在荷载比例因子即特征值问题

A要进行非线性屈曲,必须加到可以使结构发生屈曲的荷载,即荷载要比结构的实际承载力要

5

大,在进行非线性屈曲分析(极限承载力分析)时荷载逐步增加到破坏为止.计算时一定要输出每一步的计算结果,这样在post1中就会有一系列的荷载系数,而最后收敛的那个系数就是相对于极限承载力的系数,极限承载力=系数和总荷载的乘积

--------------------------------------------------------------------------------

Q如果材料一定为各向异性我用的材料为钢板该如何选用单元类型

A只要材料参数中有Ex,Ey,Ez单独输入项就能考虑各向异性

A波纹板的EXEY可以输入但是EZ该如何确定呢我自己练习了一下如果不输入EZ则计算就不能进行下去提示我材料的弹性模量要输入正值这是为什么

A如果是钢板一般认为是各向同性的厚板因轧制次数少晶体不密实存在着Z向厚度方向弱化对厚板一般地Ex=EyEz较小通过试验确定对薄板<16mm认为各向同性

A不管是什么板材料的本构关系是微小单元之间应力与应变的关系因此波纹板的几何构成是用几何模型->有限元模型计算的而不是用比较大块的板试验得出的Ex,Ey来输入你说的波纹板是不是压型钢板

A是压型钢板EXEY是通过拉伸试验确定的EX与EY不相等但是既然是考虑板壳为什么还要输入EZ况且EZ并没有试验数据

A波纹板的Ex=Ey=Ez

--------------------------------------------------------------------------------

Q怎样使梁单元把弯矩传给实体单元实体单元只有应力使它与相连的梁单元保持连续性真是件很难的事情还请各位大侠指教

A关键不是实体单元只有应力而是实体单元每个点只有三个平动自由度能不能模拟成将梁插入到实体中一小段用两个或多个实体单元节点固定梁单元根据Sevant原理忽略掉两种不同类型单元连接处的计算结果

A关键问题在于传递弯矩是为了保持变形的连续性实体单元中尽管没有弯矩但同一截面各点处的应力实际已构成了弯矩但这弯矩是多大怎样传递给相邻的梁恐怕不是简单的插入嵌固可以解决的从各结构杂志来看这一问题似乎很普遍还请各位一起深入探讨

A上述方法实质上就是施加了人工约束单元相邻近点的连线几乎是梁的切线这样就能保证转角和位移连续弯矩相等并不意味着位移和一阶导数相等这是因为两种单元的位移函数相差太大

--------------------------------------------------------------------------------

Q弹性地基梁中的弹簧(2维)在ANSYS 中应采用那一个单元?又如何*作?

A很久以前我用过弹簧单元好象是COMBIN2D你试试看

Q我在平面应变的有限元分析中,要使用初始应力这个条件,该初始应力在没加重力的时候在各点的应力是相同的.于是我在WIN2000 中相当于DOS的平台上用EDIT生成了初始应力文件,在第一个载荷步中用constsnt stress 这个选项读入初始应力并逐步调整边界条件,使计算的结果与初始应力相差不大.然后我在第二个载荷步中施加重力,进一步形成迭加应力场;在完成了上述步骤之后在第三载荷步中进行开挖,然而, 在这一步中无论挖掉多少,要么结果十分不可靠或不收敛,请众同行不惜赐教,指出错误之处.本人对您们的帮助万分感激!

模型大概如下:

模型边界为四方形,其边界取足够大,中间有圆孔需要在一定的初始应力条件下挖掉,这个初始应力条件为正应为SX=12.7MPa,SY=14.7MPa,SZ=10.2MPa,剪应力

SXY=1.7MPa,SYZ=2.1MPa,SZX=3.0MPa,在下边界及左右边界法向约束,上部边界施加稍大于SZ值的力.

A可能实施加重力时,单位不一致

--------------------------------------------------------------------------------

6

QANSYS里椭圆有没有命令可以直接画

A你可以利用ansys的局部坐标系,非常方便

gui:workplane>local coordinate systems>creat local cs>at wp origin

选择cylindrical 1,再用part1,part2,设置长短轴

这种方法在工作平面原点产生局部椭圆坐标系,在局部椭圆坐标系下产生2 keypoints,连点成线即可.

Q在椭圆坐标系下生成怎样两个点

A每一种坐标系的坐标都用x,y,z来表示,但是当你设为不同的坐标系时,x,y,z表示不同的含义

例如,在圆柱坐标下,x,y,z非别代表R,,z,输入相应的值即可.

更简单,现在直角坐标系下生成两个点,然后在生成椭圆坐标系.

可以参阅ansys基本过程手册关于坐标系的讲解.

--------------------------------------------------------------------------------

Q一个圆锥形曲面刚架屋顶怎样加风荷载

A在Ansys中施加面荷载时可以设置倾斜度首先用SFGRAD命令设置倾斜度和0点位置然后选择面或节点施加压力0点位置的压力等于施加的压力其他点的压力值等于由距0点的坐标差笛卡尔坐标系X轴或Y轴或Z轴球面坐标系R或乘于倾斜度再加上施加的压力值如水荷载倾斜度等于水的比重

Qansys默认面荷载是垂直于作用面而曲面承受的风荷载为恒定的方向那怎么办

Q如果想把上一次的分析结果作为下一次分析的初始条件该怎么做

A使用LDREAD命令

首先需要注意下面两个问题

1每一个ANSYS的实体模型的面或体都要定义对应的单元类型编号材料属性编号实常数编号单元坐标系编号这些参数在整个分析过程中保持不变而这些编号对应的属性在各个步骤中是不同的

2网格划分要满足所有步骤的要求单元类型必须兼容

步骤

1创建实体模型2

创建多个物理环境

设定一个物理环境中的单元类型材料属性实常数坐标系等

将这些参数的编号赋给实体模型的面或体

施加基本物理载荷和边界条件

设定求解选项

选择一个标题使用PHYSICS, WRITE命令将物理环境存入文件中

3清楚当前的物理环境命令是PHYSICS, CLEAR

4重复第二步准备下一个物理环境

例如可以按照下面的方法将磁场的计算结果施加到流场上

/solu 进行求解

physics,read,magnetics 读入磁场分析的物理环境

solve

finish

/solve

physics,read,fluids

7

LDREAD,FORCE,,,,2,,rmg 读入洛伦兹力

solve

更多信息请参考ansys的耦合场分析指南和命令手册

--------------------------------------------------------------------------------

Q有一个形状很不规则的刚架结构外表面是玻璃建模的时候只建了刚架风荷载怎样传递到刚架上去啊

A建模时就应该把玻璃和刚架同时建好要注意两种材料的衔接面(或LINE)和材料特性定义风载加上去后就会由玻璃传递给刚架或者简化风载成为PRESSURE直接加到刚架上

如果手上有实验数据加以比较会避免很多失误

--------------------------------------------------------------------------------

Q管路穿过厚壁,管路采用梁单元,厚壁采用体单元,在两者交接的地方怎么办?

A建模时可以考虑先建管路再建厚壁体时厚壁体由外部的平面与管路外表曲面组成这样管路和厚壁SHARE同一个曲面荷载作用下力也可以传递材料属性分配的时候先PICK UP厚壁单元定义为所选取得体单元再PICK UP管路定义为梁单元其实很多梁单元也属于体单元没有什么冲突除非管路属于SHELL或面单元时要后PICK UP此面单元定义为所要的属性

至为重要的是在MESH这两个单元面体衔接的时候要注意衔接面体对应的ELEMENTS 和NODES的COMPATIBILITY否则将会出现大的错误REFINE时也不可忽视

QA1请详细解释一下

1再建厚壁体时厚壁体由外部的平面与管路外表曲面组成这样管路和厚壁SHARE同一个曲面

2至为重要的是在MESH这两个单元面体衔接的时候要注意衔接面体对应的ELEMENTS 和NODES的COMPATIBILITY否则将会出现大的错误REFINE时也不可忽视并非英文不懂2我的情况是已经建好了厚壁体再建立管路

3这样一个问题该如何办例如

一块厚板(x1,x2,y1,y2,z1,z2)(0,1,0,1,0,0.1),一根管由

k1=x1,y1,z1(0.5,0.5,2);k2=x2,y2,z2(0.5,0.5,-2)构成求其固有频率其他条件随便管路为圆管

--------------------------------------------------------------------------------

Q模型是一个正交异性钢桥面板是个钢箱梁取了一部分来分析由主梁纵向加劲肋横隔板钢板还有沥青铺装层组成横向4800mm纵向9060mm纵向加劲肋共8个横隔板是3m一个共四个钢板层厚度14mm沥青混凝土50mm除了上面的沥青铺装层其他都是钢的我原来考虑可以使用shell63单元但我要计算铺装层与钢板层的层间剪应力似乎不好处理所以全部采用solid45建模

在中间某一荷位作用一矩形均布荷载我现在建模完成了在网格划分的时候如果按照smartsize的话可以划分但是我怎么施加荷载呢所以好像必须在荷载作用区域有整数个element如果按照smartsize的话可以划分但是我怎么施加荷载呢

这样可以press on element但是我不会自己控制网格划分在中间某一荷位作用一矩形均布荷载建模型的时候就要预先把那个部分设一个面?似乎很复杂这个面下面又有钢板又有加劲

肋

另外我即使简化成一个集中荷载来计算那也得中间有node呀

我怎么能保证我的网格在我需要的坐标值上恰好有个node我的模型并不是规则的啊

下面是梯形的加劲肋如果按map的话它会提示非法拓扑形式即使我等效成集中力我在想如果定义hard point的话可以加载到那个位置上但我好像没看到可以按坐标值定义hard point

8

A使用SWEEP或CONCATENATE命令重新MESH桥板

可以在HELP文件中找到这两种命令的用法十分有用将桥面MESH成为HEX类的ELEMENT 这样寻找加载位置会方便很多

A建议你将整个桥面划分为比较规则的四边形网格(六面体)当然网格的大小肯定与加强筋的厚度等有关这样应该能够解决你所说的问题

--------------------------------------------------------------------------------

Q一个梁上有一滑块计重力和摩擦滑块在梁上滑动求梁的响应梁为1m*0.1*0.1,滑块为0.2*0.1*0.1ex=2e11,nuxy=0.3,dens(1)=7.8e3,dens(2)=2.7e3

A1命令CONCATENATE如果你要把一个面或者体MESH成为一个四方或HEX单元这个命令十分有用你可以CONCATENATE相邻的线或者面把所要MESH的面变成一个四边形或体变成一个六面体(例如5边形CONCATENATE两条相邻的线于是ANSYS认为此时的面是四边形你就可以进行3OR4SIDES的MESH)所以无论你的面由多少边都可以进行MAPPED MESH

2SWEEP和上边的命令一样都是我所喜欢的命令只要你的体已经建成SWEEP好一个面以后整个体MESH可以一扫而过HEX的ELEMENT加载在什么地方或什么样的荷载都十分方便不需要在建模的时候就考虑加载位置

3至于选用单元具体问题具体分析但是衔接一定要注意特别是衔接面的ELEMENT要对应好以免出现错误

4自己的一点心得MESH十分重要别在这上面偷懒我曾经用FREE MESH与MAPPED MESH得到的值分别为19.2和25.9而理论值是26.3所以请小心

--------------------------------------------------------------------------------

Q我做的是一个三层的压力容器每层厚都是2mm请问我可不可以只定义一种单元类型没有shell91,99) 另外这三种材料的材料属性用什么命令分别定义给这三层另外用轴对称分析的话应该怎么做

A划分网格前,用attribute>define可定义不同的材料属性

--------------------------------------------------------------------------------

Q怎样进行螺纹的失效分析

我现在要做一个Failure analysis of a thread joint

对于螺纹失效分析问题我不知道Ansys里是如何处理的在creat 3D model时是否要建立完全体现thread的几何model,然后直接用普通的solid element进行mesh,再按照通常的分析过程分析还是在建模时将joint简化为一普通轴而采用特殊的element进行mesh然后再分析

对于前一方法用Ansys creat thread geometry比较困难至少我不知道怎么做而且mesh 时也可能会出错(thread geometry 比较复杂类似sharp corner or small area)

对于后一方法在model的时候比较简单但不知道用什么样的element才能体现tread和crew的应力集中特性

--------------------------------------------------------------------------------

Q有几个很简单的问题

1在Ansys中如何在一个2节点的单元内同时加多个集中荷载我用了一个很土的办法加了几个LoadStep

2如何使Beam3单元的内力以局部坐标输出我用List Results->ElemSolution->Nodal Force后结果仍和整体坐标系下一样

3我用一个静定结构进行分析结果输出后和实际的受力方向刚好相反受拉的成了受压但反力的方向是对的

--------------------------------------------------------------------------------

9

Q有一个圆柱桶内装有70%的水如何将桶和水的重量考虑进去

Aload>gravity..

A把水另设为体单元

--------------------------------------------------------------------------------

Q四边固定的板用体单元有一条线用梁单元从板的任意一点穿过分析其固有频率

如果不将线体联为一体分析时各算各的问如何将线体联为一体

glue无法将线\体连为一体

glue只能将同样的元素连为一体(如体和体,面和面,线和线)

--------------------------------------------------------------------------------

Q如何给模型加3g的冲击载荷

AIn short time(for example:0.0001 second),用一小time 步加载

A 3g指的是加速度冲击载荷例如在地震时的地基响应

例如

在固定端上用弹簧悬挂一根杆在固定端加一个3g的冲击载荷g指加速度求杆的响应

1之所以用弹簧悬挂是为了减振这个问题主要想知道如何加这种载荷

2类似于地震时的地基响应问题

3我在ANSYS的GUI中没有找到这个功能

4BTW在ANSYS中如何加弹簧元

--------------------------------------------------------------------------------

Q我想知道用ANSYS分析岩土洞室开挖效果好不好

我刚开始学ANSYS不知道ANSYS用于分析有裂隙岩体和松散土体方面效果好不好其精度如何望各位朋友能给予指教

Adyna做隧道开挖比较擅长而且看起来非常形象

A用ansys里的ls-dynadyna是绑在ansys上的你选用模块就可以

--------------------------------------------------------------------------------

Q请问这个错误是什么

划分网格时出现Shape testing revealed that 16 of the 30 new or modified elements violate shape warning limits.To review test results,please see the output file or issue the CHECK command.这个怎么改

A检查单元prep7>meshing>check mesh若有不合格的单元用别的办法划单元还可以设limits有时候其实这种警告对结果没有多大的影响完全根据你的经验

A如果可以进行下去可忽略该warning

一般来说温度场分析对形状测试更严格一些一般不会对问题产生太大影响可通过剖分更密的

网格来解决

--------------------------------------------------------------------------------

Q如何得到某一温度值的距离-时间关系曲线

在分析平板移动热源温度场时希望取得时间历程中熔化温度值在各时刻对应的平板上出现的位置或节点号即某一温度值的距离-时间关系曲线或者如何得到某一温度值在给定时刻出现

的位置

A我想你可以把你需要的那些温度值的点选出来设置成路径在用路径来画图

在path operation下

A谢谢你提供的建议我不知道这些温度值点出现在什么位置因此路径方法不可取也许只能采用set设置时刻,sort排序各节点温度prnsol找出温度值对应的节点此法较为繁

10

琐

--------------------------------------------------------------------------------

Q5.4版如何用extrude将2D单元拖拉出3D单元的如何进行旋转体的结构分析加离心力吗怎么加

A先用2D单元剖分面再画出拖拉路径直线修改默认的设置3D单元类型材料实常数extrude 即可

--------------------------------------------------------------------------------

Q一个后处理问题我用shell93作的壳体,节点应力输出:

PRINT S NODAL SOLUTION PER NODE

***** POST1 NODAL STRESS LISTING *****

PowerGraphics Is Currently Enabled

LOAD STEP= 1 SUBSTEP= 1

TIME= 1.0000 LOAD CASE= 0

SHELL NODAL RESULTS ARE AT TOP/BOTTOM FOR MATERIAL 1

THE FOLLOWING X,Y,Z VALUES ARE IN ELEMENT COORDINATES

NODE SX SY SZ SXY SYZ SXZ

225 0.51606E+06 0.49482E+06 0.0000 -3464.6 -10039. 38600.

225 0.50183E+06 0.53119E+06-0.48105E+06 7004.0 -10031. 38565.

3020 -61720. -0.38298E+06 0.0000 -3034.2 -9606.2 2571.7

3020 0.21884E+06 10255. 0.0000 81566. 1491.2 33824.

3020 75691. 0.31964E+06 0.0000 18182. -9597.7 2285.4

3020 -0.24616E+06 -44015. 0.0000 -82203. 1491.2 33824.

我想弄明白为什么有的点(225)输出2行,有的(3020)输出4行.如果平均化了,那应该只有2行(Top & Bottom)?

A由于你的模型中有多种材料在材料交界处的节点结果跟交界的组数有关该节点共几组材料则有几种结果节点结果是由相邻单元的高斯点结果外推出来然后平均的在PowerGraphics ON 状态下ANSYS在材料交界处的节点自动按材料区分同一组材料的单元外推点参与平均不同材料的不参与平均所以节点交界处有几种材料该节点就有几组结果

--------------------------------------------------------------------------------

Q如何根据已创建的关键点划分网格使节点恰好在关键点上

A把关键点变成硬点即可或在关键点处生成硬点

--------------------------------------------------------------------------------