你在这里
有助于二维共旋光束的制定
星期一,2009-02-23 22:11 -yawlou
我在matlab中实现了一个二维共旋转桁架公式。我得到相同的结果OpenSEES。
现在我要实现一个二维共旋转光束的公式。我遵循了克里斯菲尔德第一卷书中的简单公式。看起来这应该是桁架公式的一个相当直接的扩展。我重新推导了克里斯菲尔德所做的一切,得到了同样的结果,但还是有些地方不正确。
有人实现过简单的二维共旋转梁单元公式吗?如果你们有一个可行的二维共旋光束公式,有没有人愿意通过与我交换几封电子邮件来帮助我缩小问题范围?我们可以做一个简单的悬臂,用一个单梁元素建模。
或者,有没有人能给我指出一个二维共旋转梁单元公式的工作版本,它可以很容易地让我输出中间值来与我的matlab程序进行比较?我知道OpenSees有这样的实现,这就是为什么我知道我的程序不工作。我会查看OpenSees内部,但它非常复杂,很难看到我需要找到我的错误。
如果你愿意帮忙,请告诉我。
谢谢,
路易
![订阅“二维共旋光束公式帮助”评论](http://m.limpotrade.com/misc/feed.png)
评论
亲爱的路易先生,我有
亲爱的路易先生:
我一直在用克里斯菲尔德的二维光束旋转公式。我没有发现明显的错别字或其他什么。我同意他的书的一些版本/印刷品的印刷质量很差,特别是对于下标来说令人困惑。
你有没有试过完全按照Crisfield给出的编码元素矩阵?
如果您想比较结果,一个验证结果的来源是“Int”。j .号码。方法工程17(1981)145-153”。
如果您仍然遇到问题,那么您可以得出结论,这只是程序调试。
Kodanda
这是一些错误在我的程序,不确定什么虽然…
谢谢Kodanda的回复。
我同意,我不认为克里斯菲尔德的推导有任何错别字。我通过推导得到了和他相同的结果。
我的问题确实是一个编程错误,或者可能是我没有理解一些关键的小部分。
我试图在matlab中通过扩展我的桁架版本的二维旋转来编程梁的公式。然而不知怎么的,我的平衡木程序不工作了。现在,我一直在试图让我的结果与关闭NR迭代的OpenSees的结果相匹配。
第一个线性步骤当然是匹配的,但是第二步在几何体的旋转更新之后不匹配我的结果。所以我一直在尝试两个数据点与一个单一的光束元素。
我会继续努力。
问候,
路易
我也许能帮上忙……
嗨,更新元素有点让人困惑。如何给出元素斜率的增量?我有一个用mathematica写的代码,你可以通过电子邮件联系:
hakan.tanriover@gmail.com
祝你好运。
二维旋转梁公式(pdf文件)
感谢一路走来帮助过我的人。
给那些感兴趣的人。
我写了一篇论文,描述了二维旋转光束的公式。我的方法与Crisfield提供的方法密切相关。然而,我已经尽力展示了推导过程中的所有步骤,并仔细地解释了每一件事。给出了在计算机程序中实现的有用算法。此外,使用de Souza的直接技术克服了任意大旋转的问题。解决了代表性基准问题。
该文件可在以下地点找到:
http://people.wallawalla.edu/~louie.yaw/co_rotational%20papers.html
问候,
路易
你好,路易,我在看书
你好,
路易,我在读你关于二维旋转光束公式的论文。谢谢你的工作,这是非常有用的,因为我的知识几何非线性是非常混乱的。
我有两个问题:
-这个理论适用于最初轻微弯曲的光束吗?(我会简单地将光束离散成足够数量的直线单元。)
-它是否适用于受跟随力影响的光束?
问候,
弗雷德
弧长法
亲爱的先生,
请告诉我E Ramm的弧长法的算法。我正在使用以下类型的算法,但我没有得到正确的结果----
你{0},λ= 0
Fext
N个步数
Kg(刚度u{0})
u1=(Kg)^-1 * text '
[L D] =低密度脂蛋白(公斤)
m =侦破(D)
arc_length =规范(u1)
λ= 0.1 * arc_length /√(u1’* u1 + 1)
如果m < 0
λ=λ
结束
u =λ* u1
迭代:
公斤(u {0} + u)
弗林特(u {0} + u)
耐火的= Fint-lamda * Fext
你=发票(公斤)*耐火的”
* Fext u1 =发票(公斤)
del_lamda = - (u ' * ur) / (u ' * u1)
del_ur =你+ del_lamda * u1
u = = u + del_ur
λ=λ+ del_lamda
如果规范(即< 10 ^ 5)
结束
u {0} = {0} + u
如果λ> = 2
结束
所以请帮助我……急切地等待解决方案。
谢谢
Md Rushdie Ibne Islam
实施合作制定
尊敬的耀先生:
非常感谢您分享您关于梁和桁架的旋转公式的论文,它们对我研究这门学科非常有帮助。
我刚刚开始攻读博士学位,我的研究计划的第一步将涉及在ABAQUS中实现旋转公式,但是,对于这个主题和软件都是新的,我不确定如何进行。您是否了解将此公式应用于ABAQUS或其他软件分析的现有工作?
我从这个博客上读到,你已经编写了一个MATLAB算法来解决梁的简单问题;这将是非常有帮助的,如果你能请分享这个工作的一部分。
非常感谢您的考虑。
希望收到您的来信,
问候,
达里奥