用户登录

你在这里

有助于二维共旋光束的制定

yawlou的照片

我在matlab中实现了一个二维共旋转桁架公式。我得到相同的结果OpenSEES。

现在我要实现一个二维共旋转光束的公式。我遵循了克里斯菲尔德第一卷书中的简单公式。看起来这应该是桁架公式的一个相当直接的扩展。我重新推导了克里斯菲尔德所做的一切,得到了同样的结果,但还是有些地方不正确。

有人实现过简单的二维共旋转梁单元公式吗?如果你们有一个可行的二维共旋光束公式,有没有人愿意通过与我交换几封电子邮件来帮助我缩小问题范围?我们可以做一个简单的悬臂,用一个单梁元素建模。

或者,有没有人能给我指出一个二维共旋转梁单元公式的工作版本,它可以很容易地让我输出中间值来与我的matlab程序进行比较?我知道OpenSees有这样的实现,这就是为什么我知道我的程序不工作。我会查看OpenSees内部,但它非常复杂,很难看到我需要找到我的错误。

如果你愿意帮忙,请告诉我。

谢谢,

路易

评论

亲爱的路易先生:

我一直在用克里斯菲尔德的二维光束旋转公式。我没有发现明显的错别字或其他什么。我同意他的书的一些版本/印刷品的印刷质量很差,特别是对于下标来说令人困惑。

你有没有试过完全按照Crisfield给出的编码元素矩阵?

如果您想比较结果,一个验证结果的来源是“Int”。j .号码。方法工程17(1981)145-153”。

如果您仍然遇到问题,那么您可以得出结论,这只是程序调试。

Kodanda

yawlou的照片

谢谢Kodanda的回复。

我同意,我不认为克里斯菲尔德的推导有任何错别字。我通过推导得到了和他相同的结果。

我的问题确实是一个编程错误,或者可能是我没有理解一些关键的小部分。

我试图在matlab中通过扩展我的桁架版本的二维旋转来编程梁的公式。然而不知怎么的,我的平衡木程序不工作了。现在,我一直在试图让我的结果与关闭NR迭代的OpenSees的结果相匹配。

第一个线性步骤当然是匹配的,但是第二步在几何体的旋转更新之后不匹配我的结果。所以我一直在尝试两个数据点与一个单一的光束元素。

我会继续努力。

问候,

路易

嗨,更新元素有点让人困惑。如何给出元素斜率的增量?我有一个用mathematica写的代码,你可以通过电子邮件联系:

hakan.tanriover@gmail.com

祝你好运。

yawlou的照片

感谢一路走来帮助过我的人。

给那些感兴趣的人。

我写了一篇论文,描述了二维旋转光束的公式。我的方法与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

Dario Schiano Moriello的照片

尊敬的耀先生:

非常感谢您分享您关于梁和桁架的旋转公式的论文,它们对我研究这门学科非常有帮助。

我刚刚开始攻读博士学位,我的研究计划的第一步将涉及在ABAQUS中实现旋转公式,但是,对于这个主题和软件都是新的,我不确定如何进行。您是否了解将此公式应用于ABAQUS或其他软件分析的现有工作?

我从这个博客上读到,你已经编写了一个MATLAB算法来解决梁的简单问题;这将是非常有帮助的,如果你能请分享这个工作的一部分。

非常感谢您的考虑。
希望收到您的来信,

问候,
达里奥

订阅“二维共旋光束公式帮助”评论

最近的评论

更多的评论

辛迪加

认购辛迪加