Математическая модель реконструкции остаточных напряжений

@xd000 · Регистрация: 09.11.2022

Студворк — интернет-сервис помощи студентам

Есть функция, заданная таблично.
Надо найти z_0 такое, что delta_x=0. Это я сделал:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
import math
import numpy as np
import matplotlib.pyplot as plt
 
# Исходные данные
z = np.array([0, 0.01, 0.02, 0.04, 0.06, 0.08, 0.1, 0.14, 0.18])
delta_x = np.array([-1032, -1012.721, -968.3837, -799.7408, -580.272, -368.2734, -202.092, -31.10915, 11.3538])
 
# Построение точек
plt.plot(z, delta_x, 'o')
plt.xlabel('z')
plt.ylabel('delta_x')
plt.grid()
plt.show()
 
# Нахождение приближенного значения z_0
z_0 = np.interp(0, delta_x, z)
print('z_0 =', z_0)

Потом по формулам, которые нам дали, нужно методом перебора вычислить b, при котором f(b)=0
Моя попытка:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
def integral(b): #интеграл от 0 до H=0.18
    int_z = 0; z = 0
    for i in range(18):
        int_z += math.exp(-(z ** 2 / b ** 2))
        z += 0.01
    return int_z * 0.01
 
 
def f(b):
    return math.exp(-(z_0 ** 2 / b ** 2)) * 0.18 - integral(b)
 
 
b = 0.01; func = f(b)
 
while func * f(b) > 0:
    b += 0.01
    func = f(b)
 
b = (2 * b - 0.01) / 2
 
print("b =", b)

Считает до бесконечности, не могу найти ошибку, формулы точно переписаны правильно
Знаю, что в ответе b должно получится 0.08

@passant · 10.03.2024, 14:28

Формулы-то вы написали может быть и правильно. А вот код отлаживать вас учили?
Цикл while второй части скрипта начинается со значений b=0.01 и func = f(b)=-0.013863186024133263
В самом цикле какие бы значения func и f(b) вы не вычисляли (а они при окончании каждой итерации цикла всегда одинаковы) условие цикла (func * f(b) > 0 ) будет выполняться ВСЕГДА, потому-как слева у вас квадрат.
Вопрос - как вы себе представляете, какие значения должны иметь func и f(b) что бы приведенное выше условие НЕ выполнилось (т.е. этот цикл закончился)?

Рыжий Лис · 10.03.2024, 19:50

Что-то я этот код уже видел в другой теме…

Добавлено через 1 минуту
Так ТС этот же! Нахождение параметра

Новые блоги и статьи Все статьи Все блоги /
Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .	модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .

@passant 578 / 411 / 69 Регистрация: 09.01.2018 Сообщений: 1,363
	10.03.2024, 14:28
	Формулы-то вы написали может быть и правильно. А вот код отлаживать вас учили? Цикл while второй части скрипта начинается со значений b=0.01 и func = f(b)=-0.013863186024133263 В самом цикле какие бы значения func и f(b) вы не вычисляли (а они при окончании каждой итерации цикла всегда одинаковы) условие цикла (func * f(b) > 0 ) будет выполняться ВСЕГДА, потому-как слева у вас квадрат. Вопрос - как вы себе представляете, какие значения должны иметь func и f(b) что бы приведенное выше условие НЕ выполнилось (т.е. этот цикл закончился)? 0

Рыжий Лис Просто Лис 5973 / 3735 / 1099 Регистрация: 17.05.2012 Сообщений: 10,791 Записей в блоге: 9
	10.03.2024, 19:50
	Что-то я этот код уже видел в другой теме… Добавлено через 1 минуту Так ТС этот же! Нахождение параметра 0