非线性解析法2——程序那个.docVIP

Appendix　C　　非线性解析程序说明 平面应变状态弹性条件的应力-应变关系 ： 　　　　　　　　　(C1) 　　　　　　　　　　　(C2) 　　　　　 図.C1　6节点三角形单元 单元中任意点位移：　　　　　　　(C3) 　　　　　　　　　　(C4) 　　　　　　　　　　　　　(5) 　　　　　　　　　　　　　　(C6) 　　　　　　　　　　　(C7) 在平面应变条件下　 　　　　(C8) 　　　　　　　　　　　(C9) 　　　　　　　　　　　　　　　　　　(C10) 　　　　　　　　　　　　　　　　　　(C11) 根据虚功原理(内部) 　　　　　(C12) 为初期应力的等价节点力。 对于公式(C12)的积分，采用在三角形内的高斯5次数值积分来评价。 图.C2　数值积分点 表.C1　三角形面積積分 L1 L2 L3 重み(w) (a) 1/3 1/3 1/3 0.225 (b) ?? ?? ?? 0 (c) ?? ?? ?? 0 (d) ?? ?? ?? 0 (e) ?? ?? ?? 0 (f) ?? ?? ?? 0 (g) ?? ?? ?? 0 在这里、 ?????????????、????????????? ?????????????、????????????? 　对（(a)wa+ … (g)wg）的总和乘以三角形的面积。 外部做功 ①　节点荷载（集中力） 　　　　　　　　　　　　　　(C13) B：单元厚度。一般情况取为B=1。 ②　分布荷载（分布力） (C14) 分布荷载的等价节点力 　　　　　　　　(C14′) 对ij边进行积分时(B=1)、L3=0 　　　　　　　　　　　　　　(C15) (C16) (C17) (C18) (C19) (C20) (C21) ③物体力（重力）（B=1） 　　　　　　　　　　　　　(C22) 三角形单元的面积 　　　　　　(C23) 　　　　(C24) 　　　　(C25) 　　　　　　　　(C26) 　　　　　　　　(C27) 　　　(C28) 非线性解析 初期应力设定：　　　　　　　　　　　　　（C29） Jaky式的土压系数：　　　　　　　　　　　（C30） 　 在一维压缩条件下，不发生水平方向应变（）、根据虎克定律 　　　　　　　　　　　　　　　　　　（C31） 　　　　　　　　　　　　　　　　　　　（C32） 泊松比：　　　　　　　　　　　　　　　（C33） 解析时，土提比重取为。 各应力状态下的杨氏模量：　　　　　　　　（C34）　　 Mohr-Columb破坏准则：　　　　　　　　　　　　（C35）　 在弹性有限元计算中，杨氏模量取为E0=20Mpa。在非线性计算中，使初期杨氏模量在深度方向的平均值与弹性解析的杨氏模量相等而确定的。由此， 　　　　　　　　　　　　　　　　　　　（C36） 　　　　　　　　　　　　　　　　　　　　　　　　　（C37） 　　　　　　　　　　　　　　　　　　　　　　　　　　　（C38） 程序流程图 本程序所采用的单位为，长度：m 　　　　　　　　　　　荷载：tf 程序名：duchanchang.f 程序概略： 　1)　计算条件：平面应变条件 　2)　计算精度：双倍精度 　3)　单元：三角形6节点 　4)　非线性解法：单纯增量计算 　5） 外力：集中荷载、重力、分布荷载 　6） 初期应力：根据自重自动确立 　7)　边界条件：变位条件，荷载条件 duchanchang.f所包含的子程序的机能 　①　PLANSM：单元刚性矩阵计算 　②　ASMAT： 单元刚性矩阵写入到整体刚性矩阵 　③　DLINCB：带状矩阵的线性方程式求解 　④　DMAT: 做成单元的弹性应力-应变关系矩阵[D]（参照式(C2)） 　⑤　CINV：在做成[c] 矩阵之后，求其逆矩阵[c-1]（参照式(C5)） 　⑥　STRMX：做成[A] 矩阵（参照式(C9)） 　⑦　DINV2：求解逆矩阵 　⑧　PRST：主应力计算 duchanchang.f内所包含的子程序间的相互关系 程序内所含符号说明 NODT ：总节点数 NELT ：总单元数 NRJE ：在应力解析中的总拘束点数 NRJJE ：在应力解析中的无效点总数 NP ：集中荷载作用点数