如何使用双重非线性（几何和材料非线性）！

kanghui

好多论文进行了双重非线性分析，说可以采用增量M- N- R法求解，“双重非线性就是同时考虑几何和材料非线性，即将几何非线性迭代嵌入材料非线性的增量法中，在每级荷载增量中折减刚度不变，并用N-R方法考虑几何非线性问题。”
请问在ANSYS程序中如何实现？有没有这方面的ADPL程序可供参考？谢谢大家了

kanghui

有的书这样写到：在实际求解中采用嵌套的方法考虑双重非线性问题。结构受力全过程分析采用混合法求解，将荷载分成若干增量,给定参数,由程序控制加载步长,在各个增量荷载上进行迭代。用增量法考虑材料的非线性影响,将几何非线性嵌入材料非线性的增量法之中,在每级用修正的Newton—Raphson方法考虑几何非线性问题。其求解步骤如下:
(1)在某级荷载作用下,调用材料非线性EAI子程序迭代得EAi、EIi;
(2)把EAi、EIi代入线性刚度阵[K0],求出第一次迭代后的位移{δ}e及单元内力Ni、Mi;
(3)再由所求内力求得单元中点刚度EA′i、EI′i;
(4)求出杆件单元几何刚度矩阵[Kσ] ,从而定出切线刚度矩阵;(5)应用几何非线性迭代求出考虑几何非线性影响后的杆端位移{δ1}e、杆端力{F1}eJ;(6)再增加一级荷载,应用上述方法继续计算,直至某一断面所承受的弯矩超过其所能达到最大弯矩,输出结果。
请问这样的程序是利用ＡＮＳＹＳ的命令流编的，还是自己用其他计算机语言如：Ｃ或者ＦＯＴＲＯＮ编写的？那位高手这样做过，恳请帮忙！
EMAIL：zhadaonan@163.com, QQ:782745301   相互交流》