离散kirchhoff理论三角形单元及其在大变形和大应变分析中的应用.pdfVIP

第十届全国反应堆结构力学会议 998.11 离散kirchhoff理论三角形单元及其在大变形 和大应变分析中的应用 李朋洲 谷方毓 中‘国核动力研究设计院二所) 摘要 本文首先介绍了离散Kirchhoff理论三角形单元(DKT)的理论假设和构造方法。基于更改的 拉格朗日列式法(UpdatedLagrangianFormulation)，推导了该单元用于壳体大变形大应变静动态分析 的有限元控制方程、刚度矩阵和节点力表达式，在每一加载步 (或时间步)，不仅单元节点坐标和面 积，而且单元厚度也得到了修正。还编制了该单元的计算机程序，数值算例表明程序稳定、可靠，具 有较高精度。 关性词:DKT单元、UL法、大变形大应变 0前言 对壳体结构进行分析，低阶的三角形单元是很有效的。而其中离散的Kirchhoff理论三角 形单元 (DKT)是非常有效和可靠的，此为文献[[1]通过比较DKT单元和ADINA程序中其它 类型三角形单元而得到的结论。在文献[21中，DKT单元已被推广到大位移和大转动但是小应 变的薄壳分析中，在文献[3]该单元又被应用于一般壳体的弹塑性分析，但是应用的是合力和 合力矩表达的近似屈服条件。本文将在他们的基础上，将DKT单元应用于大位移、大转动和 大应变的壳体结构的瞬态动力分析中。 1DKT弯曲单元 由于该单元薄膜刚度矩阵与通常常应变三角形单元同，此处不再重述，主要讨论它的弯 曲特性。众所周知，在Kirchhoff薄板弯曲理论中，去构造一个位移协调的三节点弯曲单元是 很困难的，对于DKT单元此困难得到了克服，这种单元中独立地假设了:中面法线转动 几，几在单元中二次变化;横向位移沿二角形边界呈三次变化和沿三角形边界对各边法线转 角几强加线性变化规律。同时将Kirchhoff假设强加于三角形的节点及各边中点以保证好的 收敛性。 如图1所示，、、2,3为单元的角结点，4,5,6为单元各边 蜘 分别为各边的 切向和法向，而Tsi表示边ii的方向与x轴之间的夹角·该单元3个 伞 处共有9个自由度 a o z 二*T一卜、0s1。;，w2)x291Y2w3 913{ (1) 首先假设0氏在单元内部呈二次变化规律，它们可用两个角结点和3个边中点的 Bx,B，表示为: 9x=艺N,B,,B，一艺N,B, (2) 其中形函数N与6节点三角形平面单元相同。 引入下面的约束条件: (1),在角节点和各边中点强加Kirchhoff假设，即: PMkirchoffJN}Ll三角形单元及其在大变形和大应变分析中的应用 、一ll}二一a(civ 、一，”一(dw )2(在单元边界ii,假定挠度w呈三次变化规律，它完全由该边界两端的WW,Wf,wri 所决定，从而可将上式中Wsk表示成 3 nsk“一可袱一4万Ws.十可WI一4W. (4) 上式 中 、代表各边ii的中点，1。为边ii的长度，1。一叫十巾气 凡 =x，一x,，Y=戈一Y,“ (3)假定转角氏沿三角形边界线性变化，因而各边中点B*可用对应边端点的相应值 表示，即: 氏*=(B,;+0J)/2 (5) 注意到ii边上的几何关系 (见图1) }O