参考文献/References:
[1] 陈国兴. 岩土地震工程学[M]. 北京:科学出版社, . CHENG Guoxing, Geotechnical Earthquake Engineering[M]. Beijing:Science Press, .[2] LYSMER J, KULEMEYER R L. Finite dynamic model for infinite media[J]. Engineering Mechanics Division, ASCE, -877.[3] 孙海峰, 景立平, 孟宪春等. ABAQUS中动力问题边界条件的选取[J].地震工程与工程振动,):71-76. SUN Haifeng, JING Liping, MEN Xianchun, et al.The selection of boundary conditions in dynamic problems by using ABAQUS[J]. Journal of Earthquake Engineering and Engineering Vibration, ):71-76.[4] 刘晶波,吕彦东.结构一地基动力相互作用问题分析的一种直接方法[J].土木工程学报,):55-64. LIU Jingbo,LU Yandong.A direct method for analysis of dynamic soil-structure interaction[J].China Civil Engineering Journal,):55-64.(in Chinese)[5] 刘晶波,谷音,杜义欣.一致粘弹性人工边界及粘弹性边界单元[J].岩土工程学报,):. LIU Jingbo,GU Yin,,DU Yixin.Consistent viscous-spring artificial boundary and viscous-sping boundary element[J]. Chinese Journal of Geotechnical Engineering,):.(in Chinese)[6] 刘晶波,王振宇,张克峰,等.考虑土-结构相互作用大型动力机器基础三维有限元分析[J].工程力学,):34-38. LIU Jingbo,WANG Zhenyu,ZHANG Kefeng,et al.3D finite element analysis of large dynamic machne foundation considering soil-structure interaction[J]. Engineering Mechanics,):34-38.(in Chinese)[7] 刘晶波, 杜义欣, 闫秋实. 粘弹性人工边界及地震动输入在通用有限元软件中的实现[J]. 防灾减灾工程学报, ):37-42. LIU Jingbo, DU Yixing, YAN Qiushi. The implementation of viscoelastic artificial boundary and the ground motion input in the general finite element software[J]. Journal of Disaster Prevention and Mitigation Engineering, ):37-42.[8] 费康, 张建伟. ABAQUS在岩土工程中的应用[M]. 北京:中国水利水电出版社, . FEI Kang, ZHANG Jiangwei. The Application of ABAQUS in Geotechnical Engineering[M]. Beijing:China Water Power Press, .
相似文献/References:
[1]方志,陆浩亮,王龙.世界地震工程,7.　Fang Zhi,Lu Hao-liang,Wang Long.An overview on research of dynamic soil-structure interaction[J].,7.
备注/Memo:
收稿日期:;改回日期:。基金项目:国家自然科学基金资助项目（）；西安市地下铁道有限责任公司《西安地质条件下地铁地下车站抗震性能研究》科研基金项目（D4-YJ-）作者简介:尹尚之(1991-),男,硕士研究生,主要从事工程结构抗震方面研究.E-mail:
工具/Tools
统计/Statistics
摘要浏览/Viewed77
全文下载/Downloads48
评论/Comments[转载]ABAQUS常见问题汇总&-&2.0版（15-18）
．选择单元类型
云木斋开光佛珠手串奉请官网：&咨询QQ：
mr_chj : Abaqus中国培训资料，单元选择
问：我对同一有限元模型在相同的冲击载荷下，分别用这两种软件进行了计算，都选用显式算法，可结果相差比较大，请问各位热心人，是什么使他们的计算结果有如此的差异呢（网格划分一样）？
答：估计你的ABAQUS模型用的是reduced
integration单元，比如C3D8R，这种单元只能察看单元积分点上的应力，节点应力则很不准确。可以自己做一个悬臂梁线性弯曲实例试试，和解析解或者C3D20R,C3D8I的结果比较一下。详见
Getting Started with ABAQUS 4.1.2 Reduced integration。
LS-DYNA 的单元不知是不是这样。
问：请教一下各位,我在Z向加载局部的均布荷载,采用一般静态分析，而计算出来的Z向位移比理论计算值小得多,在8204节点计算为0.17971mm,而理论计算应略为0.4449mm,相差太大,从应力分布看,没什么问题。请问是什么原因?（文件比较大，分成了两部分）
检查一下你的边界条件是否和理论计算完全相同（比如理论解是无限大区域，而你的模型还不够大）。
在应力集中部位应该再细化网格。
C3D8R单元只有一个积分点，精度较差（尤其是应力结果）。用C3D8I, 得到8204节点U3=0.18.
模型中没有塑性材料和hard接触属性，可以用C3D20R。
位移很小，没必要用NLGEOM, 就可以加快速度。对称问题，可以只取1/4建模。大部分区域应力很小，网格可以划粗些。
问：什么意思永远不用C3D8？
还是在哪种情况下不要用？不明白。||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent
N2hMKU"un-?d#v
我做的轴与套的分析都用这个阿，分析结果挺好的阿
Getting Started with ABAQUS& 4.1 Element
formulation and integration& 有比较的算例。
可以自己分析一个受弯部件圆角处的应力集中问题，就可以看到C3D8比C3D8I,C3D20R差很多。
C3D8R的节点应力结果也很差，只能看积分点上的应力，所以也不适于应力集中问题。
问：非线性的接触问题从面最好不用2次单元，所以我在桩与土的接触部分partition了一个狭长的区域，用作1次单元。由于1次单元毕竟不如2次单元来的精确，所以我在其他的区域都是用2次单元。cae提示两个区域的网格不协调，提交运算后初始地应力的步都算不过去。把单元都改成一致就能算，请问向这样的问题，怎么样能实现两种单元的过渡呢？
最好不要用两种单元，在其交界面上会产生额外的应力，反而更不精确。都用一阶单元，网格划细点就行了。
问：各位前辈们好，我在做的模拟是关于钢丝拉拔成型的问题，拉拔的速度是5米/秒，模具长为0.002米，也就是成型过程的时间为0.0004秒，不知道这样的问题按静态计算还是应该按动态计算呢？
如何使得二次单元和一次单元相连接
我的模型中既包括二次单元，也包括一次单元，二者需要相连，我不知道在abaqus
cae中能否实现，如何实现？研究了一下，将稽核图形生成meshpart，编辑单元和节点是可以的，但是如何将其和原有的几何图形统一起来呢，而且单元复杂，一个一个节点的编辑，不太现实？还望各位前辈多多指点阿！版主快来啊！我迫切需要你啊！
加载时间这么短，不能忽略惯性和冲击效应，要用动态分析。
二阶单元不能用于塑性材料和接触问题:
/forum/viewthread.php?tid=647602&page=1&extra=#pid729426
最好不要用两种单元:
/forum/thread--1.html
一定要二次单元和一次单元相连的话：如果是同一个part,
先partition，再在不同的区域设不同的单元类型；如果是两个part，可以用tie。
．截面属性和材料
定义截面属性
问：在已划分好的壳体单元中，计算时程序老是提示整个壳体中间有个单元没有定义属性，很不明白是怎么回事．是不是可以将这个单元删除，然后再补上，
如果是这样，具体应如何操作，请教大家了．
solong答：
楼主的描述有点太简单，不知道到底是如何……
我是初学者，貌似遇到过类似的情况，是因为赋值材料属性的时候疏忽，cae中用shell建立模型之后单元的类型没有调整，所以cae中虽然赋值，实际上却并没有把材料的属性付给单元
问：如何在cae里定义 *SHELL GENERAL SECTION,
我估计：在CAE里用shell section，然后ABAQUS根据你的section参数和单元类型，自动选择*SHELL
GENERAL SECTION 或 *SHELL SECTION。
问：當我要畫2d的平板時
我再定義完材料和斷面為shell時
到assign要定義section時
它會出現不是solid不可以定義
在创建part时，modeling space不要用2d planar, 而要用3d.
shell都是3维的，能弯曲才叫shell, 否则就是平面应力单元了。
问：请问大家：
1。线，面模型能导入abaqus吗，有一个空间网架，12000个杆件左右，用abaqus cae建模的话很复杂，我想用auto
cad建模再导入abaqus可以么。
2。网架的12000个杆件的截面不全相同，大概有10中截面形式，分布也没有什么规律。模型建立起来后我要赋予各杆section。如果在cae下操作，这12000个杆件我得干到什么时候啊。请问大家有什么好办法么。
建议：auto
cad建模时，分别建10个子模型，每个子模型只包括截面相同的杆，导入ABAQUS后就是10个part，各赋予一次截面属性就可以了。
用OFFSET参数偏置shell的中性面
问：如圖想請問,如果我建了一個薄殼的物件,如圖
並在中間部分給厚度0.1mm ,而其他兩邊0.3mm
(1) 那在abaqus裡面實際的意義到底是前視圖A還是B ?
(2) 假設是A,那我該如何設定才能變成B那樣的厚度分佈 ?
麻煩哪位前輩可以賜教,謝謝~
guoyj8答：
通常的直接结果是A情况，这时的壳中面是所定义的参考面。如果想定义成B情况，则需要应用OFFSET参数功能定义真实的中性面。对于你的情况：厚度分别是0.1和0.3，则需要应用中间部分壳的中面做参考面，对于两边部分，则应用OFFSET功能说明这时的参考面相对于中面的相对位置----如果向上为壳的正向，则OFFSET=-0.3333。这里的负号说明定义的参考面位于真实中面的负方向，而数值0.3333是中面和参考面之间距离跟板厚的比值(0.15-0.05)/0.3=1/3=0.3333
问：本人想通过壳体单元（shell)来模拟地基上的板，地基用实体单元（solid）。由于壳体单元的节点位于壳的中面上，而实体单元的节点位于单元表面，因此在装配后，壳体的下表面位于土中，这个会在后续的计算中存在问题吗，听说可以对壳体单元进行偏置，
tingzung答：
SHELL-TO-SOLID COUPLING中，只能shell的边跟solid的面耦合在一起~
另外可以在iput档中使用offset这个指令去偏移shell的中性面~
如果在CAE中，你用tie连接shell与solid的面，
打开Adjust slave node initial
position的话，abaqus就会自动调整slave那个面，使得shell的面与solid的面接在一起 shown in
如果要再实体上加一层薄壳可以用skin这个功能
ABAQUS/CAE property-special-skin
问：在我的模拟中，土体材料采用理想弹塑性模型，塑性部分采用 mohr-coulomb
plasticity,其中的hardening里的cohension yield
stress是随着土体深度的增加而增加的，如果土体深度设为Z,则cohension yield stress=1Z ( KPa ),
本例中，土体总深度为10米，则土体顶部cohension yield stress=0，底部cohension yield
stress=10( KPa ),同时，弹性部分的young's modulus与cohension yield
stress比例保持不变，始终为500，即young's modulus / cohension yield
stress=500，请教高手这个材料性质应该如何定义，
chenyu 答：
使用场变量，然后定义材料时指明参数依赖于该变量
使用自定义材料或场变量都可以，前者复杂，因此不推荐。
G"nuNvs*|仿真|设计|有限元|虚拟仪器仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM~
可将Z坐标定义成一个场变量，然后再定义材料时，指明某个或某几个参数（此处为cohension yield stress和young's
modulus ）依赖于该场变量（此处手动修改Inp文件比较方便，仅以弹性为例）：
&R&x&| s1L*
DEPENDENCIES=1仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
E+S&y-T3v&J:]^-B
&&100.E9, 0.3, ,
& 200.E9, 0.4, , 10.
#l1C8q"pRB$|仿真|设计|有限元|虚拟仪器计算时，当场变量为(z=)0时，abaqus查表确定模量为100.E9，场变量变为10时，abaqus查表确定模量为200.E9，中间的则线性插值。
OM5x,@ j-k%V a,[
x%y4~w)U$vd#]^
因为该场变量（积分点z坐标）连续变化，所以可以使用USDFLD子程序（25.2.38 USDFLD
r@?kL%b"|仿真|设计|有限元|虚拟仪器User subroutine to redefine
field variables at a material point.SimWe仿真论坛N1w5WB)gB'JAL"A
）来连续改变该变量。该子程序中，只需将FIELD(NFIELD)赋值为用变量COORD保存的z坐标即可。&
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentQ+m2N?`
我有一个小问题，你的模型变形后（z坐标变化后），还遵从那个关系吗？如果材料性能只是取决于未变形时的z坐标，可能得定义一个状态变量STATEV，在计算开始时将z坐标保存在STATEV中（IF
TIME .EQ. ZERO THEN STATEV(1)=COORD(?)
...），然后计算过程中，对每个积分点将FIELD赋值为该状态变量(ELSEIF FIELD(1)=STATEV(1) ENDIF)。
问：一个板，20层，各层材料方向不同
请问在abaqus怎么给各层材料属性，能直接加到网格上吗？
多谢高手帮助
在帮助里搜索"composite"，可以找到
ABAQUS/CAE User's
Manual&&&&
12.8.4 Creating composite shell sections
问：請問fluid（ex.水）的材料性質要怎麼設定阿...
我看user's manual 7.12的那個例題沒有給inp檔或CAE的教學...
ABAQUS Analysis User's
Manual&&&&
7.12.1 Modeling fluid-filled cavities
ABAQUS Benchmarks Manual&&
&&2.5.1 Fluid filled rubber
问：在置换时，同一个part要在不同的步骤用到不同的材料属性，
对于同一个part给出两个单元集，对于两个单元集定义不同的材料，第一步移除一个单元集，第二步激活这个单元集，同时移除另外一个单元集，这样我觉得比较合适
关于ABAQUS中粘弹性模型
想用burgus模型之类的模拟沥青混合料的变形仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent1EVc&hG16k6j@M
目前只知道时间硬化的那个蠕变模型*`XdrS
还不知道这些模型在哪仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM)C&
参数怎么输入
有两种方式可以定义，一种是定义粘弹性，通过creep test
data定义；还有一种是可以定义蠕变模型，通过creep来定义time－harding模型。
．载荷 / 速度场
施加荷载时出现错误,求教.
在load步中，在模型上新建一个reference
point，然后建一个只包括该点的set，把荷载加到这个set上，这样对么？SimWe仿真论坛k2Dj Q
o&r9v& & &
&&&为什么在job中data
check时会出错呢？
|+s‑P(^ r9d5P+B9cNode set assembly_set-1 has not been
&f5s+P'|$n*p q#u7cW
A concentrated load has been specified on node set
assembly_set-1.&&this node set is
not active in the model
要用coupling把reference point和受力区域连接起来，或者不用reference
point，直接在受力区域上加载。
为了便于施加集中荷载和添加beam element， 我在shell 面上定义了一些datum,
但是家荷载让我指定点时却不能选这些点，请问应该如何解决呢？谢谢
不能在datum上施加荷载，要通过partition，使shell在加载位置出现一个顶点。
问：有一个旋转体,如图,节点上定义集中力载荷.
第一个和最后一个节点上的力的分量均为10N,中间节点上的力的分量为20N.单元边长度为0.4714mm.
转化为均布载荷,即从第一个到最后一个节点的中间部分上施加Pressure,大小为20/0.MPa.
(这部分斜面为45度)
最后的结果,施加压力情况下的应力值比施加集中力大了40多倍.
类似的情况,使用plane strain单元,结果完全一样.
想问问,在轴对称情况下,集中力代表什么?沿周长方向的单元长度的力?
答：附件是一个更简单的圆筒模型，内径r=10, 外径 R=20, 在step 2 中只施加均布载荷 p = 0.4,
分析后得到的底面反作用力RF总共为376.8，可见总载荷为 p * 面积，即 p * (3.14 * R * R - 3.14 * r
在step 1 中只施加集中载荷，圆筒顶面三个节点上的力大小分别为1，2，1, 分析后得到的底面反作用力RF总共为4
，可见ABAQUS没有把集中载荷乘以面积。
另外注意到，施加集中载荷得到的变形是不均匀的。因此，在轴对称问题中应尽量施加分布载荷，不应施加集中载荷。
如果把单元类型改为平面应变单元CPE4I，边界条件变为底面固支，上述step 1和step
2得到的结果就是完全一样的（不像你说得那样）, 总的RF都是4。
还有一个问题，CAX4，CPS4,
CPE4，C3D8等完全积分单元过于刚硬，得到的结果精度很差，建议使用非协调模式单元CAX4I，CPS4I,
CPE4I，C3D8I。使用非协调模式单元时要注意，在所关心部位的单元各个角都应尽量接近90度，不要有尖角或钝角，否则会影响精度。
多谢回复.问题已经解决.关键是在轴对称问题中,集中载荷是沿圆周积分所得,也就是在圆周上总的载荷,而我开始把集中载荷理解为了单位长度上的载荷,这样在轴对称问题中,由于X的值不同导致圆周长度不同,从而在集中载荷与分布载荷中产生了一个差异.
在Analysis user's manual的14.1.6 Axisymmetric solid element
library开头convection中有关于载荷的解释:
Point loads and moments, concentrated (nodal) fluxes, electrical
currents, and seepage should be given as the value integrated
around the circumference (that is, the total value on the
在上面的模型中,step1采用原始的集中载荷1,2,1.在step2中,三个节点的X值分别是10,15,20,即它们的半径.因此需要计算单位长度的集中力,分别是f/(2*pi*X),得到f1,f2,f3.两个单元边长均为5,因此两个单元面上的压力分别为(f1+f2/2)/5,(f3+f2/2)/5,结果为0.7.计算结果比较一致,如果细化网格,计算结果更准确.
问：我在做钢筋混凝土构件的动力分析，目的是想得到荷载位移曲线。首先给构件施加力荷载，但是得出的位移曲线很不理想（比试验数据小很多），所以想改用施加位移荷载，用Amplitude定义位移荷载很自然就得到位移时间曲线，可是不知道如何得到荷载时间曲线。通过修改历史输出变量可以得出一个点的RF-时间曲线，可是该曲线太乱，根本无法应用。不知哪位能不能指导一下该如何得到荷载位移曲线呢？谢谢！
假如是在某个面上施加固支边界条件，就在这个面和一个参考点之间建立coupling，去掉此面上的边界条件，在参考点上定义边界条件。参考点上的RF就等于载荷。
问：我修改了模型给悬臂梁加了加强筋之后，应力相差就不是很大了。我能不能得出下面这个结论：给构件施加集中力时，为了避免受力点应力集中和扭曲可以把该点和该点所在的面耦合？
还是用均布荷载来代替集中力已达到避免受力点应力集中和扭曲的目的？
你说的这两种方法的效果应该是差不多的。集中载荷还是均布载荷应该根据工程实际来确定。
对于塑性材料，尽量不要只把载荷施加在一个节点上，否则可能因为局部变形太大而不收敛。
请问如何对一条边施加边载荷，谢谢！
是solid的边吗？在与边重合的位置定义弹性模量很小的梁单元，把梁和那条边tie在一起 ，Slave surface
类型为Node Region, 选中梁的所有节点 (如果Slave surface 类型为surface,
选中梁的单元，会出现异常结果). 在梁上施加line load (不是pressure！).
但是我定义了梁单元后，我要creat instance这个梁单元才能够tie 梁和solid
bar。但是当我creat的时候：
& w*`)w;@@`#G@
fAn error was reported "part dimensionality does not match the
dimensionalty of the existing instances."
创建梁的 part时，要把modeling space设为3D, 而不是2D planar.
问：在LOAD里加线载荷和集中载荷（不是在模型边上，在中间），始终不知道怎么加上去，请大虾指教！
集中载荷: 通过partition，使加载位置出现一个顶点。
或者用另外一种方法： 定义一个参考点，在边和参考点之间建立coupling约束，然后把力施加在参考点上。
tie和coupling都在interaction模块，菜单constraint / create.
我要在一個圓形薄管施力斷面施ㄧ拉力我的指令對嗎?
看来pressure是不行了，可以在interaction模块中定义一个reference point,
然后菜单constraint / create, 选coupling, 类型为distributing
coupling。把受力面和参考点连接起来，力施加在参考点上，见附件。
你用的是kinematic coupling，受约束的边无法再变形。应该用distributing coupling,
受约束的边就可以变形。见附件。为了看到均匀的应力场，边界条件我也改了一下。为看到变形，我加大了载荷。得到的应力场还是不完全均匀，可能是distributing
coupling不能使拉力完全均匀地分布。
我的方法也不一定对，你再查查帮助，看shell单元的面内载荷该怎么施加，用s4r单元是否合适。
位移加载时，加载面的平均应力怎么求？
因为最后给出的加载面的应力分布很不均匀，支反力也很不均匀
如何定义加载面的平均应力啊？？
用distributing coupling把加载面和一个参考点耦合起来，把力加在参考点上。
加载面的平均应力 = 参考点上的支反力 / 加载面的面积
用三维单元模拟1/4的空间轴对称问题，结果并不轴对称？
我在part里定义了*Orientation，并把该方向定义应用到了单元的section中仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentU:[Y-A9s#yl
在边界条件我还是用的是整体坐标系。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent‑U(G`7^YO;v
另外cae中odb文件中显示出的应该是：应力、应变以*Orientation定义的局部坐标系显示，而位移意整体坐标系显示
f2l'g0{||仿真|设计|有限元|虚拟仪器模型是一个空心桩，除其自重外，还考虑在外壁上作用静水地压和负摩擦力。在两个对称面上分别给定相应位移约束，在底部给出竖向约束
。但计算结果却不轴对称，不知原因在哪里？请予指点。1fx`on/c f$A
把静水压力改为普通的pressure，结果就轴对称了：
WAIBI, P, 2.34e+06
可以找找帮助里的例子，看看静水压力的用法。
我又作了一遍，用普通压力的确没问题。用静水压力不行，主要是因为方向的问题，abaqus默认的两点（定义静水压力为零和某定值的）坐标须为3轴，我上面的模型是通过sweep建立的，因此沿圆筒轴向为2方向。所以不对，后来我用*ORIENTATION的局部坐标方向，结果还是有问题，可能是我还没弄太明白。再后来，干脆重新建模用的是EXTUSION，结果圆筒轴向变成了3轴，不用再用那些麻烦的局部坐标系。谢谢。
s9v6R['N e#A
mH‑`@||仿真|设计|有限元|虚拟仪器但为了明白，还是把我用局部坐标方向等命令放上来，请给看看仿真分析,有限元,模拟,计算,力学,航空,航
*ORIENTATION 主要用来定义材料方向：
Define a local axis system for material or element property
definition, for kinematic coupling constraints, for free directions
for inertia relief loads, or for connectors.
用于载荷和边界条件的局部坐标系要用*TRANSFORM。
Body Force /
12:05, 29.35 K)
我直接计算的重力场怎么是这样的？两种方法上表面的力怎么区别这么大？什么原因？
让模型只包含一个单元，就可以很容易看出区别：
_PickedSet17, GRAV, 10., 0., -1., 0.
10 是重力加速度，合力为 密度2000 * 一个单元的体积 0.008 * 重力加速度 10 = 160 N.
body force:
_PickedSet16, BY, -10.
10 是单位体积上的力，和密度无关，合力为& 一个单元的体积 0.008 * 单位体积上的力
10 = 0.08 N.
所以你应该用gravity。
随时间变化的载荷 /
问：关于分级加载的问题
现在想将10000kN的水平集中荷载分10级加在桩的顶部，请高手给讲解一下分级加载的原理和方法，如果能有一个例子最好，再次谢过。仿真分析,有限元,模拟,计算,
施加集中载荷10000，step的time period为1，initial increment 和 maximum
increment 都设为0.1,&
即step分为10步完成。这样time=0.1时载荷就为1000，time=0.2 时载荷就为2000 ......
问：在Amplitude中定义一个随时间线性变化的力，起始点是（0，0），最高点是（0.5，1），最低点是（1，0），magnitude是1000，却怎么也是现不了啊！大家看看，我到低错在了什么地方，是不是step中的步长的问题啊？
changyj答：
论坛上有定义正弦载荷的情况，而且我用amplitude定义的随机载荷也可以，看看帮助文献，在论坛上找找。
在模型数据部分：
*AMPLITUDE,TIME=STEP
TIME,VALUE=RELATIVE,NAME=EQ,INPUT=＊＊＊＊.inp（你定义的曲线的文件名）
在历史数据部分：
*BOUNDARY,TYPE=displacement,AMPLITUDE=EQ
我定义的是边界条件，你可以参考定义载荷曲线。
问：在load中使用amplitude但是怎么好像amplitude不起作用啊？
你amplitude中的时间点是0.25, 0.75,& step中可以把increment
size的initial 和maximum 都设为0.25。
问：請問使用*AMP tabular設定完的三角波
循環5次要怎設定阿...
应该不难想出来啊。
方法1：在*Amplitude中多给几个数据点，把5次循环的每个顶点都列出来。
方法2：在*Amplitude中只描述一个循环，在模型中创建5个相同的step.
firce答：我知道我错在哪里了。把一个周期在amp中完整的描述出来。然后加上definition=PERIODIC，指定为循环，可能会好用
aba大侠的办法恐怕好用但不易于扩展的，如果要100个周期，写起来就太麻烦了
问：目前希望的情況是block先受壓力產生位移...
再以此位置為最高點進行循環負載...
可小弟設出來的結果在進入step2的時候壓力就不見...
且循環位移會以未受力的起始位置為amp=0進行...
我有試過更改起始的amp大小，不過會無法收斂...
所以想請教前輩有沒有可行的解決辦法...
在各个step中，位移边界条件和载荷中的数值都是相对于模型最开始状态的值，和上一个step的结果无关。所以应在定义amplitude时设置initial
amplitude。注意这时要把amplitude中的B设为5，边界条件中的位移设为-1。
在同一个节点（或参考点）的同一个自由度上不能即施加载荷，又施加边界条件，因为这二者所导致的位移不同，会产生矛盾。可以去掉载荷，在step
1中加一个位移边界条件，把橡胶压到initial amplitude的位置。或者去掉位移边界条件，只保留载荷。
接触面为什么不是橡胶顶面，而是侧面？我的模型见附件。
底部角上单元扭曲太大，位移大了就不能收敛，见Getting Started with
ABAQUS&&&&
10.8 Mesh design for large distortions。我把初始位移改小了。
问：我想在一些列节点上分别加上随时间变化的力，比如说，t1时刻，分别在node1,node2,node3上加上集中力f1(t1),f2(t1),f3(t1),t2时刻分别在node1,node2,node3上加上集中力f1(t2),f2(t2),f3(t2)。其中的f1
都是我通过matlab程序算出来的，需要从外部导入到数组里，而不是随时间线性变化或是时间的函数，也就是说不能在abaqus里面定义的。时间从t1到tn,加在每个节点上的力也是变化的，对这个结构进行时程分析，请问这个在abaqus里面能否实现，
可以使用用户子程序。
或者先用CAE生成一个没有载荷的inp文件，然后用c或fortran编程，来修改inp文件，用*cload为各个节点添加载荷。
小弟初学ANSYS &
ABAQUS，不是很熟悉。现一个大变形小应变分析的例子，真诚地请各位高手帮助解决。由衷地谢意！SimWe仿真论坛4JY
q%C{)W[1]D^B)c
E8s)C9y.hu%DHSimWe仿真论坛杨氏模量是 200e9 N/m^2 （可改变），泊松比是
0.3，密度7800。三边长度N1-N2 = N1-N3 = N2-N3 = 1.2 m.
M1，M2以及M3均为三边的中点，梁均为矩形截面，尺寸为宽 b＝2.5 mm 高h＝1.25 mm.
:w1rG$ie仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,#hG
要对左图所示的一正三角形进行大变形加载（力或者位移），使得结构最后变成右图所示形状。现有一种加载方式：
e5R't pX7C.W0CSimWe仿真论坛
~:C+s;~:QZ6jBSimWe仿真论坛把这一结构变形过程分为NN步，每一步变形为N3 & N2 分别沿垂直于平面P3
&P2的方向（两个方向是在同一直线的不同方向）发生位移，位移的大小与NN有关。P3
是由当前变形中边N3－N1和边N3-N2中靠近N3端的两小段的轴线确定的平面。P2是由当前变形中边N2－N1和边N2-N3中靠近N2端的两小段的轴线确定的平面。在有限元的实现过程中，两小段的轴线可以取靠近N3或N2的两单元的轴线。
h_iw仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
sm9{ap仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM问题是一个关于结构在非保守荷载（？位移）作用下的跌代求解的过程。关键有以下两点：
c,v仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1．变形应分成NN步计算，每一步的荷载（位移的大小）可以设为一致（方向肯定不一致），也可根据当前变形的大小择取。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM0k,Y:R%nc
2．每一步的荷载（位移）是定义一个局部坐标上的，这个局部坐标是定义在当前的结构变形之上的，随变形而变化的，所以是个非保守荷载问题（？）。A7dyz
SimWe仿真论坛a$u4hP7M
程序循环可以描述为：||仿真|设计|有限元|虚拟仪器*Ot&w#G9P9pv!I&f
$O3w&ED$s2b!pSimWe仿真论坛*do &
&!&&这个DO
循环可以在ANSYS很好的实现。
3~#]/I,p,vns3kj$?|仿真|设计|有限元|虚拟仪器SimWe仿真论坛Q ^Y
根据当前变形定义局部坐标& &
CX仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
-E o~6f8D|按照局部坐标施加荷载（位移约束）LOAD , SimWe仿真论坛0v&F*F ^
[&]^7p%snp仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM**
ABAQUS可以定义局部坐标－DATUM，但不知道能否在局部坐标下？
mm&yp6U&i%x仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent加荷载？ANSYS
可定义局部坐标吗？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentOl&s+MX
||仿真|设计|有限元|虚拟仪器p t‑HXD8x E求解
|仿真|设计|有限元|虚拟仪器~s[1]{)_q_JV;c.RW.k
/pW"kS[1]|S-g;T?更新结构变形&&UPDATE&&GEOMETRY!&&&D"d
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent/w3?:W
**在每一步求解中，程序是自动完成UPDATE的。我们需要的是让程序完成一步求解后，跳出来，重新定义荷载，再求解。&&**(这个不知道ANAQUS或ANSYS是否有能够实现？)SimWe仿真论坛
_5j@%jt‑h*J@;H0Q"z[
Rl:_S& & &
u;e:Z*|仿真|设计|有限元|虚拟仪器J+Z.K8`*J
nc[!d b yOSimWe仿真论坛其实，结构的变形也可以通过施加一个parametric, non-conservative
的力来实现，力的方向与以上位移一致，也是分别垂直与当前的P2&P3平面。注意力不是一直增加的。个人认为施加位移比较好控制。
YM+M)_1}||仿真|设计|有限元|虚拟仪器||仿真|设计|有限元|虚拟仪器@;r/R't8V
R"e;c8q0R8`
问题还是一个对称的变形问题。整个变形过程对称于N1-M3轴。所以可以取一半结构（即N1-N2-M3段）来计算。施加的边界条件为N1固定，M3可以有沿N1-M3方向的位移以及绕N1-M3轴的转动，其它自由度固定。
H I1s;|仿真|设计|有限元|虚拟仪器9D&L
总之，目的是使结构能实现如图所示的变形。以上提出的加载只是一种比较可以控制的方式。如果您有什么更好的加载方式，只要能实现以上变形，请告诉我。或者您知道见过类似问题，也请不吝赐教。谢谢！3i
IK1NL1M]&Rg
SimWe仿真论坛|9c[ B(K'Z&^P‑X
问题如有进展，我将把更详细的说明承上，与大家一起提高。
尽量不要手写inp文件，出错不容易找。应该用ABAQUS/CAE建模。附上我的模型：
1。Beam.cae: 梁单元，
加力载荷: 不收敛，原因不详。
加位移载荷：
同样出现了类似屈曲的变形，我认为只有让位移载荷随时随刻都垂直于梁的一端，才不会出现“屈曲”，只要有沿着梁轴向的位移，就会“屈曲”。我的位移载荷是随便给的（从0到某个值线性增大），存在轴向位移。
2。solid-pressure.cae:
solid单元, 为减小模型规模，缩短了梁的长度。施加压力pressure,
从结果来看，pressure的方向应该是始终垂直于加载面的，但弯到一半就不收敛了。
如果施加点载荷，即使选择了follow nodal
rotation，载荷方向也不会随变形而改变，因为solid单元没有旋转自由度。如果是beam单元，follow nodal
rotation应该会起作用。
3。solid-diplacement.cae:
模型同上，但不加pressure,
而是加位移载荷（从0到某个值线性增大），同样出现了类似屈曲的变形，原因同Beam.cae。
4. 加载方式
或加位移载荷应该都有可能成功。位移载荷必须随时随刻都垂直于梁的一端，方法可以是：根据你试验的位移曲线定义一个随时间变化的amplitude,
用于定义位移。ABAQUS6.5需要给出多个amplitude上的点，ABAQUS6.6可以直接定义函数曲线。不需要手工定义循环，也不需中止分析来判断方向。
建模时要注意：
- step中的nlgeom要设为on(几何非线性)
-&&&&&&&&&
大变形应该用塑性材料（先把弹性做收敛了，再改为塑性）。塑性问题收敛会更加困难，需要在大应变处细化网格，step的minimum
increment要设得很小，比如10-9.
找出为什么使用梁单元不容易收敛了，ABAQUS/Standard分析几何非线性问题时要使用hybrid单元，例如
附件中的模型使用
B31H单元和力载荷，看到了大变形，并且载荷是始终垂直于梁的一端的，但变形到一定程度就进行不下去了（原来的模型使用B31单元，只一点点微小变形就不收敛了）。看来这样的大变形还是要施加位移载荷才行。
有关刚体的旋转
我知道刚体需要在参考点上定义好6个位移自由度（3转动+3平动），否则就欠约束。
xF&|仿真|设计|有限元|虚拟仪器但是我想让某一个转动自由度，依靠其他物体随动，就是说我事先不能约束这个自由度，这该怎么处理呢？？
问题解决了，真不好意思问了一个这么...的问题：
M*Q[1]o"|仿真|设计|有限元|虚拟仪器
P/jQ-_8P6LSimWe仿真论坛在property里定义好刚体的转动惯量就ok了
如图所示，如要在悬臂梁（3维）的一段施加纯弯矩的话，该怎样施加？施加载荷时我选择了moment，但是在作业分析中总是出错，说5个自由度不能分析等等，但是我怎么也想不通，请高手指教ABAQUS施加moment的时候应注意的问题，我在帮助文件中找了，但是没有找到！
如果我在一端施加集中力的话，那么这个会随着材料变形而一直改变方向吗？或者是一直竖直向下？谢谢赐教！
三维SOLID单元没有旋转自由度，所以不能直接施加moment.
可以先在梁的一端定义一个参考点，在梁的端面和参考点之间建立coupling约束，然后把moment施加在参考点上。
创建载荷时，如果在Edit Load对话框中选中Follow nodal
rotation，就可以使载荷的方向随着梁的变形而变化，同时要在step中设置几何非线性参数 nlgeom 为 on.
我想问在一个圆柱体加旋转速度怎么加啊，我怎么都是在面上加哦，还有两个物体相互接触 ，他们之间的自由度怎么可以耦合哦
可能可以用速度场：load模块，菜单field / Create
看不清你的薄圆筒在那里。两个面之间可以用tie.
error message：The rigid bodies with the reference nodes
contained in node set ErrNodeRefNodeNoRot1Axis have atleast one
axis about which the rotary inertia magnitude is zero, and the
rigid bodies have atleast one rotational degree of freedom which is
unconstrained. Either rotary inertia must be defined at the
reference nodes or all of the rotational degrees of freedom at the
reference nodes must be constrained.
其中提示了两种解决办法：1）定义rotary inertia ，这比较麻烦。2）all of the rotational
degrees of freedom at the reference nodes must be
constrained，但这样芯模就没法转动（不过我觉得你的芯模转动没什么意义，你的芯模和坯料之间既没有摩擦，也没有力使它们紧贴在一起，芯模无法带动坯料转动）。
有一个解决办法是不用刚体单元，而是用普通的c3d8r单元,
但让其弹性模量比正常值大1000倍（太大则可能不收敛），就相当于刚体了。
------------------
这次做得挺不错的。那几个warning都不重要。
首先检查一下你的单位对不对，钢的弹性模量应该是2.1e11 N/m2。
旋轮圆弧处单元太少，这样大的尖角会阻碍坯料的变形。圆弧上至少要有10个单元，最好改用analytical rigid。
C3D8R单元在厚度上至少要有4个单元，才能保证精度。或者用shell单元。
摩擦系数0.2，相当于没有润滑，很难发生相对滑动。
你的AMPLITUDE的含义是速度从0线性增加到给定的速度，而不是始终匀速，不知这是不是你想要的，可以参考Getting
Started with ABAQUS&&7.5 Example:
cargo crane—dynamic loading，Table 7&1 。
目前的field output的输出间隔是20，改为1就可以看到更完整的变形过程。
还不行的话，可以试试减小旋轮和芯模的速度。
问：因為我是從Hypermeh做設定到ABAQUS做求解~我是做碰撞問題的嚕~
F(ak仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM結果產生如下~~請問是那裡的問題?
‑d#o3k-H
b Z‑y)i
~[/@I@e‑|_tU
錯誤訊息~~仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
Underlying element 9487 instance part-1-1 can only have spos, sneg,
e1, e2, e3, e4, edge, or twosided as a face
identifier仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentn@#ev&V
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,y&J(R9u_U
警告訊息~~SimWe仿真论坛
y i[1]c+]2@7| I
Output request e23 is not available in a nonlinear step -- le (log
strain) will be output instead.]3Xy+{/o fuq"x
我要如何去定義速度!
^g+以及碰撞分析!要定義座標嗎?仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent6[[1]n;JQ&Ih+S6v
感覺自己越做越亂!
RI]�Zq}
z-i但很感謝你的幫忙!
ifP&p[1]~ a真希望站主~
WD1u仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM能夠給你積分!SimWe仿真论坛;ke)A(W!W-l'e
我要修改的model~
A1S4X0_&q[1]C&FC+|仿真|设计|有限元|虚拟仪器如附件!$b5vv&Sl"E
在dat文件里看到
&&*SURFACE, NAME = ballset, TYPE
***ERROR: UNDERLYING ELEMENT 9486 CAN ONLY HAVE SPOS, SNEG, E1, E2,
&&&EDGE, OR
TWOSIDED AS A FACE IDENTIFIER
说明这个surface的定义有问题，可以删掉名为 ballset 和 ironset 的 surface.
定义接触时可以用名为iron的surface。
warning先不用管，没有error就可以运行了。
你直接用ABAQUS/CAE建模多好，从hypermesh转过来会出现很多问题。
***ERROR: D.O.F. 6 IS NOT ACTIVE IN THIS MODEL
&LINE IMAGE: ,6,45000
实体单元没有旋转自由度，不能用*INITIAL CONDITIONS, TYPE = VELOCITY
来定义旋转速度场。*INITIAL CONDITIONS下面所有关于自由度4，5，6的数据行都要删掉。定义旋转速度场的方法见
/forum/viewthread.php?tid=710486 。
&***ERROR: INVALID LOAD TYPE ON ELEMENT 12620.
LOAD TYPE P IS NOT VALID FOR
&&&&&&&&&&
ELEMENT TYPE C3D4
不能在c3d4单元上用*DLOAD施加pressure。不过你在所有单元上施加这个pressure有什么意义吗，每个单元上的pressure方向都是不同的。可以定义一个面，然后在面上施加pressure,
也就是说用*DSLOAD.
--------------
定义初始平移速度： 用
*INITIAL CONDITIONS, TYPE = VELOCITY
节点编号或集合,&&1,&
&x方向速度
节点编号或集合,&&2,&
&y方向速度
节点编号或集合,&&3,&
&z方向速度& &
定义初始旋转和平移速度： 用
*INITIAL CONDITIONS, TYPE = ROTATING VELOCITY
First line:
1. Node set or node number.
2. Angular velocity about the axis defined from point a to point b,
where the coordinates of a and b are given below.
3. Global x-component of translational velocity.
4. Global y-component of translational velocity.
5. Global z-component of translational velocity.
Second line:
1. Global x-component of point a on the axis of rotation.
2. Global y-component of point a on the axis of rotation.
3. Global z-component of point a on the axis of rotation.
4. Global x-component of point b on the axis of rotation.
5. Global y-component of point b on the axis of rotation.
6. Global z-component of point b on the axis of rotation.
旋转速度只能定义在reference point上，不能定义在实体单元上。
问：6.5里如何加扭矩或扭转角。
把变形放大四十倍以上，发现变形有问题。一头大一头小，而且端面变成斜的了.
一头大一头小是你把它显示成了这样。固定端没有位移，而把受力端变形放大四十倍，自然是一头大一头小。如果你察看一下受力端面上任意两点的距离，base
distance 和deformed distance(unscaled)应该是几乎一样的。
或者试一下，把扭矩增大40倍，变形系数为1，应该不会一头大一头小。
端面变成斜的：
如果旋转轴是z，且一直保持此方向不变, 则要约束参考点的U1,U2,UR1,UR2
(这是一个经常被忽视的问题，参考点上也要定义边界条件).
另外，如果用coupling, 且不想让受力端面发生任何变形(如同刚体约束)，就要用kinematic coupling，
这样端面就不会变斜;& 如果允许受力端面发生变形，就要用distribute coupling，
这样端面会略微变斜。
问：参考点与面绑定的时候，选SURFACE与选 NODE region 有什么区别.
选SURFACE: 受约束的面是基于单元的：
*Surface, type=element
选 NODE region： 受约束的面是基于节点的：
*Surface, type=NODE
我通常都选SURFACE。
问：鼓式制动器接触应力分析。(简单模型)
工作过程：
二个制动蹄在促动力的作用下，分别围绕各自的销轴向外转动，与制动鼓接触，发生摩擦，其摩擦系数为0.3；(注：本模型中，在未制动时，制动蹄与制动鼓间有间隙)
本模型只是简单的模拟一下制动器的工作过程，看自已所设置的参数用abaqus能否算出结果。(此结果不具有实际的价值)。
材料参数在模型中已经给出。(是在文献中找到的数值)
由于本人对abaqus的应用还是刚入门阶段，对于机械力与温度场耦合问题，现在还不清楚如何处理。所以现在想简化模拟过程。暂时先不考虑温度场问题。如果有人帮我一下哪我可真是感激不尽啊!
想让高手指点一下，分析步如何设置以及边界条件，加载是否合理!
1．&&&&&&&&&&&&&&
以销轴中心对制动蹄施加一个小的转动位移，如何设置？对其施加转动还需要在旋转中心处定义一个多点约束吗？(如：制动蹄与销轴联接处的结点与销轴轴心需要定义多点约束吗？)
2．&&&&&&&&&&&&&&
在静力分析时，不可以有刚体位移，请问：一个刚体六个自由度全部约束，但其中1方向的位移设置为非零的值，在计算时为可以还以计算呢？？？(我晕!)
我的联系方式：sun_
模型挺漂亮的，看你那么勤奋，收集了那么多资料，一定能搞出点东西来。我的建议是这个样子地：
1。静态 / 动态
首先你模拟的目的是什么，比如想看在几秒内能使转动完全停止，就要考虑惯性和速度，需要用动态分析（显式或隐式），模拟整个转动过程。
如果只是想看最大的应力应变，这是研究某个时刻的状况，如同用照相机拍了个定格，是根据这一瞬间的状态来建模，用静态分析，不模拟转动过程。
另外，静态分析也可以模拟运动过程，但这代表的是一种非常缓慢的运动，缓慢得可以不考虑惯性，也没有时间的概念。例如，你可以设定某个实体的ur1等于某个弧度，其他五个自由度为0,
在静态分析的STEP中把INCREMENT设为0.1，就可以看到实体一点点转动过去。
“静态分析不能有刚体位移”，意思是不能因为缺少约束而出现不确定的刚体位移，并不是不可以让实体发生大的位移或转动，否则就没有“几何非线性”的说法了。
2。BC 和 LOAD
动态分析的话，制动蹄受到一个力，使它向外转动，鼓以一定的初始速度转动。
静态分析的话，关键是想清楚在哪里施加BC，在哪里施加载荷,
在同一个节点的同一个自由度上，不能既施加载荷，又定义BC。比如模拟接触面完全接触的状态：
制动蹄: 仍受到那个压紧的力，在销轴孔处约束所有自由度。
鼓：圆孔和圆心上的一个参考点建立coupling，根据鼓的转速，质量和刹车时间，可以计算出一个扭距（这一点我不是很确定），施加在参考点上。约束参考点的其他5个自由度。
要注意的是，静摩擦力 = 压紧力 *
摩擦系数，如果上述扭矩大于静摩擦力的力矩，模型就不再处于静力平衡状态，估计静态分析不会收敛。
2。COUPLING
实体单元的节点本身没有旋转自由度，只有u1,u2,u3，设置这些节点的ur1,ur2,ur3为0或者某个值，都不起作用（尽管CAE提供了这种操作界面）。实体的旋转是整体绕着某一个点旋转，比如上面说的参考点。所以要把一个区域和这个参考点coupling在一起，边界条件和载荷定义在参考点上。
C3D8R单元在厚度方向上至少要有4个单元，才能保证精度。静态分析可以用c3d8i单元，厚度方向2个单元就够了，且节点应力精度比C3D8R高。
目前网格太细，计算时间太长。网格不是越细越好，有限元的主要问题就是用尽量少的计算
时间得到尽量高的精度，只需在接触面和应力集中的部位细化网格。
调试参数阶段，用粗网格可以节约时间（但接触面上网格太粗可能会不收敛），调好后再对关键部位用细网格。建议你根据截面图建一个平面应力模型，调通了再改为三维实体。
如果关心危险部位的应力，在拐弯处要倒圆角。
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。