《爐石戰記》官方網站遊戲卡牌系列電子競技遊戲卡牌系列電子競技在《爐石戰記》世界大賽中，由你決定誰來為國參賽！投票結束、隊伍名單公布後，各代表隊互相較勁，看地表最強的爐石高手來自哪個國家/地區。瀏覽下方的隊伍名單，為喜愛的選手加油打氣！世界大賽即將展開！國家/地區代表隊已出爐，比賽即將展開。敬請於7月17日（PST）觀賞賽事。亞太區歐洲美洲澳洲中華人民共和國香港印尼日本南韓馬來西亞紐西蘭菲律賓新加坡台灣泰國越南奧地利白俄羅斯比利時保加利亞克羅埃西亞捷克共和國丹麥芬蘭法國德國希臘匈牙利以色列義大利哈薩克荷蘭挪威波蘭葡萄牙羅馬尼亞俄羅斯西班牙斯洛伐克瑞典瑞士土耳其烏克蘭英國阿根廷巴西加拿大智利墨西哥秘魯美國種子FroStee選擇CitizenNappa選擇GingaNinja選擇Akumaker種子Leaoh選擇OmegaZero選擇YouLove選擇Trunks種子blitzchung選擇kin0531選擇losefield選擇Shy種子Rezdan選擇DouAhou選擇Joth703選擇nexok40種子Tredsred選擇trahison選擇oyatsu選擇posesi種子DacRyvius選擇Sooni選擇Steelo選擇Surrender種子wkyew90選擇porkmypig選擇Gilgamesh選擇ChongGEr種子Mage選擇Pathra選擇Jakattack選擇TINGT0NG種子Staz選擇WaningMoon選擇Chalk選擇CaraCute種子Sequinox選擇Khaius選擇Crumpled選擇katsucurry種子tom60229選擇ahqRoger選擇reall選擇ahqShaxy種子Disdai選擇ViperPro選擇patoyao選擇Kapi種子Tuan選擇Neilyo選擇Frosty選擇C0ntrolBoy種子fateos選擇Crumoet選擇Xandi選擇SebastianJo種子NoName選擇fzmushko選擇NewMeta選擇Картошечка種子Maverick選擇Sjoesie選擇oliech選擇Chinoize種子Jengo選擇SilentStorm選擇MorGalad選擇superlubo種子Insane選擇nuke選擇Lavvek選擇MagicZYYZA種子Jarla選擇Findan選擇pokrovac選擇Faeli種子Twink選擇Crane333選擇NikolajHoej選擇Hoej種子Maur1選擇zumpp選擇Vardu選擇Janetzky種子Swidz選擇Windello選擇Kalàxz選擇Odemian種子Seiko選擇Viper選擇Sintolol選擇Bunnyhoppor種子Fenomeno選擇Dethelor選擇tholwmenos選擇SpinaLogaMan種子exigo選擇PenhaDani選擇zridex選擇SteveFrench種子Glaser選擇Falular選擇yonatan24z25選擇Hatul種子SCACC選擇Meliador選擇Pool8選擇Turna種子Snowdrop選擇Cagey選擇limonadik選擇Naiman種子Tyler選擇Mitsuhide選擇Theo選擇ThijsNL種子hunterace選擇Zorkthar選擇Kolbein選擇ntR種子dawido選擇Tuttek選擇dogmeat選擇A83650種子Zuka選擇TheViper94選擇DrBoom選擇LazyA種子Rdu選擇angeloshaggy選擇Charzuull選擇Kevineter種子Fox選擇SilverName選擇RENMEN選擇ШтанУдачи種子riku97選擇Dalesom選擇Jac0b0選擇AKAWonder種子martim選擇cekmara選擇HiImJoZo選擇zataro種子Bozzzton選擇Kenobe選擇Orange選擇Arcade種子Gusonaj選擇Nights選擇Hypno選擇TheRabbin種子ThunderUP選擇Woven選擇Fujitora選擇berserkci種子Kolento選擇Neirea選擇OldBoy選擇iNS4NE種子BoarControl選擇GreenSheep選擇Toastmonster選擇J4CKIECHAN種子Joaquin選擇PNC選擇Nalguidan選擇Pinche種子Rase選擇Neves選擇Perna選擇Lucas種子Monsanto選擇DocPwn選擇Purple選擇TheJordude種子Anguistar選擇Tempo選擇topopablo11選擇Macheen種子Yinus選擇Arreador選擇kolitas選擇Empanizado種子chosenone選擇edanco選擇Izzy選擇UchihaSaske種子justsaiyan選擇Muzzy選擇Fr0zen選擇TerrenceM/獎勵條款：無需負擔任何費用。每位參加者只能選擇三名選手。您必須在日晚上11:59（PST）以前做出選擇，並且必須在日晚上11:59（PST）以前擁有《爐石戰記》帳號才能兌換並使用所獲得的卡牌包。活動過程中獲得的卡牌包無法更換、轉讓或兌現。若與法律牴觸時則無效。注意！您必須在《爐石戰記》競技場或休閒/排名模式中贏得三場比賽才能進行投票。接收最新資訊未开启cookie无法登录使用i春秋的完整服务，请设置开启浏览器cookie
第十一届全国大学生信息安全竞赛 (创新实践能力赛-分区半决赛)
比赛开始时间： 00:00:00
主办方：教育部高等学校信息安全专业教学指导委员会
比赛模式：线上赛
第十一届全国大学生信息安全竞赛 (创新实践能力赛-总决赛)
比赛开始时间： 00:00:00
主办方：教育部高等学校信息安全专业教学指导委员会
比赛模式：线上赛
XCTF分站赛—SCTF
比赛开始时间： 09:00:00
主办方：三叶草小组（Syclover）
比赛模式：线上赛
2018年第四届全国网络空间安全技术大赛——线下决赛
比赛开始时间： 09:00:00
主办方：中共中央网络安全和信息化委员会办公室、教育部高等学校兵器类专业教学指导委员会、中国兵工学会、国家国防科技工业局、中国兵器装备集团、中国兵器工业集团、陕西省兵工学会、西安工业大学
比赛模式：线下赛
TSCTF2018“京东安全杯”第四届 北京邮电大学信息网络安全技术挑战赛
比赛开始时间： 09:00:00
主办方：北京邮电大学网络空间安全学院
比赛模式：线上赛
第二届红帽杯网络安全攻防大赛-线下赛
比赛开始时间： 00:00:00
主办方：广东省计算机信息网络安全协会
比赛模式：线下赛
2018世界智能驾驶挑战赛（WIDC）——信息安全组汽车破解挑战赛
比赛开始时间： 00:00:00
主办方：天津市人民政府、国家发展和改革委员会、科学技术部、工业和信息化部、国家互联网信息办公室、中国科学院、中国工程院、中国科学技术协会
比赛模式：线下赛
竞赛训练营
Writeup分享
分享人：An0n7m
分享人：snowleopardbin
分享人：闪闪的夜空
2018赛季网络安全竞赛排名
<td title="44.310
<td title="61.609
<td title="86.340
<td title="65.745
<td title="01.887
<td title="86.965
<td title="36.073
<td title="05.045
个人能力榜
关于CTF大本营
CTF大本营意在为网络安全爱好者提供最新、最全的比赛资讯、比赛练习题、解题分享等信息。同时，i春秋致力于打造最好的网络安全竞赛平台，目前我们已推出了5种竞赛模式，最多可支持千人同时在线竞技，具体详情请
北京五一嘉峪科技有限公司
海淀区中关村软件广场C座
京ICP证150695号
京公网安备37号待解决问题
咨询者：小明12&&城市：
&& 咨询时间： 12:15&&解答数：69
中国股市什么时候实行T+0
离您最近的顾问回答
尚未登录，可以匿名提问。也可按此
一个电话，一个垂询，一份财富
想理财？从咨询开始！
5分钟快速响应，100%回复率
推荐理财顾问HOT 月积分
手机：积分：76
手机：积分：1188
手机：积分：795
手机：积分：386
最新理财产品
最新理财文章
全国找理财顾问abaqus问答自己总结
&&&&1、q:预拉钢筋怎样施加预应力，请各位指点~~~~ q:我在文档里看到要在 inp 文件定义一个 rebar， 但是 rebar 只能用于 shell, membrane, and solid elements 。我现在想做的是一个预应力拉索，不是镶嵌在 shell, membrane, and solid 这些单 元里的，而是独立的一根拉锁。拉索单元打算用 truss,但是怎样&&&&在 truss 上使用 rebar 啊请高 手指点 还有个问题，我看到别人的 inp 文件，如下： *rebar,element=continuum,material=rebar2,name=ubar top1,1.005e-4,0.15,0.0,0.5,1 第二行第一个是 setname（top1),第二个是 rebar 的截面面积（1.005e-4)，那第三、第四、第五 是指什么（0.15,0,0.5),最后一个应该是方向，是 1 方向。 哪位高人指点下第三、四、五项分 别代表什么 a:施加预应力 *initinial conditions,type=stress,rebar elset,rebar name,所施加预应力的值 ,另 prestress hold 为保持所施加的预应力的值不变，我的 理解是防止别的构件吃掉所施加的预应力， 造成所施加预应力的损失。 使用了这个命令之后就避 免了这种损失，保证所施加的预应力都施加到了钢筋上。 a:谢谢指点，你所说的应该是把预应力加在 rebar 上面，但我发觉 truss 单元不能定义成 rebar， 其实是我多想了，truss 本来就可以当拉索，实际工程中加预应力只是为了使钢绞线拉紧，起到 张拉作用，而在 abaqus 里，truss 本身就是拉紧的，不用施加预应力 a:我知道模拟加强筋的时候需要用 rebar，但 钢筋确实可以直接用 truss 来模拟 ,而 lz 所说的预 应力其实其实只是施工时的张力而并不是真正意义上的预应力， 比如螺栓预应力之类的。 如果是 索的话可能是要施加预应力的，仅个人看法。 q:请教：做一个空间钢框架结构，梁柱用梁元，板采用壳元，打算采用 tie 命令（共用节点） ， 但不知该如何实现 a:我想可以用*equation 实现，共用节点的约束情况自己在这一命令下定义。 a:我因为用命令比较多，但是用 cae 我想一样， 在 cae 里进入命令编辑器，然后编辑就是了， 写入*equation 命令， 指定约束的自由度（这个看一下标准手册，写得很清楚） q:“response spectrum analysis（响应谱分析）与 modal dynamic analysis（模态动力分析）区 别在什么地方 如 response spectrum analysis 可以进行结构设计 但 modal dynamic analysis 是用来干什么的阿 a:就我知道的，modal dynamic analysis 应该是振型分解法做动力解析。分解为单自由度体系再 取有限个进行组合求反应。 q:abaqus 如何施加地震荷载 a:可以参考 abaqus 6.3 的例子，seismic analysis of a concrete gravity dam 可以使用： 1。*amplitude, name=amp, input=seismicdata.dat 输入地震波 2。*boundary, type=acceleration, amplitude=amp 来施加荷载。 在&abaqus example problems manual&的 2.1.15 seismic analysis of a concrete gravity dam a:这是一个相对的问题，你可以推导一下那个动力方程， 结果是：ansys 是取基础固定，解出结构相对基础的相对时程,而 abaqus 是在边界上施加加速 度时程，解出的是结构的绝对时程。因而两者相差一个所施加加速度时程。在边界上施加加速度时程用显式计算速度较快（主要取决于你的 cpu）abaqus 的 standerd 里面也是显式计算，另外它的 explict 模块也可以，具体例子可参见 seismic analysis of a concrete gravity dam 的分析； ansys 用的是隐式求解（主要取决于你的硬盘速度和 cpu） ， 并且结果暴大（应为它默认输出一大堆数据） ，可用 ansys/lsdyna 进行显式计算。 偶曾经作过共同作用的时程分析，有一点点经验。 q:我老板要求我用的是 concrete damage plasticity model 1。在做非线性分析的时候，怎么知道结构破坏了，如何得知什么时候获取了结构的极限承载力 2。如何确定 concrete compression damage 和 concrete tension damage二者似乎都无法通过 试验得到。如果说可以忽略 compression damage，那么由于 tension stiffening 的存在（下降段 相对而言变长了） ，tension damage 的定义就显得很重要了。在 abaqus 关于重力坝的例子中， tension stiffening 和 tension damage 的数据，似乎有一些关系。 3。定义 concrete tension stiffening 的时候，我看一些论文上面关于 tesion stiffening 说明的时 候，说如果 tensiong stiffening 定义得太小，收敛会很难。但是 tension stiffening 的大小在用 *concrete tension stiffening 定义如何体先呢，在前两个数据是 remaining direct stress after cracking 和 direct cracking strain. 4。tension stiffenig 与配筋率是否有关系呢如果有的话，关系是怎么样的呢 谢谢大家了。 a:我认为 damage plasticity model 主要用于混凝土受循环往复荷载的情况，concrete compression damage 和 concrete tension damage 也只在这时需要定义， tension stiffening 好像 就是定义混凝土单轴受拉本构关系，与 tension damage 应该没有关系。与配筋率好像也没有关 系。 damage plasticity model 既可以模拟循环拟静荷载的情况，也可计算结构在动力荷载（地震、风 以及冲击荷载的作用） 通过 concrete compression damage 和 concrete tension damage 考虑混 ， 凝土的拉压异性，tension stiffening 则可考虑拉伸屈服后的软还阶段本构行为，同样是混凝土损 伤全过程的一部分，它与配筋率及网格划分密度的相互关系很大程度上影响到计算的收敛性。 a:损伤主要针对循环加载问题，如果你不是循环加载的话可以不考虑损伤 ,我做过比较 a:楼上的大侠好象问题考虑的不是很全面吧.如果在单向加载时同时存在拉压,就应该考虑.而其 中的恢复系数是循环加载下考虑损伤恢复 a:我现在也在搞 concrete damage plasticity model,希望能和各位大虾交流交流 a:。 在做非线性分析的时候， 怎么知道结构破坏了， 如何得知什么时候获取了结构的极限承载力
这个问题牵涉到结构在损伤到多少时即认为是损坏.也就是通常说的怎么把结构算倒.不知道我 理解的是否正确.如果是这个问题就牵涉到损伤系数的取法.因为是整个结构,就须要把所有的损 伤系数加权平均.加权的取法也很有讲究,有人认为达到 0.7 即认为倒塌,可以参考王光远的文章. 现有有限元程序还不能完成,是个急待解决的难题. 当然是有关系的了 concrete compression damage 和 concrete tension damage 分别是受压损伤 系数和受拉损伤系数函数,在单项受力时是直接乘入混凝土刚度的.在往复荷载下要乘损伤恢复 系数，在三轴应力下还要乘加权因子。 考虑了该参数后，开裂应变就不等于塑性应变。开裂应变是考虑弹性模量折减的非弹性应变。因 此此时塑性应变不等于非弹性应变。这是关键啊 a:谢谢各位的回复, 帮助很大. 现在 riks 去求解的时候,可以获取极限何在,结果与实验结果差别 不大,但是在极限荷载下的位移与实验结果差别很大,无论用 shell 还是 solid 来模拟钢筋混凝土的 话, 位移误差都很大. 请问大家有没有碰到着个问题呢 q:钢筋混凝土梁中的箍筋如何加入用 cae 操作可以直接加入吗a:cae 中可以用 truss 单元模拟箍筋和纵筋,再采用 embed 约束实际到混凝土实体中.钢筋单元的 位移由他所在的混凝土单元的节点位移来决定.abaqus 中的帮助文件里面有详细的说明,abaqus 加钢筋的办法有几种： 1。采用 rebar layer 的办法，在 part 里面画一个面，然后在 property 里面定义一个 surface 为 rebar layer，把这个 surface 的属性赋给前面的 part 里面的那个面。然后在 interation 中 embed 中把钢筋层 embed 到混凝土实体中去。 2。采用桁架单元的办法，在 part 里面建好纵筋和箍筋的钢筋骨架，在 property 中分别赋予截 面和属性，在 interation 中的 embed 把钢筋骨架 embed 到混凝土的实体中去。 3。如果是作构件的话，第二种办法建的比较精确，而且后处理比较方便，查看钢筋单元的应力 比较直观，如果是作结构的话，第一种钢筋层的办法比较好，但是个人觉得钢筋层的办法纵筋和 箍筋的位置定义的不是很明确。 楼上高手多谢指点。你说的建钢筋骨架是指在 part 中用 wire 画吗一根一根的画我做一根混 凝土梁的话，我先创建了一个实体代表混凝土梁，然后再 create---part 画钢筋，画钢筋的坐标 系与前面的坐标系应该是相同的在 part 模块中创建两个 part,对吗这两个 part 的坐标系完全 相同，要不然在 interation 中没法 embedded.。不知我的想法对不 liuxinmei2005 兄见笑,我学 abaqus 也不到一个月,绝对不是什么高手.这是我自己想的办法,我自 己的办法是这样的,part 里面画一个纵筋,画一个箍筋,在 assembly 中逐个加入纵筋和箍筋并且定 位,形成钢筋骨架,然后再 embed.我的办法很笨的,好处就是模型中钢筋的位置非常明确,纵筋箍 筋的后处理也都很方便,这种办法对于建一个构件来说还可以,一般 10 分钟就建完了.对于建结构 来说,钢筋层可能好些,很方便,但是后处理想直观的看钢筋的应力就不方便了. 箍筋的建立的办法画图说明一下: 1.在 part 中画一个 wire, 2.赋予 wire 截面和材料属性 3.在 asemnly 中插入,平移或者旋转来布置箍筋的位置 下图就是箍筋的布置,同样的办法也可以建立纵筋,这样就形成了钢筋的骨架,这么建模比较仔细, 后处理比较方便。 对于剪切破坏的钢筋混凝土构件,箍筋的作用比较关键,必须要建立箍筋单元,对于受弯曲破坏的 钢筋混凝土构件,可以不建箍筋。 对于剪切破坏的混凝土构件， 箍筋的作用在任何一本混凝土教材里已经叙述的很清楚啦， 由于箍 筋的作用对该类型构件影响非常大，当然一定要考虑。 对于轴向受压构件，由于箍筋对侧向膨胀的约束作用在后期非常明显，所以当然也要考虑。 而对于受弯破坏的构件， 起主要作用的是受压区混凝土以及纵向配置的钢筋， 其性能受箍筋的影 响很小，从计算成本及建模的方便与否等因素考虑，可以忽略。 请问 ilxy，assembly 中的插入操作怎么实现俺在菜单中只找着了平移，旋转 所谓的&插入&只是形象的说法,并不是说 constrain 里有相应的操作,都是通过那几个基本命令实 现的.我是用点重合的命令插入钢筋的.在混凝土的相应位置定义一个 datum point,再让该点与钢 筋相应的关键点重合就好啦.比较苯的方法,供参考. 钢筋嵌入单元可以直接横穿嵌入钢管混凝土吗也就是，即嵌入了混凝土又嵌入了钢管 请问另外一个问题， 不同材料属性的两个 instant 可以用 merge geometry 在一起， 而属性不变 吗我要把型钢梁和栓钉焊接在一起，所以想用这个命令，总说有属性错误，有办法吗还是要 用 merge mesh 吗还不会用 merge mesh，请高手指点。 q:怎么在属性里面定义应力与应变的关系 a:可以直接在材料定义的地方输入数据。自己把本构关系的曲线分段，然后输入每一段的应力和塑性应变数据对。 看一下入门手册关于材料非线形那章说明就知道了。 可以在 inp 文件中直接输入， *plastic,hardening=（强化方法） &yield stress&,&plastic strain& 就可以了。 q:我以前没接触过 abaqus，有一点 ansys 的基础，但也只是刚入门，没做过具体的项目。 目前在做竖井施工对软粘土的扰动效应研究的硕士论文， 其中一部分工作是用有限元软件模拟竖 井施工的挤土效应，考虑土体自重，得到竖井周围土体的塑性区范围，以便与理论解互相验证。 应该是一个比较简单的圆孔扩张问题：轴对称建模，在土体中给定一个孔半径大小的位移，不需 要考虑接触问题。 理论解采用的是摩尔库仑屈服条件，而 ansys 里面没有这个模型，所以打算采用 abaqus。由 于只是论文的其中一部分， 不可能花大量的时间来深入学习 abaqus， 这两天在版上看了许多贴 子，发现在短时间内掌握 abaqus 并不是一件容易的事情，甚是着急，希望高手能给小弟指点 一下迷津，如何针对我这个课题做一下分析3x a:把 abaqus documentation——abaqus/cae user&#39;s manual 里面的第二个例子做一下，就了解 abaqus cae 的基本过程了。例子里面没有而你需要解决的问题是初始地应力的形成。 过程如下： （假定基本操作你已经会了） 1、建好模型，设置好材料参数，只施加重力荷载，计算一遍，点 view results 进入 visualization 模块，点 report 菜单，将 mises 应力的 s11,s22,s33,s12,s13,s23 输出至文件，计算方式可以选 择 average，将文件命名为 xxx.inp。 2、将文件中的非数据行删掉，再在每个数据之间加一个逗号（用 excel 很容易的） ，每行的最后 一个数据不用加，保存。 3、在 cae 中重新进入你原先建好的模型，在 step 模块里面再新建一个 step，进入 load 模块， 沿着圆孔径向加位移（和加重力没什么区别） ，这个步应该发生在新建的 step 里面，注意原来加 的重力荷载不要删了！ 4、进入 job 模块，在 job manage 里面点 write input，于是 inp 文件被刷新了，用文本编辑器打 开 inp 文件，找到*step 区域和*material 区域，在这两个区域之间加入*initial conditions, type=stress, input=xxx.inp 5、在开始——程序里面将 abaqus command 界面打开，假设你原先的工作文件名为 yyy.inp，打 入 abaqus job=zzz.inp input=yyy.inp，所以路径都应该一致。 6、 打开 log 文件， 不管有错误还是顺利都会有提示， 计算完成后， command 界面打入 abaqus 在 viewer database=zzz.odb 查看结果 忘了说了，abaqus 可以形成两种 inp 文件，应该用没有 part 和 assembly 的那种，这样单元和 结点的编号不会重复，*initial conditions, type=stress, input=xxx.inp 才能正确执行，参考 &cae_no_parts_input_file& q:就是一根三点弯曲的普通钢筋混凝土梁，我想看它的荷载－位移曲线，如何操作 a:利用 riks 分析，记录位移和反力，再在 xy plot operation 里面作图。 *step,inc=100, nlgeom *static,riks 0.001,1.0,,0.04,,4,3,-0.080 riks 指的是弧长法也可以利用*node print 命令，把节点位移和支座反力输出到.dat 文件，用 excel 画图，很方 便，修改也容易。 *node print 只能在 abaqus/standard 中起作用，在 abaqus/explicit 中不能使用 q:有个问题想请教大侠： 我在做一个双材料界面的斜裂纹问题，模型见图。裂纹两侧是不同材料。如果不定义局部坐标， 输出的裂纹附近的正应力和剪应力应该都是整体坐标系下的，显然这时候的 s22,s12 并不满足 连续条件。 现在就是想定义局部坐标系，并且输出局部坐标系定义下的 s22&#39;，s12&#39;，不知道 abaqus 具体怎 么样实现，请指教。 看到过帮助里提到过 orient 子程序，不知道怎么用 a:在 simwe 上坛友的帮助下 （当时只记了他的 mail:asia_,忘了他的 id， 很抱歉） ， 解决了这个问题。 1。在打开 odb 文件的时候，前面的 read_only 前面的勾去掉。 2。create coordinate system 选择一种方法建立局部坐标。 3。coordinate system manager 选择建立的 csys--&move to odb 4。重新打开 odb，main menu bar--&result--&options...--&transformation--&user-spcified--&choose your csys--&apply q:adaptive mesh 是不是只能用在 explicit 里头呀 a:如果想要在 standard 里面用，在 cae 里该如何操作 在 step 模块，在 other 菜单里设定 我昨天凌晨三点终于把弹簧加上了，幸亏我是一个人一个宿舍，否则别人会有意见的。在我学习 abaqus 的过程中，从中 simwe 论坛中获益不少，有了办法，不敢独享，特传上来。 要点：1、弹簧单元只能在 inp 文件中手工加入。2、*element，----和*spring，-----这两个 语句要挨着。3、新加的节点号和单元号，不要和已有的冲突。4、有不明白的，我愿意解释， 最好是晚上，白天我和老师在一个办公室，包教包会 qq 号：，请注明是 abaqus 同行 q:我想将 abaqus cae 中的图形导成 jpeg， bmp 文件格式。 或 找了好几天， 还是不知道怎么弄 知道的，说一下具体操作。谢谢！ a:file -& print -& to file, select tif or png or ..., you will get it. q:section 那项里没有 point section 这一项啊，不知道如何定义,望大虾们指教，谢谢 a:create section 对话框， category 一栏选 other，然后在 type 中有 point。 在定义属性的模块下，工具兰有一个定义特殊属性的菜单可以定义惯性属性 q:有哪位大侠知道怎样在 abaqus 中输入初应力文件吗： ） a:abaqus example problems manual 1.1.10 stress-free element reactivation 里面有部分相关内容，可以查阅之。 ** -------------------------------------*stepstep 1: add initial stress state & remove liner *static *dload soil, by, -20.0e3 *model change,remove liner, *cload,amplitude=relax 100,1, 5.,1, 4.3918e+04 q:如何让 abaqus 以 double precision 运行 a:job manager-&precision q:无意间发现 abaqus 里输出的都是积分点的应力，如果想得到节点处的应力，该怎么办 a:这是节点处的应力输出命令： *el print, position=nodes *el file, position=nodes 这是单元质心处的应力输出命令： *el print, position=centroidal *el file, position=centroidal q:我现在在分析 abaqus example problem manual 的 1.1.7 例子， 但是当把 inp 文件导入时 出现这个问题： warning: the following keywords/parameters are not yet supported by the input file reader: --------------------------------------------------------------------------------*elfile *elprint *initialconditions, type=hardening *nodefile *nodeprint 为什么硬化初始条件导不进去了如果导不进去又该怎么处理了请各位高手帮忙 q:请问如何实现，用户子程序里如何实现分层填土荷载 a: subroutine dload(f,kstep,kinc,time,noel,npt,layer,kspt, 1 coords,jltyp,sname) c include &#39;aba_param.inc&#39; c dimension time(2), coords (3) character*80 sname h=kstep*2 width=12 if(coords(1).le.width)thenpload=(width-coords(1))*2/3*20000 else pload=0 end if if(coords(1).le.(width-1.5*h))then nload=-(width-1.5*h-coords(1))*2/3*20000 else nload=0 end if f=pload+nload return end q:在 abaqus/explicit 计算中，请问： *ouput ,field,number interval=n,time marks=yes (or no) *restart,write,number interval=n,time marks=yes (or no) 中，这里的 n 是什麽取多少合适marks=yes 或 no 有何区别谢谢！ ！ 上述的 n 在 abaqus/cae 中怎样设置 a:n 是你取的输出结果的间隔数， n 为 4 是指每计算 4 个迭代步骤输出一次结果。 取 time mark 取 yes 就行了。 q:在 abaqus 中做钢筋混凝土的分析，钢筋怎么样加进去，能不能象 ansys 中直接画进去还是 要在 inp 文件中自己输进去！请大侠指点！ a:使用 rebar element q:abaqus 的 cae 中有没有弹簧元件， 怎么调用 a:you should use editor to edit *.inp. add spring element and the property of spring. q:请问各位高手， 如何在 viewer 里将 contour 最大值与最小值所在处以 max 及 min 文字符号表 示 a:u will see the max. & min. node # in legend. q:我在做一个例子时发现，在一个 input 文件中不能定义两种不同的材料属性。后来我在 keywords manual 中关于*material 的用法中看到了以下的文字： required parameter: name set this parameter equal to the name to be used for the material when it is referenced in the element property options. material names in the same input file must be unique. material names adhere to the naming convention for labels (see ``input syntax rules,&#39;&#39; section 1.2.1 of the abaqus/standard user&#39;s manual and the abaqus/explicit user&#39;s manual), except that they cannot begin with a number. 这是不是就是说 abaqus 在一个算例中没法定义一种以上的材料这个局限性好像有点大了 吧很多问题都需要定义多种材料属性的！我印象里好像其 cae 中还可以分别给不同的部分施 以不同的材料属性，难道反而在程序流中不能适用吗还是我的理解有偏差 万望高手能指点 迷津！a:material name is unique. you could definie many different materials by different name. q:abaqus 里做静力分析时怎么加重力加速度在 inp 文件里的 card 是什么 如： *dload elset, grav, 9.81, 0., 0., 1. a:有疑问查一些*dload 卡就清楚了 *dload 和*cload 有什么不同 *cload : concentrated load, 集中負載 *dload : distributed load, 分佈負載 q:请问如何设置 abaqus 的默认路径 a:我是这样做的： 比如，工作目录在 d:abaqusdata 从开始运行 abaqus command。 c:&d: d:&cd abaqusdata d:abaqusdata&abaqus job=**** interactive or d:abaqsudata&abaqus viewer/cae 所有工作文件都保存在 d:abaqusdata 中。 对于 command 方式，一般就是在文件的目录下运行 abaqus 就可以！生成的文件全部在运行的 文件夹下面。 而对于 cae 方式也可以通过另存为...而指定路径的。2、昨天还好好的，今天开机发现 abq command 没法使用了。 提示“不是内部或外部命令....” 解决办法： 右键我的电脑－》属性－》高级－》环境变量－》****的用户变量－》path 项添加： x:abaquscommands（x 为 abq 的安装盘） ，添加时注意与前面的项用“；”分隔。 小弟刚从帮助文档选了个率相关的用户自定义材料例子， 用命令行提交后出现如下症状， 请各位 给把把麦，谢了先 c: emp&cd elastic c: empelastic&abaqus job=umatmsst3 user=umatmst3 input file : umatmst3 abaqus error: the following file(s) could not be located: umatmst3.inp abaqus/analysis exited with error(s). c: empelastic&_ there are sereval different ways to invoke the user subroutine. 1. you can separate you input file and subroutine.for example: input file: example.inp subroutine: umat.for you can use following command to run &abaqus job=a input=example user=uamt 2. put your umat inside you input file named &example.inp&, simply run with &abaqus job = aa input = example 谢谢 pjnb 的指点，请问怎样将用户子程序夹到 inp 文件里呀 命令里如有 user 参数，就会去调用相应的 user subroutine 〖问题解决〗 〖求助〗abaqus 如何在命令行中提交多个 job. [精华] 在一个 job 分析完了以后，自动分析下一个 job. 谢谢！ 命令行问题： abaqus 后面跟不同的参数的意义，在什么地方可以查到，翻手册没有找到。 －－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－ 谢谢 wjytj abaqus analysis user&#39;s manual 3.2.2 execution procedure for abaqus/standard and abaqus/explicit 如果提交的是 inp 文件，比如在当前工作目录下有 1.inp 和 2.inp 两个 inp 文件，需要调用 damage.for 这个子程序，可以建立如下 python 文件，我将其命名为 job.py，于是就可以 在 command 下运行（当然你也可以采用其他的运行方式，mdb.jobfrominputfile 这个 命令还可以有其他的参数，详细请查看 manuals）: abaqus cae nogui=job.py 或者： abaqus cae startup=job.py 前者不用打开 cae，比较节约系统资源。 #---------------------------------------from abaqus import * from abaqusconstants import * #---------------------------------------import job # create an analysis job for the model and submit it. myjob1 myjob2 = = mdb.jobfrominputfile(name=&#39;1&#39;, mdb.jobfrominputfile(name=&#39;2&#39;, inputfilename=&#39;1&#39;, inputfilename=&#39;2&#39;,usersubroutine=&#39;damage.for&#39;)usersubroutine=&#39;damage.for&#39;) myjob1.submit() myjob1.waitforcompletion() myjob2.submit() myjob2.waitforcompletion() #---------------------------------------如果提交的是 cae 文件，则应该相应的使用 mdb.job， 而不是 inp 文件时候的 mdb.jobfrominputfile，其用法亦可参见 manuals。 至于提交大量的文件或者远程提交，要用到 queue，并需要修改环境变量，比较复杂， 我也没试过，如有必要，请查阅如下章节： abaqus analysis user&#39;s manual 3.2.2 execution procedure for abaqus/standard and abaqus/explicit [转贴] you want to run several jobs in succession without entering the abaqus command line for each job. to do this you can create a batch file as follows. in windows create a batch file (*.bat) with the following: __________________________________ call abaqus job=filename1 int call abaqus job=filename2 int ... call abaqus job=filenamex __________________________________ make sure all input files are in the same directory as the batch file. the &#39;call&#39; command will submit the first run (filename1) and once that run is complete or stops (even if the run does not complete successfully, the batch file will iterate to the next run) then the &#39;int&#39; will tell it to go to the next line and call the next job. thus, you do not need the &#39;int&#39; on the last line since there are no more jobs. i create batch files in a text editor (*.txt file) then rename the file to a batch file (*.bat). then simply type in the name of your batch file in the dos shell and hit enter and it should run each job in the order that you told it in the batch file. an easier alternative, i learned of it damage.for 这个子程序应该不是用户自己编的吧 作为用户是否不必管它的内容是什么或者在什么地方而只管引用就行了 到哪里可以找到这个子程序 晕，这是我自己编的一个子程序而已，只是个示例。 如果提交 job 无需子程序就无需 usersubroutine 这个参数， 如果需要子程序，则将 usersubroutine 的值改为你需要的子程序名. 如何在 cae 中使用用户定义子程序 cae 中在 creat job 的 job manager 中的 general 中可以指定子程序；com 下 将用户子程序，aba_param.inc 或 vaba_param_dp.inc， 以及子程序调试所得到的.obj 文件放到同一工作目录下， 然后，运行 abaqus job=your_input_file user=your_subroutine 就行了。 在线帮助文件：http://sufi.nchc.org.tw:/。 3、当用到某个用户子程序时，用户所关心的主要有两方面：一是 abaqus 提供的用户子程序的 接口参数。有些参数是 abaqus 传到用户子程序中的，例如 subroutine dload 中的 kstep， kinc，coords；有些是需要用户自己定义的，例如 f。二是 abaqus 何时调用该用户子程序， 对于不同的用户子程序 abaqus 调用的时间是不同的。 有些是在每个 step 的开始， 有的是 step 结尾，有的是在每个 increment 的开始等等。当 abaqus 调用用户子程序是，都会把当前的 step 和 increment 利用用户子程序的两个实参 kstep 和 kinc 传给用户子程序，用户可编个 小程序把它们输出到外部文件中，这样对 abaqus 何时调用该用户子程序就会有更深的了解。 （子程序中很重要的就是要知道由 abaqus 提供的那些参量的意义，如下） 首先介绍几个子程序： 一．subroutine dload(f,kstep,kinc,time,noel,npt,layer,kspt,coords, jltyp,sname) 参数： 1． f 为用户定义的是每个积分点所作用的荷载的大小； 2． kstep,kinc 为 abaqus 传到用户子程序当前的 step 和 increment 值； 3． time(1),time(2)为当前 step time 和 increment time 的值； 4． noel,npt 为积分点所在单元的编号和积分点的编号； 5． coords 为当前积分点的坐标； 6． 除 f 外，所有参数的值都是 abaqus 传到用户子程序中的。 功能： 1． 荷载可以被定义为积分点坐标、时间、单元编号和单元节点编号的函数。 2． 用户可以从其他程序的结果文件中进行相关操作来定义积分点 f 的大小。 例 1：这个例子在每个积分点施加的荷载不仅是坐标的函数，而且是随 step 变化而变化的。 subroutine dload(p,kstep,kinc,time,noel,npt,layer,kspt,coords, 1 jltyp,sname) include &#39;aba_param.inc&#39; c dimension time(2),coords(3) character*80 sname parameter (pload=100.e4) if (kstep.eq.1) then ！当 step=1 时的荷载大小 p=pload else if (kstep.eq.2) then ！当 step＝2 时的荷载大小 p=coords(1)*pload ！施加在积分点的荷载 p 是坐标的函数 else if (kstep.eq.3) then ！当 step=3 时的荷载大小 p=coords(1)**2*pload else if (kstep.eq.4) then ！当 step=4 时的荷载大小 p=coords(1)**3*pload else if (kstep.eq.5) then ！当 step=5 时的荷载大小 p=coords(1)**4*pload end if return endumat 子程序具有强大的功能，使用 umat 子程序： (1) 可以定义材料的本构关系，使用 abaqus 材料库中没有包含的材料进行计算，扩充程序功 能。 (2) 几乎可以用于力学行为分析的任何分析过程，几乎可以把用户材料属性赋予 abaqus 中的 任何单元； (3) 必须在 umat 中提供材料本构模型的雅可比（jacobian）矩阵，即应力增量对应变增量的变 化率。 (4) 可以和用户子程序“usdfld”联合使用，通过“usdfld”重新定义单元每一物质点上传递到 umat 中场变量的数值。 由于主程序与 umat 之间存在数据传递，甚至共用一些变量，因此必须遵守有关 umat 的书写 格式，umat 中常用的变量在文件开头予以定义，通常格式为： subroutine umat(stress,statev,ddsdde,sse,spd,scd, 1 rpl,ddsddt,drplde,drpldt, 2 stran,dstran,time,dtime,temp,dtemp,predef,dpred,cmname, 3 ndi,nshr,ntens,nstatv,props,nprops,coords,drot,pnewdt, 4 celent,dfgrd0,dfgrd1,noel,npt,layer,kspt,kstep,kinc) include &#39;aba_param.inc&#39; character*80 cmname dimension stress(ntens),statev(nstatv), 1 ddsdde(ntens,ntens),ddsddt(ntens),drplde(ntens), 2 stran(ntens),dstran(ntens),time(2),predef(1),dpred(1), 3 props(nprops),coords(3),drot(3,3),dfgrd0(3,3),dfgrd1(3,3) user coding to define ddsdde, stress, statev, sse, spd, scd and, if necessary, rpl, ddsddt, drplde, drpldt, pnewdt return end umat 中的应力矩阵、应变矩阵以及矩阵 ddsdde ，ddsddt ，drplde 等，都是直接分量 存储在前，剪切分量存储在后。直接分量有 ndi 个，剪切分量有 nshr 个。各分量之间的顺序 根据单元自由度的不同有一些差异，所以编写 umat 时要考虑到所使用单元的类别。下面对 umat 中用到的一些变量进行说明： ddsdde（ ntens， ntens） 是一个 ntens 维的方阵， 称作雅可比矩阵， ， 是应力的增量， 是应变的增量， ddsdde i， （ j）表示增量步结束时第 j 个应变分量的改变引起的第 i 个应力分量的变化。通常雅可比是一个 对称矩阵，除非在“*user material”语句中加入了“unsymm”参数。 stress （ntens） 应力张量矩阵，对应 ndi 个直接分量和 nshr 个剪切分量。在增量步的开始，应力张量矩阵中 的数值通过 umat 和主程序之间的接口传递到 umat 中， 在增量步的结束 umat 将对应力张量 矩阵更新。对于包含刚体转动的有限应变问题，一个增量步调用 umat 之前就已经对应力张量 的进行了刚体转动，因此在 umat 中只需处理应力张量的共旋部分。umat 中应力张量的度量 为柯西（真实）应力。 statev （nstatev） 用于存储状态变量的矩阵，在增量步开始时将数值传递到 umat 中。也可在子程序 usdfld 或uexpan 中先更新数据，然后增量步开始时将更新后的数据传递到 umat 中。在增量步的结束 必须更新状态变量矩阵中的数据。和应力张量矩阵不同的是：对于有限应变问题，除了材料本构 行为引起的数据更新以外，状 态变量矩阵中的任何矢量或者张量都必须通过旋转来考虑材料的刚体运动。 状态变量矩阵的维数， 等于关键字“*depvar”定义的数值。 状态变量矩阵的维数通过 abaqus 输 入文件中的关键字“*depvar”定义，关键字下面数据行的数值即为状态变量矩阵的维数。 材料常数的个数，等于关键字“*user material”中“constants”常数设定的值。 props （nprops） 材料常数矩阵，矩阵中元素的数值对应于关键字“*user material”下面的数据行。 sse ， spd ， scd 分别定义每一增量步的弹性应变能，塑性耗散和蠕变耗散。它们对计算结果没有影响，仅仅作为 能量输出。 其他变量： stran( ntens) ：应变矩阵； dstran( ntens) ：应变增量矩阵； dtime ：增量步的时间增量； ndi ：直接应力分量的个数； nshr ：剪切应力分量的个数； ntens ：总应力分量的个数， ntens ndi nshr = + 。 使用 umat 时需要注意单元的沙漏控制刚度和横向剪切刚度。通常减缩积分单元的沙漏控制刚 度和板、壳、梁单元的横向剪切刚度是通过材料属性中的弹性性质定义的。这些刚度基于材料初 始剪切模量的值， 通常在材料定义中通过“*elastic”选项定义。 但是使用 umat 的时候， abaqus 对程序输入文件进行预处理的时候得不到剪切模量的数值。所以这时候用户必须使用 “*hourglass stiffness” 选 项 来 定 义 具 有 沙 漏 模 式 的 单 元 的 沙 漏 控 制 刚 度 ， 使 用 “*transverse shear stiffness”选项来定义板、壳、梁单元的横向剪切刚度。 几个关于子程序的问题及相应解答 q: 本人在用 umat 作本构模型时， *static, 1,500,0..1 此时要求的增量步很多，即每次增量要很小， *static 1,500 时，在弹性向塑性过度时，出现错误，增量过大，出现尖点. a: you can try as follows: *step,extrapolation=no,inc=2000000 *static 0.001,500.0,0.。 q: 在 abaqus 中，如果采用 umat,利用自己的本构，如何让 abaqus 明白这种材料的弹塑性应 变，也就是说，如何让程序返回弹性应变与塑性应变，好在 output 中输出，我曾想用最笨 地方法，在 uvarm 中定义输出，利用 getvrm 获取材料点的值，但无法获取增量应力，材料常 数等，研究了帮助中的例子，umatmst3.inp,umatmst3.for,他采用 mises j2 流动理论，我 在 output history 显示他已进入塑性状态，但他的 pe 仍然为 0！ ！ a: 用 uvar( )勉强成功 。q: 偶在 umat 中调用求主应力函数 call sprinc(stress,ps,lstr,ndi,nshr) 后，存储主应力得数组 ps 中 各个主应力排列顺序是什么 ps1&ps2&ps3
ps1&ps2&ps3
ps1&ps3&ps2
a: 第二个 。 个人觉得： umat 实现自己的本构没有固定的方法，对于不同的本构有可能必须采用不同的方法。这要 靠自己不断地摸索。有可能一种方法对于简单加载问题还行，但有可能对于复杂问题并不 收敛。最重要一点，就是 umat 中采用的算法必须 consistent.再就是 ddsdde 必须正确， （如 果采用 back_euler 方法等一些算法，ddsdde 错误有时不影响结果（对于简单加载问题没有 影响，能收敛，，但对于复杂问题不收敛。 ） uptonow，你这个算法对于 mises,hill,j2，j2d 等一类的屈服函数是正确的，但具体的本构 还要灵活运用，这我也正学习，正在摸索。 有时，umat 需要很强的有限元基础，并且对采用的本构要很熟悉，不要在一颗树上吊死才 好。首先要确认自己的 umat 没有错误，如果没有，但就是不收敛（在不断减小加载步长的 情况下，当然最好对步长不敏感，特别是对于粘弹性，粘塑性，内变量一类的材料，有的 本构取决于背应力的计算） 那就应该考虑换一种算法。 。 一点体会，请大家探讨。 q: abaqus-uamt 的老问题，缺少 &#39;aba_param.inc&#39;文件
a: 在 cvf6.5 调试时，显示 缺少 &#39;aba_param.inc&#39;文件！ 这个没有任何关系的，这个错误将在 abaqus 调用 umat 的时候自动会找到，仅仅有这个错误 将没有任何影响的。也就是说，abaqus 中调用的时候，实际并不存在这个错误。ft，忘了 应该将 说一句了， 你把 aba_param.inc.dp 或 aba_param.inc.sp 拷到你的程序工作空间后 ， 把 aba_param.inc.dp 或 aba_param.inc.sp 的后缀.sp 或.dp 去掉， 即将 aba_param.in c.dp 或 aba_param.inc.sp 改名为 aba_param.inc。 呵呵，他的意思是在 visual fortran 中调试其子程序，我觉得这是一个好办法，我当时也 是这麽办的，毕竟在 abaqus 中调试是非常麻烦的，只有当你的 umat 没有语法或者明显的逻 辑错误，你在 abaqus 中调试才能事半功倍。 q：uvarm 可以输出到哪里.odb 可以么 另外那个 strav就是自己定义用于 umat 的那个数组里的数可以输出到.odb 里么 a：在 umat 中，statev 是不能用在 output 中的，statev 只是作为一个解的 状态变量，说来惭愧，我是在 umat 中定义 peeq,(peeq)的求解一般在 弹塑性力学书上有（等效塑性应变） ，用 write()写入一个临时文件， （ （切记：这个文件 unit 号不要与 abaqus 中的重合，因为他有一些系统默认 的文件号，） ） 然后在 uvarm 中读取，以 uvarm 输出，因为 uvarm 可以以 odb 的形式输出， 支持 output,field,output,history q：在本版看了一个一维固结的例子，其中含有用户子程序，如下subroutine ufield(field,kfield,nsecpt,kstep,kinc,time,node, coords,temp,dtemp) include &#39;aba_param.inc&#39; dimension field(nsecpt),time(2),coords(3),temp(nsecpt), dtemp(nsecpt) kfield=1 field(1)=coords(2) return end 我在 vf6.5 中进行调试，提示找不到 aba_param.inc 请问大侠这如何解决。 还有我因为是初次接触用户子程序，我查阅了本版所有的相关贴子，都讲的不太详细，我将问题 总结一下， 大侠们能不能详细的讲解一下， 1 子程序格式(程序后缀是.f; .f90; ..obj) 2 cae 中如何调用，command 下如何调用 3 若有多个子程序同时存在，如何处理 4 我对 vf 不是很熟，是否可以用 vc，c＋＋编写子程序 a: 若要在 vf 中调试，那么应该根据需要把 site 文件夹中的 aba_param_dp.inc（双精度）或 aba_param_sp.inc（单精度）拷到相应的位置，并改名为 aba_param.inc 即可。 1。我试过，.for 格是应该是不可以的，至少 6.2 和 6.3 版本应该是不行，其他的没用过，没有 发言权。 在 abaqus 中，运行 abaqus j=jobname user=username 时，默认的用户子程序后缀名是.for （.f,.f90 应该都不行的，手册上也有讲过） ，只有在 username.for 文件没有找到的情况下，才会 去搜索 username.obj，如果两者都没有，就会报错误信息。 如果 username 包括扩展名 for 或 obj， 那么就根据各自的扩展名 abaqus 会自动选择进行操作。 2。 中在 creat job 的 job manager 中的 general 中可以指定子程序； cae command 下用命令:abaqus j=jobname user=userfilename （无后缀） ； 3。将其写在一个文件中即可，然后用一个总的子程序调用（具体参见手册） 4。据说 6.4 的将可以，6.3 的你可以尝试着将 vc，c＋＋程序编译为 obj 文件，没试过。 在你的工作目录下应该已经存在 ufield.obj 和 uvarm.obj 这两个文件（这两个文件应该是你分别 单独调试 ufield.for 和 uvarm.for 时自动编译生成的，你可以将他们删掉试试看） ，但是由于 你的 for 文件中已经有了 uvarm 和 ufield 这两个 subroutine， 显然会造成重复定义， 请查实。 q: 假定采用 mises 屈服准则。 1 在调用 umat 之前，abaqus 传递给 umat 本次增量开始时的应力 sigma(0)，总应变 e,应变 增 量 delta(e)。 状态变量保存：弹性应变，塑性应变，等效塑性应变。（请问，状态变量保存的弹性应变＋塑性应变是否等于 abaqus 传递给 umat 的总应变 ） 2，然后在 umat 中利用上述的 abaqus 传递的量和状态变量得到 ddsdde 矩阵，然后返回给 abaqus，abaqus 根据 delta(sigma)=ddsdde*delta(e)，并且得到本次增量结束时的应力 sigma(1)=sigma(0)+delta(sigma)=sigma(0)+ddsdde*delta(e) 3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下次调 用 umat 请问手册上 umat 必须更新应力，可是根据上述我的理解好像是 abaqus 根据 umat 提供的 本 次增量的 ddsdde 在 abaqus 中更新，请问到底是怎么回事谢谢！ a: 1 在调用 umat 之前， abaqus 传递给 umat 本次增量开始时的应力 sigma(0)， 总应变 e,应.. : 量 delta(e)。 : 状态变量保存：弹性应变，塑性应变，等效塑性应变。 2，然后在 umat 中利用上述的 abaqus 传递的量和状态变量得到 ddsdde 矩阵，然后返回给 : abaqus，abaqus 根据 delta(sigma)=ddsdde*delta(e)，并且得到本次增量结束时的应力 3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下 以上有些 answers 本人并没有亲自证实，如有问题请及时更正！ . 市面上最多用户 cae 软件---可能是 ansys 有关市面上各种 cae 软件用户究竟有多少 很难评估，但可以肯定不管什么用户(合法或非法用 户)，我个人统计、调研(包括网上论坛、文章发表等)可以肯定 ansys 软件的用户最多。其主要 原因是: 1. ansys 软件自带的前后处理功能较强; 2. ansys 有教育版&＝2000 节点(流体)和 &=1000 节点(固体)、大学版(节点 1600 和 3200)以及商用版(无限制)；3. ansys 软件涉及的面 较广(应力场、温度场、流场和电磁场、优化设计、拓扑优化设计、随机有限元等)。ansys 软件 的高校教育计划做的比较早、比较好。 2. abaqus 用户较少的原因 abaqus 软件进入中国，在我的记忆中，也大约在 90 年代左右，它本身没有前后处理，国内其 前后处理主要借助于 patran 软件，这也是发展其用户致命的弱点，到 90 年代中后期，才开始其 前后处理软件的开发，也就是现在的 abaqus-cae/viewer,但是有某些原来的参数引入不进其 cae 中，些参数主要是涉及到岩土工程中的某些参数，abaqus5.8.14 的命令行可以全部引入， 但其 cae 中不行，后来我一直跟踪 abqus6.0~6.3.3 某些参数还是没有解决，只有回归老版本。 但对于机械行业、板筋成形方面一点没有问题。abaqus 在北美的石油行业的研究院所和石油专 业的院下均有，他们主要是“租”软件，而不是“买断”软件，只要每年交租金，那么均可以免费升 级。而我国主要是“买断”某一版本，要升级就在交升级费。 本身 abaqus 软件的开发一开始就是基于高度非线性问题，其理论性较强，专业性较强，要求 用户背景知识的起点较高，加上“以前”又没有自己的前后处理功能，其微机版的问世也比其它软 件晚，有没有高校培养计划……因此，导致其用户较少。 3. 单从在结构方面的应用 abaqus、marc、adina 和 ansys 的功能 1).如果做接触问题，选择软件的顺序为 abaqus、adina、marc 和 ansys，接触问题本身就是一 个高度非线性问题，前三者本身就是基于高度非线性问题而开发的，从建立接触对(因为接触对 中按材料硬度可分硬－硬、硬－软、软－软，如果相同硬度，那么那个接触体谁大、谁小， 那个是凸面、那个是凹面等－&来确定谁是接触面、谁是目标面等考虑)的方便程度和收敛程度为 以上顺序。 2).如果对结构要做结构优化设计或拓扑优化设计，那么 ansys 最强，ansys 软件中直接有优化 设计模块，是单目标优化设计，设计变量有结构尺寸变量和状态变量(如某些地方的某种应力不 能超过某一值，或某一变形不能超过多少)，优化结构变量写入 apdl 程序中，如果对 apdl 程 序不是很熟悉，那么可以通过 ansys 软件界面菜单完成建模和目标变量和设计变量设置，然后 把所有操作过程写入*.log 或*.lgw 文件中，它们是文本文件，以 apdl 程序保存的，用记事本 等调出此*.log 文件进行整理，整理出循环迭代结构，另存文件名，在菜单中执行优化模块时， 直接调此文件，一次性优化出结果。其它几个软件中没有结构优化设计模块，但也可以通过自己 编写个小程序，用 marc、adina 和 abaqus 对结构进行优化设计，但首先要熟悉如何取某节点 或某单元的结果数据，使其在设计范围内寻求最优。 3).如果从界面菜单上建模方面来讲，目前 adina、abaqus 与 ansys 旗鼓相当，marc 最弱，甚 至前两者超过 ansys 软件的建模，adina-m 和 abaqus/cae 的建模方式是基于现代 cad 的建模 方式(如类似 pro/e、ug、solidwork 等，其蒙皮技术、复杂曲面扫描技术远强于 ansys)。 4). 如果从编程序建模，那么 ansys 最强，因为它有自己的 apdl 程序语言，所有结构尺寸都可 以参数化， 这也是其率先开发结构优化设计和拓扑优化设计模块的基础。 marc 也有一个 python, 但很不好用。adina 可以在 adina-in 准备文本模型文件，但不能设置变量参数，可以通过文本 编辑处理模型数据。abaqus 与 adina 一样，可以编辑输入模型文件参数。 5). 如果从结构网格划分的方便程度来讲(这里不指自由网格划分)，设置网格线、面、体的分段 数和质量较好的映射网格方面，这几个软件的排序是 abaqus、ansys、adina 和 marc。 6). 从用于教学方面，只有 ansys(&2000 节点)和 adina 有教学版(900 节点)。 4、请问各位大侠 我想将计算出来的结点的温度值直接输出 成 ascii，能直接编辑，如何在 output 中实现 signalat wrote: *node file 这样生成的文件不是文本文件，只能供第三方后处理软件打开。 文本文件用*node print 或者*el print。 对，使用*node print 或者*el print，结果会写入*.dat 文件，可以用文本文件打开或者直接 改后缀为.txt。 我在 abaqus/cae 里的 command window 里写： abaqus job=****** 然后回车，总是说语法错误 还有，我看帮助里说可以设定内存，我写了这样的一句，也不好用 abaqus job=****** memory= 大家教教我吧，谢谢了！ 开始--〉运行--〉abaqus--&abaqus command&abaqus interactive job=tnnd 我觉得你应该先看看你的 inp 文件是不是在 abaqus 的 temp 文件夹中,如果在的话就应该没问 题。命令是对的。 abaqus 的 temp 文件夹在哪里啊我没找到。 ekulaoxue wrote: 我觉得你应该先看看你的 inp 文件是不是在 abaqus 的 temp 文件夹中,如果在的话就应该没问 题。命令是对的。一句话害死人 是 abaqus 的工作目录，安装的时候会指定，安装完了也可以指定，自己查以前的帖子吧。 右键你的 abaqus 快捷方式应该就可以看到的吧。 可以在 cae 下边直接弄吧 ，为啥都用 command 啊 temp 文件夹如果你安装的时候没有指定路径，就应该在 c 盘的根目录下。 我在使用自定义单元 uel 时, 在 inp 文件中把一个单元集定义为 rigid body. 之后在运行 inp 文件时, dat 文件中出现了如下的提示: ***warning: only a uel property option may be used with element 700. check that no other property refers to this element. ***warning: 14 elements have incorrect property definitions. the elements have been identified in element set warnelemincorrectproperty. 请诸位指教一二! “element set warnelemincorrectproperty” 是何含义 wjytj wrote:单元属性定义有问题，单元 700 只能够赋予 uel 属性，请检查是否重复赋予了其 它的属性。 这么简单的英语不会看不明白吧 wjytj wrote: 单元属性定义有问题，单元 700 只能够赋予 uel 属性，请检查是否重复赋予了其它的属性。 这么简单的英语不会看不明白吧 谢谢你的答复! 我将 14 个 uel 单元定义为一个 rigid body. abaqus 运行中突然停电，重新启动能否进行 比如说，在运行到 step 5 的 incriment＝125，停电，restart 可不可以从 incriment＝120 启动每一步计算结果默认都是保存的。运行过程中并没有随时 点击保存。 主要是我的程序要运行了 1100 多步，时间是 12 个小时，想重新运行时间浪费太多。 完全可以，在 restart 文件中有中断前的重启动文件就行！ 你这样说，我觉得有问题，至少突然停电，还有*.lck 文件存在，有这个文件存在，可以重新启 动 warning: the database has been opened with readonly flag on. it will remain readonly.the job input file &s-restart.inp& has been submitted for analysis. error in job s-restart: error in unknown rdm operation attempt to access a package before allocating it job s-restart: solver input file processor aborted due to errors. error in job s-restart: solver input file processor exited with an error - please see the data file for possible error messages if the file exists. job s-restart aborted due to errors. “ rdm” operation 是什么意思 你看就出现了问题，我不知道改如何处理。谢谢指点。 可以在 dat,log,msg 中找到中断的增量步，假如为 step5 中的 inc=125 中断，创建一个 重启动 文件，加入 *restart,read,step=5,inc=125 如果原文件是 1.inp，重启动文件是 2.inp，则在 command 里运行 abaqus job=2 oldjob=1 显然，可以从任意个 step 任意个 inc 开始重启。5、本人在进行热传导分析中,由于外界原因停电,分析中断.现在想继续分析可是却找不到*.res 文件.所以我怀疑默认情况下不会输出*.res 文件文件,不知道这样理解对不对.我的这个问题还有 没有别的方法可以重新启动分析呢by default, no restart information is written for an abaqus/standard analysis and restart information is written only at the beginning and end of each step for an abaqus/explicit analysis. you can use the step module to change the frequency at which restart information is written. 6、在关于 abaqus 的 user subroutine 学习材料中, 我看到这样的一句话: all abaqus common block names start with the letter “c,” so name conflicts can easily be avoided. 我想问一下 abaqus common block 是何含义 请指教! fortran 语言中的数据块公用区, 用来在程序单元间传递和共享数据. 公用区有&有名公用区& 和&无名公用区&之分, abaqus 中已定义的有名公用区的名字都以&c&开头. 现在明白了! 比如说我看到在一个程序中有 include ‘common-files.inc’ 这么一句话, 在后面有 一个 common_block.inc 的程序, 其中部分内容为 common/spline /igc(25),spgc(6,5) common/ntube /n,m,nunitcell,acc,a,totatom,tl common/cpot /b1,b2,b3,beta1,beta2,beta3, & a1,alpha,q1,dmin,dmax 这个 common_block.inc 就是在程序单元间传递和共享数据! 我只是提前没有想到 abaqus 的 inp 输入文件也可以调用 common_block.inc 以实现在程序单元间传递和共享数据.利用 common 语句可以在程序单元间传递和共享数据并且节约内存, 我知道这一用法, 但是没 有好好用过, 现在我明白了! 8、我首先安装了 fortran6.5, 然后按照 abaqus 6.4 的安装指南，建立了一个个所需的环境 变量，一步一步在进行安装，但是在安装过程中出现如下提示： there was a problem constructing the value for the abaquslm_licence_file parameter to be entered into your abaqus _v6.env file. please enter the value for this variable exactly as it should appear in the environmental file. 我输入 27000@&my computer name& (27000@bc029219)之后，安装能够继续进行，但是在最 后一步 verification 时，与 user defined element 相关的几个部分不能通过。 后续如果利用 user defined element 功能时，出现了问题, 要么是 df.exe not found in path, 要 么是 cannot open input file &oldnames.lib&, 总之 abaqus error: problem during linking abaqus/standard user subroutines, abaqus/analysis exited with errors. 还是环境变量的原因，有没有建立 lm_license_file 这个环境变量 9、我是 abaqus 初学者，但要用到用户单元，我先把安装文件夹中的例子程序 uelnonli.f uelnonli.inp 编译、链接，生成了结果文件，可是我打开结果文件，用户单元根本没生成！不知 道是怎么回事 我又把 uelriksx.f uelriksx.inp 运行，结果还是没生成用户单元！请大家帮帮忙啊！ 还有两个问题： 1.用户定义单元能在 cae 界面中调用吗我没找到！ 2.cae 中能直接定义节点坐标吗如何定义 2.cae 中应该不能直接定义节点坐标，cae 的建模方式和其他软件不太一样，part 中的模型实 际上只是相对坐标，最后还要在 assembly 中加以定位...10、在运行 abauqs 6.4 时，经常 abaquslm.exe 出现问题，提示程序不能进行初始化，关闭提 示窗口，后又自动打开，最后出现内存资源不足，且 lisence 也提示连接不上，不知为何请帮 助解决此问题。 很可能是端口被占用了 你找个能看端口的软件看看 把 lmtools 设为 start server at power up 去掉, 然后重新启动，把那个 abaquslm,lmtools 删掉， 重新把 abaqus 光盘里的 crack 里的 abaquslm,lmtools 复制过来，重新设定一下，每次用的时候 手动 start server 我以前也老出现窗口不能初始化提示，而且发现进程里开了无数个 lmtools,最后资源耗尽，机器 差不多被拖死 11、 q ：a）畸变网格怎么划分 我目前采用的是 icem 划分网格的方法，在 abaqus cae 中如何实现不知道，直接修改 inp 文件肯 定可以，但是 inp 文件修改真的是觉得比较麻烦。 b） 、断裂参量如果计算 ki 和 j 怎么计算cae 中有办法实现 a：首先你必须有 abaqus65 在 65 版前，abaqus 的 cae 不支持断裂特性计算装好后 在的定义接触、连接的那部分模块可以定义裂纹前缘、或者体内部的隐含裂纹 进入后，直接在菜单里面可以找到 crack 子菜单，定义奇异单元、隐含裂纹都在那里 这样 mesh 完后，会得到奇异单元，但是要注意，裂尖要用三角形或三棱柱的单元，这样系统才 会智能的给你产生所需奇异单元 定义好裂纹轮廓后 进入 step 设定那里 直接定义历史的输出，可以找到定义 jk g 的选项。这样定义后就行了 注意塑性材料不能定义 k 什么的，弹性材料不能定义。 12、 from eng-tips q: i need to draw contours of stress ratio(mises/pressure).abaqus can draw individual stress component like s11, s22. but i need to draw the ratio of the stress components.does any one know how to implement this into abaqus a1 don&#39;t know whether abaqus can do this or not, but this is possible in other softwares which allow the user to specify a result entity as an equation that can include any of the standard stress components and displacements plus geometric position. for instance roshaz can import abaqus ascii fil files from solid models and display expressions like:vm/(s11+s22+s33) or vm/(p1+p2+p3) etc..... where vm is the von mises stress and s11,s22 and s33 are the three direct stresses in global axes, p1,p2 and p3 are principal stresses. equations can have multiple nested brackets and contain standard geometric functions like sin, cos, sqrt etc.. a2: the input pressure at a face will be one of the stress components if the face is orthogonal to one of the global directions. using python script you could create a new variable of vm/s11 for example, and plot that. a3: an alternative to the previous suggestions is to define user output variables for the material(s) in question, e.g. *material, name=mat1 *elastic , type= iso
, 0.30 ** etc ** and within the same material properties block: *user output variables 4, ** *material, name=mat2 *elastic , type= iso
, 0.30 ** etc ** and within the same material properties block: *user output variables 4, the above would assign space for 4 user output variables at each element integration point for these two materials. you then write a short uvarm user subroutine (dead easy: see below). this allows you to call the getvrm routine at each element integration station and obtain values of any element integration point variable (eg any stress component). in your uvarm subroutine you then calculate the 4 user output variables. for example: subroutine uvarm( + uvar,direct,t,time,dtime,cmname,orname, + nuvarm,noel,npt,nlayer,kspt,kstep,kinc,ndi,nshr) c check manual to ensure that above is ok! c include, dimension and character statements as c described in the manual - check these are correct! include &#39;aba_param.inc&#39; character*8 cmname,orname,flgray(15) dimension uvar(nuvarm),direct(3,3),t(3,3),time(2) dimension array(15),jarray(15) c call getvrm (&#39;s&#39;,array,jarray,flgray,jrcd) s11=array(1) s22=array(2) s33=array(3) s12=array(4) call getvrm (&#39;sinv&#39;,array,jarray,flgray,jrcd) smises=array(1) spress=array(3) uvar(1) = s11/s33 uvar(2) = s22/s33 uvar(3) = s12/s33 uvar(4) = -spress/smises c remember to check for zero divides! return end the user output variables are stored along with other element variables at the integration stations. user output variables can be printed to the .dat file or plotted like any other element variable, e.g.*el print, s11, s22, uvarm1, uvarm2, uvarm3, uvarm4 13、请问 abaqus 中剪应力的正负符号如何确定 和一般材料力学中应该并无区别：使单元面顺时针旋转为正，逆时针旋转为负 14、 我最近在算一个题目,一个大筒体上有三个接管端面固定,大筒体两端加载扭距,但是不知道如 何加载.请各位高手能赐教 1. 将大筒体两端要施加扭矩的节点分别定义为两个 nset:left, right. 2. 分别在大筒体两端的圆心处定义两个 reference node: rp-left,rp-right. 3. 用如下命令将两个节点集绕 3 轴旋转的自由度与参考点耦合起来,其他自由度度是否耦合根据 具体问题而定: *kinematic coupling, ref node=rp-left left, 6,6 *kinematic coupling, ref node=rp-right right, 6,6 4. 在两个参考点上施加绕 3 轴旋转的弯矩. 友情提醒:reference node 也有自由度,注意相应的边界条件. 15、例如一个简支梁,怎么加支座塌陷位移,然后计算结构变形 当然可以拉，就在第二个 step 中的边界条件里加，可以，一般约束条件位移为 0，在 step 里给 约束的支座一个位移就可以了，当然没问题 ，你还可以设定速度加速度等。 16、不好意思，目前想将弹塑性和粘弹性结合起来（abaqus 里的粘弹性模型不能和塑性结合） ， 粘弹性部分没问题，可对弹塑性的算法不清楚（好像要用什么 return-mapping 算法)， 。不知有 没有高手指点（我愿意合作）一下塑性模型 umat 该怎么编写，最好给出个实例，简单的双线 性模型甚至完全塑性模型都行。 abaqus verification manual 4.1.17 umat and uhyper ，umatmst3.f 这个看明白就会了. 16、材料的疲劳性能一般以单轴应力－循环次数的形式表示（s－n 曲线。此处不考虑基于断裂 力学的疲劳理论） ，应力随时间的变化也很有规律，如正弦波、方波或脉冲等。除此之外，平均 应力对疲劳性能的影响也很少考虑 （也即 r=smin/smax!=-1 的影响） 。但实际的应力状态多是多轴应力，应力变化规律性较差， 并且 r!=-1。 如何将实际的应力 （应力变化无规律， 多轴， r!=-1)和实验室测得的材料疲劳性能(应 力变换有规律，单轴，r=1) 对应起来，就构成了疲劳分析的基础和依据。 （1）平均应力影响的处理 如果有不同 r 值下的 s-n 曲线， 一般采用插值方法确定未知 r 值下的 s-n 曲线。 如果只有 r=-1 的 s-n 曲线，可采用如下的公式计算等效的应力（就是将 r！＝－1 的单轴应力转换为 r=-1 时的 单轴应力，即等效应力） ： （sa/se)+(sm/su)^n=1 ^为指数运算符。 其中，sa 为半应力幅值，se 为欲求的等效应力，sm 为平均应力，su 和 n 不同的取值，构 成不同的理论： theory su n -----------------------------------------------------------------soderberg yield stress (sy) 1 goodman ultimate tensile stress (su) 1 gerber ultimate tensile stress (su) 2morrow true fracture stress (sf) 1 ----------------------------------------------------------------（2）多轴应力转换为单轴应力 这个转换其实就是采用何种应力（或分量） 。只能有以下选择： ） von-mises 等效应力；最大剪应力；最大主应力；或某一应力分量（sx,syz 等等） 。有时也 采用带符号的 mises 应力，其大小不变，符号取最大主应力的符号，好处是可以考虑拉或压的影 响（反映在平均应力或 r 上） 。同强度理论类似，von-mises 等效应力和最大剪应力转换适用于 延展性较好的材料，最大主应力转换用于脆性材料。 （3）无规律应力的处理 本质上是从无规律的高高低低的等效单轴应力--时间曲线中提取出一系列的简单应力循环 （sa,sm）以及对应的次数。有很多种方法可以完成此计数和统计工作，其中又分为路径相关方 法和路径无关方法。用途 最广的雨流法(rain flow counting method)就是一种路径相关方法。其算法和原理可见“downing, s., socie, d. (1982) simplified rainflowcounting algorithms. int j fatigue,4, 31–40“。哪位高手 把这个方法做成一个类吧：，好不好 ） 经过雨流法的处理后，无规律的应力--时间 曲线转化为一系列的简单循环（sa,sm 和 ni,ni 为该循环的次数） 。这样就可以应用损伤累计理论 （miner 准则）计算分析了：sum(ni/ni) ni 为该应力循环对应的寿命（考虑 sa,sm，见上） 。 根据此和可以衡量一定循环次数后的安全系数，或者一定复杂应力循环相应的寿命等等。 目前商品化的疲劳分析软件多基于以上流程。 同时应当指出的是，疲劳分析是一个经验型的分析，还没有成熟完备的理论。 17、abaqus 的 remesh 中遇到的 shape=deformed 问题 查了一下精华区，wjytj 版主曾经回答过类似的问题，但当时提供的附件现在都没有了，所以只 好再次请教： 使用/standard 中的 remesh 功能，首先需要提取已变形网格的形状，按照 manual，使用下面的 命令：mdb.models[&#39;model-1&#39;].partfromodb(name=&#39;a&#39;,filename=&#39;job-1.odb&#39;)没有问题，能够提 取出未变形的网格，但一旦增加参数 shape=deformed 或 undeformed 后都提示出错信息 “there is no variable named deformed&。琢磨了一下午，也没找出问题，还请版主和各位高手 指教。 先给出 inp 文件或是 odb 文件。 没问题，附件是旧网格的 inp 文件，还请版主指教：我是在 cae 底部的命令窗口中直接输入提 取变形网格命令的，没有使用 script 文件，不加 shape 参数成功，加了就不行，很是奇怪。 将下面一段红色文字拷贝到写字板中，并另存为一个 rpy 文件，比如 job-1.rpy 注意将其中的 odb 路径该为你的 odb 文件的路径，我已在相应的地方用中文标出。 # # abaqus/cae version 6.4-1 replay file # internal version: -11.18.28 46457 # run by administrator on sat may 15 10:11:32 2004 ## from driverutils import executeoncaegraphicsstartup # executeoncaegraphicsstartup() #: executing &oncaegraphicsstartup()& in the site directory ... from abaqus import * from abaqusconstants import * session.viewport(name=&#39;viewport: 1&#39;, origin=(0.0, 0.0), width=185., height=128.) session.viewports[&#39;viewport: 1&#39;].makecurrent() session.viewports[&#39;viewport: 1&#39;].maximize() from driverutils import executeoncaestartup executeoncaestartup() #: executing &oncaestartup()& in the site directory ... mdb() #: a new model database has been created. #: the model &model-1& has been created. session.viewports[&#39;viewport: 1&#39;].setvalues(displayedobject=none) import part import regiontoolset import displaygroupmdbtoolset as dgm #: number of assemblies: 1 #: number of assembly instances: 0 #: number of part instances: 1 #: number of meshes: 1 #: number of element sets: 5 #: number of node sets: 4 #: number of steps: 1 # 注意这里的路径该为你的 odb 文件的路径 odb = session.openodb(&#39;f:/workstudio/abaqus/job1.odb&#39;) # 注意这里的路径该为你的 odb 文件的路径 #: modeling space is 2d planar, type is deformable. p = mdb.models[&#39;model-1&#39;].partfromodb(name=&#39;part-1-1&#39;, instance=&#39;part-1-1&#39;, odb=odb, shape=deformed, step=0, frame=10) p = mdb.models[&#39;model-1&#39;].parts[&#39;part-1-1&#39;] session.viewports[&#39;viewport: 1&#39;].setvalues(displayedobject=p) #: the part &part-1-1& has been imported from the mesh of part instance &part-1-1& on the output database. q = mdb.models[&#39;model-1&#39;].part2dgeomfrom2dmesh(name=&#39;b&#39;,part=p,featureangle=0) 然后在 command 下运行：abaqus cae replay=job1.rpy 即可得到你要的变形后的 orphan mesh 和相应的 part. 需要注意的是,由于比例的原因,odb 中显示的变形比实际变形要大得多,所以看起来得到的 part 变形很小. 变形后的 orphan mesh: 感谢 wjytj 版主，问题解决。参考您的 rpy 文件和 scripting manual，原因可能如下：我最初在 command line interface 下按照 example manual 中的格式直接键入了 openfromodb 命令行， 正确的做法应该在此之前首先键入 from abaqus import * from abaqusconstants import * odb=session.openodb(&#39;your path&#39;) 再次表示感谢。 18、.py 和.rpy 应该都是 python script 的文件吧二者有何不同呢 .rpy 文件时记录下的 cae 操作时对应的 python 语言命令，就是所谓的 replay 文件， 运行方式是: abaqus cae replay=your_replay_file.rpy .py 文件是你自己写的 python 语言或相应的宏文件，运行方法是： abaqus cae startup=you_script_file.py 或 abaqus cae nogui=you_script_file.py 按道理，将.rpy 文件改名为.py 文件也应该可以运行。 19、ansys 的模型能不能导入到 abaqus 中去 以前走过 ansys&--hypermesh&--abaqus 感觉不错 可以！先把 ansys 中的模型用*.igs 文件导出，在 hypermesh 中导入然后导出*.inp 或者*.sat 数据丢失比较严重,一些实体之间的接触关系没有了,装配体成了一个实体. 如果是单个零件,问题不是很严重. 不过,要处理好两个模型的坐标对应关系. 20、 〖求助〗手册里的一个例子怎么算不出来呢 目前在学用户子程序 usdfld，以下是手册 25.2.38 usdfld 中的一个例子 example: damaged elasticity model。 可我在计算时，.dat 文件提示 ***warning: all cards before the first keyword card are ignored ***error: no step definition was found input 文件里明明有 step 定义啊不知有没有人碰到过，版主给瞅瞅，是个什么问题，谢谢啦!! input file *heading damaged elasticity model with user subroutine usdfld *element, type=t2d2, elset=one 1, 1, 2 *node 1, 0., 0. 2, 10., 0. *solid section, elset=one, material=elastic1. *material, name=elastic *elastic, dependencies=1 ** table of modulus values decreasing as a function ** of field variable 1. 2000., 0.3, 0., 0.00 1500., 0.3, 0., 0.01 1200., 0.3, 0., 0.02 1000., 0.3, 0., 0.04 *user defined field *depvar 1 *boundary 1, 1, 2 2, 2 *step *static 0.1, 1.0, 0.0, 0.1 *cload 2, 1, 20. *end step *step *static 0.1, 1.0, 0.0, 0.1 *cload 2, 1, 0. *end step *step, inc=20 *static 0.1, 2.0, 0.0, 0.1 *cload 2, 1, 40. *end step user subroutine subroutine usdfld(field,statev,pnewdt,direct,t,celent, 1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer, 2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo, 3 laccfla) c include &#39;aba_param.inc&#39; c character*80 cmname,ornamecharacter*3 flgray(15) dimension field(nfield),statev(nstatv),direct(3,3), 1 t(3,3),time(2) dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*) c c absolute value of current strain: call getvrm(&#39;e&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, matlayo,laccfla) eps = abs( array(1) ) c maximum value of strain up to this point in time: call getvrm(&#39;sdv&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, matlayo,laccfla) epsmax = array(1) c use the maximum strain as a field variable field(1) = max( eps , epsmax ) c store the maximum strain as a solution dependent state variable statev(1) = field(1) c if error, write comment to .dat file: if(jrcd.ne.0)then write(6,*) &#39;request error in usdfld for element number &#39;, 1 noel,&#39;integration point number &#39;,npt endif c return end give a try *heading damaged elasticity model with user subroutine usdfld *element, type=t2d2, elset=one 1, 1, 2 *node 1, 0., 0. 2, 10., 0. *solid section, elset=one, material=elastic 1. *material, name=elastic *elastic, dependencies=1 ** table of modulus values decreasing as a function ** of field variable 1. 2000., 0.3, 0., 0.00 1500., 0.3, 0., 0.01 1200., 0.3, 0., 0.02 1000., 0.3, 0., 0.04*user defined field subroutine usdfld(field,statev,pnewdt,direct,t,celent, 1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer, 2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo, 3 laccfla) c include &#39;aba_param.inc&#39; c character*80 cmname,orname character*3 flgray(15) dimension field(nfield),statev(nstatv),direct(3,3), 1 t(3,3),time(2) dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*) c c absolute value of current strain: call getvrm(&#39;e&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, matlayo,laccfla) eps = abs( array(1) ) c maximum value of strain up to this point in time: call getvrm(&#39;sdv&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, matlayo,laccfla) epsmax = array(1) c use the maximum strain as a field variable field(1) = max( eps , epsmax ) c store the maximum strain as a solution dependent state variable statev(1) = field(1) c if error, write comment to .dat file: if(jrcd.ne.0)then write(6,*) &#39;request error in usdfld for element number &#39;, 1 noel,&#39;integration point number &#39;,npt endif c return end *depvar 1 *boundary 1, 1, 2 2, 2 *step *static 0.1, 1.0, 0.0, 0.1 *cload 2, 1, 20.*end step *step *static 0.1, 1.0, 0.0, 0.1 *cload 2, 1, 0. *end step *step, inc=20 *static 0.1, 2.0, 0.0, 0.1 *cload 2, 1, 40. *end step 谢谢关注。但是好像不行哎。.dat 文件提示： ***error: unknown keyword &nodepar&. the keyword may be misspelled, obsolete, or invalid. ***note: due to an input error the analysis pre-processor has been unable to interpret some data. subsequent errors may be caused by this omission ***error: unknown keyword &userdefinedfieldpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &depvarpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &boundarypar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &steppar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &staticpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &cloadpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &endsteppar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &steppar&. the keyword may be misspelled, obsolete, or invalid.***error: unknown keyword &staticpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &cloadpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &endsteppar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &staticpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &cloadpar&. the keyword may be misspelled, obsolete, or invalid. ***error: unknown keyword &endsteppar&. the keyword may be misspelled, obsolete, or invalid. ***error: in keyword *element, file &compare.inp&, line 4: non integer in element connectivity for element 1 ***error: value cannot be converted to an integer: 1par ***error: value cannot be converted to an integer: 1par ***warning: error(s) during part, instance and assembly processing. analysis is continued without the information. end processing part, instance, and assembly information *********************************************************** ***note: the part map file &d:abaqususdfldcompare.prt& is corrupt, does not exist, or cannot be opened. ***error: error during reading part map the program has discovered 19 fatal errors ** execution is terminated ** 兄弟，abaqus 手册中的子程序有点小问题，两个 call 语句丢了两个续行符而已， 很小的小问题而已，估计是 abaqus 的雇员打了个盹 我在下面已经标出来了，拷贝为 for 文件就没问题了。subroutine usdfld(field,statev,pnewdt,direct,t,celent, 1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer, 2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo, 3 laccfla) c include &#39;aba_param.inc&#39; c character*80 cmname,orname character*3 flgray(15) dimension field(nfield),statev(nstatv),direct(3,3), 1 t(3,3),time(2) dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*) c c absolute value of current strain: call getvrm(&#39;e&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, 1 matlayo,laccfla) eps = abs( array(1) ) c maximum value of strain up to this point in time: call getvrm(&#39;sdv&#39;,array,jarray,flgray,jrcd,jmac,jmatyp, 1 matlayo,laccfla) epsmax = array(1) c use the maximum strain as a field variable field(1) = max( eps , epsmax ) c store the maximum strain as a solution dependent state variable statev(1) = field(1) c if error, write comment to .dat file: if(jrcd.ne.0)then write(6,*) &#39;request error in usdfld for element number &#39;, 1 noel,&#39;integration point number &#39;,npt endif c return end谢谢老大，老大说的这个接续号，我在第一次试算时已经发现并改了过来，我在帖贴子时忘了标 出来。 如果子程序有误， abaqus 会提示用户子程序编译错误的， 应该不会}